![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
The Gompertz inverse exponential (GoIE) distribution using the Gompertz generalized family of distributions was derived and introduced in this article. Some basic statistical properties of the model were derived and discussed in minute details. The model parameters were estimated using the maximum likelihood estimation method. Real-life applications were provided and the GoIE distribution provides better fits than the Gompertz exponential, Gompertz Weibull and Gompertz Lomax distributions.
PUBLIC INTEREST STATEMENT
Extending existing probability distributions became necessary in recent years because of the need for models that can adequately fit or describe real-life events. The three-parameter probability model that was introduced in this article is however an extension of the one-parameter inverse exponential distribution. The extra two parameters are shape parameters due to the Gompertz family of distribution that was used. The role of these two extra parameters is to induce skewness into the inverse exponential distribution. As a consequence, the Gompertz inverse exponential distribution appears better than the Gompertz exponential, Gompertz Weibull and Gompertz Lomax distributions when applied to real-life datasets.
1. Introduction
The inverse exponential distribution is a special case of the inverse Weibull distribution; it has been introduced as far back as 1982 by Keller & Kamath (Citation1982) and is capable of modelling datasets with inverted bathtub failure rate. It is a modification of the well-known exponential distribution and has an advantage of not having a constant failure rate. Because of this quality, it can be applied to describe real-life events in engineering, medicine and biology (especially, events with inverted bathtub failure rates). More details about the usefulness of the inverse exponential distribution have been discussed by Bakoban & Abu-Zinadah (Citation2017), Lin, Duran, & Lewis (Citation1989), Oguntunde & Adejumo (Citation2015), Oguntunde, Babatunde, & Ogunmola (Citation2014) and many others.
Recent studies in this line of research involve extending exiting probability distributions with the aim of increasing their modelling capability. Some attempts in the literature to increase the modelling capacity of the inverse exponential distribution include the works of Oguntunde, Adejumo, & Owoloko (Citation2017a, Citation2017b, Citation2017c, Citation2017d) and Singh & Goel (Citation2015); these works used different families of distributions to extend the inverse exponential distribution; a list of these families of distributions which include the beta generalized family of distribution (the first family of distribution developed from the logit of a random variable) and many others can be found in Alizadeh, Cordeiro, Bastos Pinho, & Ghosh (Citation2017); Cordeiro, Alizadeh, Nascimeto, & Rasekhi (Citation2016); Oguntunde, Adejumo, Okagbue, & Rastogi (Citation2016); Oguntunde, Khaleel, Ahmed, Adejumo, & Odetunmibi (Citation2017) and Owoloko, Oguntunde, & Adejumo (Citation2015).
The interest of the article is, however, on the Gompertz family of distributions because it is relatively new and flexible. Unlike the beta generalized family of distribution, it does not include complex functions and thus is simpler to work with. It has been used in earlier works by Alizadeh et al. (Citation2017) and Abdal-Hameed, Khaleel, Abdullah, Oguntunde, & Adejumo (Citation2018), and its densities are as follows:
and
where and
are extra shape parameters and it was developed using the following transformation:
is the probability density function (pdf) of the Gompertz distribution and t is a random variable.
and
are the cumulative distribution function (cdf) and pdf of the baseline distribution, respectively. In our case, the baseline distribution is the inverse exponential distribution defined as:
and
where is regarded as the scale parameter.
2. The Gompertz inverse exponential distribution
The cdf of the Gompertz inverse exponential (GoIE) distribution is derived by substituting Equation (4) into Equation (1) to have:
Also, the corresponding pdf is derived by substituting Equations (4) and (5) into Equation (2) to have:
where and
are shape parameters and
is the scale parameter
The pdf and cdf of the GoIE are represented in Figures and , respectively, using different parameter values.
2.1. Some basic statistical properties of the Gompertz inverse exponential distribution
Some of the basic statistical properties of the GoIE are obtained as follows:
First, the reliability function is obtained by using the relation:
So, the reliability function of the GoIE distribution is:
This is represented graphically in Figure .
The failure rate of the GoIE distribution is obtained by dividing the pdf in Equation (7) by the reliability function in Equation (8) to have:
This is represented graphically in Figure .
The reversed hazard function of the GoIE distribution is obtained by dividing the pdf in Equation (7) by the cdf in Equation (6) to have:
Also, the odds function is obtained by dividing the cdf in Equation (6) by the reliability function in Equation (8) to have:
2.2. Quantile function and median
The quantile function is obtained from the relation:
So, the quantile function of the GoIE distribution is derived as:
where . In other words, random samples from the GoIE distribution can be generated using:
The median of the GoIE distribution can conveniently be derived by making the substitution of in Equation (12) to have:
The corresponding first quartile and third quartile can also be obtained by making the substitution of and
, respectively, into Equation (12).
2.3. Order statistics
The pdf of the order statistic for a random sample of size
from a distribution function
and an associated pdf
is given by:
where and
are the pdf and cdf of the GoIE, respectively. The pdf of the jth order statistics for a random sample of size
from the GoIE distribution is, however, given as follows:
So, the pdf of minimum order statistics is obtained by substituting in Equation (15) to have:
While the corresponding pdf of maximum order statistics is obtained by making the substitution of in Equation (15) as:
2.4. Estimation
The method of maximum likelihood estimation is used to estimate the parameters of the GoIE distribution. For a random sample distributed according to the cdf of the GoIE distribution, the log-likelihood function is obtained as:
The log-likelihood function is obtained as:
Differentiate with respect to parameters
and
and equate the results to zero, and solving the resulting simultaneous equations gives the parameter estimates. The solution may not be obtained in closed form; so, software can be used to obtain the estimates numerically; example of such software include R, MAPLE and so on.
3. Model validation and application
The GoIE distribution is applied to two datasets and comparisons are made with the Gompertz exponential, Gompertz Weibull and Gompertz Lomax distributions. R software was used for the analysis and the criteria used for model selection are Akaike information criteria (AIC), consistent Akaike information criteria (CAIC), Bayesian information criteria (BIC), negative log-likelihood (NLL) and Hannan and Quinn information criteria (HQIC).
The criteria for selecting the distribution with the best fit depends on the values of the AIC, CAIC, BIC, NLL and HQIC, and lower values of this criteria indicate a better fit.
First Illustration: A data relating to the strengths of 1.5cm glass fibres which has been analysed previously by Oguntunde et al. (Citation2017) and Bourguignon, Silva, & Cordeiro (Citation2014), Merovci, Khaleel, Ibrahim, & Shitan (Citation2016) Smith & Naylor (Citation1987) was used. The observations are:
0.55, 0.74, 0.77, 0.81, 0.84, 1.24, 0.93, 1.04, 1.11, 1.13, 1.30, 1.25, 1.27, 1.28, 1.29, 1.48, 1.36, 1.39, 1.42, 1.48, 1.51, 1.49, 1.49, 1.50, 1.50, 1.55, 1.52, 1.53, 1.54, 1.55, 1.61, 1.58, 1.59, 1.60, 1.61, 1.63, 1.61, 1.61, 1.62, 1.62, 1.67, 1.64, 1.66, 1.66, 1.66, 1.70, 1.68, 1.68, 1.69, 1.70, 1.78, 1.73, 1.76, 1.76, 1.77, 1.89, 1.81, 1.82, 1.84, 1.84, 2.00, 2.01, 2.24
The result of the analysis is displayed in Table .
Table 1. Table of result
The GoIE is adjudged the distribution with the best fit because it has the lowest AIC, CAIC, BIC, NLL and HQIC values.
The histogram of the data with the competing distributions is displayed in Figure .
The empirical cdf of the competing distributions with respect to the dataset used is displayed in Figure .
Second Illustration: A sport data which is downloadable at www.stat.auckland.ac.nz/~lee/330/datasets.dir/sport.data was used and the observations are:
19.75, 21.30, 19.88, 23.66, 17.64, 15.58, 19.99, 22.43, 17.95, 15.07, 28.83, 18.08, 23.30, 17.71, 18.77, 19.83, 25.16, 18.04, 21.79, 22.25, 16.25, 16.38, 19.35, 19.20, 17.89, 12.20, 23.70, 24.69, 16.58, 21.47, 20.12, 17.51, 23.70, 22.39, 20.43, 11.29, 25.26, 19.39, 19.63, 23.11, 16.86, 21.32, 26.57, 17.93, 24.97, 22.62, 15.01, 18.14, 26.78, 17.22, 26.50, 23.01, 30.10, 13.93, 26.65, 35.52, 15.59, 19.61, 14.52, 11.47, 17.71, 18.48, 11.22, 13.61, 12.78, 11.85, 13.35, 11.77, 11.07, 21.30, 20.10, 24.88, 19.26, 19.51, 23.01, 8.07, 11.05, 12.39, 15.95, 9.91, 16.20, 9.02, 14.26, 10.48, 11.64, 12.16, 10.53, 10.15, 10.74, 20.86, 19.64, 17.07, 15.31, 11.07, 12.92, 8.45, 10.16, 12.55, 9.10, 13.46, 8.47, 7.68, 6.16, 8.56, 6.86, 9.40, 9.17, 8.54, 9.20, 11.72, 8.44, 7.19, 6.46, 9.00, 12.61, 9.03, 6.96, 10.05, 9.56, 9.36, 10.81, 8.61, 9.53, 7.42, 9.79, 8.97, 7.49, 11.95, 7.35, 7.16, 8.77, 9.56, 14.53, 8.51, 10.64, 7.06, 8.87, 7.88, 9.20, 7.19, 6.06, 5.63, 6.59, 9.50, 13.97, 11.66, 6.43, 6.99, 6.00, 6.56, 6.03, 6.33, 6.82, 6.20, 5.93, 5.80, 6.56, 6.76, 7.22, 8.51, 7.72, 19.94, 13.91, 6.10, 7.52, 9.56, 6.06, 7.35, 6.00, 6.92, 6.33, 5.90, 8.84, 8.94, 6.53, 9.40, 8.18, 17.41, 18.08, 9.86, 7.29, 18.72, 10.12, 19.17, 17.24, 9.89, 13.06, 8.84, 8.87, 14.69, 8.64, 14.98, 7.82, 8.97, 11.63, 13.49, 10.25, 11.79, 10.05, 8.51, 11.50, 6.26
The result of the analysis is displayed in Table .
Table 2. Table of result
The GoIE is also adjudged the distribution with the best fit because it has the lowest AIC, CAIC, BIC, NLL and HQIC values.
The histogram of the second data with the competing distributions is displayed in Figure .
The empirical cdf of the competing distributions with respect to the second dataset is displayed in Figure .
4. Discussion
The results from Tables and show that the GoIE has the lowest values for all the criteria used. So, the GoIE distribution is considered the best model to fit the datasets among the other distributions used. This is also evident in the plots provided in Figure –. In Figures and , the red curve which represents GoIE distribution has the highest peak compared to others. Also in Figures and , the red line is closer to the empirical cdf of the datasets than others.
5. Conclusion
The GoIE distribution has been successfully introduced in this paper and some of its basic statistical properties have been obtained. The pdf of the model and its failure rate have unimodal shapes; this means that the model would be useful to fit real-life events with unimodal failure rates. The model is tractable and flexible and shows high modelling capability as it performs better than the Gompertz exponential, Gompertz Weibull and Gompertz Lomax distributions; this was judged based on the AIC, CAIC, BIC, NLL and HQIC values of these distributions. The GoIE distribution is of no doubt a competitive model, and it is hoped that it would be of use in fields like engineering, biology and medicine. Some other statistical properties of the distribution (which were not considered in this paper) can be explored and simulation study can also be conducted.
Acknowledgements
The authors appreciate Covenant University for providing an enabling environment for this research. The efforts of the reviewers are greatly appreciated.
Additional information
Funding
Notes on contributors
Pelumi E. Oguntunde
Dr. Pelumi E. Oguntunde is a Statistician and a faculty in the Department of Mathematics, Covenant University, Nigeria. He has Bachelor of Science (B.Sc.) and Master of Science (M.Sc.) degrees in Statistics from University of Ilorin and University of Ibadan respectively. He earned his Ph.D. degree in Industrial Mathematics (Statistics Option) from Covenant University, Nigeria. His area of specialization is Mathematical Statistics. He teaches Statistics courses and he has made noticeable contributions to several scholarly journals especially in the area of probability distributions.
References
- Abdal-Hameed, M., Khaleel, M. A., Abdullah, Z. M., Oguntunde, P. E., & Adejumo, A. O. (2018). Parameter estimation and reliability, hazard functions of Gompertz Burr Type XII distribution. Tikrit Journal for Administration and Economics Sciences, 1(41 part 2), 381–400.
- Alizadeh, M., Cordeiro, G. M., Bastos Pinho, L. G., & Ghosh, I. (2017). The Gompertz-g family of distributions. Journal of Statistical Theory and Practice, 11(1), 179–207. doi:10.1080/15598608.2016.1267668
- Bakoban, R. A., & Abu-Zinadah, H. H. (2017). The beta generalized inverted exponential distribution with real data applications. REVSTAT-Statistical Journal, 15(1), 65–88.
- Bourguignon, M., Silva, R. B., & Cordeiro, G. M. (2014). The Weibull-G family of probability distributions. Journal of Data Science, 12, 53–68.
- Cordeiro, G. M., Alizadeh, M., Nascimeto, A. D. C., & Rasekhi, M. (2016). The exponentiated Gompertz generated family of distributions: Properties and applications. Chilean Journal of Statistics, 7(2), 29–50.
- Keller, A. Z., & Kamath, A. R. (1982). Reliability analysis of CNC machine tools. Reliability Engineering, 3, 449–473. doi:10.1016/0143-8174(82)90036-1
- Lin, C. T., Duran, B. S., & Lewis, T. O. (1989). Inverted gamma as life distribution. Microelectron Reliability, 29(4), 619–626. doi:10.1016/0026-2714(89)90352-1
- Merovci, F., Khaleel, M. A., Ibrahim, N. A., & Shitan, M. (2016). The beta type X distribution: Properties with applications. SpringerPlus, 5, 697. doi:10.1186/s40064-016-2271-9
- Oguntunde, P. E., & Adejumo, A. O. (2015). The transmuted inverse exponential distribution. International Journal of Advanced Statistics and Probability, 3(1), 1–7. doi:10.14419/ijasp.v3i1.3684
- Oguntunde, P. E., Adejumo, A. O., Okagbue, H. I., & Rastogi, M. K. (2016). Statistical properties and applications of a new Lindley exponential distribution. Gazi University Journal of Science, 29(4), 831–838.
- Oguntunde, P. E., Adejumo, A. O., & Owoloko, E. A. (2017a). Application of Kumaraswamy inverse exponential distribution to real lifetime data. International Journal of Applied Mathematics and Statistics, 56(5), 34–47.
- Oguntunde, P. E., Adejumo, A. O., & Owoloko, E. A. (2017b, July 5–7). On the flexibility of the transmuted inverse exponential distribution. Lecture Notes on engineering and computer science: Proceeding of the World Congress on Engineering (pp. 123–126). London, UK.
- Oguntunde, P. E., Adejumo, A. O., & Owoloko, E. A. (2017c, July 5–7). On the exponentiated generalized inverse exponential distribution. Lecture Notes on engineering and computer science: Proceeding of the World Congress on Engineering (pp. 80–83). London, UK.
- Oguntunde, P. E., Adejumo, A. O., & Owoloko, E. A. (2017d, July 5–7). The Weibull-inverted exponential distribution: A generalization of the inverse exponential distribution. Lecture Notes on engineering and computer science: Proceeding of the World Congress on Engineering (pp. 16–19). London, UK.
- Oguntunde, P. E., Babatunde, O. S., & Ogunmola, A. O. (2014). Theoretical analysis of the Kumaraswamy-inverse exponential distribution. International Journal of Statistics and Applications, 4(2), 113–116.
- Oguntunde, P. E., Khaleel, M. A., Ahmed, M. T., Adejumo, A. O., & Odetunmibi, O. A. 2017. A new generalization of the Lomax distribution with increasing, decreasing and constant failure rate. Article ID 6043169 Modelling and Simulation in Engineering, 2017, 6. doi:10.1155/2017/6043169
- Owoloko, E. A., Oguntunde, P. E., & Adejumo, A. O. (2015). Performance rating of the transmuted exponential distribution. An Analytical Approach, SpringerPlus, 4, 818.
- Singh, B., & Goel, R. (2015). The beta inverted exponential distribution: Properties and applications. International Journal of Applied Sciences and Mathematics, 2(5), 132–141.
- Smith, R. I., & Naylor, J. C. (1987). A comparison of maximum likelihood and Bayesian estimators for the three-parameter Weibull distribution. Applied Statistics, 36, 258–369. doi:10.2307/2347795