Matrix variate Macdonald distribution

Abstract

In this article, we generalize the univariate Macdonald distribution to the matrix case and give its derivation using matrix variate gamma distribution. We study several properties such as cumulative distribution function, marginal distribution of submatrix, triangular factorization, moment generating function, and expected values of several functions of the Macdonald matrix. Some of these results are expressed in terms of special functions of matrix arguments and zonal polynomials.     

Bimatrix Variate Beta Type IV Distribution: Relation to Wilks's Statistic and Bimatrix Variate Kummer-Beta Type IV Distribution

In this article, the bimatrix variate beta Type IV distribution is derived from independent Wishart distributed matrix variables. We explore specific properties of this distribution which is then used to derive the exact expressions of the densities of the product and ratio of two dependent Wilks's statistics and to define the bimatrix Kummer-beta Type IV distribution.     

On matrix variate skew-normal distributions

Two new families of matrix variate distributions are introduced. They are based on matrix normal distribution and yet can be used to model data involving skewness. The properties of the two families are investigated. Among others, the marginals, conditionals, stochastic representation, linear and quadratic forms are studied.     

A new class of matrix variate elliptically contoured distributions

Summary A new class of matrix variate elliptically contoured distributions is defined. Properties of this class of distributions are studied. Examples of distributions which belong to this class are also presented. Rose-Hulman Institute of Technology Research supported by the FRC Major Grant, Bowling Grant, Bowling Green State University.     

The power series skew normal class of distributions

A class of power series skew normal distributions is introduced by generalizing the geometric skew normal distribution of Kundu. Various mathematical properties are derived and estimation addressed by the method of maximum likelihood. The data application of Kundu [Sankhyā B, 76, 2014, 167–189] is revisited and the proposed class is shown to provide a better fit.     

A new method for generating distributions with an application to exponential distribution

A new method has been proposed to introduce an extra parameter to a family of distributions for more flexibility. A special case has been considered in detail, namely one-parameter exponential distribution. Various properties of the proposed distribution, including explicit expressions for the moments, quantiles, mode, moment-generating function, mean residual lifetime, stochastic orders, order statistics, and expression of the entropies, are derived. The maximum likelihood estimators of unknown parameters cannot be obtained in explicit forms, and they have to be obtained by solving non linear equations only. Further, we consider an extension of the two-parameter exponential distribution also, mainly for data analysis purposes. Two datasets have been analyzed to show how the proposed models work in practice.     

Inference for two-stage adaptive treatment strategies using mixture distributions

Summary.  Treatment of complex diseases such as cancer, leukaemia, acquired immune deficiency syndrome and depression usually follows complex treatment regimes consisting of time varying multiple courses of the same or different treatments. The goal is to achieve the largest overall benefit defined by a common end point such as survival. Adaptive treatment strategy refers to a sequence of treatments that are applied at different stages of therapy based on the individual's history of covariates and intermediate responses to the earlier treatments. However, in many cases treatment assignment depends only on intermediate response and prior treatments. Clinical trials are often designed to compare two or more adaptive treatment strategies. A common approach that is used in these trials is sequential randomization. Patients are randomized on entry into available first-stage treatments and then on the basis of the response to the initial treatments are randomized to second-stage treatments, and so on. The analysis often ignores this feature of randomization and frequently conducts separate analysis for each stage. Recent literature suggested several semiparametric and Bayesian methods for inference related to adaptive treatment strategies from sequentially randomized trials. We develop a parametric approach using mixture distributions to model the survival times under different adaptive treatment strategies. We show that the estimators proposed are asymptotically unbiased and can be easily implemented by using existing routines in statistical software packages.     

On random variate generation for the generalized hyperbolic secant distributions

We give random variate generators for the generalized hyperbolic secant distribution and related families such as Morris's skewed generalized hyperbolic secant family and a family introduced by Laha and Lukacs. The rejection method generators are uniformly fast over the parameter space and are based upon a complex function representation of the distributions due to Harkness and HarknessThe author's research was sponsored by NSERC Grant A3456 and FCAR Grant 90-ER-0291.     

