FINITE MIXTURE OF CERTAIN DISTRIBUTIONS

ABSTRACT

A two parameter extended form of standard gamma function is suggested which provide extra flexibility to the density function over positive range. A finite mixture of beta distribution is defined by using the suggested extended form. The shape of density function of extended gamma distribution and also that of finite mixture of beta distribution for various values of the parameters are shown. Inverted distribution of extended gamma and that of finite mixture of beta distribution are given.     

Discrete Linnik Weibull distribution

Abstract

We introduce a new family of distributions using truncated discrete Linnik distribution. This family is a rich family of distributions which includes many important families of distributions such as Marshall–Olkin family of distributions, family of distributions generated through truncated negative binomial distribution, family of distributions generated through truncated discrete Mittag–Leffler distribution etc. Some properties of the new family of distributions are derived. A particular case of the family, a five parameter generalization of Weibull distribution, namely discrete Linnik Weibull distribution is given special attention. This distribution is a generalization of many distributions, such as extended exponentiated Weibull, exponentiated Weibull, Weibull truncated negative binomial, generalized exponential truncated negative binomial, Marshall-Olkin extended Weibull, Marshall–Olkin generalized exponential, exponential truncated negative binomial, Marshall–Olkin exponential and generalized exponential. The shape properties, moments, median, distribution of order statistics, stochastic ordering and stress–strength properties of the new generalized Weibull distribution are derived. The unknown parameters of the distribution are estimated using maximum likelihood method. The discrete Linnik Weibull distribution is fitted to a survival time data set and it is shown that the distribution is more appropriate than other competitive models.     

The beta Weibull-geometric distribution

A new five-parameter distribution called the beta Weibull-geometric (BWG) distribution is proposed. The new distribution is generated from the logit of a beta random variable and includes the Weibull-geometric distribution of Barreto-Souza et al. [The Weibull-geometric distribution, J. Stat. Comput. Simul. 81 (2011), pp. 645–657], beta Weibull (BW), beta exponential, exponentiated Weibull, and some other lifetime distributions as special cases. A comprehensive mathematical treatment of this distribution is provided. The density function can be expressed as an infinite mixture of BW densities and then we derive some mathematical properties of the new distribution from the corresponding properties of the BW distribution. The density function of the order statistics and also estimation of the stress–strength parameter are obtained using two general expressions. To estimate the model parameters, we use the maximum likelihood method and the asymptotic distribution of the estimators is also discussed. The capacity of the new distribution are examined by various tools, using two real data sets.     

The beta exponentiated Weibull distribution

The Weibull distribution is one of the most important distributions in reliability. For the first time, we introduce the beta exponentiated Weibull distribution which extends recent models by Lee et al. [Beta-Weibull distribution: some properties and applications to censored data, J. Mod. Appl. Statist. Meth. 6 (2007), pp. 173–186] and Barreto-Souza et al. [The beta generalized exponential distribution, J. Statist. Comput. Simul. 80 (2010), pp. 159–172]. The new distribution is an important competitive model to the Weibull, exponentiated exponential, exponentiated Weibull, beta exponential and beta Weibull distributions since it contains all these models as special cases. We demonstrate that the density of the new distribution can be expressed as a linear combination of Weibull densities. We provide the moments and two closed-form expressions for the moment-generating function. Explicit expressions are derived for the mean deviations, Bonferroni and Lorenz curves, reliability and entropies. The density of the order statistics can also be expressed as a linear combination of Weibull densities. We obtain the moments of the order statistics. The expected information matrix is derived. We define a log-beta exponentiated Weibull regression model to analyse censored data. The estimation of the parameters is approached by the method of maximum likelihood. The usefulness of the new distribution to analyse positive data is illustrated in two real data sets.     

Statistical analysis of the Lomax–Logarithmic distribution

