Gumbel-Exponentiated Weibull {Logistic} Lifetime Distribution and Its Applications ()
1. Introduction
Modeling and analyzing lifetime data have received tremendous attention in recent years due to the changing nature of lifetime data structures. Data, which arise from several real-life problems, are assumed to follow specific probability distributions defined by a probability density function (pdf) and the shape of the hazard rate function (hrf) which could be monotonically increasing or decreasing, increasing-decreasing-increasing, decreasing-increasing-decreasing, and bathtub shaped or upside-down bathtub shaped. Several lifetime distributions have been used in literature to model these data sets with varying characteristics, including the Weibull distribution, exponential, gamma, Rayleigh, Lindley, Pareto distributions and their generalizations. It is worthy of note that most generalizations have followed several known methods for generating new families of distribution. Some of the continuous univariate generators which have received great deal of attention in recent decade include: method based on differential equations introduced by [1], method of generation based on transformations by [2], method based on quantiles proposed by [3] and [4], methods based on generating skewed distributions developed by [5], method based on adding parameters and generators proposed by [6], method based on compounding univariate continuous distributions and family of discrete distributions suggested by [7], methods based on generators proposed by [8], composite method by [9], and the transformed-transformer (T-X) approached introduced by [10] which was later modified as T-R{Y} by [11] and [12].
[6] introduced an extension of the Weibull distribution called the exponentiated Weibull distribution by adding additional shape parameter to this distribution and since then many works have been carried out by researchers to make the exponentiated weibull distribution more flexible. These efforts include: the Exponentiated Weibull family by [13], extended exponentiated Weibull lifetime distribution developed by [14], exponentiated Weibull-Weibull distribution studied by [15], while [16] suggested a new generalization of the exponentiated modified Weibull known as beta-exponentiated modified Weibull. For other modifications of the Weibull distribution, see [17] and [18]. The exponentiated Weibull distribution and its extensions have a great deal of usage and current generalization will be able to handle more intricate lifetime modeling situations.
This paper therefore focuses on deriving a new exponentiated Weibull lifetime distribution that will be capable of modeling both unimodal and bimodal datasets. Gumbel-Weibull, Gumbel-Rayleigh and Gumbel-exponential distributions form the special cases of the proposed distribution. The proposed distribution is flexible and can play important role in reliability analysis, because the hazard function can assume several shapes. The kurtosis of the Gumbel-Exponentiated Weibull distribution will achieve higher flexibility when compared with the baseline distribution.
The remainder of this paper is organized as follows: the cdf and pdf of the new generalized exponentiated Weibull distribution is derived in Section 2. The statistical and reliability properties of the proposed model are presented in Section 3. Estimation of parameters using the maximum likelihood approach is discussed in Section 4. Simulation of result is presented in Section 5. Section 6 presents the application of real-life dataset and conclusion is Section 7.
2. The Gumbel-Exponentiated Weibull {Logistic} Distribution (GEWLD)
Let T be a random variable following Gumbel distribution with the cumulative distribution and density function given by
and
,
respectively and b is a scale parameter and s is a location parameter. Suppose R is an exponentiated Weibull distribution with cumulative distribution function (cdf) and probability density function (pdf) given by
pdf is given by,
where
is a scale parameter and
are the shape parameters (see [6] ).
Let the density function
and quantile function
of a Y random variable following a standard logistic distribution be given respectively by
We introduce the 5-parameter GEWLD in this section, by adopting the T-R{Y} family of distributions generation approach introduced by [11]. Here we define the cdf of the T-R{Y} as
(1)
The corresponding pdf associated with (1) is
(2)
where
where,
,
and
are the cdfs of Gumbel, exponentiated Weibull and standard Logistic random variables with their corresponding cdfs
,
and
and the quantile function of Y random variable
,
, respectively.
Suppose
are the parameters of the proposed model, where
,
,
,
,
,
. Applying (1) and (2) above, then the cdf of the Gumbel-exponentiated Weibull {Logistic} distribution (GEWLD) is obtained as
(3)
Let
be a random variable with cdf (3), then the pdf of X is
(4)
2.1. Shape of the Density
The study of the first and second derivative of the cdf (3) helps us to understand the main features of the density shape. The process to obtain the first derivative of
is provided below.
Let
. Then
,
and
. Rewrite (4) into a function in k. For
, we have,
Taking logarithm of both sides of the equation, we have
Taking derivatives, we obtain
The first derivative of
of the Gumbel-Exponentiated Weibull {Logistic} distribution after replacing k by
accordingly, is:
The
does not exist when
. If there are modes for GEWLD it will satisfy
.
Let,
, since it is a scale parameter and does not affect the shape of the density, then
(5)
Multiplying the first part of the Equation (5) by
, we obtain
(6)
So, the modes of
are the roots of Equation (6). Another possible root of the density is
and there may be other roots. If we have
as a root of (6) then it corresponds to the local maximum of
for all
and
for all
. It corresponds to local maximum of
for all
and
for
. The solution tallies to the point of inflexion if either
for all
or
for all
.
The plots in Figure 1 give (a) bimodal density shapes with fixed values of
and varying values of α, b and s, (b) left-skewed density shapes
Figure 1. Plots of the GEWLD pdf for
and some values of α, b and s.
with fixed values of
and varying values of β, b, and s, (c) reversed-J, unimodal right-skewed density shapes with fixed values of
and varying values of α, β and b.
The Survival function which is used in engineering to model time-to-event surviving beyond time x in reliability analysis is presented in Equation (6)
(7)
It is easy to verify that the hazard rate function (hrf) is of GEWLD is given by
(8)
From the graph of the hazard rate function (hrf) plots in Figure 2. below, the GEWLD displays decreasing, upside-down bathtub shape, bathtub shape and increasing hazard rate with time at different values of the parameters.
Figure 2. Hazard rate function plots of GEWLD for selected parameter value.
2.2. Submodels of the GEWL Distribution
The following results in the proposed distribution provide some immediate characterization of the Gumbel-Exponentiated Weibull {Logistic} distribution (GEWLD).
1) If the shape parameter
, the GEWLD reduces to the four parameter Gumbel-Weibull (GW) distribution introduced by [19] with cdf and pdf
with
and
, and α, b are the shape parameters, λ and s are the scale parameters.
2) If the parameter
, then the GEWLD reduces to the Gumbel exponential distribution with cdf given as
where
and
and λ, s are the scale parameters and b is the shape parameter.
3) If
, then the GEWLD reduces to Gumbel-Rayleigh distribution with cdf and pdf given respectively as
and
where
and
, b is the shape parameter, b and s are the scale parameters.
3. Some structural Properties
3.1. Expansions for the Cumulative and Density Function of the GEWLD
Very important mixture expression of the cumulative distribution function (cdf) and probability density function pdf of the GEWLD are presented here. By using some series expansions with obtain the cdf of GEWL distribution which is given by
(9)
where
we also obtain the expression for the pdf of GEWL distribution as a linear combination of Gumbel and exponentiated Weibull densities as
, (10)
where
3.2. Moments
The rth ordinary moment of
can be expressed using (10) as
For any positive integer n, the first four moments of GEWLD are
3.3. Incomplete Moments
The incomplete moment of a distribution is very useful and answers many questions in medicine, economics, reliability analysis, demography and insurance when dealing with the Bonferroni and Lorenz curves.
The sth incomplete moment
, is defined as
Set
,
Then
(11)
where
, is the lower incomplete Gamma function. Besides, the sth conditional moments of X, say
, is defined as
Hence by (3.42), we obtain
, (12)
where
, is the upper incomplete gamma function.
3.4. Generating Function
The moment generating function for GEWL distribution is given by
(13)
3.5. Shannon Entropy
The Shannon Entropy as a measure of randomness or uncertainty [11] has provided important solution to several studies such as anomalous diffusion, DNA sequencing, daily temperature randomness, study of content of information signals, etc. According to [18], the Shannon’s entropy for the GEWLD is will be written as
But T is a Gumbel, Y is logistic, and R is Exponentiated Weibull distribution random variables respectively, hence
in equation is the Shannon entropy of the Gumbel distribution given by
, where
is the Euler Mascheroni constant.
is approximately equal to 0.5772.
But
Recall that
which implies that
Therefore
Substituting
and
in
we have
(14)
3.6. Quantiles
The pth quantile
of the Gumbel-Exponentiated Weibull{Logistic} distribution is obtained by direct reversal of the cdf given in (3)
(15)
The first three quantiles of GEWLD can be obtained by using
above for
and
. Generally, to generate a random variate “x” from GEWL from the uniform random number “u” we can use the formula
(16)
3.7. Order Statistics
Let
be a random sample GEWL distribution, Suppose
denote the kth order statistic. The probability density function of the
can be expressed as
(17)
where
Given the general expression of the cdf (9) and pdf (10) of the GEWL distribution, we get the pdf of the kth order statistic for the GEWL distribution as
The rth moment of the kth order statistic of the Gumbel-exponentiated Weibull {Logistic} distribution is given by
Let
(18)
4. Estimation
The model parameters of the GEWL distribution can be estimated using the method of maximum likelihood. Suppose
is a random independent sample size n from the GEWLD distribution with parameter vector
, then log-likelihood function is expressed as
(19)
Solving the log-likelihood function analytically is complex though it can be maximized numerically by using global iterative optimization methods available with software like R and Mathematica. Taking partial derivatives of the log-likelihood function with respect to the parameters we obtain the follows:
(20)
(21)
(22)
(23)
(24)
Solving the equations
simultaneously gives the maximum likelihood estimate (MLE)
of
.
Inverting the Fisher information matrix
which can be gotten using the second partial derivatives of the log-likelihood function w. r. t each parameter
, the asymptotic variance-covariance matrix of the MLE parameters is obtained.
As
,
.
5. Simulation Study
In this section, assessment of the performance of the MLEs of the GEWLD parameters given in (20), (21), (22), (23) and (24) based on simulation study with respect to sample size n is carried out. The simulation is repeated for
times each with sample size
and 700 and arbitrarily, the fixed choice of parameter values is
,
,
,
,
. The evaluation of estimates is based on mean estimates (ME), average bias (AVB), roots mean square error (RMSE), average width (AW) and coverage probability (CP) using the R-software and result presented in Table 1. The results in Table 1 indicate that the estimates are stable and the values of the mean estimates approached true values as the sample size increases. Moreover, Table 1 indicates that the values of the RMSE, AVB and AW decreased as the sample size increases and that the coverage probabilities (CP) of the confidence intervals are quite close to the nominal level of 95%.
6. Applications
In this section, empirical illustrations of the GEWLD using real-life datasets are
Table 1. Results of Monte Carlo simulations for the GEWLD M, AVB, RMSE, AW and CP.
analyzed. All the datasets used are based on complete observations in different lifetime situations showing the usefulness of the GEWLD. The goodness-of-fit criterion, Akaike information criterion, Kolmogorov-Smirnov (K-S) statistic and p-value are used to compare the GEWLD with other competing models.
The first dataset represents the Breaking Stress of Carbon Fibers of 50 mm Length (GPa). The data was obtained from [20]. The dataset is unimodal and is approximately symmetric (Skewness = −0.1285 and excess kurtosis = 0.1261). Six distributions are used to fit the dataset as shown in Table 2 namely: the proposed Gumbel-exponentiated Weibull {logistic} distribution (GEWLD), the Weibull distribution (WD) [21], the Gumbel distribution (GD) [22], the exponentiated Weibull distribution (EWD) [6], the beta exponential distribution
(standard error of estimate in parenthesis)
Table 2. Maximum likelihood fit of the breaking stress of carbon Fibers of 50 mm Length (GPa).
(BED) [23] and the beta generalized exponential distribution (BGED) [24]. Results in Table 2 clearly show that the proposed GEWLD not only provided a good fit to the dataset but also outperformed the other distribution in fitting the data set given that its p-value of the K-S statistic is highest and very close to unity. Again, the p-value of the K-S statistic for all the fitted distributions is greater than the nominal 0.05 level of significance indicating that all the distributions fitted the data considerably well and that the proposed GEWLD presented the best fit. From Figure 3, the pdfs and cdfs plots suggest good adjustments to the dataset and that GEWLD fits almost symmetric and unimodal dataset very well.
The data set represents Kevlar 49/epoxy strands failure times data (pressure at 70%). The dataset was used in the work of [19] to compare the fit of the Gumbel-Weibull distribution (GWD) with that of the EWD and the beta-normal distribution (BND) [25]. The data is multimodal, platykurtic, and approximately symmetric. (Skewness = 0.0938, excess kurtosis = −0.9154). Four distributions are used to fit the dataset namely: the proposed GEWLD, the EWD, the GWD and the BND. The results of the maximum likelihood fit of all the distributions are contained in Table 3. Results in Table 3 clearly show that the proposed GEWLD not only provided a good fit to the bimodal dataset but also outperformed the other distribution in fitting the data set given that its p-value of the K-S statistic is highest and very close to unity. Again, the p-value of the K-S statistic for all the fitted distributions is greater than the nominal 0.05 level of significance indicating that all the distributions fitted the data considerably well and that the proposed GEWLD presented the best fit. The estimated pdfs and cdfs are presented in Figure 4 and the figure indicated GEWLD provides a better than the competing models.
The third dataset represents Kevlar 49/epoxy strands failure times data (pressure at 90%). The dataset was collected by [26] and obtained from [27]. The dataset is highly skewed to the right with a reverse-J shape. The data is unimodal and leptokurtic (Skewness = 2.9573, excess kurtosis = 13.3798) (Figure 5).
Figure 3. Fitted pdfs (3a) and cdfs (3b) of GEWLD and other completive models for dataset 1.
Figure 4. Fitted pdfs (4a) and cdfs (4b) for the bimodal Kevlar 49/epoxy strands failure times data (pressure at 70%).
Table 3. Maximum likelihood fit of the Kevlar 49/epoxy strands failure times data (pressure at 70%).
(standard error of estimate in parenthesis).
Figure 5. Fitted pdfs (5a) and cdfs (5b) for the right-skewed and heavy-tailed Kevlar 49/epoxy strands failure times data (pressure at 90%).
7. Conclusion
This work has introduced a new generalized exponentiated Weibull univariate continuous distribution called the Gumbel-exponentiated Weibull {Logistic} distribution (GEWLD). The probability density function showed that the distribution is capable of modeling bimodal lifetime dataset. Several mathematical properties of the distribution were derived. Certain characterization results were also presented in the work. Using the method of maximum likelihood estimation method, the estimates of the model parameters were obtained. The new model provides adequate fits when compared with other models competing very well in terms of AICs but best when the smallest K-S statistics and highest P-values are considered.