Moments and quadratic forms of matrix variate skew normal distributions

ABSTRACT

In 2007, Domínguez-Molina et al. obtained the moment generating function (mgf) of the matrix variate closed skew normal distribution. In this paper, we use their mgf to obtain the first two moments and some additional properties of quadratic forms for the matrix variate skew normal distributions. The quadratic forms are particularly interesting because they are essentially correlation tests that introduce a new type of orthogonality condition.     

Regular mles for nonregular distributions

Distributions whose extremity values of the support depend on unknown pa¬rameters are usually known as nonregular distributions. In most cases, the MLEs for these parameters cannot be obtained by differentiation. Familiar examples are the uniform distribution on the interval (0,0) and the truncated exponential distribution with truncation parameter 0. However, there exist distributions whose extremity points of the support depend on unknown pa¬rameters, which nevertheless are regular in the sense that the MLEs can be obtained by differentiation. This note provides a method of constructing such nonregular distributions with regular MLEs.     

On a multivariate Lindley distribution

This article proposes a multivariate extension of the generalized Lindley distribution introduced by Abouammoh et al. (2015 Abouammoh, A.M., Alshangiti Arwa, M., Ragab, I.E. (2015). Reliability estimation of the generalized Lindley distributions. J. Stat. Comput. Simul. 85(17):3662–3678.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). The proposed model is based on a probabilistic construction in which several Lindley random variables are connected by a common shock. Many statistical properties of this new distribution are explored. In particular, explicit forms of product moments, moment-generating function, and conditional moments are derived. Explicit expressions of moment-based estimators of the underlying parameters of the proposed model are established. Estimation using the maximum likelihood procedure is also investigated. Simulations attesting to the quality of the proposed estimators are presented. Application of the proposed model in reliability is also discussed.     

On compounded geometric distributions and their applications

Here, we introduce two-parameter compounded geometric distributions with monotone failure rates. These distributions are derived by compounding geometric distribution and zero-truncated Poisson distribution. Some statistical and reliability properties of the distributions are investigated. Parameters of the proposed distributions are estimated by the maximum likelihood method as well as through the minimum distance method of estimation. Performance of the estimates by both the methods of estimation is compared based on Monte Carlo simulations. An illustration with Air Crash casualties demonstrates that the distributions can be considered as a suitable model under several real situations.     

Maximum likelihood estimation of dirichlet distributions

Global concavity of the likelihood function is proved by means of an inequality involving the trigamma function. The computation of maximum likelihood estimates is discussed.     

A new lifetime model with different types of failure rate

A new class of lifetime distributions, which can exhibit with upside-down bathtub-shaped, bathtub-shaped, decreasing, and increasing failure rates, is introduced. The new distribution is constructed by compounding generalized Weibull and logarithmic distributions, leading to improvement on the lifetime distribution considered in Dimitrakopoulou et al. (2007 Dimitrakopoulou, T., K. Adamidis, and S. Loukas. 2007. A lifetime distribution with an upside-down bathtub-shaped hazard function. IEEE Transactions on Reliability 56:308–11.[Crossref], [Web of Science ®] , [Google Scholar]) by having no restriction on the shape parameter and extending the result studied by Tahmasbi and Rezaei (2008 Tahmasbi, R., and S. Rezaei. 2008. A two-parameter lifetime distribution with decreasing failure rate. Computational Statistics and Data Analysis 52:3889–901.[Crossref], [Web of Science ®] , [Google Scholar]) in the general form. The proposed model includes the exponential–logarithmic and Weibull–logarithmic distributions as special cases. Various statistical properties of the proposed class are discussed. Furthermore, estimation via the maximum likelihood method and the Fisher information matrix are discussed. Applications to real data demonstrate that the new class of distributions is more flexible than other recently proposed classes.     

An extension of the generalized linear failure rate distribution

