On the exact distribution of the MLEs based on Type-II progressively hybrid censored data from exponential distributions

ABSTRACT

Distributions of the maximum likelihood estimators (MLEs) in Type-II (progressive) hybrid censoring based on two-parameter exponential distributions have been obtained using a moment generating function approach. Although resulting in explicit expressions, the representations are complicated alternating sums. Using the spacings-based approach of Cramer and Balakrishnan [On some exact distributional results based on Type-I progressively hybrid censored data from exponential distributions. Statist Methodol. 2013;10:128–150], we derive simple expressions for the exact density and distribution functions of the MLEs in terms of B-spline functions. These representations can be easily implemented on a computer and provide an efficient method to compute density and distribution functions as well as moments of Type-II (progressively) hybrid censored order statistics.     

Exact likelihood inference based on an unified hybrid censored sample from the exponential distribution

Recently, exact confidence bounds and exact likelihood inference have been developed based on hybrid censored samples by Chen and Bhattacharyya [Chen, S. and Bhattacharyya, G.K. (1998). Exact confidence bounds for an exponential parameter under hybrid censoring. Communications in Statistics—Theory and Methods, 17, 1857–1870.], Childs et al. [Childs, A., Chandrasekar, B., Balakrishnan, N. and Kundu, D. (2003). Exact likelihood inference based on Type-I and Type-II hybrid censored samples from the exponential distribution. Annals of the Institute of Statistical Mathematics, 55, 319–330.], and Chandrasekar et al. [Chandrasekar, B., Childs, A. and Balakrishnan, N. (2004). Exact likelihood inference for the exponential distribution under generalized Type-I and Type-II hybrid censoring. Naval Research Logistics, 51, 994–1004.] for the case of the exponential distribution. In this article, we propose an unified hybrid censoring scheme (HCS) which includes many cases considered earlier as special cases. We then derive the exact distribution of the maximum likelihood estimator as well as exact confidence intervals for the mean of the exponential distribution under this general unified HCS. Finally, we present some examples to illustrate all the methods of inference developed here.     

On the existence and uniqueness of the maximum likelihood estimates of the parameters of Birnbaum–Saunders distribution based on Type-I,Type-II and hybrid censored samples

The Birnbaum–Saunders (BS) distribution is a positively skewed distribution and is a common model for analysing lifetime data. In this paper, we discuss the existence and uniqueness of the maximum likelihood estimates (MLEs) of the parameters of BS distribution based on Type-I, Type-II and hybrid censored samples. The line of proof is based on the monotonicity property of the likelihood function. We then describe the numerical iterative procedure for determining the MLEs of the parameters, and point out briefly some recently developed simple methods of estimation in the case of Type-II censoring. Some graphical illustrations of the approach are given for three real data from the reliability literature. Finally, for illustrative purpose, we also present an example in which the MLEs do not exist.     

Exact Bayesian variable sampling plans for the exponential distribution with progressive hybrid censoring

From the exact distribution of the maximum likelihood estimator of the average lifetime based on progressive hybrid exponential censored sample, we derive an explicit expression for the Bayes risk of a sampling plan when a quadratic loss function is used. The simulated annealing algorithm is then used to determine the optimal sampling plan. Some optimal Bayes solutions under progressive hybrid and ordinary hybrid censoring schemes are presented to illustrate the effectiveness of the proposed method.     

Robust Bayesian analysis for exponential parameters under generalized Type-II progressive hybrid censoring

The Type-II progressive hybrid censoring scheme has received wide attention, but it has a disadvantage in that long time may be required to complete the life test. The generalized progressive Type-II hybrid censoring scheme has recently been proposed to solve this problem. Under the censoring scheme, the time on test does not exceed a predetermined time. In this paper, we propose a robust Bayesian approach based on a hierarchical structure when the generalized progressive Type-II hybrid censored sample has a two-parameter exponential distribution. For unknown parameters, marginal posterior distributions are provided in closed forms, and their statistical properties are discussed. To examine the robustness of the proposed method, Monte Carlo simulations are conducted and a real data set is analyzed. Further, the quality and adequacy of the proposed model are evaluated in an analysis based on the real data.     

