Research Article  Open Access
Kyeongjun Lee, JungIn Seo, "Different Approaches to Estimation of the Gompertz Distribution under the Progressive TypeII Censoring Scheme", Journal of Probability and Statistics, vol. 2020, Article ID 3541946, 7 pages, 2020. https://doi.org/10.1155/2020/3541946
Different Approaches to Estimation of the Gompertz Distribution under the Progressive TypeII Censoring Scheme
Abstract
This paper provides an estimation method for an unknown parameter by extending weighted leastsquared and pivotbased methods to the Gompertz distribution with the shape and scale parameters under the progressive TypeII censoring scheme, which induces a consistent estimator and an unbiased estimator of the scale parameter. In addition, a way to deal with a nuisance parameter is provided in the pivotbased approach. For evaluation and comparison, the Monte Carlo simulations are conducted, and real data are analyzed.
1. Introduction
The probability density function (PDF) and cumulative distribution function (CDF) of the random variable with the Gompertz distribution are, respectively, given bywhere is the shape parameter and is the scale parameter. This distribution was first introduced by Gompertz [1] and often used to model human mortality. For this distribution, some studies have been done on the inference method based on the pivotal quantity to estimate confidence intervals (CIs) of the parameters of interest more efficiently. Chen [2] provided exact CIs under the TypeII censoring scheme. Wu et al. [3] provided exact CIs under the progressive TypeII censoring scheme. Studies of various estimation methods for this distribution have also been conducted by several researchers. Dey et al. [4] proposed different methods to estimate the PDF and CDF and compared the estimation methods based on the Monte Carlo simulations. Dey et al. [5] provided various mathematical and statistical properties and compared different estimation methods from both frequentist and Bayesian point of view. Moala and Dey [6] provided Bayesian analysis methods under the objective and subjective priors including Jeffreys prior, maximal data information prior, Singpurwalla’s prior, and elicited prior.
In addition to this distribution, inferences based on the pivotal quantity have been studied for many distributions because the pivotalbased approach provides exact CIs even for small samples, as well as more efficient estimators than the maximum likelihood estimators (MLEs) in terms of bias. Wu [7] studied on a bathtubshaped lifetime distribution under the progressive TypeII censoring scheme. Wang et al. [8] provided a generalized pivotbased method based on the progressive TypeII censored data from the reverse hazard distributions. Under the same censoring scheme, Seo and Kang [9, 10] provided closed forms of exact CIs for the scale parameter for each with and without the nuisance parameter such as the location parameter in a half logistic distribution.
Recently, a new estimation method based on a weighted regression framework has been proposed under some censoring schemes. Lu and Tao [11] applied the method based on a regression framework to the Pareto distribution where there is no censoring. Seo et al. [12] extended the idea of Lu and Tao [11] to the progressive TypeII censoring scheme to estimate unknown parameters of the Pareto distribution. Seo et al. [13] provided closed forms for unknown parameters of a twoparameter Rayleigh distribution using the same approach under the progressive TypeII censoring scheme.
This paper focuses on point estimation using the weighted regression framework and pivotbased methods based on the progressive TypeII censored data from the Gompertz distribution with the PDF (1). The remainder of the paper is structured as follows. Section 2 provides the weighted regression framework and pivotbased estimation methods corresponding to the progressive TypeII censored data from the Gompertz distribution. Section 3 evaluates and validates the proposed method through the Monte Carlo simulations and analyzes real data for purposes of illustration. Section 4 concludes the paper.
2. Estimation
This section gives a brief description for the progressive TypeII censoring scheme that is the generalization of the TypeII censoring scheme and that is one of the most popular censoring schemes and provides different approaches on estimation for unknown parameters of the Gompertz distribution with the PDF (1). For brevity’s sake, the following notations are used throughout this paper: Exp (1): the standard exponential distribution : the chisquared distribution with degrees of freedom : the inverse gamma distribution with the shape parameter α and the scale parameter β : the uniform distribution on the interval (0, 1)
By Balakrishnan and Aggarwala [14], the progressive TypeII censoring scheme is described as follows. Let be the number of failures and denotes the censoring scheme, where is the number of removals at the th censoring time. In addition, suppose that the number of failures and the censoring scheme are fixed in advance. At the occurrence of the first failure, units are randomly withdrawn (or censored) from the surviving units. Subsequently, following the second observed failure, units are randomly censored from the surviving units, and the process continues in this way until the th failure is observed. Finally, at the time of the th observed failure, all remaining units are censored from the test. As mentioned earlier, this scheme includes as a special case the conventional TypeII right censoring scenario that is defined when . In addition, the scenarios and induce the complete sample situation.
Suppose that is a progressive TypeII censored sample with the censoring scheme from the Gompertz distribution. Then, the corresponding likelihood function is given bywhere
The MLEs and can be found by maximizing the logarithm of the likelihood function (2) for and . The following provides approaches based on the regression framework and pivotal quantity which lead to a consistent estimator and an unbiased estimator of for known , respectively.
Let
Then, are progressive TypeII censored order statistics from with the meanand the variancegiven by Theorem 7.2.1 by Balakrishnan and Cramer [15]. From that fact, consider the following linear regression model:where is the error term with the mean 0. For known , the regression model (7) provides a leastsquare estimator of as
By minimizing the following quantity for :
However, the approach gives same weight on each point, and it is not proper because the variances of do not satisfy the condition of being constant. As an alternative, Lu and Tao [11] considered weights that are inversely proportional to the corresponding variances. That is, a point with a low variance will be given a higher weight and a point with a higher variance will be given a low weight. Here, the weights corresponding to the progressive TypeII censoring scheme are defined as . Then, an estimator of based on the weighted regression framework is given by
By minimizing the following quantity with the weighted square term for :
Theorem 2.1. For known , the weighted squareleast estimator is a consistent estimator.
Proof. LetThen, the estimator can be written asHere, both and converge in the mean to 0 becausewhich implies convergence in probability [16]. In addition, does not converge to 0 as by the following inequality:Therefore, the fraction term in (13) converges in probability to 1, which completes the proof.
For unknown , its estimator can be derived by minimizing the quantity (11) with for and denoted as .
The pivotal quantity provided in Wu et al. [3], , can be used to derive another estimator of for known . According to Lemma 1 by Seo and Kang [9], converges to one in probability as because it has , and the estimator of can be obtained as from the equation . In addition, is an unbiased estimator because it has .
Note that the estimator depends on the parameter as in the case of . To deal with it, a pivotal quantity based on the work of Wang et al. [8] is provided in the following lemma.
Lemma 2.2. A quantityhas .
Proof. By Wang et al. [8], a quantity based on the quantity (3) is given byand it induces the pivotal quantity that is order statistics from . Then, from the fact that has , the following quantity:has . This completes the proof.
By Lemma 1 in the work of Seo and Kang [9] and Lemma 2.2, converges to one in probability as , which leads to an equation that has a unique solution for because the argument of the log term in (16) can be written asand the term in (19) is an increasing function of (resp. decreasing) for (resp. ). The unique solution is denoted as .
3. Application
This section assesses and compares the estimation methods provided in Section 2 through the Monte Carlo simulations and real data analysis.
3.1. Simulation Study
For evaluation and comparison, the mean squared errors (MSEs) and biases of the provided estimators are reported in Table 1. The progressive TypeII censored samples are generated from the Gompertz distribution with and under the following scenarios Scheme I: Scheme II: Scheme III: Scheme IV: by using the algorithm of Balakrishnan and Aggarwala [14]. All MSEs and biases are computed based on the generated 1,000 progressive TypeII censored datasets.