In this paper we introduce a three-parameter lifetime distribution following the Marshall and Olkin [New method for adding a parameter to a family of distributions with application to the exponential and Weibull families. Biometrika. 1997;84(3):641–652] approach. The proposed distribution is a compound of the Lomax and Logarithmic distributions (LLD). We provide a comprehensive study of the mathematical properties of the LLD. In particular, the density function, the shape of the hazard rate function, a general expansion for moments, the density of the rth order statistics, and the mean and median deviations of the LLD are derived and studied in detail. The maximum likelihood estimators of the three unknown parameters of LLD are obtained. The asymptotic confidence intervals for the parameters are also obtained based on asymptotic variance–covariance matrix. Finally, a real data set is analysed to show the potential of the new proposed distribution.     

Empirical approximations for Hoeffding’s test of bivariate independence using two Weibull extensions

The sampling distributions are generally unavailable in exact form and are approximated either in terms of the asymptotic distributions, or their correction using expansions such as Edgeworth, Laguerre or Cornish–Fisher; or by using transformations analogous to that of Wilson and Hilferty. However, when theoretical routes are intractable, in this electronic age, the sampling distributions can be reasonably approximated using empirical methods. The point is illustrated using the null distribution of Hoeffding’s test of bivariate independence which is important because of its consistency against all dependence alternatives. For constructing the approximations we employ two Weibull extensions, the generalized Weibull and the exponentiated Weibull families, which contain a rich variety of density shapes and tail lengths, and have their distribution functions and quantile functions available in closed form, making them convenient for obtaining the necessary percentiles and p-values. Both approximations are seen to be excellent in terms of accuracy, but that based on the generalized Weibull is more portable.     

The beta-Weibull geometric distribution

We propose a new distribution, the so-called beta-Weibull geometric distribution, whose failure rate function can be decreasing, increasing or an upside-down bathtub. This distribution contains special sub-models the exponential geometric [K. Adamidis and S. Loukas, A lifetime distribution with decreasing failure rate, Statist. Probab. Lett. 39 (1998), pp. 35–42], beta exponential [S. Nadarajah and S. Kotz, The exponentiated type distributions, Acta Appl. Math. 92 (2006), pp. 97–111; The beta exponential distribution, Reliab. Eng. Syst. Saf. 91 (2006), pp. 689–697], Weibull geometric [W. Barreto-Souza, A.L. de Morais, and G.M. Cordeiro, The Weibull-geometric distribution, J. Stat. Comput. Simul. 81 (2011), pp. 645–657], generalized exponential geometric [R.B. Silva, W. Barreto-Souza, and G.M. Cordeiro, A new distribution with decreasing, increasing and upside-down bathtub failure rate, Comput. Statist. Data Anal. 54 (2010), pp. 935–944; G.O. Silva, E.M.M. Ortega, and G.M. Cordeiro, The beta modified Weibull distribution, Lifetime Data Anal. 16 (2010), pp. 409–430] and beta Weibull [S. Nadarajah, G.M. Cordeiro, and E.M.M. Ortega, General results for the Kumaraswamy-G distribution, J. Stat. Comput. Simul. (2011). DOI: 10.1080/00949655.2011.562504] distributions, among others. The density function can be expressed as a mixture of Weibull density functions. We derive expansions for the moments, generating function, mean deviations and Rénvy entropy. The parameters of the proposed model are estimated by maximum likelihood. The model fitting using envelops was conducted. The proposed distribution gives a good fit to the ozone level data in New York.     

General results for the Kumaraswamy-G distribution

