Abstract
In a life testing experiment, the successive failure times at putting units under test are recorded under consideration that the advance fixed experiment is terminated at time T. This type of censoring scheme, called conventional Type-I censoring scheme, is used widely. In this paper, we assume that the items’ failure times are independent and distributed with exponential lifetime distribution with parameter . The estimator with maximum likelihood method is obtained in an exact form and its distribution is also obtained with the unknown parameter. We proposed the exact confidence interval, for estimators and , asymptotic confidence intervals, confidence interval under likelihood ratio test, and finally, two bootstrap confidence intervals. Under the Bayesian approach, the unknown parameter is estimated and the corresponding credible interval is obtained considering the prior information formulated with the inverted gamma distribution. The Monte Carlo simulation study is used to compare different methods. Finally, for illustrative purposes, the real dataset is used and analyzed.
1. Introduction
Suppose that identical items are placed under a life testing experiment. Then, independent and identically distributed (i.i.d.) lifetime of each item is distributed with exponential lifetime distribution with parameter . Hence, the lifetime random variable has the probability density function (PDF) with parameter presented by
The test continues to advance until fixed time is reached to terminate the experiment. The experiment is done under consideration that the failed items through the test are not replaced. This type of censoring scheme is known as conventional Type-I censoring scheme which is quite important in reliability acceptance test in MIL-STD-781C [1].
The Type-I censoring scheme was introduced by Epstein [2], and the estimator of two-sided confidence interval for the parameter was proposed by Epstein [3] without any formal proof. The proposition presented by Epstein was modified slightly by Fairbanks et al. [4] to present a simple new set of confidence intervals. The maximum likelihood estimator (MLE) of the parameter and the corresponding exact distribution of the estimator as well as one-sided confidence interval were obtained by Chen and Bhattacharya [5]. Under consideration that, Bayesian approach with inverted gamma prior distribution, the Bayesian point and corresponding two-sided credible interval of the mean lifetime estimator is considered by Draper and Guttman [6]. Recently, the maximum likelihood estimator under right-censored linear regression data is presnted by Yu [7] and for mean ranked acceptance sampling plan under exponential distribution is presented by Hussain et al. [8]. The problem of building two-sided confidence interval with respect to the exact distribution of was presented by Chen and Bhattacharya [5].
The problem of the exact distribution of under Type-I censoring scheme is more serious to obtain. And, several books on reliability or survival analysis such as, Bain [9], Barlow and Proschan, [10] and Lawless [11] do not talk about this problem. It might be worth mentioning that, specialy as an application in the analysis repairable system, the problem of random sampling from a truncated exponential distribution. This is discussed in detail given by Bain and Engelhardt [12] (Chapter 9) as well as distribution of the sum of the observations taken from a truncated exponential distribution. As the same conditional density of the sum of observed failure times given number of failures in a time-censored Weibull process, it turns out in the repairable systems framework. In natural setting this density is given of sampling from the truncated exponential distribution see, Bain, Engelhardt and Wright [13], see also Bain and Weeks [14] in this respect.
In statistics, we are using a sample of data to estimate an interval of interested parameter value. The interval estimation in the paper is discussed in different forms. Under frequentist method, different confidence intervals are proposed, but under Bayesian approach, the credible interval is proposed. Also, there is a common form of interval estimation, such as fiducial intervals, tolerance intervals, and prediction intervals. In the literature, there are different novel situations as a guide on how interval estimates are formulated. Both credible intervals and confidence intervals are different but have a similar standing. Confidence intervals can be applied in more situations in parametric and non-parametric models other than credible interval. The problem of testing the performance of interval estimation procedures involves approximations of various kinds, and there is a need to check that the actual performance of a procedure is close to what is claimed. The coverage probability and interval length are the key concepts associated with interval estimators. We know that a higher coverage is obtained under longer interval length and lower coverage is obtained under shorter interval length. In statistics, we face a number of interval estimators of the population parameters, and a decision has to be made on what the “best” method of estimation is. The key concepts are that the better coverage probability goes with smaller interval length and vice versa, and it is useful to have some practical way of combining these measures.
Under consideration the Type-I censoring scheme, the exact and asymptotic distributions of are obtained. Also, based on the two values and , the exact and the corresponding asymptotic confidence intervals of are derived. Also, based on the likelihood ratio test (LRT), the confidence interval, the Bayes credible interval for inverted gamma prior, and the exact expressions for and are derived from considering the unknown parameter. Finally, different methods are compared with Monte Carlo simulations, and the real-life dataset is applied to all different methods.
This paper is organized as follows. The problem and the MLE of the unknown parameter with its distribution are formulated in Section 2. Confidence intervals as well as credible intervals are presented in Section 3. Different methods are compared through the numerical experiments in Section 4. The real dataset is analyzed in Section 5. Finally, in Section 6, we draw the conclusions from our work.
2. The MLE and Its Distribution
In this section, we constructed the exact conditional moment generating function that can be used to build exact conditional confidence intervals of the parameter under the assumption of exponential distribution.
2.1. MLE
Suppose that denote the lifetimes of the identical items put under test, where is distributed with exponential distribution with parameter and PDF given in (1). The ordered lifetimes of these items is denoted by . Suppose that, the number of items that fails upto and including the prefixed time point is denoted by . So, the observations under conventional Type-I censoring scheme are defined bywhere
Let us definewhere defines the total time on test. Therefore, based on the exponential Type-I censored data, the likelihood function of can be written as
From the likelihood function (5) and for , the MLE of exists and is given byand that of does not exist for .
2.2. Distribution of the MLE
Here, the distribution of show by the conditional distribution of for and in other words obtain the distribution of given , and , then unconditional distribution is obtained. It can be done, as given by Hoem [15] by using the convolution property of the uniform distribution or as suggested by Chen and Bhattacharya [5] by the conditional moment generating function (MGF) approach. Here we are using the MGF approach. The conditional MGF of given is given bywhere and . The proof of (7) is given in Appendix A.
Note that, the conditional distribution of the sum of observations given () is conditional distributed as Type-I censoring observation taken from one-parameter exponential distribution is the same as the sum of observations in a random sample of size , taken from truncated exponential distribution Bain [9] or Chen and Bhattacharya [5]. Then, the MLE distributed as a function obtained by taking the expectation relative to (truncated above 0) of the conditional distribution of given .
Now we are going to calculate the value of the mean and variance of for Type-I censoring scheme. It is clear that the problem of obtaining the first two moments may be complicated. However, from the conditional MGF, they can be obtained as follows.where
Similarly,
It is interesting to observe that the bias is positive unless, and as T ⟶ ∞, i.e., when the Type-I censoring scheme becomes complete sample life testing experiment. In that case which is nothing but an expected value of complete life testing experiment. The second term on the right-hand side (8) can be written in recursive relation (see, for example, Govindrajulu [16]).
3. Different Confidence and Credible Intervals
In this section, the different confidence intervals are proposed as follows.(1)Based on, the exact distribution of .(2)Based on, the asymptotic distribution of and .(3)Based on, the likelihood ratio test.(4)Bootstrap-p and bootstrap-t confidence intervals.
We have also proposed one Bayes credible interval of based on the inverted gamma prior on .
3.1. Exact Confidence Interval
The exact distribution of was used by Chen and Bhattacharya [5] to present the one-sided confidence interval of the parameter . Also, the two-sided confidence interval of can be obtained easily based on the exact distribution of , in the conventional Type-I censoring scheme under the consideration that is an increasing function of as it was originally proposed by Chen and Bhattacharya [14]. Using Chen and Bhattacharya [14], we observe thatwhere present the chi-square PDF with degrees of freedom, and
Under the consideration of Type-I censoring (11), we observe the cumulative distribution function (CDF) of . Hence, the symmetric confidence interval under Type-I censoring is formulated by choosing and such that
3.2. Asymptotic Confidence Interval
The limiting properties of the MLEs are used to obtain the asymptotic distribution of . Then, the asymptotic distribution of , is asymptotically normally distributed with mean and variance where
The asymptotic confidence interval based on the asymptotic distribution of can be easily calculated. It is observed that
The proof of (15) is given in Appendix B and
Hence, using (15) and (16), we can obtain (14). Therefore, for , the 100 asymptotic confidence interval for is
As reported in Meeker et al. [17], the confidence interval based on the asymptotic theory of In is superior than . Then, the approximate 100 confidence interval under of is given by
Therefore, the 100 approximate confidence interval for for becomeswhich can be easily calculated by using (15) and (16).
3.3. Confidence Interval Based on the LRT
From Meeker and Escobar [11], we observe that the confidence interval based on the LRT is often superior than the confidence interval based on the asymptotic distribution of the MLE. So, we propose, as given by Lawless [8], the confidence interval based on the LRT for constructing confidence intervals for the gamma parameters. But similar methods can be easily adopted here. Under testing of hypothesis problem, firstly considerand the likelihood ratio statistic presented bywhere is the likelihood function (see, for example, (5)). For a large sample size, he the distribution of is approximated under as with one degree of freedom. Then, the exact distribution of is independent of for scale parameter . Therefore, we constructed the 100 confidence interval, of , bywhere the upper 100 percentile point of the , distributed with one degree of freedom. Since, by definition, , is a unimodal function and therefore for a given , is a unique interval.
3.4. Bootstrap Confidence Intervals
In this section, based on the bootstrap technique, we propose two confidence intervals, namely, (a) the percentile bootstrap (Boot-p) confidence interval proposed by Efron [18] and (b) the bootstrap-t (Boot-t) confidence interval proposed by Hall [19]. Different authors have shown interest in bootstrap techniques in the problem of building the bootstrap confidence intervals (see, for example, Almarashi and Abd-Elmougod [20] and Abd-Elmougod and Mahmoud [21]). The following algorithms described the steps needed to construct the Boot-up and Boot-t confidence intervals of .
Boot-p method:(1)Obtain and the MLE of as discussed before.(2)Generate a bootstrap sample using and . Obtain and the MLE of based on the bootstrap sample.(3)Step 2 is repeated, NBOOT times.(4)Let be the CDF of . Define for a given . The approximate 100 confidence interval of is given by
The following method may be used to construct the Boot-t confidence interval of .
Boot-t method:(1)Obtain and the MLE of as discussed before.(2)By using and , the bootstrap sample is generated. Obtain and the MLE of based on the bootstrap sample. Also, compute using (fish-info) for .(3)The statistic value is determined by .(4)Repeat Step 2 and Step 3, NBOOT times.(5)Let be the CDF of Define for a given . So, the approximate 100 confidence interval of is given by
3.5. Bayesian Credible Interval
In this section, we provide the Bayesian analysis of the above-mentioned problem. In the context of exponential lifetimes, may be reasonably modelled by the inverted gamma prior. Following the approach of Draper and Guttman [6], it is assumed that has an inverted gamma prior with parameters and with the following PDF:
When , one obtains the non-informative prior on The posterior PDF of based on the above inverted gamma prior is
From post-theta, it is clear that the posterior distribution is also an inverted gamma distribution, and we can easily obtain the Bayes estimate of under the squared error loss function as
Interestingly, the Bayes estimate based on the non-informative prior coincides with the MLE. Hence, the posterior distribution of is used to obtain credible interval of Note that the posterior distribution of is a distribution having a positive integer degrees of freedom. Hence, 100 credible interval of is presented byfor . For the non-integer value of , gamma distribution can be applied to formulate credible interval of . Then, the non-informative prior can be used if no prior information is available to construct a credible interval for . If some prior information is available, then some positive values of and may be used to construct a credible interval for (Table 1).
4. Numerical Experiments
Since the performances of the different methods cannot be compared theoretically, we use Monte Carlo simulations to compare different methods for different sample sizes and for different censoring times (). Pentium IV processor is used for all the computations, and we use the random number generator RAN2 of Press et al. [1].
We consider different and values, and in all the cases, we have considered . For each case, we have computed the 95% confidence intervals using all the five methods proposed here. For comparison purposes, we have also computed the 95% credible interval of based on the non-informative prior as discussed before. For each dataset, we have computed the length of the confidence interval of based on (a) the exact distribution of (Table 1), (b) the asymptotic distribution of (Table 2), (c) the asymptotic distribution of (Table 3), (d) the likelihood ratio test (Table 4), (e) Boot-p procedure (Table 5), and (f) Boot-t procedure (Table 6) as discussed in Section 3. We have also computed 95% Bayes credible interval as discussed in Section 3 based on the non-informative prior (Table 7). For the asymptotic confidence interval, if the lower limit is negative, it is replaced by zero. We have repeated the experiment 1000 times and reported the average lengths of the confidence intervals and the coverage percentages. In each box, the first quantity represents the average length of the confidence/credible interval, the second quantity within brackets represents the coverage percentages, and the third quantity within the brackets if any, represents the non-existence percentage of the MLE.
Some of the points are quite clear from the above experiments. It is observed that for fixed , as increases, the length of the intervals based on all the above-mentioned procedures decreases. It indicates the consistency property of the MLE. Similarly, for fixed , as increases, the length of the intervals decreases. All the methods more or less maintain the coverage percentages at the nominal level (95% in this case) unless and or is very small.
It is also observed that the exact confidence interval does not work very well for small values of and . It is also observed that the Bayes credible interval works quite well for small values of and . By comparing Table 5 with Table 7, at = 50, we observed that the coverage percentages of Boot-p serve well and are closer to 95% than Bayes credible intervals (4 best out of 7 choices). Also, the widths of Boot-p are almost all shorter than those of Bayes credible intervals. In fact, most of the times, it is better than the other methods considered here. The two terms of their average lengths and coverage percentages are used to compare different confidence intervals. It is observed that the confidence intervals based on the asymptotic distribution of and and Boot-p confidence interval perform quite well, unless and are very small. Since the Boot-t method is involved numerically and the confidence intervals based on the asymptotic distributions are larger than the Bayes credible interval, we recommend to use the Bayes credible interval in all the cases. It may be mentioned that the exact confidence intervals are difficult to obtain. One needs to solve two non-linear equations to compute the exact confidence intervals. Therefore, it may be avoided. We have further reported the bias values (Table 8) for different and ; as expected, the biases tend to zero, as and increase.
5. Data Analysis
In this section, we have analyzed one dataset taken from Bain [9] and applied different methods discussed so far. So, the following dataset is used (see Table 9).
Dataset. From exponential population, suppose 20 items are put on a life testing experiment and the experiment is continued for 150 hours, which is prefixed before the starting of the experiment. During that period, 13 items have failed at the following hours: 3, 19, 23, 26, 37, 38, 41, 45, 58, 84, 90, 109, and 138.
In this case, , , and . Further
Hence,
We have reported different 95% confidence and credible intervals of in Table 9.
6. Conclusion
In this section, we have considered the Type-I censoring scheme in case of exponential failure distributions. We have proposed two-sided exact confidence interval of based on the exact distribution of the conditional MLE of . We have compared different confidence intervals, namely, (a) exact confidence interval, (b) confidence interval based on the asymptotic distribution of , (c) confidence interval based on the asymptotic distribution of , (d) confidence interval based on the LRT, (e) Boot-p confidence interval, and (f) Boot-t along with the credible interval based on the non-informative inverted gamma prior, through extensive computer simulations. We have also analyzed one real dataset for illustrative purposes. It is observed that the Bayes credible interval with non-informative prior works quite well in terms of the length of the interval and coverage percentages. It is also observed that in this case, the Boot-t method does not work at all.
When the data are observed under Type-I censoring scheme, another important point is that at which is typically ignored, the MLE does not exist. One may look at the interesting discussion by Bain and Engelhardt [14] in this aspect.
Appendix
A. Proof of Equation (7)
The conditional MGF of given is given by
Nowwhere are the order statistics of a random sample of size from the distribution with PDFwhere . Let ; then,
Again,
Hence, we get equation (7) as
B. Proof of Equation (15)
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This project was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, under grant no. KEP-PhD-75-130-42. The authors, therefore, acknowledge with thanks the DSR for technical and financial support.