A new distribution with decreasing, increasing and upside-down bathtub failure rate
Introduction
When modeling monotone hazard rates, the Weibull distribution may be an initial choice because of its negatively and positively skewed density shapes. However, the Weibull distribution does not provide a reasonable parametric fit for modeling phenomenon with non-monotone failure rates, such as the upside-down bathtub failure rates, which are common in reliability and biological studies. For example, such failure rates curves can be observed in the course of a disease whose mortality reaches a peak after some finite period and then declines gradually. The lifetime models that present upside-down bathtub shaped failure rates are very useful in survival analysis. Recently, probability distributions with such properties have been investigated in Shen et al. (2009), Navarro et al. (2009), Carrasco et al. (2009), Barreto-Souza et al. (2009), Ghitany et al. (2009) and Bebbington et al. (2008).
Adamidis and Loukas (1998) introduced the exponential geometric (EG) distribution whose cumulative distribution function (cdf) is given by for and . In the same way, Kus (2007) introduced the exponential-Poisson distribution. Gupta and Kundu, 1999, Gupta and Kundu, 2001a, Gupta and Kundu, 2001b, Gupta and Kundu, 2007 proposed another lifetime distribution which referred to the generalized exponential (GE) (also called the exponentiated exponential) distribution, and investigated several of its mathematical properties. See, also, Raqab and Ahsanullah (2001). The GE distribution can be defined by elevating the CDF of the exponential distribution with a power .
A further exponentiated type distribution has been introduced and studied in the literature. The exponentiated Weibull (EW) distribution was proposed by Mudholkar and Srivastava (1993) to extend the GE distribution. This distribution was also studied by Mudholkar et al. (1995), Mudholkar and Hutson (1996) and Nassar and Eissa (2003). Nadarajah and Kotz (2006a) introduced four more exponentiated type distributions: the exponentiated gamma, exponentiated Weibull, exponentiated Gumbel and exponentiated Fréchet distributions by generalizing the gamma, Weibull, Gumbel and Fréchet distributions in the same way that the GE distribution extends the exponential distribution. In a recent paper, Barreto-Souza and Cribari-Neto (2009) introduced the generalized exponential-Poisson distribution which extends the exponential-Poisson distribution in the same way that the GE distribution extends the exponential distribution.
Following the same idea of the GE distribution, the CDF of the generalized exponential geometric (GEG) distribution is defined from Eq. (1) where , and .
We now give a characterization for the GEG distribution. We consider integer and a random sample from the EG distribution. For an interpretation of the EG distribution, see Adamidis and Loukas (1998). The definition leads to the GEG distribution which can be used to model the maximum lifetime of a random sample from the EG distribution. Further, it also has a desirable physical interpretation. If there are components in a parallel system and the lifetimes of the components are independent and identically distributed having the GEG distribution, then the system lifetime also has the GEG distribution.
The rest of the paper is organized as follows. In Section 2, we present the probability density function (PDF) and failure rate function and provide plots of such functions for selected parameter values. In Section 3, we obtain the moment generating and characteristic functions. We also give the moments of the order statistics. The Rényi entropy is derived in Section 4. Maximum likelihood estimation of the parameters and the expected information matrix are discussed in Section 5. Section 6 deals with the estimation of the stress–strength parameter. An application of the GEG model to real data is illustrated in Section 7. Concluding remarks are given in Section 8.
Section snippets
Density function, failure rate and order statistics
The density of the GEG distribution corresponding to the CDF (2) is where , and . The parameters and are scale and shape parameters, respectively. A random variable with density function (3) is denoted by .
Clearly, the EG distribution is obtained from (3) when . The GE distribution comes as the limiting distribution (the limit is defined in terms of the convergence in distribution) of the GEG distribution when
Moment generating function
If a random variable has the distribution, Nadarajah and Kotz (2006b) showed that its moment generating function (mgf) and characteristic function (cf) are and respectively, where is a complex number.
Combining the expansion (5) and these results, we obtain the MGF and CF of the GEG distribution and respectively.
Rényi entropy
The entropy of a random variable is a measure of uncertainty variation. The Rényi entropy is defined by , where and . Let be the GEG density (3). We have
If we expand as in Eq. (4), we obtain where Hence, can be written as an infinite linear combination of BE
Estimation and Fisher information matrix
We examine maximum likelihood estimation and inference for the GEG distribution. Let be a random sample from and let be the vector of the model parameters. The log-likelihood function for reduces to
The score vector , has components easily calculated by differentiating (10). We obtain
Estimation of the stress–strength parameter
In the context of reliability, the stress–strength model describes the life of a component which has a random strength subjected to a random stress . The component fails at the instant that the stress applied to it exceeds the strength, and the component will function satisfactorily whenever . Hence, is a measure of component reliability. It has many applications, especially in engineering concepts such as structures, deterioration of rocket motors, static fatigue of ceramic
Application
We present an application of the GEG distribution to a real data set. The data set is given by Birnbaum and Saunders (1969) on the fatigue life of 6061-T6 aluminum coupons cut parallel to the direction of rolling and oscillated at 18 cycles per second. The data set consists of 101 observations with maximum stress per cycle 31,000 psi. The data are: 70, 90, 96, 97, 99, 100, 103, 104, 104, 105, 107, 108, 108, 108, 109, 109, 112, 112, 113, 114, 114, 114, 116, 119, 120, 120, 120, 121, 121, 123,
Concluding remarks
The two-parameter exponential geometric (EG) distribution introduced by Adamidis and Loukas (1998) is generalized by introducing an extra parameter, thus creating the generalized EG model (GEG) with a broader class of hazard functions. This is achieved by (the well known technique) raising the cumulative distribution function of the EG to the power of the extra parameter. A detailed study on the probabilistic characteristics of the distribution is presented. The new model includes the
Acknowledgements
We gratefully acknowledge grants from CAPES and CNPq (Brazil). The authors are also grateful to an associate editor and two referees for helpful comments and suggestions.
References (26)
- et al.
On a generalization of the exponential-geometric distribution
Statistics & Probability Letters
(2005) - et al.
A lifetime distribution with decreasing failure rate
Statistics & Probability Letters
(1998) - et al.
Estimating the turning point of a bathtub-shaped failure distribution
Journal of Statistical Planning and Inference
(2008) - et al.
Generalized exponential distribution: Existing results and some recent developments
Journal of Statistical Planning and Inference
(2007) A new lifetime distribution
Computational Statistics & Data Analysis
(2007)- et al.
The beta exponential distribution
Reliability Engineering and System Safety
(2006) - et al.
A generalization of the exponential-Poisson distribution
Statistics and Probability Letters
(2009) - et al.
A beta generalized exponential distribution
Journal of Statistical Computational and Simulation
(2009) - et al.
Estimation for a family of life distributions with applications to fatigue
Journal of Applied Probability
(1969) - et al.
A generalized modified Weibull distribution for lifetime modeling
Computational Statistics & Data Analysis
(2009)
The estimation of reliability from stress strength relationships
Technometrics
The estimation of in gamma case
Communications in Statistics—Simulation and Computation
Bathtub and related failure rate characterizations
Journal of the American Statistical Association
Cited by (68)
Exponentiated Generalized General Class of Inverted Distributions: Estimation and Prediction
2023, G Families of Probability Distributions: Theory and PracticesZero Truncated Poisson - Pareto Distribution: Application and Estimation Methods
2023, WSEAS Transactions on MathematicsClassical and Bayesian estimations of improved Weibull–Weibull distribution for complete and censored failure times data
2022, Applied Stochastic Models in Business and IndustryA NEW FAMILY OF COMPOUND EXPONENTIATED LOGARITHMIC DISTRIBUTIONS WITH APPLICATIONS TO LIFETIME DATA
2022, Mathematica Slovaca