For any continuous baseline G distribution [G.M. Cordeiro and M. de Castro, A new family of generalized distributions, J. Statist. Comput. Simul. 81 (2011), pp. 883–898], proposed a new generalized distribution (denoted here with the prefix ‘Kw-G’ (Kumaraswamy-G)) with two extra positive parameters. They studied some of its mathematical properties and presented special sub-models. We derive a simple representation for the Kw-G density function as a linear combination of exponentiated-G distributions. Some new distributions are proposed as sub-models of this family, for example, the Kw-Chen [Z.A. Chen, A new two-parameter lifetime distribution with bathtub shape or increasing failure rate function, Statist. Probab. Lett. 49 (2000), pp. 155–161], Kw-XTG [M. Xie, Y. Tang, and T.N. Goh, A modified Weibull extension with bathtub failure rate function, Reliab. Eng. System Safety 76 (2002), pp. 279–285] and Kw-Flexible Weibull [M. Bebbington, C.D. Lai, and R. Zitikis, A flexible Weibull extension, Reliab. Eng. System Safety 92 (2007), pp. 719–726]. New properties of the Kw-G distribution are derived which include asymptotes, shapes, moments, moment generating function, mean deviations, Bonferroni and Lorenz curves, reliability, Rényi entropy and Shannon entropy. New properties of the order statistics are investigated. We discuss the estimation of the parameters by maximum likelihood. We provide two applications to real data sets and discuss a bivariate extension of the Kw-G distribution.     

A modified Kolmogorov-Smirnov test for the inverse gaussian density with unknown parameters

The Kolmogorov-Smirnov (KS) test is an empirical distribution function (EDF) based goodness-of-fit test that requires the underlying hypothesized density to be continuous and completely specified. When the parameters are unknown and must be estimated from the data, standard tables of the KS test statistic are not valid. Approximate upper tail percentage points of the KS statistic for the inverse Gaussian (IG) distribution with unknown parameters are tabled in this paper.

A study of the power of the KS test for the IG distribution indicates that the test is able todiscriminate between the IG distribution and distributions such as the uniform and exponentialdistributions that are very different in shape, but is relatively unable to discriminate between the IG distribution and distributions that are similar in shape such as the lognormal and Weibull distributions. In modeling settings the former distinction is typically more important to make than the latter distinction.     

The skew generalized t distribution as the scale mixture of a skew exponential power distribution and its applications in robust estimation

In this paper, we consider the family of skew generalized t (SGT) distributions originally introduced by Theodossiou [P. Theodossiou, Financial data and the skewed generalized t distribution, Manage. Sci. Part 1 44 (12) ( 1998), pp. 1650–1661] as a skew extension of the generalized t (GT) distribution. The SGT distribution family warrants special attention, because it encompasses distributions having both heavy tails and skewness, and many of the widely used distributions such as Student's t, normal, Hansen's skew t, exponential power, and skew exponential power (SEP) distributions are included as limiting or special cases in the SGT family. We show that the SGT distribution can be obtained as the scale mixture of the SEP and generalized gamma distributions. We investigate several properties of the SGT distribution and consider the maximum likelihood estimation of the location, scale, and skewness parameters under the assumption that the shape parameters are known. We show that if the shape parameters are estimated along with the location, scale, and skewness parameters, the influence function for the maximum likelihood estimators becomes unbounded. We obtain the necessary conditions to ensure the uniqueness of the maximum likelihood estimators for the location, scale, and skewness parameters, with known shape parameters. We provide a simple iterative re-weighting algorithm to compute the maximum likelihood estimates for the location, scale, and skewness parameters and show that this simple algorithm can be identified as an EM-type algorithm. We finally present two applications of the SGT distributions in robust estimation.     

Mixture Models and the Krätzel Integral Transform

Continuous mixture Weibull models arise in many areas of sciences such as reliability studies, communications theory, etc. Due to its wide applicability, we introduce a class of continuous mixture Weibull models which is a combination of Weibull and generalized gamma distributions. Some characteristics of the distribution are obtained. It is seen that Krätzel integral enters into the model naturally, and then the model can be called as a Krätzel density. Applications of the density function related to fading channels and ultrasonic backscatter signals modeling are discussed. A real data analysis is given to illustrate the use of this distribution.     

General results for the beta Weibull distribution