In this article, we introduce a new extension of the generalized linear failure rate (GLFR) distributions. It includes some well-known lifetime distributions such as extension of generalized exponential and GLFR distributions as special sub-models. In addition, it can have a constant, decreasing, increasing, upside-down bathtub (unimodal), and bathtub-shaped hazard rate function (hrf) depending on its parameters. We provide some of its statistical properties such as moments, quantiles, skewness, kurtosis, hrf, and reversible hrf. The maximum likelihood estimation of the parameters is also discussed. At the end, a real dataset is given to illustrate the usefulness of this new distribution in analyzing lifetime data.     

Compounded inverse Weibull distributions: Properties,inference and applications

ABSTRACT

In this paper two probability distributions are analyzed which are formed by compounding inverse Weibull with zero-truncated Poisson and geometric distributions. The distributions can be used to model lifetime of series system where the lifetimes follow inverse Weibull distribution and the subgroup size being random follows either geometric or zero-truncated Poisson distribution. Some of the important statistical and reliability properties of each of the distributions are derived. The distributions are found to exhibit both monotone and non-monotone failure rates. The parameters of the distributions are estimated using the expectation-maximization algorithm and the method of minimum distance estimation. The potentials of the distributions are explored through three real life data sets and are compared with similar compounded distributions, viz. Weibull-geometric, Weibull-Poisson, exponential-geometric and exponential-Poisson distributions.     

Testing the equality means of several log-normal distributions

Researches propose various methods for comparing the means of two log-normal distributions. Some of these methods have been recently extended to test the equality means of several log-normal populations. Investigations show that none of the established methods is satisfactory. In this article, we provide three methods based on the computational approach test, which is a parametric bootstrap approach, for testing the means of several log-normal distributions. Further, we compare our methods with the existing methods through Monte Carlo simulation. The numerical results show that the Type I errors of these procedures are satisfactory regardless of the sample size, number of populations, and the true parameters. Finally, we explain the considered methods by real examples.     

On MLE of a nonlinear discriminant function from a mixture of two Gompertz distributions based on small sample size

The property of identifiability is an important consideration on estimating the parameters in a mixture of distributions. Also classification of a random variable based on a mixture can be meaning fully discussed only if the class of all finite mixtures is identifiable. The problem of identifiability of finite mixture of Gompertz distributions is studied. A procedure is presented for finding maximum likelihood estimates of the parameters of a mixture of two Gompertz distributions, using classified and unclassified observations. Based on small sample size, estimation of a nonlinear discriminant function is considered. Throughout simulation experiments, the performance of the corresponding estimated nonlinear discriminant function is investigated.     

On length and area-biased Maxwell distributions

This article addresses the various properties and different methods of estimation of the unknown parameter of length and area-biased Maxwell distributions. Although, our main focus is on estimation from both frequentist and Bayesian point of view, yet, various mathematical and statistical properties of length and area-biased Maxwell distributions (such as moments, moment-generating function (mgf), hazard rate function, mean residual lifetime function, residual lifetime function, reversed residual life function, conditional moments and conditional mgf, stochastic ordering, and measures of uncertainty) are derived. We briefly describe different frequentist approaches, namely, maximum likelihood estimator, moments estimator, least-square and weighted least-square estimators, maximum product of spacings estimator and compare them using extensive numerical simulations. Next we consider Bayes estimation under different types of loss function (symmetric and asymmetric loss functions) using inverted gamma prior for the scale parameter. Furthermore, Bayes estimators and their respective posterior risks are computed and compared using Markov chain Monte Carlo (MCMC) algorithm. Also, bootstrap confidence intervals using frequentist approaches are provided to compare with Bayes credible intervals. Finally, a real dataset has been analyzed for illustrative purposes.     

Parametric estimation under a class of multivariate distributions

In this paper the parameters of some members of a class of multivariate distributions, which was constructed by AL-Hussaini and Ateya (2003), are estimated by using the maximum likelihood and Bayes methods.