Research Article | Open Access
Masood Anwar, Jawaria Zahoor, "The Half-Logistic Lomax Distribution for Lifetime Modeling", Journal of Probability and Statistics, vol. 2018, Article ID 3152807, 12 pages, 2018. https://doi.org/10.1155/2018/3152807
The Half-Logistic Lomax Distribution for Lifetime Modeling
We introduce a new two-parameter lifetime distribution called the half-logistic Lomax (HLL) distribution. The proposed distribution is obtained by compounding half-logistic and Lomax distributions. We derive some mathematical properties of the proposed distribution such as the survival and hazard rate function, quantile function, mode, median, moments and moment generating functions, mean deviations from mean and median, mean residual life function, order statistics, and entropies. The estimation of parameters is performed by maximum likelihood and the formulas for the elements of the Fisher information matrix are provided. A simulation study is run to assess the performance of maximum-likelihood estimators (MLEs). The flexibility and potentiality of the proposed model are illustrated by means of real and simulated data sets.
The commonly used lifetime distributions (exponential, gamma, Weibull, Lomax, lognormal, etc.) have a limited range of behavior and do not provide adequate fit to complex data sets in different sciences. Generalizations of these distributions offer more flexibility and provide reasonable parametric fits to complex data sets. Motivated by the various applications of Lomax and half-logistic distributions in areas of income and wealth inequality, firm size, size of cities, queuing problems, actuarial science, medical and biological sciences, and engineering, we propose a two-parameter continuous lifetime distribution by compounding the half-logistic and the Lomax distribution called half-logistic Lomax (HLL) distribution.
The Lomax  (or Pareto Type-II) distribution was introduced to model business failure data. For more detail about the Lomax distribution, we refer the readers to Rady et al. , Tahir et al. , and the references therein. In literature, there are several generalizations of the Lomax distribution. Abdul-Moniem  developed the exponentiated Lomax distribution, and Al-Awadhi and Ghitany  introduced the discrete Poisson–Lomax distribution by using the Lomax distribution as a mixing distribution for the Poisson parameter. Asgharzadeh et al.  proposed the Pareto Poisson–Lindely distribution, and Cordeiro et al.  investigated the gamma-Lomax distribution and studied its properties. Ghitany et al.  and Gupta et al.  considered the Marshal–Olkin approach and extended the Lomax distribution, and Lemonte and Cordeiro  proposed and studied the McDonald-Lomax, the beta Lomax, and the Kumaraswamy Lomax distributions. Other models constitute flexible family of distributions in terms of the variates of shapes and hazard functions; see, for example, Al-Zahrani and Sagor , El-Bassiouny et al. , Rady et al. , Kilany , and Tahir et al. . These generalizations of the Lomax distribution are considered to be useful life distribution models.
The cumulative distribution function (cdf) of Lomax distribution is given by where is a shape parameter and is a scale parameter. The probability density function (pdf) corresponding to (1) is
Cordeiro et al.  define the cdf of the new type I half-logistic-G (TIHL-) family of distributions by where is the baseline cdf depending on a parameter vector and an additional shape parameter . As a special case, if , then the TIHL- is the half-logistic- (HL-) distribution with cdf
The corresponding pdf to (4) is given by
This paper aims to provide a new lifetime model with a minimum number of parameters by compounding the half-logistic and the Lomax distribution called half-logistic Lomax (HLL) distribution. The proposed distribution is heavy-tailed and has a decreasing or upside-down bathtub (or unimodal) shaped hazard rates depending on its parameters. Upside-down bathtub shaped hazard rates are common in reliability, engineering, and survival analysis. The HLL distribution can also be applied in engineering as the Lomax  distribution and can be a useful alternative to other well-known densities in lifetime applications. It is interesting to note that the HLL distribution is a special case of Marshall-Olkin–Lomax distribution introduced by Ghitany et al. . We obtain some mathematical properties of the proposed distribution and parameters of the model are estimated by the maximum-likelihood estimation method.
The rest of this paper is organized as follows. In Section 2, we introduce the half-logistic Lomax distribution and provide plots of its density function. In Section 3, we investigate various mathematical properties of the HLL distribution including survival and hazard rate function, quantile function, moments, mean residual life function, mean deviation from the mean and the mean deviation from the median, entropies, and order statistics. In Section 4, estimation of parameters is given by MLE method and the asymptotic distribution of the estimators is studied via Fisher’s information matrix. Simulation results on the behavior of the MLEs are presented in Section 5. A real data application is conducted in Section 6. Finally, in Section 7, we conclude that the HLL distribution is the best model as compared to other competing models.
2. The HLL Distribution
We define the half-logistic-Lomax (HLL) density function by inserting (1) and (2) into (5). So, we obtain The corresponding cumulative density function (cdf) follows from (1) and (4) and is given byHereafter, we will denote a random variable having pdf (6) by . The limit of the HLL density (6) as is 0 and the limit as is . Figure 1 depicts some of the possible shapes of density (6) for selected parameter values. The mode of density (6) is obtained from solving , which is given by
3. Mathematical Properties
This section describes mathematical properties of the HLL distribution.
3.1. The Survivor and Hazard Rate Functions
The limit of as is ; that is, is bounded from below and continuous in its parameters. The limit of as is 0. According to Glaser , we determine the parameter intervals for which the hazard rate function of the HLL distribution is decreasing or upside-down bathtub. Let
Then its first derivative is given by
If , then for all . Then the hazard rate function is decreasing. By Glaser’s theorem , it is sufficient to show that there exists such that for all . implies that and for all , which implies unimodal shape of the hazard rate function. For , for all and for all . Thus, the hazard rate function is upside-down bathtub-shaped for . Figure 2 illustrates some of the possible shapes of hazard rate function for selected values of the parameters.
3.2. Quantile and Random Number Generation
The cumulative distribution function is given by (7). Inverting , we obtain
Equation (12) can be used to simulate the HLL variable.
The following theorem gives moments of the HLL distribution and its mean moments and cumulants.
Theorem 1. The th moment about the origin of is given by
Proof. We haveSetting , yields
Corollary 2. The first four moments of are given, respectively, by
Proof. Putting in (13) yields the desired result.
The mean moments and cumulants of can be obtained by the following relation;Here . We have , , , and so on.
The skewness and kurtosis follow from the second, third, and fourth cumulants.
The variance of , denoted by , is given byThe th descending factorial moment of is where is the Sterling number of the first kind which counts the number of ways to permute a list of items into cycles.
Theorem 3. The moment generating function of is given by
Proof. We have Setting , yieldsAfter simplification, we get
3.4. Mean Residual Life (MRL) Function
The MRL function is important in reliability, survival analysis, actuarial sciences, economics, and social sciences for characterizing lifetime distributions. It also plays an important role in repair and replacement strategies and summarizes the entire residual life function.
Theorem 4. The MRL function of is given by
Proof. The MRL function is defined asTherefore, we haveSetting , yieldsAfter simplification, we get (25).
3.5. Mean Deviations
The mean deviations about mean and median are, respectively, defined by
Theorem 5. The mean deviations about mean and mean deviation about median of are given by
Proof. Let .
By inserting pdf of the HLL into (29) and setting , , we get Using the series representations , where is positive integer, we getSetting into (29) and after some manipulations we get the desired result.
Similarly, the measure can be obtained.
The entropy of a random variable is a measure of variation of uncertainty. Two popular entropies are the Rényi and Shannon entropies.
The Rényi entropy is defined as The Rényi entropy for the HLL distribution is given bySetting , yields Thus, we haveAnother entropy measure defined by is known as the Shannon entropy and plays a similar role as the kurtosis measure in comparing the shapes of densities and measuring tail heaviness. We have By using the binomial series expansion and (2.727) of Gradshteyn and Ryzhik , it follows that Hence, the Shannon entropy can be expressed in the form
3.7. Order Statistics
Let be a random sample of size from a distribution with pdf and cdf . Let denote the corresponding order statistics. Then the pdf and cdf of are where is the complete beta function.
Theorem 6. Let and be the pdf and cdf of . Then the pdf and cdf of are
4. Estimation and Asymptotic Distribution
We consider the maximum-likelihood estimation of the parameters , . The distributional properties of , are obtained using Fisher information matrix.
4.1. Maximum-Likelihood Estimation
Let be a random sample of size from distribution. Then the corresponding likelihood function is The log-likelihood function is given byCalculating the first-order partial derivatives of (46) with respect to , and equating them to zero, we get the following nonlinear equations: To find out the maximum-likelihood estimates (MLEs) of , , we have to solve the above nonlinear equations. Apparently, there is no closed form solution in , . We have to use a numerical technique method, such as Newton-Raphson method, to obtain the solution.
4.2. Asymptotic Distribution
The normal approximation for the large-sample for the MLE can be used to compute confidence intervals (CIs). For finding the observed information matrix of , , we compute the second-order partial derivatives of (46) in the appendix.
The variance-covariance matrix may be approximated by , where is the observed information matrix. Since involves the parameters , , we replace the parameters by the corresponding MLEs in order to obtain an estimate of , which is denoted by , where , when , is replaced by , . Using these results, an approximate CI for , , respectively, is given by where is the upper percentile of the standard normal distribution.
5. Simulation Study
The assessment of the maximum-likelihood estimates is based on a simulation study. The following steps were followed:(1)Generate 5,000 samples of size from (6). The HLL variables are generated using given in (12), where .(2)Compute the MLEs for the 5,000 samples, say for (3)Compute the mean square errors (MSE) for each parameter.
We repeat these steps for with , . The comparison is based on MSEs. Figure 3 plots the MSEs of the MLEs of two parameters. The assessment based on this simulation study is that the MSEs for each parameter decrease to zero with increasing sample size.
6.1. Cancer Data
We have considered a real data set, which records the remission times (in months) of a random sample of 128 bladder cancer patients reported in . We have fitted the HLL distribution to the data set using MLE and compared the proposed distribution with the following distributions:(i)The Lomax distribution, introduced by Lomax , with pdf is(ii)The McLomax distribution, introduced by Lemonte and Cordeiro , with pdf is(iii)The beta Lomax (BLomax) distribution is a submodel of McLomax for with pdf(iv)The Kumaraswamy Lomax (KwLomax) distribution is a submodel of McLomax for , with pdf(v)The Exponentiated Standard Lomax (ESLomax) distribution, considered by Gupta et al. , with pdf is(vi)The Lomax-Logarithmic (LL) distribution, introduced by Al-Zahrani and Sagor , with pdf isFor identification of the shape of the hazard function of the cancer data set, we use a graphical method based on the Total Time on Test (TTT) plot . Let be a random variable with nonnegative values, which represents the survival time. The empirical TTT curve is constructed by plotting , against , where and are the order statistics of the sample. If the empirical TTT transform is straight diagonal, hrf is constant; if it is convex, hrf is decreasing; if it is concave, hrf is increasing; if it is convex and then concave, hrf is bathtub; and if it is concave and then convex, hrf is upside-down bathtub . Figure 4 shows the empirical TTT plot for the data set. The shape of the data set is unimodal, as its TTT plot is concave and then convex. The model selection is carried out using Akaike information criterion (AIC), the Bayesian information criterion (BIC), the consistent Akaike information criterion (CAIC), and the Hannan-Quinn information criterion (HQIC). where denotes the log-likelihood function evaluated at the MLEs, is the number of model parameters, and is the sample size. The model with the lowest values for these statistics could be chosen as the best model to fit the data. First, we describe the data set in Table 1. Then, we report the maximum-likelihood estimators (and the corresponding standard errors in parentheses) of the parameters and values of AIC, BIC, HQIC, and CAIC in Table 2. In Table 2, we provide some results reported in [10, 11]. We conclude that our new model HLL has the smallest AIC, BIC, HQIC, and CAIC values among all fitted models, and so it could be chosen as the best model. These results are obtained using the AdequacyModel package version 1.0.8 available for the programming language R .
6.2. Simulated Data
Here, we used a simulation study to check the flexibility of the proposed distribution. We have generated a sample of size from the six distributions and compared the fit of the HLL distribution with competing models.
The details of generated samples with the specified values of the parameters from the six models are as follows:(1)Lomax distribution with and (see Table 3)(2)McLomax distribution with parameters , , , , and (see Table 4)(3)BLomax distribution with parameters , , , and (see Table 5)(4)KwLomax distribution with parameters , , , and (see Table 6)(5)ESLomax distribution with parameters and (see Table 7)(6)LL distribution with parameters , , and (see Table 8)
The comparison is based upon the formal goodness-of-fit tests given above. The values of best fitted model are highlighted in the tables. Based on the goodness-of-fit tests, we conclude that the HLL distribution has better fit than the other competing models.
7. Concluding Remarks
In this paper, we introduce a two-parameter HLL distribution by compounding the half-logistic and the Lomax distributions. The shape of hazard function of the new compounding distribution can be monotonically decreasing or upside-down bathtub (unimodal). Some mathematical and statistical properties of the new model including moments and moment generating functions, mean deviations from mean and median, quantile function, mean residual life function, mode, median, order statistics, and entropies are studied. We estimate the model parameters by maximum likelihood and determined the observed information matrix. We present a simulation study to illustrate the performance of MLEs. The flexibility and potentiality of the proposed model are illustrated by means of a real data set. We hope that the HLL distribution may attract wider range of applications in areas such as engineering, survival and lifetime data, economics, meteorology, hydrology, and others.
Elements of Observed Information Matrix
The elements of observed information matrix () are given by
The authors acknowledge that the abstract of the manuscript had been presented in a conference paper under the link http://archive.stat.uconn.edu/lida17/program.pdf.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
The authors have equally contributed to this work. All authors read and approved the final manuscript.
- K. S. Lomax, “Business failures: another example of the analysis of failure data,” Journal of the American Statistical Association, vol. 49, no. 268, pp. 847–852, 1954.
- E.-H. A. Rady, W. A. Hassanein, and T. A. Elhaddad, “The power Lomax distribution with an application to bladder cancer data,” SpringerPlus, vol. 5, 1838 pages, 2016.
- M. H. Tahir, G. M. Cordeiro, M. Mansoor, and M. Zubair, “The Weibull-Lomax distribution: properties and applications,” Hacettepe Journal of Mathematics and Statistics, vol. 44, no. 2, pp. 455–474, 2015.
- I. B. Abdul-Moniem, “Recurrence relations for moments of lower generalized order statistics from exponentiated Lomax distribution and its characterization,” International Journal of Mathematical Archive, vol. 3, pp. 2144–2150, 2012.
- S. A. Al-Awadhi and M. E. Ghitany, “Statistical properties of Poisson-Lomax distribution and its application to repeated accidents data,” Journal of Applied Statistical Science, vol. 10, no. 4, pp. 365–372, 2001.
- A. Asgharzadeh, H. S. Bakouch, and L. Esmaeili, “Pareto Poisson-Lindley distribution with applications,” Journal of Applied Statistics, vol. 40, no. 8, pp. 1717–1734, 2013.
- G. M. Cordeiro, E. M. Ortega, and B. z. Popović, “The gamma-Lomax distribution,” Journal of Statistical Computation and Simulation, vol. 85, no. 2, pp. 305–319, 2015.
- M. E. Ghitany, F. A. Al-Awadhi, and L. A. Alkhalfan, “Marshall-Olkin extended Lomax distribution and its application to censored data,” Communications in Statistics—Theory and Methods, vol. 36, no. 9-12, pp. 1855–1866, 2007.
- R. C. Gupta, M. E. Ghitany, and D. K. Al-Mutairi, “Estimation of reliability from Marshall-Olkin extended Lomax distributions,” Journal of Statistical Computation and Simulation, vol. 80, no. 7-8, pp. 937–947, 2010.
- A. J. Lemonte and G. M. Cordeiro, “An extended Lomax distribution,” Statistics, vol. 47, no. 4, pp. 800–816, 2013.
- B. Al-Zahrani and H. Sagor, “Statistical analysis of the Lomax-logarithmic distribution,” Journal of Statistical Computation and Simulation, vol. 85, no. 9, pp. 1883–1901, 2015.
- A. H. El-Bassiouny, N. F. Abdo, and H. S. Shahen, “Exponential lomax distribution,” International Journal of Computer Applications, vol. 121, no. 13, pp. 24–29, 2015.
- N. M. Kilany, “Weighted Lomax distribution,” SpringerPlus, vol. 5, no. 1, article no. 1862, 2016.
- G. M. Cordeiro, M. Alizadeh, and P. R. D. Marinho, “The type I half-logistic family of distributions,” Journal of Statistical Computation and Simulation, vol. 86, no. 4, pp. 707–728, 2016.
- R. E. Glaser, “Bathtub and related failure rate characterizations,” Journal of the American Statistical Association, vol. 75, no. 371, pp. 667–672, 1980.
- I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series, and Products, Academic Press, New York, NY, USA, 2007.
- E. T. Lee and J. W. Wang, Statistical Methods for Survival Data Analysis, John Wiley and Sons, New York, NY, USA, 3rd edition, 2003.
- R. C. Gupta, P. L. Gupta, and R. D. Gupta, “Modeling failure time data by Lehman alternatives,” Communications in Statistics—Theory and Methods, vol. 27, no. 4, pp. 887–904, 1998.
- M. V. Aarset, “How to identify a bathtub hazard rate,” IEEE Transactions on Reliability, vol. 36, no. 1, pp. 106–108, 1987.
- R Core Team, R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Austria, 2015, https://www.R-project.org/.
Copyright © 2018 Masood Anwar and Jawaria Zahoor. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.