In this paper, we study some mathematical properties of the beta Weibull (BW) distribution, which is a quite flexible model in analysing positive data. It contains the Weibull, exponentiated exponential, exponentiated Weibull and beta exponential distributions as special sub-models. We demonstrate that the BW density can be expressed as a mixture of Weibull densities. We provide their moments and two closed-form expressions for their moment-generating function. We examine the asymptotic distributions of the extreme values. Explicit expressions are derived for the mean deviations, Bonferroni and Lorenz curves, reliability and two entropies. The density of the BW-order statistics is a mixture of Weibull densities and two closed-form expressions are derived for their moments. The estimation of the parameters is approached by two methods: moments and maximum likelihood. We compare the performances of the estimates obtained from both the methods by simulation. The expected information matrix is derived. For the first time, we introduce a log-BW regression model to analyse censored data. The usefulness of the BW distribution is illustrated in the analysis of three real data sets.     

Graphical Tests for the Assumption of Gamma and Inverse Gaussian Frailty Distributions

The common choices of frailty distribution in lifetime data models include the Gamma and Inverse Gaussian distributions. We present diagnostic plots for these distributions when frailty operates in a proportional hazards framework. Firstly, we present plots based on the form of the unconditional survival function when the baseline hazard is assumed to be Weibull. Secondly, we base a plot on a closure property that applies for any baseline hazard, namely, that the frailty distribution among survivors at time t has the same form as the original distribution, with the same shape parameter but different scale parameter. We estimate the shape parameter at different values of t and examine whether it is constant, that is, whether plotted values form a straight line parallel to the time axis. We provide simulation results assuming Weibull baseline hazard and an example to illustrate the methods.     

On a Class of Distributions Defined by the Relationship Between Their Density and Distribution Functions

Knowledge concerning the family of univariate continuous distributions with density function f and distribution function F defined through the relation f(x) = F ^α(x)(1 ? F(x))^β, α, β ? , is reviewed and modestly extended. Symmetry, modality, tail behavior, order statistics, shape properties based on the mode, L-moments, and—for the first time—transformations between members of the family are the general properties considered. Fully tractable special cases include all the complementary beta distributions (including uniform, power law and cosine distributions), the logistic, exponential and Pareto distributions, the Student t distribution on 2 degrees of freedom and, newly, the distribution corresponding to α = β = 5/2. The logistic distribution is central to some of the developments of the article.     

Properties and Inference for a New Class of Skew-t Distributions

In this article we introduce a new generalization of skew-t distributions, which contains the standard skew-t distribution, as a special case. This new class of distributions is an adequate model for modeling some dataset rather than the standard skew-t distributions. This kind of distributions can be represented as a scale-shape mixture of the extended skew-normal distributions. The main properties of this family of distributions are studied and a recurrence relation for the cumulative distribution functions (cdf) of them is presented. We derive the distribution of the order statistics from the trivariate exchangeable t-distribution in terms of our distribution and then an exact expression for the cdf of order statistics is derived. Likelihood inference for this distribution is also examined. The method is illustrated with a numerical example via a simulation study.     

On the Exact and Near-Exact Distributions of the Product of Generalized Gamma Random Variables and the Generalized Variance

In this article, the authors first obtain the exact distribution of the logarithm of the product of independent generalized Gamma r.v.’s (random variables) in the form of a Generalized Integer Gamma distribution of infinite depth, where all the rate and shape parameters are well identified. Then, by a routine transformation, simple and manageable expressions for the exact distribution of the product of independent generalized Gamma r.v.’s are derived. The method used also enables us to obtain quite easily very accurate, manageable and simple near-exact distributions in the form of Generalized Near-Integer Gamma distributions. Numerical studies are carried out to assess the precision of different approximations to the exact distribution and they show the high accuracy of the approximations provided by the near-exact distributions. As particular cases of the exact distributions obtained we have the distribution of the product of independent Gamma, Weibull, Frechet, Maxwell-Boltzman, Half-Normal, Rayleigh, and Exponential distributions, as well as the exact distribution of the generalized variance, the exact distribution of discriminants or Vandermonde determinants and the exact distribution of any linear combination of generalized Gumbel distributions, as well as yet the distribution of the product of any power of the absolute value of independent Normal r.v.’s.     

