Research Article  Open Access
Masood Anwar, Jawaria Zahoor, "The HalfLogistic 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 HalfLogistic Lomax Distribution for Lifetime Modeling
Abstract
We introduce a new twoparameter lifetime distribution called the halflogistic Lomax (HLL) distribution. The proposed distribution is obtained by compounding halflogistic 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 maximumlikelihood estimators (MLEs). The flexibility and potentiality of the proposed model are illustrated by means of real and simulated data sets.
1. Introduction
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 halflogistic 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 twoparameter continuous lifetime distribution by compounding the halflogistic and the Lomax distribution called halflogistic Lomax (HLL) distribution.
The Lomax [1] (or Pareto TypeII) distribution was introduced to model business failure data. For more detail about the Lomax distribution, we refer the readers to Rady et al. [2], Tahir et al. [3], and the references therein. In literature, there are several generalizations of the Lomax distribution. AbdulMoniem [4] developed the exponentiated Lomax distribution, and AlAwadhi and Ghitany [5] introduced the discrete Poisson–Lomax distribution by using the Lomax distribution as a mixing distribution for the Poisson parameter. Asgharzadeh et al. [6] proposed the Pareto Poisson–Lindely distribution, and Cordeiro et al. [7] investigated the gammaLomax distribution and studied its properties. Ghitany et al. [8] and Gupta et al. [9] considered the Marshal–Olkin approach and extended the Lomax distribution, and Lemonte and Cordeiro [10] proposed and studied the McDonaldLomax, 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, AlZahrani and Sagor [11], ElBassiouny et al. [12], Rady et al. [2], Kilany [13], and Tahir et al. [3]. 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. [14] define the cdf of the new type I halflogisticG (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 halflogistic (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 halflogistic and the Lomax distribution called halflogistic Lomax (HLL) distribution. The proposed distribution is heavytailed and has a decreasing or upsidedown bathtub (or unimodal) shaped hazard rates depending on its parameters. Upsidedown bathtub shaped hazard rates are common in reliability, engineering, and survival analysis. The HLL distribution can also be applied in engineering as the Lomax [1] distribution and can be a useful alternative to other wellknown densities in lifetime applications. It is interesting to note that the HLL distribution is a special case of MarshallOlkin–Lomax distribution introduced by Ghitany et al. [8]. We obtain some mathematical properties of the proposed distribution and parameters of the model are estimated by the maximumlikelihood estimation method.
The rest of this paper is organized as follows. In Section 2, we introduce the halflogistic 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 halflogisticLomax (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
Using (6) and (7), the survival function (sf) and hazard rate function (hrf) of are, respectively, given by
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 [15], we determine the parameter intervals for which the hazard rate function of the HLL distribution is decreasing or upsidedown bathtub. Let
Then its first derivative is given by
If , then for all . Then the hazard rate function is decreasing. By Glaser’s theorem [15], 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 upsidedown bathtubshaped for . Figure 2 illustrates some of the possible shapes of hazard rate function for selected values of the parameters.
(a)
(b)
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.
3.3. Moments
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.
3.6. Entropies
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 [16], 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
Proof. Inserting (6) and (7) into (44), we get the result.
4. Estimation and Asymptotic Distribution
We consider the maximumlikelihood estimation of the parameters , . The distributional properties of , are obtained using Fisher information matrix.
4.1. MaximumLikelihood Estimation
Let be a random sample of size from distribution. Then the corresponding likelihood function is The loglikelihood function is given byCalculating the firstorder partial derivatives of (46) with respect to , and equating them to zero, we get the following nonlinear equations: To find out the maximumlikelihood 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 NewtonRaphson method, to obtain the solution.
4.2. Asymptotic Distribution
The normal approximation for the largesample for the MLE can be used to compute confidence intervals (CIs). For finding the observed information matrix of , , we compute the secondorder partial derivatives of (46) in the appendix.
The variancecovariance 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 maximumlikelihood 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. Application
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 [17]. 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 [1], with pdf is(ii)The McLomax distribution, introduced by Lemonte and Cordeiro [10], 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. [18], with pdf is(vi)The LomaxLogarithmic (LL) distribution, introduced by AlZahrani and Sagor [11], 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 [19]. 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 upsidedown bathtub [19]. 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 HannanQuinn information criterion (HQIC). where denotes the loglikelihood 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 maximumlikelihood 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 [20].


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 goodnessoffit tests given above. The values of best fitted model are highlighted in the tables. Based on the goodnessoffit tests, we conclude that the HLL distribution has better fit than the other competing models.
7. Concluding Remarks
In this paper, we introduce a twoparameter HLL distribution by compounding the halflogistic and the Lomax distributions. The shape of hazard function of the new compounding distribution can be monotonically decreasing or upsidedown 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.
Appendix
Elements of Observed Information Matrix
The elements of observed information matrix () are given by
Disclosure
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.
Authors’ Contributions
The authors have equally contributed to this work. All authors read and approved the final manuscript.
References
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 M. H. Tahir, G. M. Cordeiro, M. Mansoor, and M. Zubair, “The WeibullLomax distribution: properties and applications,” Hacettepe Journal of Mathematics and Statistics, vol. 44, no. 2, pp. 455–474, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 I. B. AbdulMoniem, “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. View at: Google Scholar
 S. A. AlAwadhi and M. E. Ghitany, “Statistical properties of PoissonLomax distribution and its application to repeated accidents data,” Journal of Applied Statistical Science, vol. 10, no. 4, pp. 365–372, 2001. View at: Google Scholar  MathSciNet
 A. Asgharzadeh, H. S. Bakouch, and L. Esmaeili, “Pareto PoissonLindley distribution with applications,” Journal of Applied Statistics, vol. 40, no. 8, pp. 1717–1734, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 G. M. Cordeiro, E. M. Ortega, and B. z. Popović, “The gammaLomax distribution,” Journal of Statistical Computation and Simulation, vol. 85, no. 2, pp. 305–319, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 M. E. Ghitany, F. A. AlAwadhi, and L. A. Alkhalfan, “MarshallOlkin extended Lomax distribution and its application to censored data,” Communications in Statistics—Theory and Methods, vol. 36, no. 912, pp. 1855–1866, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 R. C. Gupta, M. E. Ghitany, and D. K. AlMutairi, “Estimation of reliability from MarshallOlkin extended Lomax distributions,” Journal of Statistical Computation and Simulation, vol. 80, no. 78, pp. 937–947, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 A. J. Lemonte and G. M. Cordeiro, “An extended Lomax distribution,” Statistics, vol. 47, no. 4, pp. 800–816, 2013. View at: Publisher Site  Google Scholar
 B. AlZahrani and H. Sagor, “Statistical analysis of the Lomaxlogarithmic distribution,” Journal of Statistical Computation and Simulation, vol. 85, no. 9, pp. 1883–1901, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 A. H. ElBassiouny, N. F. Abdo, and H. S. Shahen, “Exponential lomax distribution,” International Journal of Computer Applications, vol. 121, no. 13, pp. 24–29, 2015. View at: Publisher Site  Google Scholar
 N. M. Kilany, “Weighted Lomax distribution,” SpringerPlus, vol. 5, no. 1, article no. 1862, 2016. View at: Publisher Site  Google Scholar
 G. M. Cordeiro, M. Alizadeh, and P. R. D. Marinho, “The type I halflogistic family of distributions,” Journal of Statistical Computation and Simulation, vol. 86, no. 4, pp. 707–728, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 R. E. Glaser, “Bathtub and related failure rate characterizations,” Journal of the American Statistical Association, vol. 75, no. 371, pp. 667–672, 1980. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 M. V. Aarset, “How to identify a bathtub hazard rate,” IEEE Transactions on Reliability, vol. 36, no. 1, pp. 106–108, 1987. View at: Publisher Site  Google Scholar
 R Core Team, R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Austria, 2015, https://www.Rproject.org/.
Copyright
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.