Abstract
We propose a new lifetime model called the exponentiated power generalized Weibull (EPGW) distribution, which is obtained from the exponentiated family applied to the power generalized Weibull (PGW) distribution. It can also be derived from a power transform on an exponentiated Nadarajah-Haghighi random variable. Since several structural properties of the PGW distribution have not been studied, they can be obtained from those of the EPGW distribution. The model is very flexible for modeling all common types of hazard rate functions. It is a very competitive model to the well-known Weibull, exponentiated exponential and exponentiated Weibull distributions, among others. We also give a physical motivation for the new distribution if the power parameter is an integer. Some of its mathematical properties are investigated. We discuss estimation of the model parameters by maximum likelihood and provide two applications to real data. A simulation study is performed in order to examine the accuracy of the maximum likelihood estimators of the model parameters.
Key words
Exponential distribution; lifetime data; Nadarajah-Haghighi distribution; power generalized Weibull distribution; survival function
1-INTRODUCTION
There has been an increased interest in defining new continuous distributions by adding shape parameters to an existing baseline model. One of the most widely-accepted methods on this parameter induction is the exponentiated-G (exp-G) class. Let G(y) and g(y)be the baseline cumulative distribution function (cdf) and the probability density function (pdf) of a random variable W, respectively. We obtain the exp-G cdf by raising G(y) to a positive power shape parameter to the baseline model. Thus, a random variable Y has an exp-G distribution if its cdf is given by
for , where represents the additional parameter. The corresponding pdf is given by
Tahir and Nadarajah (2015)TAHIR MH AND NADARAJAH S. 2015. Parameter induction in continuous univariate distributions: Well-established G families. An Acad Bras Cienc 87: 539-568. traced this approach back to the first half of the nineteenth century and found twenty-eight different exp-G models published in recent literature. Most of these models are motivated by their usefulness in exploring tail properties and also for improving the goodness-of-fit in comparison with their baselines. Another current reason for introducing exp-G distributions is their applications in lifetime data analysis.
Thus, the classical lifetime distributions have received a great deal of attention as baselines on the exp-G class, among other generated families. Using the exponential lifetime model as baseline,Gupta et al. (1998)GUPTA R, GUPTA P AND GUPTA R. 1998. Modeling failure time data by Lehman alternatives. Commun Stat Theory Methods 27: 887-904. pioneered the exponentiated exponential (EE) distribution.
The exponentiated Weibull (EW) distribution was introduced by Mudholkar and Srivastava (1993).
Another model that has been considered for modeling lifetime data is the Nadarajah-Haghighi (NH) distribution pioneered byNadarajah and Haghighi (2011).NADARAJAH S AND HAGHIGHI F. 2011. An extension of the exponential distribution. Statistics 45: 543-558.. The NH model is a generalization of the exponential distribution with cdf given by (for z > 0)
where λ and α are the scale and shape parameters, respectively. If Z has the cdf(1), we write NH(). The pdf of Z is given by
The motivations for studying the NH model are: the relationship between the pdf(2) and its hazard rate function (hrf), the ability (or inability) to model data with mode fixed at zero and the fact that it can be interpreted as a truncated Weibull distribution. Further details and general properties can be found inNadarajah and Haghighi (2011)NADARAJAH S AND HAGHIGHI F. 2011. An extension of the exponential distribution. Statistics 45: 543-558.. The exponentiated Nadarajah-Haghighi (ENH) model was proposed byLemonte (2013)LEMONTE AJ. 2013. A new exponential-type distribution with constant, decreasing, increasing, upside-down bathtub and bathtubshaped failure rate function. Comput Stat Data An 62: 149-170..
The exponential, NH and Weibull distributions are all special cases of the power generalized Weibull (PGW) distribution proposed byBagdonavicius and Nikulin (2002)BAGDONAVICIUS V AND NIKULIN M. 2002. Accelerated life models: modeling and statistical analysis. Boca Raton: Chapman and Hall/CRC, 360 p. in the context of accelerated failure time models. The original PGW cdf is given by (for )
where and . By setting and , the cdf, pdf and hrf of this distribution reduce to
and
respectively. Dimitrakopoulou et al. (2007)DIMITRAKOPOULOU T, ADAMIDIS K AND LOUKAS S. 2007. A lifetime distribution with an upside-down bathtub-shaped hazard function. IEEE T Reliab 121: 308-311. presented a lifetime model that has the same formulation as that one in(3) with motivation in competing risk scenario. Lai (2013)LAI CD. 2013. Constructions and applications of lifetime distributions. Appl Stoch Model Bus 29: 127-140 described the PGW among the Weibull generalizations that are often required to prescribe the nonmonotonic nature of the empirical hazard rates.
Nikulin and Haghighi (2006)NIKULIN M AND HAGHIGHI F. 2006. A chi-squared test for the generalized power Weibull family for the head-and-neck cancer censored data. J Math Sci 133: 1333-1341. introduced a chi-square statistic for testing the validity of the PGW distribution and presented an application to censored survival times of cancer patients.Nikulin and Haghighi (2009)NIKULIN M AND HAGHIGHI F. 2009. On the power generalized Weibull family. Metron 67: 75-86. presented shape analysis for the PGW pdf and hrf. They also obtained a series representation for the sth moment of this distribution, but only for integer values of . They do not provide a general expression for the PGW ordinary moments. We also note that there is a lack of other structural properties of the PGW distribution like incomplete moments, skewness, mean deviations, Bonferroni and Lorenz curves and Rényi entropy.
In this paper, we use the concept of exponentiated distributions for introducing a new four-parameter Weibull-type family, so-called the exponentiated power generalized Weibull(EPGW) distribution. The proposed distribution is obtained by considering the PGW model as baseline in the exp-G family. Thus, the EPGW cdf and pdf are given by (for t > 0)
and
respectively. Here, λ is the scale parameter and γ, α and β are shape parameters. Henceforth, we denote by T a random variable having pdf(6), say EPGW(). Identifiability is a property which a model must satisfy for precise inference to be possible, which refers to whether the unknown parameters in the model can be uniquely estimated. Equation(5) is clearly identifiable.
The hrf of T is given by
By inverting (5), we obtain an explicit expression for the quantile function (qf) of the EPGW distribution, say Q(u) given by
Its median M follows by setting in(8). The simulation of the EPGW random variable is straightforward. If U(0,1), then the random variable T=Q(U) follows the EPGW distribution given by (6).
Some motivations for introducing the EPGW distribution are:
-
The new distribution is quite flexible because it contains several well-known lifetime distributions as special models, see Table I. This feature is also suitable for testing the goodness-of-fit of these distributions.
-
The current distribution can also be derived from a power transform on an ENH random variable. If ENH(), the cdf of Y is given by
-
where α > 0, β > 0 and λ > 0. Consider the transformation , where γ > 0. Thus, the cdf of T has the form given by (5). A similar approach was addressed byGomes et al. (2008).GOMES O, COMBES C AND DUSSAUCHOY A. 2008. Parameter estimation of the generalized gamma distribution. Math Comput Simul 79: 963-995.. They proposed a new method of estimation for the generalized gamma distribution through the power transformation , where X is a generalized gamma random variable and W has the gamma distribution.
-
Once several structural properties of the PGW distribution have not been studied, they shall be obtained from those of the EPGW distribution.
-
By pioneering a PGW generalization of the exp-G family, it is also possible to obtain several properties of other generated families based on linear combinations from those of the EPGW distribution. For example, for the beta-G famil (Eugene et al. 2002EUGENE N, LEE C AND FAMOYE F. 2002. Beta-normal distribution and its applications. Commun Stat Theory Methods 31: 497-512., the density function can be expressed as a linear combination of exp-G pdfs for any baseline G. Similar results can also be demonstrated for the Kumaraswamy-G introduced byCordeiro and Castro (2011)CORDEIRO GM AND CASTRO MD. 2011. A new family of generalized distributions. J Stat Comput Sim 81: 883-89, among several others generated families of distributions.
-
Let β > 0 be an integer. Thus, F(t) given in(5) represents the cdf of the maximum value on a β-variate random sample from the PGW distribution, say: . In other words, the EPGW distribution can be used to model the maximum lifetime of a random sample from the PGW distribution with size β. Further, as part of the exp-G family, the EPGW distribution has the following physical interpretation. Consider a parallel system consisting of β = n components, which means that the system works if at least one of the n-components works. If the lifetime distributions of the components are independent and identically distributed PGW random variables, then the lifetime distribution of the system becomes the EPGW cdf with power parameter β = n.
-
The EPGW may provide consistently ‘better fits’ than other Weibull generalizations including its special models. This fact is shown by fitting the proposed distribution to two data sets in Section 13. The applications illustrate that the EPGW distribution can also be very competitive to other widely known lifetime models.
The paper is outlined as follows. Some mathematical properties of the new distribution are provided in Sections 3-10 They include ordinary and incomplete moments, mean deviations about the mean and the median, Bonferroni and Lorenz curves, Rényi entropy, reliability and order statistics. In Section 11, we present the maximum likelihood method to estimate the model parameters. In Section 12, a simulation study evaluates the performance of the maximum likelihood estimators (MLEs). Applications to two real data sets are presented in Section 13. Section 14 offers some concluding remarks.
2-DENSITY AND HAZARD SHAPES
Note that the pdf (6) can be expressed in terms of the cdf and pdf given in (3) and (4), respectively, in the form . Thus, the multiplicative factor is greater (smaller) than one for and for larger values of t, and the opposite occurs for smaller values of t. The inclusion of the extra shape parameter β provides greater flexibility in terms of skewness and kurtosis of the new distribution. The pdf (6) can take various forms depending on the values of the shape parameters α,β and γ. It is easy to verify that
and
Seting , we can rewrite the EPGW pdf as
Differentiating twice with respect to z, we obtain
Note that implies that . Thus, we can verify that for , , and , . This implies that the EPGW pdf is log-convex. Further, for , , and , , which implies that the EPGW pdf is log-concave. Figure 1 displays plots of the pdf(6) for some parameter values. It illustrates the flexibility of the EPGW density, which allows modeling skewed and asymmetrical data.
Analogously, the EPGW hrf can be rewritten as
The critical point are obtained from
For , and the hrf is constant. For , and , and the hrf is decreasing. There may be more than one root to this equation.
Figure 2 provides plots of the hrf(7) for some parameter values. Figure 2 reveals that the EPGW distribution can have decreasing, increasing, upside-down bathtub and bathtub-shaped hazard functions. This feature makes the new distribution very attractive to model lifetime data. For example, according toNadarajah et al. (2011)NADARAJAH S, BAKOUCH H AND TAHMASBI R. 2011. A generalized Lindley distribution. Sankhya Ser B 73: 331-359. most empirical life systems have bathtub shapes for their hrfs.
3-MOMENTS
The sth ordinary moment of T is obtained as, with f(t) from(6). For illustrative purposes, we provide a small numerical study by computing the first six moments for some scenarios. Each one considers a different parametrization for γ and β, with fixed α = 1.5 and λ = 1.
These results are presented in Table II. All computations are obtained using R software, which have numerical integration routines with great precision. Based on these values, we can note that, for fixed γ, the additional parameter β has large impact on the moments of T. Note that the moments increase as β increases. For β fixed and ≤ 1, the moments decreases when γ increases.
The sth moment of Tcan also be determined from equation(8). After some algebra, we can write
where is an integral to be evaluated numerically.
Using the binomial expansion since the inverse of the denominator of (6) can be expressed as
Further, we can rewrite μ'u as
We consider the integral
Setting , we have
Hence, after some algebra, we obtain
The most general case of the binomial theorem is the power series identity
where is a binomial coefficient and ν is a real number. This power series converges for an integer, or . This general form is given byGraham (1994).GRAHAM RL. 1994. Concrete mathematics: a foundation for computer science, 2nd ed., Boston: Addison-Wesley, 672 p.. By using (11) in equation(10), since , it follows from (9) that
where denotes the complementary incomplete gamma function, which is defined for all real numbers except the negative integers. Equation(12) is the main result of this section.
4-SKEWNESS
The central moments and cumulants of can be expressed recursively from equation (12) as
respectively, where . Thus, , etc. The skewness and kurtosis can be determined from the third and fourth bookized cumulants.
The MacGillivray (1986)MACGILLIVRAY H. 1986. Skewness and asymmetry: measures and orderings. Ann Stat 14: 994-1011. skewness function of is given by
where , is the qf defined in (8),
and
It is based on quantiles and can illustrate the effects of the shape parameters α, β and γ on the skewness of T. Plots of for some parameter values are displayed in Figure 3. These plots reveal that when the parameters and increase, the function converges to zero. The closer is to the horizontal line , the density becomes more symmetrical. The quantity does not depend on the parameter since it is a scale parameter.
5-INCOMPLETE MOMENTS
The th incomplete moment of , say , follows as
An alternative expression for takes the form
6-MEAN DEVIATIONS
The mean deviations about the mean and about the median of can be expressed as
respectively, where , is the median, is easily determined from(5) and is the first incomplete moment. Hence, we can write
Alternatively, we can determine as
7-BONFERRONI AND LORENZ CURVES
Applications of the previous results to the Bonferroni and Lorenz curves are important in several fields such as economics, demography, insurance and medicine. They are defined, for a given probability , by and , respectively, where follows from ((8)). The Gini concentration () is defined as the area between the curve and the straight line. Hence,
An alternative expression is , where . The quantity is given by
This integral can be easily evaluated numerically in softwares such as R and Ox, among others. An alternative expression for takes the form
For , we can prove that this expression reduces to that one obtained byLemonte (2013)LEMONTE AJ. 2013. A new exponential-type distribution with constant, decreasing, increasing, upside-down bathtub and bathtubshaped failure rate function. Comput Stat Data An 62: 149-170..
8-ENTROPY
The entropy of a random variable is a measure of variation of the uncertainty and has been used in many fields. Several measures of entropy have been studied in the literature. However, we consider the most popular entropy measure: the Rényi entropy of a random variable with pdf defined by
for and The Rényi entropy of can be expressed as
where . The above integral can be evaluated numerically. By expanding the inverse of the denominator using the binomial expansion, we obtain
Again, by using the binomial expansion, can be expressed as
For , the last expression agrees with the result byLemonte (2013)LEMONTE AJ. 2013. A new exponential-type distribution with constant, decreasing, increasing, upside-down bathtub and bathtubshaped failure rate function. Comput Stat Data An 62: 149-170..
9-STRESS-STRENGTH RELIABILITY
The stress-strength reliability, defined as , is a measure that describes the life of a component with a random strength , which is subjected to a random stress . The failure occurs if the stress applied to the component exceeds the strength, i.e. , otherwise it will function satisfactorily. Clearly, this measure is very useful in engineering context such as deterioration of rocket motors and the aging of concrete pressure vessels.
Let and be two independent random variables with EPGW and EPGW distributions, respectively. We shall obtain the stress–strength parameter in the form
Thus, by taking , the above integral can be reduced to
10-ORDER STATISTICS
Let be a random sample from the EPGW distribution. Let denote the th order statistic. The probability density function of is
By inserting(5) and(6) in(13) and after some algebra, we obtain
Thus, we can write
where
is the EPGW density function with scale parameter and shape parameters , and . Equation(14) is the main result of this section. Based on this, we can obtain some structural properties of using similar procedures as presented in the previous sections.
11-MAXIMUM LIKELIHOOD ESTIMATION
This section addresses the estimation of the unknown parameters of the EPGW distribution by the method of maximum likelihood. Let be a random sample of size from the EPGW distribution. Let be the parameter vector of interest. The -likelihood function for based on this sample is
The components of the score vector are given by
and
Setting the above equations to zero, , and solving them simultaneously yields the MLEs of the four parameters. These equations cannot be solved analytically. We have to use iterative techniques such as the quasi-Newton BFGS and Newton-Raphson algorithms. The initial values for the parameters are important but are not hard to obtain from fitting special EPGW sub-models.
Note that, for fixed , and , the MLE of is given by
Thus, it is easily observed that fixed on ,
-
when and/or
-
when and/or
-
when and , for some
-
when and , .
This behavior anticipates that estimates for smaller and/or may require improved estimation procedures.
By replacing by in equation(15) and letting , the profile log-likelihood function for can be expressed as
We assume that the book regularity conditions for hold.
They are not restrictive and hold for the models cited in this paper. The score vector corresponding to (16), , has the components
and
Solving the equations in simultaneously yields the MLEs of , and . The MLE of is just . The maximization of the profile log-likelihood might be simpler since it involves only three parameters. Lemonte (2013)LEMONTE AJ. 2013. A new exponential-type distribution with constant, decreasing, increasing, upside-down bathtub and bathtubshaped failure rate function. Comput Stat Data An 62: 149-170. noted a similar result for the ENH model but mentioned that some of the properties that hold for a genuine likelihood do not hold for its profile version.
For interval estimation of the components of , we can adopt the observed information matrix , whose elements can be obtained from the authors upon request. The multivariate normal distribution can be used to construct approximate confidence intervals for the model parameters.
12-SIMULATION STUDY
Here, a Monte Carlo simulation experiment is performed in order to examine the accuracy of the MLEs of the model parameters. The simulations are carried out by generating observations from the EPGW distribution using the inverse transformation method for different parameter combinations. The number of observations is set at n = 100, 300 and 500 > and the number of replications at 10,000 For maximizing the log-likelihood function, we use the Optim function with analytical derivatives in R. From the results of the simulations given in Table Table III, we can verify that the root mean squared errors (RMSEs) of the MLEs of , , and decay toward zero when the sample size n increases, as expected. The mean estimates of the parameters tend to be closer to the true parameter values when n increases.
13-APPLICATIONS
In this section, we present two applications to illustrate the flexibility of the EPGW distribution. They indicate the potentiality of the new distribution for modeling positive data. The first data set represents the remission times (in months) of 128 patients with bladder cancer (Lee and WangLEE E AND WANG J. 2003. Statistical Methods for Survival Data Analysis, 3rd ed., New York: Wiley, 534 p.2003). The second one corresponds to the 101 observations representing the stress-rupture life of kevlar 49/epoxy strands that are subjected to constant sustained pressure at the 90% stress level until all had failed. Then, we obtain a complete data set with exact failure times. This data set was studied byAndrews and Herzberg (1985).ANDREWS DF AND HERZBERG AM. 1985. Data: A Collection of Problems from Many Fields for the Student and Research Worker. New York: Springer Series in Statistics, 442 p.. Table IV gives a descriptive summary of the samples. Note that both data sets present positive skewness and that the remission times show higher variance.
We fit the EPGW distribution(6) to these data sets and also present a comparative study with the fits of some embedded and not embedded models. One of these models is the Kumaraswamy Weibull (Kw-W) distribution, whose pdf is given by
where , , and . Another model is the beta Weibull (BW) distribution, whose pdf is given by
where , , and . The Beta-Fréchet (BFr) distribution, whose pdf is given by
where , , and . We also consider the Marshall-Olkin Nadarajah-Haghighi (MONH) model, whose pdf is given by
where and . The EW distribution, whose pdf is given by
where and are shape parameters and is a scale parameter. This distribution is quite flexible because its hrf presents the classic five forms (constant, decreasing, increasing, upside-down bathtub and bathtub-shaped). The Weibull model is a special case of the EW model when .
The ENH distribution can also have the same shapes for the hrf and therefore can be an interesting alternative to the EW distribution in modeling positive data. The ENH density is given by (6) when . Further, for , we have as a special model the NH distribution given by(2). We also consider the PGW model, whose pdf is given in (4), which arises from the EPGW model when .
Xie et al. (2002)XIE M, TANG Y AND GOH T. 2002. A modified Weibull extension with bathtub-shaped failure rate function. Reliab Eng and Syst Safe 76: 279-285. proposed a modified Weibull (MW) density given by
where , and . For , it becomes the Chen distribution (ChenCHEN Z. 2000. A new two-parameter lifetime distribution with bathtub shape or increasing failure rate function. Stat Probabil Lett 49: 155-161. 2000). The MW and Chen distributions can have increasing or bathtub-shaped failure rate. An extension of the Weibull model proposed byBebbington et al. (2007)BEBBINGTON M, LAI CD AND ZITIKIS R. 2007. A flexible Weibull extension. Reliab Eng and Syst Safe 92: 719-726. has pdf given by
where and . We shall use the same terminology byLemonte (2013)LEMONTE AJ. 2013. A new exponential-type distribution with constant, decreasing, increasing, upside-down bathtub and bathtubshaped failure rate function. Comput Stat Data An 62: 149-170. for this distribution, i.e., denote the flexible Weibull (FW) density. The FW model can have increasing or modified bathtub-shaped failure rate.
We use the simulated-annealing method for maximizing the log-likelihood function of the models in the two applications. The MLEs and goodness-of-fit statistics are evaluated using the AdequacyModel script in R software. Tables Table V and Table VI list the MLEs and the corresponding book errors (SEs) in parentheses of the unknown parameters for the fitted models to remission times data (first data set) and stress-rupture failure times (second data set), respectively.
The MLEs of the model parameters for the remission times data and the corresponding SEs in parentheses.
The MLEs of the model parameters for the stress-rupture data and the corresponding SEs in parentheses.
In applications there is qualitative information about the failure rate shape, which can help for selecting some models. Thus, a device called the total time on test (TTT) plot is useful. The TTT plot is obtained by plotting
against , where and are the order statistics of the sample.
The figures in Tables V and VI indicate that the MLEs of the EPGW model are precise for both data sets. Figures 4 and 5 provide the TTT plots and plots of the hrf for the EPGW fitted model for the remission times and stress-rupture times data sets, respectively. They reveal that the EPGW hrf has decreasing and decreasing-increasing-decreasing shapes, respectively. This fact is in agreement with the TTT plot based on each data set.
Chen and Balakrishnan (1995)CHEN G AND BALAKRISHNAN N. 1995. A general purpose approximate goodness-of-fit test. J Qual Technol 27: 154-161. constructed the corrected Cramér-von Mises and Anderson-Darling statistics. We adopt these statistics, where we have a random sample with empirical distribution function , and require to test if the sample comes from a special distribution. The Cramér-von Mises () and Anderson-Darling () statistics are given by
and
respectively, where is the empirical distribution function, is the postulated distribution function evaluated at the MLE of . Note that the statistics and are given by the differences of and . Thus, the lower are these statistics, we have more evidence that generates the sample. The details to evaluate the statistics and are given byChen and Balakrishnan (1995).CHEN G AND BALAKRISHNAN N. 1995. A general purpose approximate goodness-of-fit test. J Qual Technol 27: 154-161..
The , and Kolmogorov-Smirnov (KS) statistics for these models are given in Tables Table VII and Table VIII for both data sets. We emphasize that the EPGW model fits the remission times and stress-rupture failure data better than the other models according to these statistics. They indicate that the EPGW distribution yields the best fits in both applications.
More information is provided by the histogram of the data and some fitted density functions for both data sets given in Figure 6. Clearly, in both applications, the new distribution provides a closer fit to the histogram than the other competitive models. The fitted cdfs of these models are also displayed in Figure 7. Finally, we can conclude in the two situations that the EPGW distribution is quite competitive to other well-known and widely used distributions such as the Kw-W, EW and Weibull models.
Histogram and estimated densities of the (a) EPGW, Kw-W and PGW models for the remission times data; (b) EPGW, PGW and MW models for the stress-rupture data.
Estimated and empirical cdfs for (a) EPGW, Kw-W and PGW models for the remission times data; (b) EPGW, PGW and MW models for the stress-rupture data.
CONCLUSIONS
In this paper, we introduce the exponentiated power generalized Weibull (EPGW) model to extend the Weibull distribution. It has a power parameter and its hazard rate function allows constant, decreasing, increasing, upside-down bathtub or bathtub-shaped forms. The proposed distribution contains as special models several well-known lifetime distributions. It can also be derived from a power transform on an exponentiated Nadarajah-Haghighi random variable. Several structural properties of the power generalized Weibull (PGW) distribution have not been studied. However, they can be determined from those of the EPGW distribution. It can also be useful to obtain the properties for other generated families under the PGW baseline. We give a physical motivation for introducing the new distribution if the power parameter is an integer. We obtain some mathematical properties of the EPGW distribution, estimate the model parameters by maximum likelihood and prove empirically its flexibility in two applications to real data. In fact, the new distribution yields a good adjustment in both applications. We note that the EPGW distribution is quite competitive with other lifetime models and can be used effectively to provide better fits than other usual lifetime distributions.
REFERENCES
- ANDREWS DF AND HERZBERG AM. 1985. Data: A Collection of Problems from Many Fields for the Student and Research Worker. New York: Springer Series in Statistics, 442 p.
- BAGDONAVICIUS V AND NIKULIN M. 2002. Accelerated life models: modeling and statistical analysis. Boca Raton: Chapman and Hall/CRC, 360 p.
- BEBBINGTON M, LAI CD AND ZITIKIS R. 2007. A flexible Weibull extension. Reliab Eng and Syst Safe 92: 719-726.
- CHEN G AND BALAKRISHNAN N. 1995. A general purpose approximate goodness-of-fit test. J Qual Technol 27: 154-161.
- CHEN Z. 2000. A new two-parameter lifetime distribution with bathtub shape or increasing failure rate function. Stat Probabil Lett 49: 155-161.
- CORDEIRO GM AND CASTRO MD. 2011. A new family of generalized distributions. J Stat Comput Sim 81: 883-89
- DIMITRAKOPOULOU T, ADAMIDIS K AND LOUKAS S. 2007. A lifetime distribution with an upside-down bathtub-shaped hazard function. IEEE T Reliab 121: 308-311.
- EUGENE N, LEE C AND FAMOYE F. 2002. Beta-normal distribution and its applications. Commun Stat Theory Methods 31: 497-512.
- GOMES O, COMBES C AND DUSSAUCHOY A. 2008. Parameter estimation of the generalized gamma distribution. Math Comput Simul 79: 963-995.
- GRAHAM RL. 1994. Concrete mathematics: a foundation for computer science, 2nd ed., Boston: Addison-Wesley, 672 p.
- GUPTA R, GUPTA P AND GUPTA R. 1998. Modeling failure time data by Lehman alternatives. Commun Stat Theory Methods 27: 887-904.
- LAI CD. 2013. Constructions and applications of lifetime distributions. Appl Stoch Model Bus 29: 127-140
- LEE E AND WANG J. 2003. Statistical Methods for Survival Data Analysis, 3rd ed., New York: Wiley, 534 p.
- LEMONTE AJ. 2013. A new exponential-type distribution with constant, decreasing, increasing, upside-down bathtub and bathtubshaped failure rate function. Comput Stat Data An 62: 149-170.
- MACGILLIVRAY H. 1986. Skewness and asymmetry: measures and orderings. Ann Stat 14: 994-1011.
- NADARAJAH S, BAKOUCH H AND TAHMASBI R. 2011. A generalized Lindley distribution. Sankhya Ser B 73: 331-359.
- NADARAJAH S AND HAGHIGHI F. 2011. An extension of the exponential distribution. Statistics 45: 543-558.
- NIKULIN M AND HAGHIGHI F. 2006. A chi-squared test for the generalized power Weibull family for the head-and-neck cancer censored data. J Math Sci 133: 1333-1341.
- NIKULIN M AND HAGHIGHI F. 2009. On the power generalized Weibull family. Metron 67: 75-86.
- TAHIR MH AND NADARAJAH S. 2015. Parameter induction in continuous univariate distributions: Well-established G families. An Acad Bras Cienc 87: 539-568.
- XIE M, TANG Y AND GOH T. 2002. A modified Weibull extension with bathtub-shaped failure rate function. Reliab Eng and Syst Safe 76: 279-285.
Publication Dates
-
Publication in this collection
Jul-Sep 2018
History
-
Received
20 June 2017 -
Accepted
01 Nov 2017