Survival distributions arising from two families and generated by transformations

Two families of distributions are introduced and studied within the framework of parametric survival analysis. The families are derived from a general linear form by specifying a function of the survival function with certain restrictions. Distributions within each family are generated by transformations of the survival time variable subject to certain restrictions. Two specific transformations were selected and, thus, four distributions are identified for further study. The distributions have one scale and two shape parameters and include as special cases the exponential, Weibull, log-logistic and Gompertz distributions. One of the new distributions, the modified Weibull, is studied in some detail.

The distributions are developed with an emphasis on those features that data analysts find especially useful for survivorship studies, A wide variety of hazard shapes are available. The survival, density and hazard functions may be written in simple algebraic forms. Parameter estimation is demonstrated using the least squares and maximum likelihood methods. Graphical techniques to assess goodness of fit are demonstrated. The models may be extended to include concmitant information.     

Bivariate inverse Weibull distribution

ABSTRACT

Recently it is observed that the inverse Weibull (IW) distribution can be used quite effectively to analyse lifetime data in one dimension. The main aim of this paper is to define a bivariate inverse Weibull (BIW) distribution so that the marginals have IW distributions. It is observed that the joint probability density function and the joint cumulative distribution function can be expressed in compact forms. Several properties of this distribution such as marginals, conditional distributions and product moments have been discussed. We obtained the maximum likelihood estimates for the unknown parameters of this distribution and their approximate variance– covariance matrix. We perform some simulations to see the performances of the maximum likelihood estimators. One data set has been re-analysed and it is observed that the bivariate IW distribution provides a better fit than the bivariate exponential distribution.     

Reliability estimation in stress–strength models: an MCMC approach

In this paper, we consider the estimation of the stress–strength parameter R=P(Y<X) when X and Y are independent and both are modified Weibull distributions with the common two shape parameters but different scale parameters. The Markov Chain Monte Carlo sampling method is used for posterior inference of the reliability of the stress–strength model. The maximum-likelihood estimator of R and its asymptotic distribution are obtained. Based on the asymptotic distribution, the confidence interval of R can be obtained using the delta method. We also propose a bootstrap confidence interval of R. The Bayesian estimators with balanced loss function, using informative and non-informative priors, are derived. Different methods and the corresponding confidence intervals are compared using Monte Carlo simulations.     

Order statistics from trivariate normal and -distributions in terms of generalized skew-normal and skew- distributions

We consider here a generalization of the skew-normal distribution, GSN(λ₁,λ₂,ρ), defined through a standard bivariate normal distribution with correlation ρ, which is a special case of the unified multivariate skew-normal distribution studied recently by Arellano-Valle and Azzalini [2006. On the unification of families of skew-normal distributions. Scand. J. Statist. 33, 561–574]. We then present some simple and useful properties of this distribution and also derive its moment generating function in an explicit form. Next, we show that distributions of order statistics from the trivariate normal distribution are mixtures of these generalized skew-normal distributions; thence, using the established properties of the generalized skew-normal distribution, we derive the moment generating functions of order statistics, and also present expressions for means and variances of these order statistics.Next, we introduce a generalized skew-t_ν distribution, which is a special case of the unified multivariate skew-elliptical distribution presented by Arellano-Valle and Azzalini [2006. On the unification of families of skew-normal distributions. Scand. J. Statist. 33, 561–574] and is in fact a three-parameter generalization of Azzalini and Capitanio's [2003. Distributions generated by perturbation of symmetry with emphasis on a multivariate skew t distribution. J. Roy. Statist. Soc. Ser. B 65, 367–389] univariate skew-t_ν form. We then use the relationship between the generalized skew-normal and skew-t_ν distributions to discuss some properties of generalized skew-t_ν as well as distributions of order statistics from bivariate and trivariate t_ν distributions. We show that these distributions of order statistics are indeed mixtures of generalized skew-t_ν distributions, and then use this property to derive explicit expressions for means and variances of these order statistics.