From Table 1, it can be seen that the weighted leastsquare estimator is more efficient compared to the MLE in terms of the bias, but the pivotbased estimator shows the best performance in terms of both the MSE and bias. Even for , the pivotbased estimator generally shows better results compared to the MLE and the weighted leastsquare estimator . Finally, MSEs of all estimators decreases as expected with the increasing number of failures for a fixed sample size.
3.2. Real Data
Chen [2] and Lee [17] analyzed a real data set that represents tumorfree time in days of 30 rats fed with unsaturated diet. For illustration purposes, a progressive TypeII censored data are generated from the tumorfree time data after dividing each data point by 100. The censoring scheme and corresponding censored data are presented in Table 2. Prior to analysis, the goodnessoffit test based on the first moment is first conducted. To avoid computational complex, the follow empirical first moment is employed:where are samples from the marginal distribution . The first moment (3.1) is evaluated at all the estimators provided in Section 2. Figure 1 reports the results through the boxscatter plot and the correlation coefficient (r) between the generated progressive TypeII censored data and the corresponding empirical first moments, which indicates that the assumption that the censored tumorfree time data have the Gompertz distribution which is reasonable for all estimates. The estimation results for and are reported in Table 3, which shows that the weighted leastsquare and pivotbased estimators have values similar to those of the MLEs under the considered schemes.

(a)
(b)
(c)