Bayesian Analysis of Different Hybrid and Progressive Life Tests

Type-I and Type-II censoring schemes are the widely used censoring schemes available for life testing experiments. A mixture of Type-I and Type-II censoring schemes is known as a hybrid censoring scheme. Different hybrid censoring schemes have been introduced in recent years. In the last few years, a progressive censoring scheme has also received considerable attention. In this article, we mainly consider the Bayesian inference of the unknown parameters of two-parameter exponential distribution under different hybrid and progressive censoring schemes. It is observed that in general the Bayes estimate and the associated credible interval of any function of the unknown parameters, cannot be obtained in explicit form. We propose to use the Monte Carlo sampling procedure to compute the Bayes estimate and also to construct the associated credible interval. Monte Carlo Simulation experiments have been performed to see the effectiveness of the proposed method in case of Type-I hybrid censored samples. The performances are quite satisfactory. One data analysis has been performed for illustrative purposes.     

On progressive hybrid censored exponential distribution

In this paper, we introduce a new adaptive Type-I progressive hybrid censoring scheme, which has some advantages over the progressive hybrid censoring schemes already discussed in the literature. Based on an adaptive Type-I progressively hybrid censored sample, we derive the exact distribution of the maximum-likelihood estimator (MLE) of the mean lifetime of an exponential distribution as well as confidence intervals for the failure rate using exact distribution, asymptotic distribution, and three parametric bootstrap resampling methods. Furthermore, we provide computational formula for the expected number of failures and investigate the performance of the point and interval estimation for the failure rate in this case. An alternative simple form for the distribution of the MLE under adaptive Type-II progressive hybrid censoring scheme proposed by Ng et al. [Statistical analysis of exponential lifetimes under an adaptive Type-II progressive censoring scheme, Naval Res. Logist. 56 (2009), pp. 687–698] is obtained. Finally, from the exact distribution of the MLE, we establish the explicit expression for the Bayes risk of a sampling plan under adaptive Type-II progressive hybrid censoring scheme when a general loss function is used, and present some optimal Bayes solutions under four different progressive hybrid censoring schemes to illustrate the effectiveness of the proposed method.     

Generalized inverted exponential distribution under hybrid censoring

The hybrid censoring scheme is a mixture of Type-I and Type-II censoring schemes. Based on hybrid censored samples, we first derive the maximum likelihood estimators of the unknown parameters and the expected Fisher’s information matrix of the generalized inverted exponential distribution (GIED). Monte Carlo simulations are performed to study the performance of the maximum likelihood estimators. Next we consider Bayes estimation under the squared error loss function. These Bayes estimates are evaluated by applying Lindley’s approximation method, the importance sampling procedure and Metropolis–Hastings algorithm. The importance sampling technique is used to compute the highest posterior density credible intervals. Two data sets are analyzed for illustrative purposes. Finally, we discuss a method of obtaining the optimum hybrid censoring scheme.     

Nonparametric tests for stratified additive hazards model based on current status data

Stratified regression models are commonly employed when study subjects may come from possibly different strata such as different medical centers, and for the situation, one common question of interest is to test the existence of the stratum effect. To address this, there exists some literature on the testing of the stratum effects under the framework of the proportional hazards model when one observes right-censored data or interval-censored data. In this paper, we consider the situation under the additive hazards model when one faces current status data, for which there does not seem to exist an established test procedure. The asymptotic distributions of the proposed test procedure are provided. Also a simulation study is performed to evaluate the performance of the proposed method and indicates that it works well for practical situations. The approach is applied to a set of real current status data from a tumorigenicity study.     

Quantile regression based on a weighted approach under semi-competing risks data

