Alpha Power Lomax Distribution: Properties and Application
Y. Murat Bulut1,*, Fatma Zehra Doğru2 and Olcay Arslan3
1Department of Statistics, Eskisehir Osmangazi University, 26040 Eskisehir/Turkey
2Department of Statistics, Giresun University, 28200 Giresun/Turkey
3Department of Statistics, Ankara University, 06100 Ankara/Turkey
E-mail: ymbulut@ogu.edu.tr; fatma.dogru@giresun.edu.tr; oarslan@ankara.edu.tr
*Corresponding Author
Received 22 July 2020; Accepted 07 January 2021; Publication 08 March 2021
This study offers a newly proposed distribution called alpha power Lomax (APL) distribution as a new extension of the Lomax distribution using the alpha power transformation (APT) method. Some distributional properties of newly defined distribution such as density function, moments, hazard and survival functions, orders statistics etc. are investigated. Parameters of the APL distribution are estimated with the help of the maximum likelihood (ML) estimation method. The applicability of the APL distribution is conducted through a simulation study and a real data example.
Keywords: Lomax distribution, alpha power transformation, maximum likelihood estimation.
Lomax [1] proposed Lomax (or Pareto distribution of the second kind) distribution. The proposed distribution aims to model the business failure data set. Also, Hassan and Al-Ghamdi [2] used this distribution to model reliability data set and life testing. To model income and wealth data sets, Harris [3, 4] used this distribution. Lomax distribution is used to model data set from receiver operating characteristic (ROC) curves analysis in the work of Campbell and Ratnaparkhi [5]. Balakrishnan and Ahsanullah [6] investigate recurrence relations between the moments of record values. Bryson [7] proposed that the Lomax distribution is a heavy-tailed alternative to the exponential distribution. The reader can be found more detail about Lomax distribution in the book of Johnson et al. [8].
In recent years, applications in various fields of sciences indicate that classical distribution is not enough to model data sets. So, it is necessary to expand some popular distributions to model real-life data sets. Many authors extend classical distributions to apply in many fields. For example, El-Bassiouny et al. [9] proposed exponential Lomax distribution. Gamma-Lomax distribution is proposed by Cordeiro et al. [10]. Ghitany et al. [11] extended the Lomax distribution using the Marshall-Olkin distribution. Tahir et al. [12] proposed the Weibull-Lomax distribution. McDonald Lomax distribution includes several distributions as sub-models such as beta Lomax, Kumaraswamy Lomax, exponentiated Lomax, and exponentiated standart Lomax distributions, was introduced by Lemonte and Cordeiro [13]. Discrete Poisson-Lomax distribution is introduced by Al-Awadhi and Ghitany [14]. Rady et al. [15] proposed the power Lomax distribution using the power transformation method. Al-Marzouki [16] introduced exponentiated power Lomax distribution as an alternative lifetime distribution.
To obtain more flexible distribution than the usual ones, adding extra parameters to a well-established distribution family is a useful tool. In literature, there are a lot of methods to obtain a more flexible distribution than usual. One of them is proposed by Mahdavi and Kundu [17]. This method is the alpha power transformation (APT) method.
In this work, a new distribution using the APT method based on the Lomax distribution is proposed. Alpha power Weibull distribution using the APT method is defined by Nassar et al. [18]. Dey et al. [19] extended the generalized exponential distribution using the APT method to model ozone data. Dey et al. [20] introduced a Logarithmic family of distribution using the APT method. Also, alpha power transformed inverse Lindley distribution using the APT method is established by Dey et al. [21]. After work of Mahdavi and Kundu [17], Nassar et al. [22] extend APT class to the Marshall-Olkin alpha power family of distributions.
The rest of the paper is designed as follows. Section 2 defines a newly proposed distribution called APL distribution and investigates some distributional properties of this distribution. Section 3 estimates parameters of the APL distribution via the MLE method. Section 4 provides a simulation study and a real data example to demonstrate the applicability of the APL distribution. Section 5 is devoted to conclusions.
Let the random variable is said to have alpha power Lomax (APL) distribution with the parameters , and , the of is
(1) |
and corresponding cumulative distribution function is
(2) |
In this subsection, moments and moment generating function (mgf) of the APL distribution are obtained. The moments of can be obtained as follows
(3) |
To obtain moments, we use the following series representation
(4) |
Also, the mgf of APL distribution is
(5) |
where is the Beta function.
quantile function of APL distribution can be obtained as
(6) |
Hereby, the median of APL distribution can be obtained as
(7) |
Let be the order statistics, then the pdf of is given as
(8) |
In this subsection, we give hazard rate and survival functions of APL distribution as follows
(9) | ||
(10) |
This section gives the parameter estimation of the APL distribution via the ML estimation method. Suppose that be a random sample taken from the APL distribution with the pdf given in (1). We can write the likelihood function as
(11) |
Then, the log-likelihood () function corresponding to (3) can be obtained as
(12) |
Taking the first derivative of relating to the parameters , and , then setting these equations to zero, the following estimating equations can be obtained
(13) | ||
(14) | ||
(15) |
To obtain the ML estimators of parameters, we have to solve equation systems given in (13)–(15), simultaneously.
In this subsection a small simulation study to illustrate the efficiency of the APL distribution’s ML estimators is given. We generate data from the APL distribution using the quantile function given in Eq. (6). Estimates, bias and MSE (mean squared error) values are computed. We use following formulations to calculate bias and MSE values
(16) |
where
(17) | ||
(18) |
Table 1 Estimates of parameters, bias and MSE for
MLE | |||
, , | |||
Estimate | 3.0165 | 1.0876 | 1.9901 |
Bias | 0.0165 | 0.0876 | 0.0099 |
MSE | 0.1976 | 0.2867 | 0.0845 |
, , | |||
Estimate | 9.9991 | 5.0000 | 1.0252 |
Bias | 0.0009 | 0.0000 | 0.0252 |
MSE | 0.0055 | 0.0549 | 0.1892 |
, , | |||
Estimate | 4.9944 | 10.0068 | 1.0149 |
Bias | 0.0056 | 0.0068 | 0.0149 |
MSE | 0.0243 | 0.0427 | 0.1921 |
, , | |||
Estimate | 2.4980 | 1.0841 | 1.0546 |
Bias | 0.0934 | 0.0841 | 0.0546 |
MSE | 0.3515 | 0.2992 | 0.1554 |
Table 2 Estimates of parameters, bias and MSE for
MLE | |||
, , | |||
Estimate | 2.9970 | 1.0648 | 1.9853 |
Bias | 0.0030 | 0.0648 | 0.0147 |
MSE | 0.0068 | 0.1653 | 0.0393 |
, , | |||
Estimate | 9.9997 | 4.9994 | 1.0140 |
Bias | 0.0003 | 0.0006 | 0.0140 |
MSE | 0.0030 | 0.0318 | 0.1163 |
, , | |||
Estimate | 4.9975 | 10.0042 | 0.9836 |
Bias | 0.0025 | 0.0042 | 0.0164 |
MSE | 0.0079 | 0.0178 | 0.1265 |
, , | |||
Estimate | 2.0934 | 1.0841 | 1.0546 |
Bias | 0.0934 | 0.0841 | 0.0546 |
MSE | 0.3515 | 0.2992 | 0.1554 |
Table 3 Estimates of parameters, bias and MSE for
MLE | |||
, , | |||
Estimate | 2.9976 | 1.0530 | 1.9879 |
Bias | 0.0024 | 0.0530 | 0.0121 |
MSE | 0.0059 | 0.1426 | 0.0338 |
, , | |||
Estimate | 9.9992 | 5.0054 | 0.9909 |
Bias | 0.0008 | 0.0054 | 0.0091 |
MSE | 0.0028 | 0.0288 | 0.0988 |
, , | |||
Estimate | 4.9997 | 9.9991 | 1.0229 |
Bias | 0.0003 | 0.0009 | 0.0229 |
MSE | 0.0068 | 0.0161 | 0.1208 |
, , | |||
Estimate | 2.0542 | 1.0668 | 1.0379 |
Bias | 0.0542 | 0.0668 | 0.0379 |
MSE | 0.2148 | 0.2553 | 0.1121 |
In the simulation design, we get the parameter values as , , and also we set sample sizes as and . Tables 1–4 include the simulation results. In the tables, we give the true values and estimates of parameters, and the values of bias and MSE. From the simulation results, we observe that the parameters of the APL distribution can be calculated with accuracy. Also, when the sample sizes are getting bigger, the values of MSE are getting smaller for the parameter estimates.
Table 4 Estimates of parameters, bias and MSE for
MLE | |||
, , | |||
Estimate | 2.9982 | 1.0402 | 1.9914 |
Bias | 0.0018 | 0.0402 | 0.0086 |
MSE | 0.0046 | 0.1107 | 0.0251 |
, , | |||
Estimate | 9.9997 | 4.9998 | 1.0102 |
Bias | 0.0003 | 0.0002 | 0.0102 |
MSE | 0.0025 | 0.0266 | 0.0951 |
, , | |||
Estimate | 4.9986 | 10.0026 | 0.9890 |
Bias | 0.0014 | 0.0026 | 0.0110 |
MSE | 0.0054 | 0.0126 | 0.0937 |
, , | |||
Estimate | 2.0228 | 1.0480 | 1.0107 |
Bias | 0.0228 | 0.0480 | 0.0107 |
MSE | 0.1116 | 0.1781 | 0.0726 |
Table 5 MLEs and AIC and BIC values
Distributions | Estimates | logL | AIC | BIC |
Lomax | -413.832 | 831.67 | 837.37 | |
McLomax | -409.91 | 829.82 | 844.09 | |
BLomax | -411.743 | 831.486 | 842.89 | |
KwLomax | -409.94 | 827.88 | 839.29 | |
ExpLomax | -414.978 | 835.956 | 844.512 | |
G-Lomax | -410.081 | 826.162 | 834.718 | |
TE-Lomax | -410.434 | 828.868 | 840.276 | |
WLomax | -410.811 | 829.622 | 841.03 | |
Power Lomax | -409.74 | 825.48 | 834.036 | |
Alpha Power Lomax | -409.3853 | 824.7707 | 833.3268 | |
In this subsection, a real data example is given to show the superiority of the newly defined distribution over the other Lomax distribution extensions. Remission times data of bladder cancer patients, which have been used by Lee and Wang [23], are used. Proposed distribution is fitted to the dataset by using MLE method. Also, we compare newly proposed distribution with McLomax [13], Exponentiated Lomax [13], Beta Lomax [13], Kumaraswamy Lomax [13], Transmuted Exponentiated Lomax [24], Gamma-Lomax [10], Weighted Lomax [25], Exponential Lomax [9] and Power Lomax [15] distributions. The Akaike information criterion (AIC) and Bayesian information criterion (BIC) are used to compare models. AIC and BIC indicate that the best model is Alpha Power Lomax distribution to bladder cancer data set. MLEs and the measures AIC and BIC are given in Table 5. AIC and BIC values are computed as follows:
(19) | |
(20) |
where is the parameter number and is the sample size. We use optim package in R [26] to estimate parameters.
Until this time, Power Lomax distribution defined by Rady et al. [15] is the best model to bladder cancer data. But, when we look at AIC and BIC values, APL distribution modeled better than the other extensions of Lomax distribution. The histogram and fitted distributions are given in Figure 1. As a result of real data example, we can say that APL distribution is a very competitive distribution to model lifetime data.
In this study, APL distribution has been defined and the statistical properties of the APL have been studied. Also, real data and simulation studies are conducted to show practicability of the distribution. The real data example shows that APL distribution has less AIC and BIC values rather than the other extensions of the Lomax distribution.
[1] K. S. Lomax, Business failures: Another example of the analysis of failure data, Journal of the American Statistical Association 49 (268) (1954) 847–852.
[2] A. Hassan, A. Al-Ghamdi, Optimum step stress accelerated life testing for lomax distributions, Journal of Applied Sciences Research (5) (2009) 2153–2164.
[3] C. M. Harris, The pareto distribution as a queue service descipline, Operations Research 16 (2) (1968) 307–313.
[4] A. B. Atkinson, A. J. Harrison, Distribution of Personal Wealth in Britain, Cambridge University Press.
[5] G. Campbell, M. V. Ratnaparkhi, An application of lomax distributions in receiver operating characteristic (roc) curve analysis, Communications in Statistics–Theory and Methods 22 (6) (1993) 1681–1697.
[6] N. Balakrishnan, M. Ahsanullah, Relations for single and product moments of record values from lomax distribution, Sankhya: The Indian Journal of Statistics, Series B 56 (2) (1994) 140–146.
[7] M. C. Bryson, Heavy-tailed distribution: properties and tests, Technometrics 16 (1) (1974) 61–68.
[8] N. L. Johnson, S. Kotz, N. Balakrishnan, Contineous Univariate Distributions Vol. 1, 2nd Edition, Wiley, New York, 1994.
[9] A. El-Bassiouny, N. Abdo, H. Shahen, Exponential lomax distribution, International Journal of Computer Applications 121 (13) (2015) 24–29.
[10] G. M. Cordeiro, E. M. Ortega, B. V. Popović, The gamma-lomax distribution, Journal of Statistical Computation and Simulation 85 (2) (2015) 305–319.
[11] M. E. Ghitany, F. A. Al-Awadhi, L. A. Alkhalfan, Marshall–olkin extended lomax distribution and its application to censored data, Communications in Statistics - Theory and Methods 36 (10) (2007) 1855–1866.
[12] M. H. Tahir, G. M. Cordeiro, M. Mansoor, M. Zubair, The weibull-lomax distribution: properties and applications, Hacettepe Journal of Mathematics and Statistics 44 (2) (2015) 455–474.
[13] A. Lemonte, G. Cordeiro, An extended lomax distribution, Statistics 47 (4) (2013) 800–816.
[14] S. A. Al-Awadhi, M. E. Ghitany, Statistical properties of poisson–lomax distribution and its application to repeated accidents data, Journal of Applied Statistical Science 10 (4) (2001) 365–372.
[15] E.-H. A. Rady, W. A. Hassanein, T. A. Elhaddad, The power lomax distribution with an application to bladder cancer data, SpringerPlus 5 (1) (2016) 1838.
[16] S. Al-Marzouki, A new generalization of power lomax distribution, International Journal of MAthematics and its Applications 7 (1) (2018) 59–68.
[17] A. Mahdavi, D. Kundu, A new method for generating distributions with an application to exponential distribution, Communications in Statistics-Theory and Methods 46 (13) (2017) 6543–6557.
[18] M. Nassar, A. Alzaatreh, M. Mead, O. Abo-Kasem, Alpha power weibull distribution: Properties and applications, Communications in Statistics - Theory and Methods 46 (20) (2017) 10236–10252.
[19] S. Dey, A. Alzaatreh, C. Zhang, D. Kumar, A new extension of generalized exponential distribution with application to ozone data, Ozone: Science & Engineering 39 (4) (2017) 273–285.
[20] S. Dey, M. Nassar, D. Kumar, logarithmic transformed family of distributions with applications, Annals of Data Science 4 (4) (2017) 457–482.
[21] S. Dey, M. Nassar, D. Kumar, Alpha power transformed inverse lindley distribution: A distribution with an upside-down bathtub-shaped hazard function, Journal of Computational and Applied Mathematics 348 (2019) 130–145.
[22] M. Nassar, D. Kumar, S. Dey, G. M. Cordeiro, A. Z. Afify, The marshall–olkin alpha power family of distributions with applications, Journal of Computational and Applied Mathematics 351 (2019) 41–53.
[23] E. T. Lee, J. W. Wang, Statistical Methods for Survival Data Analysis, 3rd Edition, Wiley, New York, 2003.
[24] S. Ashour, M. Eltehiwy, Transmuted exponentiated lomax distribution, Australian Journal of Basic and Applied Sciences 7 (7) (2013) 658–667.
[25] N. M. Kilany, Weighted lomax distribution, SpringerPlus 5 (1) (2016) 1862.
[26] R Core Team, R: A language and environment for statistical computing (2016). URL https://www.R-project.org/
Y. Murat Bulut received his B.Sc. in Mathematics from Çukurova University, Turkey, M.Sc., and Ph.D. degrees from Eskişehir Osmangazi University, Turkey. He is an assistant professor in the Department of Statistics at Eskişehir Osmangazi University, Turkey, since 2016. His research areas are robust statistics, distribution theory, and generalized linear models.
Fatma Zehra Doğru received her B.Sc. and M.Sc. degrees in Statistics and Computer Sciences from the Karadeniz Technical University in 2009 and 2011; and a Ph.D. degree in Statistics from the Ankara University in 2015. She is currently working as an associate professor in the department of statistics at Giresun University, Turkey. Her research interests cover mixture models, mixture regression models, multivariate mixture models, robust statistics, statistical inference, and modeling.
Olcay Arslan is a full Professor of Statistics at the Department of Statistics, Ankara University, Ankara, Turkey. She received her Ph.D. in Statistics from the University of Leeds, Leeds, UK, in 1993. Prior to joining Ankara University, she was working as a full professor at the department of Statistics, Cukurova University, Adana, Turkey. She was visiting scholar at Rutgers University and visiting professor at St. Cloud State University, St. Cloud, USA. Her principal research areas are robust statistical analysis, multivariate analysis, regression analysis, EM algorithm and distribution theory. She also works on model selection, scale mixtures, mean-variance (location-scatter) mixtures, finite mixtures and linear mixed models. She has published over 75 research papers in international refereed journals and also several book chapters in scientific books. She has been involved in organizing many scientific events devoted to statistics including ICORS 2008 (International conference on robust statistics, 2008) and IC-SMHD-2016 (International conference on information complexity and statistical modeling in high dimensions with application, 2016). She is an area editor of Hacettepe Journal of Mathematics and Statistics (HJMS) and members of the editorial board of several scientific journals on statistics and probability. She is also member of the ICORS steering committee.
Journal of Reliability and Statistical Studies, Vol. 14, Issue 1 (2021), 17–32.
doi: 10.13052/jrss0974-8024.1412
© 2021 River Publishers