Admissibility Estimation of Pareto Distribution Under Entropy Loss Function Based on Progressive Type-II Censored Sample
Guobing Fan
Department of Basic Subjects, Hunan University of Finance and Economics, Changsha, China
Email address:
To cite this article:
Guobing Fan. Admissibility Estimation of Pareto Distribution Under Entropy Loss Function Based on Progressive Type-II Censored Sample. Pure and Applied Mathematics Journal. Vol. 5, No. 6, 2016, pp. 186-191. doi: 10.11648/j.pamj.20160506.13
Received: October 5, 2016; Accepted: October 14, 2016; Published: November 7, 2016
Abstract: The aim of this paper is to study the estimation of Pareto distribution on the basis of progressive type-II censored sample. First, the maximum likelihood estimator (MLE) is derived. Then the Bayes estimator of the unknown parameter of Pareto distribution is derived on the basis of Gamma prior distribution under entropy loss function. Further the empirical Bayes estimator also obtained by using maximum likelihood on the basis of Bayes estimator. Finally, the admissibility of a class of inverse linear estimators are discussed under suitable conditions.
Keywords: Admissibility, Bayes and Empirical Bayes Estimators, Progressive Type-II Censored Sample, Entropy Loss Function
1. Introduction
In practical life testing experiments, considering with the time limitation and/or other restrictions (such as cost, material resources, etc.), the experimenter may not always be in a position to observe the life times of all the products in a lifetime test [1, 2]. Censored samples often arise in practice. Progressive Type II censored sampling test is one of the most important method of obtaining data in lifetime researches. The statistical inference studies for various distributions when sample belongs to progressive censoring have attracted many authors’ attention. For example, Ng et al. [3] computed the expected Fisher information and the asymptotic variance-covariance matrix of the ML estimates based on progressively Type-II censored sample from Weibull distribution They also discussed the construction of progressively censored reliability sampling plans. Soliman et al. [4] investigated the point and interval estimations for the modified Weibull distribution based on progressively type-II censored sample. Yang [5] derived the maximum likelihood estimation of Weibull distribution under Type II progressive censoring with random removals, where the number of units removed at each failure time follows a binomial distribution.Wu [6] considered the estimation problem of the two-parameter bathtub-shaped lifetime distribution based a progressively type-II censored sample. Cho et al. [7] discussed the Bayes estimation of the entropy of a two-parameter Weibull distribution based on the generalized progressively censored sample. Bhattacharya et al.[8] proposed an optimum life-testing plans under Type-II progressive censoring scheme using variable neighborhood search algorithm. Laumen and Cramer [9] discussed the likelihood inference and statistical test procedure for the lifetime performance index in the presence of progressive censoring. Khorram and Farahani [10] considered the maximum likelihood estimation and bayes estimation of parameters of weighted exponential distribution based on progressively Type-II censored sample.
Many statistical inference problems have been discussed for various lifetime distributions, such as exponential distribution, Weibull distribution, etc. Wu and Chang [11] pointed out that Pareto distribution can be regarded as a suitable alternative distribution in modeling product’s lifetime. Statistical inference about Pareto distribution receives great attention by authors in recent yeasrs. For example, Raqab et al. [12] derived the best linear unbiased predictors, maximum likelihood predictors and approximate maximum likelihood predictors of times to failure of units censored from Pareto distribution. Kulldorff and Vannman [13] obtained the best linear unbiased estimates based on the complete sample and the asymptotically best linear unbiased estimates based on a few selected order statistics. Fu et al. [14] discussed the Bayesian estimation of Pareto distributions under progressive Type-II censoring on the basis of several types of noninformative priors, i.e. Jeffreys prior, two reference priors and two general forms of second order probability matching prior. Saldaña-Zepeda et al. [15] proposed a goodness of fit test procedure for the Pareto distribution, when the observations drawn from Type II right censoring.
Assume that the repair time follows the Pareto distribution with the following probability density function (pdf) and cumulative distribution function (cdf) respectively:
(1)
(2)
This paper will discuss the Bayes estimation of the parameter of Pareto distribution under entropy loss function on the basis of progressive type-II censored sample. The admissibility of estimators is an important topic [16-19]. Thus this paper will also study the admissibility and inadmissibility of a class of inverse linear estimators under suitable conditions. The remains of this paper are organized as follows. The MLE and Bayesian estimators of the parameter are obtained in Section 2. In Section 3, the admissibility and inadmissibility of estimators with inverse linear form are discussed. A conclusion is finally made in Section 4.
2. Preliminary Knowledge
The progressively Type II censoring scheme can be described as follows [20].
First, the experimenter places n units or individual on test.
(i) When the first failure is observed at , then randomly select surviving unites and remove them.
(ii) When the i-th failure unit is observed at , then randomly select surviving unites and remove them. i = 2, 3,...,m.
(iii) This experiment terminates when the m-th failure unit is observed at and of surviving units are all removed. Here the censoring numers are all pre-fixed.
Remark 1. Note that if , , This censored scheme reduces to conventional type II right censoring scheme. Also note that if ,, the progressively type II censoring scheme reduces to the complete sample case.
This paper we always assume that is a progressively Type-II censored sample from a life test on n items whose lifetimes follows Pareto distribution with pdf shown in (1), is the corresponding observation of and denote the corresponding numbers of units removed from the test.
2.1. Maximum Likelihood Estimation
The likelihood function of under given progressively type-II censored sample is (Balakrishnan and Aggrwala [21]):
(3)
Here is the observation of ,, andare given respectively by (1) and (2), and.
Substituting (1) and (2) into (3), the likelihood function is given by
(4)
Where is the observation of ,
The natural logarithm of likelihood function is given by
(5)
Then the MLE of can is the solution of the following equation
(6)
Thus the MLE of can be easily solved as
(7)
2.2. Bayes Estimation
This subsection will discuss the Bayes estimation of the parameter of Pareto distribution (1) based on progressively Type-II censored sample under the following entropy loss function:
(8)
Lemma 1. Let be a progressively type-II censored sample from a life test on items whose lifetimes have Pareto distribution (1), and denote the corresponding numbers of units removed from the test. Then under the entropy loss function (8), the unique Bayes estimator of is (Wang et al. [22]):
(9)
Suppose the prior distribution of is Gamma distributions, , with the following pdf
(10)
where and .
It is easy to proved that the posterior distribution of is also a Gamma distribution , i.e. ,where and from Lemma 1, the unique Bayes estimator can be obtained as
(11)
Where
The Bayes estimator can be rewritten as
(12)
which has of the form
Remark 2 The Bayes estimator in Eq. (11) is seen to depend on the prior parameter . Then When the parameter is unknown, we can use the empirical Bayes approach to estimate it. From (1) and (10), the marginal pdf of can be derived as follows
Then based on , we can get MLE ofas
.
Substituting for in the Bayes estimator (12), the empirical Bayes estimator of can obtained
(13)
3. Admissibility and Inadmissibility of (cT+d)^{-1}
In the following discussion, we always supposes that is a progressively Type-II censored sample from a life test on items whose lifetimes have Pareto distribution (1), and denote the corresponding numbers of units removed from the test.
Note that the MLE, Bayes estimator and empirical estimators obtained in Section 2 are all the special cases of a class of inverse linear estimators of the form . In the rest of this section, the admissibility of these estimators will be discussed on the basis of their risks under the entropy loss function (8). Let ,.
For later use, we need the following result shown in Lemma 2.
Lemma 2. The statistics has Gamma distribution .
Proof. Let ,, then we can easily proved that is a progressively type II censored sample distributed with the standard exponential distribution. Considering with the following transformations
(14)
Viveros & Balakrishnan [23] showed that the generalized spacings are all independent and identically distributed (i.i.d.) as standard exponential with mean 1.
Then we have , that is
Then it is easy to prove that the distribution of statistic is .
Theorem 1 The inverse linear estimator is admissible, provided and .
Proof. From the Bayes estimator in Eq. (12), we see that when and , the coefficient of is between 0 and , and the constant is strictly bigger than 0. This proves that is admissible for the case and . For the case , the estimator is admissible since it is the unique estimator for which when .
Theorem 2 The inverse linear estimator is admissible, provided and .
Proof. We first consider the case and . Let has the prior distribution with pdf
If is a nondegenerate convex subset of , then it is obvious that there exists a such that for some and all .
In fact,
When , we have.
Then
Let ,
Then the conclusion has been proved.
The Bayes estimator with respect to prior distribution under the entropy loss function (8) can be derived as in (12) as
Let is the jiont probability density function of .
Then
The difference of the Bayes risks with respect to and is
Here , and ,
For any prior Gamma distribution of , we have
Thus ，
Therefore, is admissible for any by Blyth's lemma. This proves that the inverse linear estimatoris admissible for and .
For the case and, is the limit of Bayes estimator relative to the gamma prior , as. It is easy to verify that the Bayes risk difference converges to zero as. Then we prove the admissibility of .
Theorem 3 Let the parameter space beand the action space be. The inverse linear estimator is inadmissible under the entropy loss function (8) whenever one of the following conditions holds:
(i) or ;
(ii)
(iii) and
Proof. For the case (i), because takes on negative values with positive probability. Therefore is dominated by the estimator . Then (i) is proved.
For the case (ii), the risk function of the estimatoris
Using Lemma 2, we have
,,
where .
The derivative of the risk function with respect to is
Then where , and, where.
Thus the risk function of is minimized at the value . Hence the estimator is dominated by the estimator . Thus (ii) is proved.
For the case (iii). Because
Under the condition (iii), , then
,
Therefore, is minimized at. Hence is dominated by in this case.
Remark 3. The MLE and Bayes estimator are inadmissible by using Theorem 3. They are both dominated by the generalized Bayes estimator .
4. Conclusions
This paper considers the estimation of the unknown parameter of Pareto distribution based on progressively type II censored samples. The MLE, Bayes estimator and empirical Bayes estimators are obtained. These estimators all belong to a class of inverse linear estimators with the form , where . The admissibility and inadmissibility of are discussed. As a result, the MLE and empirical Bayes estimator are inadmissible.
Acknowledgement
This study is partially supported by Natural Science Foundation of Hunan Province (No. 2016JJ4012 and No. 2015JJ3030) and Foundation of Hunan Educational Committee (No.15C0228).
References