In this article, we investigate the quantile regression analysis for semi-competing risks data in which a non-terminal event may be dependently censored by a terminal event. Due to the dependent censoring, the estimation of quantile regression coefficients on the non-terminal event becomes difficult. In order to handle this problem, we assume Archimedean Copula to specify the dependence of the non-terminal event and the terminal event. Portnoy [Censored regression quantiles. J Amer Statist Assoc. 2003;98:1001–1012] considered the quantile regression model under right-censoring data. We extend his approach to construct a weight function, and then impose the weight function to estimate the quantile regression parameter for the non-terminal event under semi-competing risks data. We also prove the consistency and asymptotic properties for the proposed estimator. According to the simulation studies, the performance of our proposed method is good. We also apply our suggested approach to analyse a real data.     

Robust Bayesian estimation of a two-parameter exponential distribution under generalized Type-I progressive hybrid censoring

A generalized Type-I progressive hybrid censoring scheme was proposed recently to overcome the limitations of the progressive hybrid censoring scheme. In this article, we provide a robust Bayesian method to estimate the unknown parameters of the two-parameter exponential distribution of a generalized Type-I progressive hybrid censored sample. For each parameter, we derive the marginal posterior density functions and the corresponding Bayesian estimators under the squared error loss function. To assess the proposed method, Monte Carlo simulations are performed using a real dataset.     

Exact distribution of the MLEs of the parameters and of the quantiles of two-parameter exponential distribution under hybrid censoring

Epstein [Truncated life tests in the exponential case, Ann. Math. Statist. 25 (1954), pp. 555–564] introduced a hybrid censoring scheme (called Type-I hybrid censoring) and Chen and Bhattacharyya [Exact confidence bounds for an exponential parameter under hybrid censoring, Comm. Statist. Theory Methods 17 (1988), pp. 1857–1870] derived the exact distribution of the maximum-likelihood estimator (MLE) of the mean of a scaled exponential distribution based on a Type-I hybrid censored sample. Childs et al. [Exact likelihood inference based on Type-I and Type-II hybrid censored samples from the exponential distribution, Ann. Inst. Statist. Math. 55 (2003), pp. 319–330] provided an alternate simpler expression for this distribution, and also developed analogous results for another hybrid censoring scheme (called Type-II hybrid censoring). The purpose of this paper is to derive the exact bivariate distribution of the MLE of the parameter vector of a two-parameter exponential model based on hybrid censored samples. The marginal distributions are derived and exact confidence bounds for the parameters are obtained. The results are also used to derive the exact distribution of the MLE of the pth quantile, as well as the corresponding confidence bounds. These exact confidence intervals are then compared with parametric bootstrap confidence intervals in terms of coverage probabilities. Finally, we present some numerical examples to illustrate the methods of inference developed here.     

Prediction of future failures times based on Type-I hybrid censored samples of random sample sizes

In this article, we consider the prediction of future failure times based on Type-I hybrid censored samples. Point predictors and prediction intervals using different procedures are discussed for a general model. The exponential and Rayleigh distributions are used as illustrative examples to show the most simplified forms of the so obtained predictors as well as prediction intervals. Intensive simulation study and a real life dataset are presented to illustrate our findings and results.     

Quantile residual life regression based on semi-competing risks data

This paper investigates the quantile residual life regression based on semi-competing risk data. Because the terminal event time dependently censors the non-terminal event time, the inference on the non-terminal event time is not available without extra assumption. Therefore, we assume that the non-terminal event time and the terminal event time follow an Archimedean copula. Then, we apply the inverse probability weight technique to construct an estimating equation of quantile residual life regression coefficients. But, the estimating equation may not be continuous in coefficients. Thus, we apply the generalized solution approach to overcome this problem. Since the variance estimation of the proposed estimator is difficult to obtain, we use the bootstrap resampling method to estimate it. From simulations, it shows the performance of the proposed method is well. Finally, we analyze the Bone Marrow Transplant data for illustrations.     

Estimation under inverse Weibull distribution based on Balakrishnan’s unified hybrid censored scheme

