|Intelligent Automation & Soft Computing |
Parameter Estimation of Alpha Power Inverted Topp-Leone Distribution with Applications
1High Institute for Management Sciences, Belqas, 35511, Egypt
2Faculty of Graduate Studies for Statistical Research, Cairo University, Giza, 12613, Egypt
3Faculty of Business Administration, Delta University for Science and Technology, Mansoura, 11152, Egypt
4Faculty of Commerce, Mansoura University, Mansoura, 35516, Egypt
*Corresponding Author: Ehab M. Almetwally. Email: firstname.lastname@example.org
Received: 03 February 2021; Accepted: 21 March 2021
Abstract: We introduce a new two-parameter lifetime model, referred to alpha power transformed inverted Topp-Leone, derived by combining the alpha power transformation-G family with the inverted Topp-Leone distribution. Structural properties of the proposed distribution are implemented like; quantile function, residual and reversed residual life, Rényi entropy measure, moments and incomplete moments. The maximum likelihood, weighted least squares, maximum product of spacing, and Bayesian methods of estimation are considered. A simulation study is worked out to evaluate the restricted sample properties of the proposed distribution. Numerical results showed that the Bayesian estimates give more accurate results than the corresponding other estimates in the majority of the situations. The flexibility of the suggested model is demonstrated given some applications related to reliability, medicine, and engineering. A real data set is used to illustrate the potentiality of the alpha power transformed inverted Topp-Leone distribution compared to inverted Topp-Leone, inverse Weibull, alpha power inverse Weibull, inverse Lomax, alpha power inverse Lomax, inverse exponential, and alpha power exponential distributions. Criteria measures and their results showed that the suggested distribution is the best candidate for the considered data sets. The alpha power transformed inverted Topp-Leone distribution operates well for lifetime modeling.
Keywords: Inverted Topp-Leone; moments; maximum likelihood; maximum product spacing; weighted least squares; Bayesian estimation; MCMC
In recent times, probability distributions play a significant role in modeling naturally occurring phenomena. In fact, the statistics literature contains hundreds of continuous univariate distributions and their successful applications. However, there still remain many real-world phenomena involving data, which do not follow any of the traditional probability distributions. So, several attempts are introduced by many researchers to provide more flexibility to a family of distributions. Mahdavi et al.  introduced the alpha power transformation (AP) technique by adding an extra shape parameter to well- known baseline distributions. The cumulative distribution function (CDF) of the AP method is defined by:
Relevant works have been provided based on the AP method, for instance; AP Weibull distribution , AP generalized exponential distribution , AP extended exponential distribution , AP Lindley distribution , AP power Lindley distribution , AP inverse Lindley distribution , AP exponentiated Lomax distribution  and AP inverse Lomax distribution .
The inverted distributions were suggested in the literature using the inverse transformation of probability distributions. These distributions display different features in the behavior of the density and hazard rate shapes. They allow applicability to the phenomenon in many fields such as; biological sciences, life testing problems, survey sampling, and engineering sciences. Inverted distributions and their applications were discussed by several authors (see [10–18]).
In Hassan et al. , the CDF and the probability density function (PDF) of the inverted Topp-Leone (ITL) distribution with shape parameter is defined as follows:
In this paper, we propose a new two-parameter related to the ITL distribution depending on the AP family. We call it alpha power inverted Topp Leone (APITL) distribution. The basic motivations to introduce the APITL model are (i) Generalizing a new useful version of the ITL distribution based on the APT method along with deriving its statistical properties, (ii) Providing flexible PDF with right-skewed and uni-modal shapes, (iii) Modeling decreasing, increasing, upside-down hazard rate function (HRF), and (iv) Introducing some real applications in some areas.
This paper is constructed as follows. Section 2 describes the APITL distribution. Section 3 gives some structural properties of the APITL distribution. Section 4 gives the maximum likelihood (ML), the weighted least squares (WLS), the maximum product of spacing (MPS), and Bayesian estimators. Section 5 examines the effectiveness of the proposed estimates through a numerical illustration. Data analyses and some concluding remarks are employed, consequently, in Sections 6 and 7.
2 Description of the Model
In this section, based on the AP family we introduce a new probability distribution related to the ITL distribution. We define the PDF, CDF, HRF and cumulative HRF of the APITL distribution.
A random variable X is said to have the APITL distribution when we substitute the CDF (Eq. (2)) and PDF (Eq. (3)) in CDF (Eq. (1)). The CDF of a random variable X has the APITL distribution with parameters and denoted by X ~ APITL is defined by:
The PDF related to Eq. (4) is given by:
Descriptive PDF plots of the APITL distribution for some choices of parameters are represented in Fig. 1. It can be seen that the PDF of APITL distribution is uni-modal as well as possesses a long tail right-skewed.
The reliability function and the HRF of X are given by:
An illustration of the HRF plots for the APITL distribution, for some choices of and is represented in Fig. 2. It describes the HRF plots of the APITL distribution which can be decreasing, increasing, and upside down shaped.
3 Distributional Properties
Here, we give some statistical properties.
3.1 Quantile Function
The APITL distribution is simulated by inverting CDF Eq. (4) as follows:
The uth quantile for the APITL random variable is obtained by solving F(x) = u for x as follows:
where, The 25th, 50th, and 75th percentiles for the random variable X is obtained by setting u = 0.25, 0.5 and 0.75 in Eq. (9). The Bowley’s skewness, depends on quartiles, is defined as follows:
where Q(.) is the APITL quantile function. The Moor’s kurtosis is given as:
Skewness and kurtosis plots of the APITL model, based on the quantiles, are exhibited in Fig. 3.
Moments of a probability distribution are crucial to deduce its properties such as measures of central tendency, dispersion, skewness, and kurtosis. The ordinary rth moment of the APITL distribution is derived. The rth moment of the APITL distribution is obtained from Eq. (5) as follows:
Since the power series representation is written as:
Using the power series expansion Eq. (13) in Eq. (12), then we get
Using the binomial expansion in Eq. (14) then takes the form
After using binomial expansion, then the rth moment is given by:
The first four moments, for r = 1, 2, 3 and 4, of the APITL distribution are obtained from Eq. (16). The rth central moment of X is given by:
Values of mean variance skewness and kurtosis of the APITL distribution for certain values of parameters are given in Tab. 1.
We concluded from Tab. 1 that, as the value of increases and for fixed value of then values of and are decreasing. As the value of increases and for fixed value of , then values of and are increasing, while the values of and decrease. The distribution is positively skewed and leptokurtic.
3.3 The Probability Weighted Moments
The class of probability-weighted moments (PWMs), denoted by s and q are positive integers, of a random variable X is defined as follows:
Substituting Eq. (4), and Eq. (5) in Eq. (18), further using expansion Eq. (13) and binomial expansion in Eq. (18), then the PWM of the APITL distribution is obtained as follows:
3.4 Residual and Reversed Residual Life
The mth moment of the residual life (RLe), say uniquely determines F(x). The mth moment of the residual life of X is defined by:
Therefore, the mth moment of RLe for APITL distribution is obtained by substituting PDF Eq. (5) in Eq. (20), then employing binomial expansions and Eq. (13) as follows:
where and is the incomplete beta function. For m = 1, in Eq. (21) we obtain the mean RLe of the APITL distribution. Also, mth moment of the reversed RLe (RRLe) for the APITL distribution is given by:
The mean of RRLe serves as the waiting time elapsed since the failure of an item on the condition that this failure had occurred.
3.5 Rényi Entropy
The entropy of a random variable measures the amounts of information (or uncertainty) contained in a random observation; i.e., large value of entropy indicates higher uncertainty in the data. Rényi entropy of X, for and is defined by:
Substituting PDF Eq. (5) in Eq. (23), we obtain Rényi entropy of the APTIL model as follows:
From Eq. (24), and after simplification, the Rényi entropy of the APITL distribution will be:
4 Parameter Estimation of the APITL Model
In this section, we deal with parameter estimators of the APITL distribution based on ML, WLS, MPS, and Bayesian estimation methods.
4.1 ML Estimators
Let X1,…, Xn e the observed values from the APITL distribution with parameters and The likelihood function, say of the APITL distribution is expressed as:
Then the log-likelihood function, say of the APITL distribution is given as:
Therefore, the ML equations are given by:
Solving the non-linear equations and numerically using optimization algorism as conjugate-gradient optimization, we get the ML estimators of and
4.2 WLS Estimators
Let X1,…, Xn be a simple random sample from the APITL distribution and let X(1)< X(2)<…<X(n) be the associated order statistics. The WLS estimators of and are attained by minimizing the following:
Furthermore, the WLS estimators followed by solving the following nonlinear equations:
4.3 MPS Estimators
The MPS method is an alternative procedure of the ML method which provides a parameter estimate of continuous distribution. The MPS estimators of and are attained by maximizing the following:
Solving the non-linear equations and via numerical technique, we obtain the MPS estimators of and To obtain the MPS estimators, we differentiate natural logarithm of the product spacing function of APITL distribution partially for and
Equating Eqs. (33) and (34) with zero and using optimization algorism (conjugate-gradient or Newton-Raphson optimization) we get the solution.
4.4 Bayesian Estimators
Here, we get the Bayesian estimator of the APITL parameters. The Bayesian estimator is regarded under symmetric (squared error loss function (SELF)) which is defined as Also, the Bayesian estimator is considered under asymmetric (Linear exponential (LINEX) loss function) which is expressed as where h reflects the direction and degree of asymmetry. Assuming that the prior distribution of and , denoted by has an independent gamma distribution. The joint gamma prior density of and can be written as:
To elicit the hyper-parameters of the informative priors, the ML estimator for and is obtained by equating the estimates and their variances by the inverse of Fisher information matrix of and
From Eqs. (26) and (35), the joint posterior of the APITL distribution with parameters and is
To obtain the Bayesian estimators, we can use the Markov Chain Monte Carlo (MCMC) approach. A useful sub-class of the MCMC techniques is the Gibbs sampling and more general Metropolis within Gibbs samplers. The Metropolis-Hastings (MH) algorithm along with the Gibbs sampling are the two most popular examples of the MCMC method. We use the MH within Gibbs sampling steps to generate random samples from conditional posterior densities of α and λ as follows:
The Bayesian estimators are obtained via SELF and LINEX loss function (for more information see [19–23]).
5 Simulation Study
A Monte-Carlo simulation study was conducted to evaluate and compare the behavior of the different estimates based on mean square errors (MSEs) and biases. Generate 10000 random samples of sizes n = 50, 100, 150 and 200 from APITL distribution. Different actual parameter values were considered.
We calculated the ML estimate (MLE), WLS estimate (WLSE), MPS estimate (MPSE), and Bayesian estimate (BE) of and . Then, the biases and MSEs of the different estimates were determined. Simulated results were scheduled in Tabs. 2 and 3 and we noticed the following:
1. The bias and MSE for all estimates decrease as n increases (see Tabs. 2 and 3).
2. As values of near to one and for fixed values, the biases and MSEs of and estimates increase.
3. For a fixed value of as well as the value of increases, the biases and MSEs of estimates are increasing, in approximately most of the cases.
4. For a fixed value of as well as the value of increases, the biases, and MSEs for all estimates increase, in approximately most of the situations.
5. The measures of WLSEs are better than MLEs and MPSEs with decreasing sample size.
6. The measures of MPSEs are preferable to MLEs and WLSEs with sample sizes.
7. The BEs under the LINEX loss function is preferable to the other estimates.
6 Real Data Illustration
Here, we fit the APITL distribution under three real data taken from fields of survival times of medicine, engineering, and reliability. The APITL model is compared to other competitive models as, ITL, inverse Weibull (IW), alpha power IW (APIW), inverse Lomax (ILo), alpha power ILo (APILo), inverse exponential (IEx), and alpha power inverse exponential (APIEx) distributions.
Tabs. 4–6 provide values of Akaike information criterion (AIC), corrected AIC (CAIC), Hannan-Quinn information criterion (HQIC), and Kolmogorov- Smirnov (KS) statistic along with its P-value for all fitted models for three real data. In addition, these tables contain the MLEs and standard errors (SEs) (appear in parentheses) of the parameters for the considered models. We compared the fits of the APITL model with the ITL, IW, APIW, ILo, APILo, IEx, and APIEx models (see Tabs. 4–6). The fitted APITL PDF, CDF, PP-plot and QQ-plot of the three real data were displayed in Figs. 4–6, respectively. These figures indicated that the APITL distribution has the smallest values of AIC, CAIC, HQIC, KS and the largest P-value among all fitted models.
The first (I) set of data was studied in Bjerkedal . It represents the rth survival times (in days) of 72 guinea pigs infected with virulent tubercle bacilli. These data were analyzed in Refs. [25,26]. Tab. 4 listed values of MLEs and statistic measures for data I. Fig. 4 provided estimated PDF, CDF, PP-plot, and QQ-plot of APITL distribution for data I.
The second (II) set of data represents the active repair times (hr) for an airborne communication transceiver . The data are as follows: 0.50, 0.60, 0.60, 0.70, 0.70, 0.70, 0.80, 0.80, 1.00, 1.00, 1.00, 1.00, 1.10, 1.30, 1.50, 1.50, 1.50, 1.50, 2.00, 2.00, 2.20, 2.50, 2.70, 3.00, 3.00, 3.30, 4.00, 4.00, 4.50, 4.70, 5.00, 5.40, 5.40, 7.00, 7.50, 8.80, 9.00, 10.20, 22.00 and 24.50. The MLEs and statistic measures for data II listed in Tab. 5. Fig. 5 provided estimated PDF, CDF, PP-plot, and QQ-plot of the APITL model for data II.
The third (III) set of data was studied in Aarset . It refers to 30 failure times of air-conditioning system of an airplane. Tab. 6 listed values of MLEs and statistic measures for data III. Fig. 6 provided estimated PDF, CDF, PP-plot and QQ-plot of the APITL for data III.
According to tables and Figs. 4–6, we observed that the APITL distribution provides well overall the fitted model and consequently could be selected as the more suitable model than other models.
Furthermore, the suggested methods of estimation (see Section 4) for the APITL parameters were considered based on the three data. Tab. 7 displayed different estimates of the APITL parameters for the three data sets. In these data, we cannot use the MPS method because there are equal observations in the data and consequently the spacing followed by the product will be zeros. For more information about this method see [29,30].
The convergence of the MCMC estimation of and are shown in Figs. 7–9 for all datasets, respectively.
7 Concluding Remarks
We proposed and studied the alpha power transformed inverted Topp-Leone distribution. Some structural properties of the APITL distribution were provided. Bayesian and non-Bayesian methods of estimation were considered. We obtained the ML, WLS, and MPS estimators of the population parameters. The Bayesian estimator was deduced under LINEX and SELF. The Monte Carlo simulation study was worked out to assess the behavior of estimates. Generally, we concluded that the Bayes estimates are preferable to the corresponding other estimates in approximately most of the situations. We proved empirically that the APITL model reveals its superiority over other competitive models for different real data.
Acknowledgement: Thank a lot for every one help full.
Funding Statement: The authors received no specific funding for this study.
Conflicts of Interest: The authors declare that no conflicts of interest to report regarding the present study.
|This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.|