4. Conclusions
This paper provides approaches based on the weighted regression framework and pivotal quantity to estimate unknown parameters of the Gompertz distribution with the PDF (1) under the progressive TypeII censoring scheme. The proposed methods are comparatively concise and easy to perceive compared with the existing methods such as the maximum likelihood method. In addition, it was proved that the pivotbased estimators are superior to the MLEs and weighted leastsquare estimators in terms of the MSE and bias. Although the result for the goodnessoffit test does not have statistically significant differences, when considered with the simulation results, it is highly suggested the use of the estimation method based on the pivotal quantity for case where progressive TypeII censored data from the Gompertz distribution are observed.
Data Availability
The censored data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Research Foundation of Korea (NRF) Grant Funded by the Korea Government (Ministry of Education) (no. 2019R1I1A3A01062838).
References
 B. Gompertz, “On the nature of the function expressive of the law of human mortality and on the new mode of determining the value of life contingencies,” Philosophical Transactions of the Royal Society of London, vol. 115, pp. 513–580, 1825. View at: Google Scholar
 Z. Chen, “Parameter estimation of the Gompertz population,” Biometrical Journal, vol. 39, no. 1, pp. 117–124, 1997. View at: Publisher Site  Google Scholar
 S. J. Wu, C. T. Chang, and T. R. Tsai, “Point and interval estimations for the Gompertz distribution under progressive typeII censoring,” Metron, vol. 61, no. 3, pp. 403–418, 2003. View at: Google Scholar
 S. Dey, T. Kayal, and Y. M. Tripathi, “Evaluation and comparison of estimators in the Gompertz distribution,” Annals of Data Science, vol. 5, no. 2, pp. 235–258, 2018. View at: Publisher Site  Google Scholar
 S. Dey, F. A. Moala, and D. Kumar, “Statistical properties and different methods of estimation of Gompertz distribution with application,” Journal of Statistics and Management Systems, vol. 21, no. 5, pp. 839–876, 2018. View at: Publisher Site  Google Scholar
 F. A. Moala and S. Dey, “Objective and subjective prior distributions for the Gompertz distribution,” Anais da Academia Brasileira de Ciências, vol. 90, no. 3, pp. 2643–2661, 2018. View at: Publisher Site  Google Scholar
 S. J. Wu, “Estimation of the twoparameter bathtubshaped lifetime distribution with progressive censoring,” Journal of Applied Statistics, vol. 35, no. 10, pp. 1139–1150, 2008. View at: Publisher Site  Google Scholar
 B. X. Wang, K. Yu, and M. C. Jones, “Inference under progressively type II rightcensored sampling for certain lifetime distributions,” Technometrics, vol. 52, no. 4, pp. 453–460, 2010. View at: Publisher Site  Google Scholar
 J. I. Seo and S. B. Kang, “Pivotal inference for the scaled half logistic distribution based on progressively TypeII censored samples,” Statistics & Probability Letters, vol. 104, pp. 109–116, 2015. View at: Publisher Site  Google Scholar
 J. I. Seo and S. B. Kang, “Inference for the twoparameter halflogistic distribution using pivotal quantities under progressively typeII censoring schemes,” Communications in Statistics—Simulation and Computation, vol. 46, no. 7, pp. 5462–5478, 2017. View at: Publisher Site  Google Scholar
 H. L. Lu and S. H. Tao, “The estimation of Pareto distribution by a weighted least square method,” Quality & Quantity, vol. 41, no. 6, pp. 913–926, 2007. View at: Publisher Site  Google Scholar
 J. I. Seo, S. B. Kang, and H. Y. Kim, “New approach for analysis of progressive typeII censored data from the Pareto distribution,” Communications for Statistical Applications and Methods, vol. 25, no. 5, pp. 569–575, 2018. View at: Publisher Site  Google Scholar
 J. I. Seo, B. G. Seo, and S. B. Kang, “Estimation on a twoparameter Rayleigh distribution under the progressive typeII censoring scheme: comparative study,” Communications for Statistical Applications and Methods, vol. 26, no. 2, pp. 91–102, 2019. View at: Publisher Site  Google Scholar
 N. Balakrishnan and R. Aggarwala, Progressive Censoring—Theory, Methods, and Applications, Birkhäuser, Boston, MA, USA, 2000.
 N. Balakrishnan and E. Caramer, The Art of Progressive Censoring, Springer, New York, NY, USA, 2014.
 A. F. Karr, Probability, SpringerVerlag, New York, NY, USA, 1993.
 E. T. Lee, Statistical Methods for Survival Data Analysis, Wiley, New York, NY, USA, 2nd edition, 1992.
Copyright
Copyright © 2020 Kyeongjun Lee and JungIn Seo. 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.