In this article, point and interval estimations of the parameters α and β of the inverse Weibull distribution (IWD) have been studied based on Balakrishnan’s unified hybrid censoring scheme (UHCS), see Balakrishnan et al. In point estimation, the maximum likelihood (ML) and Bayes (B) methods have been used. The Bayes estimates have been computed based on squared error loss (SEL) function and Linex loss function and using Markov Chain Monte Carlo (MCMC) algorithm. In interval estimation, a (1 ? τ) × 100% approximate, bootstrap-p, credible and highest posterior density (HPD) confidence intervals (CI′s) for the parameters α and β have been introduced. Based on Monte Carlo simulation, Bayes estimates have been compared with their corresponding maximum likelihood estimates by computing the mean squared errors (MSE′s) of all estimators. Finally, point and interval estimations of all parameters have been studied based on a real data set as an illustrative example.     

Hybrid censoring schemes with exponential failure distribution

The mixture of Type I and Type I1 censoring schemes, called the hybrid censoring, is quite important in life–testing experiments. Epstein(1954, 1960) introduced this testing scheme and proposed a two–sided confidence interval to estimate the mean lifetime, θ, when the underlying lifetime distribution is assumed to be exponential. There are some two–sided confidence intervals and credible intervals proposed by Fairbanks et al. (1982) and Draper and Guttman (1987) respectively. In this paper we obtain the exact two–sided confidence interval of θ following the approach of Chen and Bhattacharya (1988). We also obtain the asymptotic confidence intervals in the Hybrid censoring case. It is important to observe that the results for Type I and Type II censoring schemes can be obtained as particular cases of the Hybrid censoring scheme. We analyze one data set and compare different methods by Monte Carlo simulations.     

Fisher information in progressive hybrid censoring schemes

The hybrid censoring scheme, which is a mixture of Type-I and Type-II censoring schemes, has been extended to the case of progressive censoring schemes by Kundu and Joarder [Analysis of Type-II progressively hybrid censored data, Comput. Stat. Data Anal. 50 (2006), pp. 2509–2528] and Childs et al. [Exact likelihood inference for an exponential parameter under progressive hybrid censoring schemes, in Statistical Models and Methods for Biomedical and Technical Systems, F. Vonta, M. Nikulin, N. Limnios, and C. Huber-Carol, eds., Birkhäuser, Boston, MA, 2007, pp. 323–334]. In this paper, we derive a simple expression for the Fisher information contained in Type-I and Type-II progressively hybrid censored data. An illustrative example is provided applicable to a scaled-exponential distribution to demonstrate our methodologies.     

Exact likelihood inference for two exponential populations based on a joint generalized Type-I hybrid censored sample

Following the work of Chen and Bhattacharyya [Exact confidence bounds for an exponential parameter under hybrid censoring. Comm Statist Theory Methods. 1988;17:1857–1870], several results have been developed regarding the exact likelihood inference of exponential parameters based on different forms of censored samples. In this paper, the conditional maximum likelihood estimators (MLEs) of two exponential mean parameters are derived under joint generalized Type-I hybrid censoring on the two samples. The moment generating functions (MGFs) and the exact densities of the conditional MLEs are obtained, using which exact confidence intervals are then developed for the model parameters. We also derive the means, variances, and mean squared errors of these estimates. An efficient computational method is developed based on the joint MGF. Finally, an example is presented to illustrate the methods of inference developed here.     

Testing the validity of the exponential model for hybrid Type-I censored data

The exponential distribution has been used in life-testing and reliability studies. In this article, we first express the entropy of Type-I hybrid censoring scheme in terms of hazard function and provide an estimate of the entropy of Type-I hybrid censored data. Then, we construct a goodness-of-fit test statistic based on Kullback–Leibler information for Type-I hybrid censored data. The test statistic is used to test for exponentiality. A Monte Carlo simulation is conducted to obtain the power of the proposed test against various alternatives. Finally, a data example is presented for illustrative purpose.