Next Article in Journal
A Comprehensive Comparison of the Performance of Metaheuristic Algorithms in Neural Network Training for Nonlinear System Identification
Previous Article in Journal
On the Use of Morpho-Syntactic Description Tags in Neural Machine Translation with Small and Large Training Corpora
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical Inference and Optimal Design of Accelerated Life Testing for the Chen Distribution under Progressive Type-II Censoring

School of Mathematics and Statistics, Beijing Jiaotong University, Beijing 100044, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(9), 1609; https://doi.org/10.3390/math10091609
Submission received: 9 April 2022 / Revised: 3 May 2022 / Accepted: 5 May 2022 / Published: 9 May 2022

Abstract

:
This paper discusses statistical inference and optimal design of constant-stress accelerated life testing for the Chen distribution under progressive Type-II censoring. The scale parameter of the life distribution is assumed to be a logarithmic linear function of the stress level. The maximum likelihood estimates of the parameters are obtained. Then, the observed Fisher information matrix is derived and utilized to construct asymptotic confidence intervals. Meanwhile, the parametric bootstrap methods are provided for the interval estimation. In addition, the Bayes estimates under the squared error loss function are obtained by applying the Tierney and Kadane technique and Lindley’s approximation. As for the optimal design, D- and A-optimality criteria are considered to determine the optimal transformed stress level. Finally, the simulation is carried out to demonstrate the proposed estimation techniques and the optimal criteria, and a real data set is discussed.

1. Introduction

In factory production, life testing is of great significance to evaluate product quality. Conventional life testing is to mainly observe the failure time information of all test units under normal usage in order to reveal the life characteristics. However, products are produced more and more reliably currently, so it is difficult or even impossible for them to fail under normal usage. To reduce the cost and time, it is necessary to make the test units fail more quickly. A common way is to perform life testing at higher stress levels than normal, which is called accelerated life testing (ALT). Then, the failure time information obtained is analyzed to extrapolate the life characteristics under the normal stress level. One type of classifying ALT is the way stress changes with time. ALTs are mainly classified into constant-stress ALT (CSALT), step-stress ALT (SSALT) and progressive-stress ALT (PSALT). The different focus of the three types is the relationship between stress loading and testing time. For CSALT, the stress is fixed with time; for SSALT, the stress is increased step by step at specific moments; for PSALT, stress is increased continuously with time. Another type of classifying ALTs is by the number of stress levels. They are mainly divided into simple ALT and multiple ALT. The simple ALT has only two stress levels while the multiple ALT has more than two stress levels.
Some of the earlier works on ALT are mainly based on McCool [1], Nelson [2] and Miller and Nelson [3]. Gradually, the study of ALT has attracted the interest of more and more researchers. For the studies of CSALT, Ismail [4] discussed the estimation of parameters of Weibull distribution in the constant-stress partial ALT under hybrid censoring. Statistical inferences on the acceleration factor and the distribution parameters were made utilizing maximum likelihood and percentile bootstrap. Mohie EI-Din et al. [5] derived the optimal plans for k-level CSALT for Lindley failure data and compared the optimal plans calculated by two optimal criteria with the traditional plan calculated by the asymptotic variance of maximum likelihood estimation. In addition, many researchers studied SSALT. Mohie El-Din et al. [6] explored parametric inference for the extension of exponential distribution on SSALT under progressive Type-II censoring. They assumed the cumulative exposure model and obtained the point estimates and interval estimates such as Bayes estimation and parametric bootstrap intervals. Chandra and Khan [7] studied the optimum test plan and parameter estimation for a Lomax distribution for 3-step SSALT under modified progressive Type-I censoring. Interval and point estimates of the model parameters were obtained, and the Bayes estimates are obtained based on Gibbs sampling. Regarding the studies of PSALT, Abdel-Hamid and Abushal [8] applied PSALT when the product lifetime under design stress followed the exponential distribution. Bayes estimates were obtained and compared with maximum likelihood estimates.
Obtaining complete failure time data is ideal in life testing. However, considering the time cost and the possibility of unexpected product failure in practice, researchers increasingly focus on the application of the censoring method. The two basic censoring methods are Type-I and Type-II censoring. A life test is terminated at a fixed time and the surviving test units are withdrawn in Type-I censoring, whereas a life test is terminated when a specified number of units fail in Type-II censoring. Considering that the surviving units can only be withdrawn at the final termination point of the life test, progressive censoring has been more and more popular which enables to withdraw units at various points in the life test. For example, the progressive Type-II censoring is commonly applied in various fields. For more studies on progressive censoring, one can refer to Balakrishnan and Aggarwala [9] and Balakrishnan [10]. Jaheen et al. [11] investigated Bayes inference of constant-partial ALT when the lifetime of test units follows the generalized exponential distribution. Under progressive Type-II censoring, the maximum likelihood and Bayes methods of estimation were utilized for estimating acceleration factor and parameters. Basak and Balakrishnan [12] discussed the problem of predicting survival times of units for the exponential distribution in the simple step-stress testing under progressive Type-II censoring. The maximum likelihood predictors and the conditional median predictors were derived. Abd El-Raheem [13] considered optimal design and estimation of multiple CSALT for generalized half-normal distribution under progressive Type-II censoring. Using various estimation methods, point estimation, interval estimation and predictive estimation were obtained and the optimal transformed stress level was derived using the optimal criteria. Algarni et al. [14] considered E-Bayesian estimation of the scale parameter, hazard and reliability rate functions of the Chen distribution under Type-I censoring. On the basis of the balanced squared error loss function, the gamma distribution was used as a conjugate prior for the scale parameter and three different distributions were discussed for the hyper-parameters.
Chen [15] proposed a two-parameter lifetime distribution with an increasing or bathtub-shaped hazard function. It has been considered by many researchers because the bathtub-shaped hazard function is a suitable conceptual model for a number of machinery and electronics; see Hjorth [16] and Olkin [17]. Due to the various good properties of the Chen distribution, many authors have discussed its statistical inference. Sarhan et al. [18] obtained the Bayes estimation of parameters with gamma priors for the Chen distribution by different approximations methods and compared the maximum likelihood estimates with the approximate Bayesian estimates under non-informative prior. Rastogi and Tripathi [19] used hybrid censoring to estimate unknown parameters of the Chen distribution. Maximum likelihood estimates were obtained by the EM algorithm and Bayes estimates were evaluated by the Tierney and Kadane method and Lindley method. Asymptotic and parametric bootstrap intervals were also proposed. Ahmed [20] derived the Bayes estimators of the parameters, reliability, and hazard function for the Chen distribution under both the balanced linear-exponential loss and balanced-squared error loss utilizing the MCMC method. Moreover, to construct the asymptotic confidence intervals, delta and parametric bootstrap methods were applied to find the variance.
Considering that ALT is relatively rarely used for the Chen distribution which has a bathtub-shaped hazard function and that optimal design studies are less concerned with solving for optimal transformed stress levels, we apply the CSALT model to test units whose lifetimes follow the Chen distribution in the presence of progressive Type-II censoring in this paper. To better estimate the lifetime characteristics of products, we discuss a variety of parameter and interval estimation methods. The maximum likelihood estimates of the parameters for the Chen distribution are obtained under progressive censoring. The Bayes estimates of the parameters are derived by both the Tierney and Kadane technique and Lindley’s approximation. We also discuss three types of confidence intervals which are the asymptotic confidence interval, the bootstrap percentile confidence interval, and the bootstrap-t confidence interval. In addition, to better estimate the life parameters of products under the normal stress level, two optimal criteria are utilized to determine and specify the optimal transformed stress level of CSALT. To measure the performance of the estimates, we conduct some simulations.
The rest of this paper is organized as follows. In Section 2, we discussed some properties of the Chen distribution. In Section 3, we propose the progressive censoring scheme and necessary model assumptions. Then, the maximum likelihood estimates are obtained in Section 4. In Section 5, we derive the Fisher information matrix and the variance–covariance matrix. The asymptotic confidence intervals are also constructed. In Section 6, we present the parametric bootstrap confidence intervals. In Section 7, we discuss the Bayesian estimation and utilize the Tierney and Kadane technique and Lindley’s approximation to obtain the approximate estimates. In Section 8, we apply two optimal criteria to determine the optimal transformed stress levels. In Section 9, a numerical study is conducted to illustrate the methods proposed and examine their performance. In addition, a real data set is discussed in Section 10. Finally, we make conclusions.

2. The Chen Distribution

The probability density function (pdf), the cumulative distribution function (cdf), and the hazard function (hf) for the Chen distribution are expressed respectively as follows:
f x ; α , β = α β x β 1 exp α 1 e x β + x β , x > 0 , α > 0 , β > 0 ,
F x ; α , β = 1 exp α 1 e x β , x > 0 , α > 0 , β > 0 ,
h x ; α , β = f x ; α , β 1 F x ; α , β = α β x β 1 e x β , x > 0 , α > 0 , β > 0 ,
where α is the scale parameter and β is the shape parameter. When β is less than one, the hazard function has a bathtub shape, while the shape is increasing when β 1 . Different plots of pdf and hf are shown in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6. The shape of pdf or hf varies for different parameter values. Assume that β is fixed. If 0 < β < 1 : when 0 < α < 1 , pdf is monotonically decreasing or decreasing, then increasing and finally decreasing, and the tail is thinner with a larger value of α ; when α 1 , pdf is monotonically decreasing, and the decrease is sharper with a larger value of α . If β = 1 : when 0 < α < 1 , the pdf is unimodal; when α 1 , the pdf is monotonically decreasing and the decrease is sharper with a larger value of α . If β > 1 : the pdf is unimodal no matter what α value is taken, and the peak value is larger with a larger value of α . From Figure 5 and Figure 6, it can be seen that the hazard rate is larger with a larger value of α if 0 < β < 1 ; the rising rate of hf is higher with a larger value of α if β 1 .
When β is fixed at some value less than 1, the relationship between the mean of the Chen distribution and α is shown in Figure 7. The mean decreases when the value of α increases. Additionally, letting the cdf of the Chen distribution equal 1 2 , we can solve for the median x m = ln 1 + 1 α ln 2 1 β . Similarly, when β is fixed and α increases, the median decreases.
In addition, the Chen distribution has a great deal of flexibility. Let X be a nonnegative random variable and follow the Chen distribution. Following Ahmed et al. [21], Y = e X β 1 follows the exponential distribution and Y = e X β 1 1 θ follows the Weibull distribution. Furthermore, the Chen distribution reduces to Gompertz distribution when β = 1 .

3. Model Assumptions

For CSALT, assume a total of n units under the test and k 2 stress levels, s 1 < < s k . At each stress level s i , i = 1 , 2 , , k , n i identical units are assigned for the test, such that i = 1 k n i = n . The lifetime of each unit X i j , j = 1 , 2 , , n i , follows the Chen distribution. s 0 is the normal stress level. The progressive Type-II censoring scheme in CSALT is briefly described as follows: For i = 1 , 2 , , k , when the first failure occurs at time t i 1 : m i : n i , R i 1 units are randomly withdrawn from the remaining n i 1 units whose lifetime X i j > t i 1 : m i : n i . Then, when the second failure occurs at time t i 2 : m i : n i , R i 2 units are randomly withdrawn from the remaining n i 2 R i 1 units whose lifetime X i j > t i 2 : m i : n i . Similarly, when the m i -th 1 m i n i failure occurs at a random time of t i m i : m i : n i , all remaining R i m i = n i m i j = 1 m i 1 R i j units are withdrawn and the test is terminated. The censoring scheme is shown in Figure 8. Under stress level s i , i = 1 , 2 , k , R i 1 , R i 2 , , R i m i is the progressive censoring scheme and t i j : m i : n i is the observed progressive censored data that t i 1 : m i : n i < < t i m i : m i : n i . Obviously, the sample is complete when R i 1 , R i 2 , , R i m i = 0 , 0 , , 0 and the sample is Type-II censored when R i 1 , R i 2 , , R i m i = 0 , , 0 , n i m i .
The following model assumptions are used throughout the whole paper:
(1)
Under stress level s i , i = 1 , 2 , , k , the lifetimes of the units follow
f i x ; α i , β = α i β x β 1 exp α i 1 e x β + x β , x > 0 , α i > 0 , β > 0 .
(2)
The relationship between the stress level s i and the lifetime distribution parameter α i is
ln α i = b 0 + b φ i , i = 0 , 1 , , k ,
where φ i = φ s i is an incremental function of s i . b 0 and b > 0 are unknown parameters.
To simplify the expression, according to the model (5), the scale parameter α i can be presented as
α i = α 0 exp b φ i φ 0 = α 0 λ h i , i = 0 , 1 , , k ,
where α 0 is the scale parameter for the Chen distribution under normal stress level s 0 , φ 0 = φ s 0 , λ = exp b φ 1 φ 0 = α 1 α 0 > 1 and
h i = φ i φ 0 φ 1 φ 0 , i = 1 , 2 , , k ,
which satisfy h k > h k 1 > > h 1 = 1 . Since the Jacobian determinant from b 0 , b , β to α 0 , λ , β is 1 / λ α 0 φ 1 φ 0 > 0 , the transformation from the model parameter b 0 , b , β to α 0 , λ , β is mapped one-to-one.

4. Maximum Likelihood Estimation

Based on the above assumptions and notations, when the stress level is s i , the likelihood function of parameters α 0 , λ and β is given by
L i α 0 , λ , β = C i j = 1 m i α i β t i j β 1 exp α i 1 e t i j β + t i j β exp α i 1 e t i j β R i j ,
where C i = n i n i 1 R i 1 n i m i + 1 j = 1 m i 1 R i j and t i j = t i j : m i : n i . Then, the likelihood function is given by
L α 0 , λ , β = i = 1 k L i α 0 , λ , β = i = 1 k C i j = 1 m i α i β t i j β 1 exp α i 1 e t i j β + t i j β exp α i 1 e t i j β R i j .
Correspondingly, the log-likelihood function ignoring the constants is given by
l α 0 , λ , β = ln L α 0 , λ , β = ln β i = 1 k m i + i = 1 k m i ln α i + β 1 i = 1 k j = 1 m i ln t i j + i = 1 k j = 1 m i α i R i j + 1 1 e t i j β + i = 1 k j = 1 m i t i j β .
The partial derivatives of l α 0 , λ , β with respect to α 0 , λ and β are given by
l α 0 = 1 α 0 i = 1 k m i + i = 1 k j = 1 m i λ h i R i j + 1 1 e t i j β ,
l λ = 1 λ i = 1 k m i h i + i = 1 k j = 1 m i α 0 h i λ h i 1 R i j + 1 1 e t i j β ,
l β = 1 β i = 1 k m i + i = 1 k j = 1 m i t i j β + 1 ln t i j i = 1 k j = 1 m i α i e t i j β t i j β ln t i j R i j + 1 .
We can solve the maximum likelihood estimates (MLE) α ^ 0 , λ ^ , β ^ for α 0 , λ and β by making the functions (10)–(12) equal to 0 simultaneously. Unfortunately, explicit solutions cannot be obtained and numerical methods such as the Newton–Raphson method can be utilized to calculate the estimates. According to MacDonald [22], for a constrained optimization problem, one can enforce the constraint by means of a transformation or use of a constrained optimizer. Considering the three parameters are constrained ( α > 0 , λ > 1 , β > 0 ), we use the optim function in R software with method L-BFGS-B and default settings. L-BFGS-B is an application of the quasi-Newton method.

5. Fisher Information Matrix

Let Θ = α 0 , λ , β and T be the observed failure data. For a large-size sample, the inverse of the Fisher information matrix (FIM) can be regarded as the approximation of the variance–covariance matrix. First, FIM is given as follows by evaluating the negative of the expected values of the second-order derivatives and mixed partial derivatives of log-likelihood function in (9) with respect to α 0 , λ and β ,
I T Θ = E 2 l Θ Θ 2 = E 2 l Θ α 0 2 E 2 l Θ α 0 λ E 2 l Θ α 0 β E 2 l Θ λ α 0 E 2 l Θ λ 2 E 2 l Θ λ β E 2 l Θ β α 0 E 2 l Θ β λ E 2 l Θ β 2 .
Due to the complexity of the pdf and the integral formula, in practice, the observed FIM which utilizes the MLEs of parameters Θ ^ = α ^ 0 , λ ^ , β ^ is used to approximate the exact values,
O T Θ = 2 l Θ α 0 2 2 l Θ α 0 λ 2 l Θ α 0 β 2 l Θ λ α 0 2 l Θ λ 2 2 l Θ λ β 2 l Θ β α 0 2 l Θ β λ 2 l Θ β 2 α 0 = α ^ 0 , λ = λ ^ , β = β ^ ,
where
2 l Θ α 0 2 = i = 1 k m i α 0 2 ,
2 l Θ λ 2 = i = 1 k m i h i λ 2 i = 1 k j = 1 m i α 0 λ h i 2 h i h i 1 R i j + 1 1 e t i j β ,
2 l Θ β 2 = i = 1 k m i β 2 + i = 1 k j = 1 m i α i R i j + 1 ln 2 t i j e t i j β t i j β t i j β + 1 i = 1 k j = 1 m i t i j β ln 2 t i j ,
2 l Θ α 0 λ = 2 l Θ λ α 0 = i = 1 k j = 1 m i h i λ h i 1 R i j + 1 1 e t i j β ,
2 l Θ α 0 β = 2 l Θ β α 0 = i = 1 k j = 1 m i λ h i R i j + 1 e t i j β t i j β ln t i j ,
2 l Θ λ β = 2 l Θ β λ = i = 1 k j = 1 m i α 0 h i λ h i 1 R i j + 1 e t i j β t i j β ln t i j .
According to the observed FIM in (14), the asymptotic variance–covariance matrix of Θ ^ can be constructed as O T 1 Θ . Then, on the basis of the respective asymptotic variances and the asymptotic normality of Θ ^ , the 100 ( 1 τ ) % asymptotic confidence intervals can be given by
θ ^ l , θ ^ u = θ ^ ± Z 1 τ 2 V a r θ ^ ,
where θ is α 0 , λ or β , and Z r is the 100 r -th percentile of the standard normal distribution.

6. Parametric Bootstrap Intervals

The asymptotic confidence interval based on the normal approximation which we mentioned in Section 5 has a better performance for the large sample size. In this section, two parametric bootstrap methods are discussed to find more accurate approximate confidence intervals for unknown parameters α 0 , λ and β when the sample size is not large. The algorithms are as follows.

6.1. The Bootstrap Percentile Confidence Interval

-
Step 1: Utilize the progressive censored sample to compute the MLE Θ ^ = α ^ 0 , λ ^ , β ^ .
-
Step 2: Generate a progressive censored sample by regarding Θ ^ as the parameter values for the Chen distribution.
-
Step 3: Compute the MLE Θ ^ * = α ^ 0 * , λ ^ * , β ^ * using the samples from Step 2.
-
Step 4: Repeat Step 2 and 3, B times.
-
Step 5: Estimate the cdf of θ * ^ by F 1 ^ x = P θ * ^ x where θ is α 0 , λ or β and let θ ^ B o o t p = F 1 ^ 1 x be the inverse of F 1 ^ for a given x 0 < x < 1 .
By following the above steps, the 100 ( 1 τ ) % bootstrap percentile confidence interval of θ is given by
θ ^ l , θ ^ u = θ ^ B o o t p τ 2 , θ ^ B o o t p 1 τ 2 .

6.2. The Bootstrap-t Confidence Interval

-
Steps 1 and 2 are the same as those in the algorithm of the bootstrap percentile confidence interval.
-
Step 3: Compute the MLE θ * ^ using the progressive censored samples from Step 2 where θ is α 0 , λ or β . Then, obtain the statistic
T * = θ * ^ θ ^ V a ^ r θ * ^ .
-
Step 4: Repeat Step 2 and 3, B times.
-
Step 5: Estimate the cdf of T * by F 2 ^ x = P T * x and let θ ^ B o o t t = θ ^ + V a ^ r θ * ^ F 2 ^ 1 x for a given x 0 < x < 1 .
By following the above steps, the 100 ( 1 τ ) % bootstrap-t confidence interval of θ is given by
θ ^ l , θ ^ u = θ ^ B o o t t τ 2 , θ ^ B o o t t 1 τ 2 .

7. Bayesian Estimation

The aim of this section is to obtain the Bayes estimates (BE) of unknown parameters α 0 , λ and β , utilizing progressive censored samples of CSALT. First, considering α 0 > 0 , λ > 1 and β > 0 , the prior distributions of α 0 and β are assumed to be the gamma distributions and the prior distribution of λ is assumed to be the left truncated gamma distribution as
π 1 α 0 α 0 a 1 1 e b 1 α 0 , α 0 > 0 , a 1 > 0 , b 1 > 0 ,
π 2 λ λ a 2 1 e b 2 λ 1 u a 2 1 e b 2 u d u , λ > 1 , a 2 > 0 , b 2 > 0 ,
π 3 β β a 3 1 e b 3 β , β > 0 , a 3 > 0 , b 3 > 0 ,
where all the hyperparameters a i and b i , ( i = 1 , 2 , 3 ) are considered to be known and non-negative. When the hyperparameters all tend to zero, the BEs are obtained in the case of non-informative prior. Because of the independence of α 0 , λ and β , their joint prior distribution is given by
π α 0 , λ , β α 0 a 1 1 λ a 2 1 β a 3 1 e b 1 α 0 + b 2 λ + b 3 β 1 u a 2 1 e b 2 u d u .
Denote the observed failure data as T. Then using both the likelihood function (8) and the joint prior distribution (25), the posterior distribution is given by
π * α 0 , λ , β T π α 0 , λ , β L T α 0 , λ , β 0 1 0 π α 0 , λ , β L T α 0 , λ , β d α 0 d λ d β .
Let g α 0 , λ , β be an arbitrary function of α 0 , λ and β . The BE of g α 0 , λ , β under squared error loss function is given by
g B E ^ α 0 , λ , β = E g α 0 , λ , β T = 0 1 0 g α 0 , λ , β π α 0 , λ , β L T α 0 , λ , β d α 0 d λ d β 0 1 0 π α 0 , λ , β L T α 0 , λ , β d α 0 d λ d β .
For the ease of calculation, the BE of g α 0 , λ , β is rewritten as
g B E ^ α 0 , λ , β = 0 1 0 g α 0 , λ , β e l α 0 , λ , β + ρ α 0 , λ , β d α 0 d λ d β 0 1 0 e l α 0 , λ , β + ρ α 0 , λ , β d α 0 d λ d β ,
where l α 0 , λ , β = ln L α 0 , λ , β and ρ α 0 , λ , β = ln π α 0 , λ , β .
Considering the complexity of the likelihood function (8) and the ratio of two multiple integrals, g B E ^ α 0 , λ , β cannot be given by an analytical form. Therefore, it is necessary to employ some approximation methods to obtain the approximate posterior expectations. Among others, the methods suggested by Tierney and Kadane [23] and Lindley [24] are commonly used and perform well. In the next subsections, the Tierney and Kadane technique and Lindley’s approximation are applied.

7.1. Tierney and Kadane Technique

By applying the Tierney and Kadane technique, we can approximate (28) into an explicit expression. First, two functions ψ α 0 , λ , β and ψ θ * α 0 , λ , β are defined as follows:
ψ α 0 , λ , β = l α 0 , λ , β + ρ α 0 , λ , β n ,
ψ θ * α 0 , λ , β = ψ α 0 , λ , β + ln g α 0 , λ , β n .
Then, (28) can be approximated as
g T K ^ α 0 , λ , β = Σ θ * Σ exp n ψ θ * α ^ 0 ψ * , λ ^ ψ * , β ^ ψ * ψ α ^ 0 ψ , λ ^ ψ , β ^ ψ ,
where α ^ 0 ψ , λ ^ ψ , β ^ ψ and α ^ 0 ψ * , λ ^ ψ * , β ^ ψ * maximize ψ α 0 , λ , β and ψ θ * α 0 , λ , β , respectively, and Σ and Σ θ * are the determinants of negative inverse Hessian of ψ α 0 , λ , β and ψ θ * α 0 , λ , β at α ^ 0 ψ , λ ^ ψ , β ^ ψ and α ^ 0 ψ * , λ ^ ψ * , β ^ ψ * , respectively.
From the (9) and (25), ψ α 0 , λ , β ignoring the constants can be expressed as
ψ α 0 , λ , β = 1 n ln α 0 a 1 1 + i = 1 k m i + ln λ a 2 1 + i = 1 k m i h i + ln β a 3 1 + i = 1 k m i + β 1 i = 1 k j = 1 m i ln t i j + i = 1 k j = 1 m i α 0 λ h i R i j + 1 1 e t i j β + i = 1 k j = 1 m i t i j β b 1 α 0 + b 2 λ + b 3 β .
To obtain α ^ 0 ψ , λ ^ ψ , β ^ ψ , the derivatives of ψ are derived
ψ α 0 = 1 n 1 α 0 a 1 1 + i = 1 k m i + i = 1 k j = 1 m i λ h i R i j + 1 1 e t i j β b 1 ,
ψ λ = 1 n 1 λ a 2 1 + i = 1 k m i h i + i = 1 k j = 1 m i α 0 h i λ h i R i j + 1 1 e t i j β b 2 ,
ψ β = 1 n 1 β a 3 1 + i = 1 k m i i = 1 k j = 1 m i α 0 λ h i e t i j β t i j β ln t i j R i j + 1 + i = 1 k j = 1 m i ln t i j + i = 1 k j = 1 m i t i j β ln t i j b 3 .
To calculate Σ , the second-order derivatives of ψ are derived
2 ψ α 0 2 = 1 n 1 α 0 2 a 1 1 + i = 1 k m i ,
2 ψ α 0 λ = 1 n i = 1 k j = 1 m i h i λ h i 1 R i j + 1 1 e t i j β ,
2 ψ α 0 β = 1 n i = 1 k j = 1 m i λ h i e t i j β t i j β ln t i j R i j + 1 ,
2 ψ λ 2 = 1 n 1 λ 2 a 2 1 + i = 1 k m i h i + i = 1 k j = 1 m i α 0 h i h i 1 λ h i 2 R i j + 1 1 e t i j β ,
2 ψ λ β = 1 n i = 1 k j = 1 m i α 0 h i λ h i 1 e t i j β t i j β ln t i j R i j + 1 ,
2 ψ β 2 = 1 n 1 β 2 a 3 1 + i = 1 k m i + i = 1 k j = 1 m i t i j β ln t i j 2 i = 1 k j = 1 m i α 0 λ h i ln t i j R i j + 1 e t i j β t i j β 2 ln t i j + e t i j β t i j β ln t i j .
Take solving α ^ 0 T K as an example. Let g α 0 , λ , β = α 0 . Then, the ψ α 0 * is given by
ψ α 0 * α 0 , λ , β = ψ + 1 n ln α 0 .
Correspondingly, the derivatives and second-order derivatives of ψ α 0 * are as follows, respectively,
ψ α 0 * α 0 = ψ α 0 + 1 n α 0 , ψ α 0 * λ = ψ λ , ψ α 0 * β = ψ β ,
2 ψ α 0 * α 0 2 = 2 ψ α 0 2 1 n α 0 2 , 2 ψ α 0 * λ 2 = 2 ψ λ 2 , 2 ψ α 0 * β 2 = 2 ψ β 2 , 2 ψ α 0 * α 0 λ = 2 ψ α 0 λ , 2 ψ α 0 * α 0 β = 2 ψ α 0 β , 2 ψ α 0 * λ β = 2 ψ λ β .
From the above derivation, we obtain the estimates α ^ 0 T K by
α ^ 0 T K = Σ α 0 * Σ exp n ψ α 0 * α ^ 0 ψ * , λ ^ ψ * , β ^ ψ * ψ α ^ 0 ψ , λ ^ ψ , β ^ ψ .
Similarly, by making g α 0 , λ , β = λ and g α 0 , λ , β = β , respectively, the estimates λ ^ T K and β ^ T K can be obtained by
λ ^ T K = Σ λ * Σ exp n ψ λ * α ^ 0 ψ * , λ ^ ψ * , β ^ ψ * ψ α ^ 0 ψ , λ ^ ψ , β ^ ψ ,
β ^ T K = Σ β * Σ exp n ψ β * α ^ 0 ψ * , λ ^ ψ * , β ^ ψ * ψ α ^ 0 ψ , λ ^ ψ , β ^ ψ .

7.2. Lindley’s Approximation

Jung and Chung [25] mentioned that the Lindley’s approximation can be utilized to the model with three parameters. Denote Θ = θ 1 , θ 2 , θ 3 = α 0 , λ , β . By applying the Lindley’s approximation, (28) becomes
g L A ^ α 0 , λ , β g Θ + 1 2 i = 1 3 j = 1 3 g i j Θ + 2 g i Θ ρ j σ i j + 1 2 i = 1 3 j = 1 3 k = 1 3 l = 1 3 l i j k g l Θ σ i j σ k l ,
where g i = g θ i , g i j = 2 g θ i θ j , ρ i = ρ θ i , l i j k = 3 l θ i θ j θ k and σ i j is the i , j th element of the matrix O T 1 Θ .
By making g Θ = α 0 , g Θ = λ and g Θ = β , respectively, the estimates α ^ 0 L A , λ ^ L A and β ^ L A can be calculated by
α ^ 0 L A = α ^ 0 + ρ 1 σ 11 + ρ 2 σ 12 + ρ 3 σ 13 + 1 2 i = 1 3 j = 1 3 k = 1 3 l i j k σ i j σ k 1 ,
λ ^ L A = λ ^ + ρ 1 σ 21 + ρ 2 σ 22 + ρ 3 σ 23 + 1 2 i = 1 3 j = 1 3 k = 1 3 l i j k σ i j σ k 2 ,
β ^ L A = β ^ + ρ 1 σ 31 + ρ 2 σ 32 + ρ 3 σ 33 + 1 2 i = 1 3 j = 1 3 k = 1 3 l i j k σ i j σ k 3 .

8. Optimal Design

In addition to the study of parameter estimation, this paper also includes the problem of selecting the optimal transformed stress level h i * i = 1 , , k of CSALT for the Chen distribution under progressive Type-II censoring. In CSALT, test units are assigned to run under different stress levels to better estimate the life characteristics of the units under the normal stress level. For different combinations of stress levels, the performance of parameter estimation will be different. Therefore, to improve the accuracy of product life characteristic estimation, the selection of stress levels is very important.
According to Nelson [2], the best test plan has two stress levels to produce the most precise estimation of life distribution parameters (minimum variance) under design stress levels. Therefore, we consider and discuss the case k = 2 when solving the optimal transformed stress level. Since the minimum transformed stress level is constant at h 1 = 1 , we only need to solve the optimal stress level of h 2 . Here, we investigate this problem according to the following two optimal criteria.

8.1. D-Optimality

D-optimality is very common in designing ALT, which maximizes the determinant of FIM. The determinant of FIM is known to be proportional to the determinant of the inverse of the asymptotic variance–covariance matrix, which is proportional to the total volume of the joint confidence interval of the parameter α 0 , λ , β . To estimate the life parameters of the unit under the normal stress level with greater accuracy, when the progressive censoring scheme is determined, the smaller the total volume of the joint confidence region of the parameter α 0 , λ , β , the better the estimation. Therefore, maximizing the determinant of the FIM within a suitable range can achieve this aim. Here, we continue to use the observed FIM O T Θ derived from the preceding section. The criterion can be expressed as follows
max O T α 0 , λ , β .

8.2. A-Optimality

Similarly, according to the relationship between the asymptotic variance–covariance matrix and the joint confidence region of the parameter α 0 , λ , β , we can also minimize the trace of the inverse of FIM to obtain the optimal solution. It provides a global indicator to measure the overall change of MLE. We still use observed FIM O T Θ to obtain the optimal transformed stress to estimate the life parameters of the unit under the normal stress level more accurately when the progressive censoring scheme is determined. The criterion can be expressed as follows
min t r O T 1 α 0 , λ , β .

9. Simulation Study

In this section, a numerical study is carried out with the Monte Carlo simulation using R software to assess the performance of the methods for point and interval estimation discussed in the preceding sections. Initially, we set the true values of the parameters to α 0 = 0.45 , λ = 2.0 and β = 0.7 . We also set the CSALT with two stress levels that h = 1 , 8.0 and k = 2 . Under both stress levels s 1 and s 2 , suppose that the number of assigned test units n i , ( i = 1 , 2 ), the number of failed test units m i , and the censoring schemes R i do not depend on i. We set three different sample sizes (30, 15), (45, 25) and (60, 35), and consider three different censoring schemes in each sample size. Among the censoring schemes, ( 0 * 7 , 15, 0 * 7 ) denotes that no unit is withdrawn for the first 7 times, then 15 units are withdrawn at one time, and no unit is withdrawn for the last 7 times. It is similar to other notations. By utilizing the algorithm from Balakrishnan and Sandhu [26], the progressive Type-II censored samples are generated.
For all point and interval estimation methods, the simulation is repeated 5000 times. Then, we take the expected values (EV) as the final point estimates and take the mean values of lower bounds (LB) and upper bounds (UB) as the final interval estimates. Additionally, we calculate the mean square error (MSE) for point estimation and the average interval length (AIL) and the coverage probability (COVP) for interval estimation to compare the performance of the estimation from various methods.
According to the log-likelihood function (9), we use the progressively censored samples to obtain the MLEs of the parameters in Table 1. In computing BEs, we set a 1 , a 2 , a 3 = 0.1 , 0.1 , 0.1 and b 1 , b 2 , b 3 = 0.2 , 0.2 , 0.2 as the values of hyperparameters for the prior distributions. Then, we apply the Tierney and Kadane technique and Lindley’s approximation to obtain BEs that are shown in Table 2 and Table 3, respectively. Further, Table 4 shows the corresponding 95 % asymptotic confidence intervals based on the asymptotic property of MLEs and the approximate variance of the parameters in the matrix O T 1 Θ . For the interval estimation, we also apply the parametric bootstrap method to find the 95 % bootstrap percentile confidence interval and the 95 % bootstrap-t confidence interval that are shown in Table 5 and Table 6, where we set B = 1000 . Finally, utilizing D- and A- optimal criteria, we obtain the average optimal transformed stress level for each progressive censoring scheme that are shown in Table 7 and Table 8.
From Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7 and Table 8, the following conclusions can be drawn.
(1)
According to Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6, the parameter estimation under different progressive censoring schemes all roughly follows the pattern of better performance with increasing sample size. Different censoring schemes have different performance in different estimation methods. As the censoring scheme of withdrawing one surviving unit for each unit observed to fail from the beginning of the test (for example, 1 * 15 ), the performance in point estimation is poor, especially when the sample size is relatively small and the estimation bias is relatively large. However, this censoring scheme has no obvious difference from other schemes in interval estimation.
(2)
According to the results of the point estimates in Table 1, Table 2 and Table 3, the estimated values are all close to the true values and, roughly, the MSE decreases as the sample size increases, so all the point estimation methods mentioned are valid. Based on the MSEs of the parameter estimators, the three methods perform similarly when the progressive censoring scheme is determined. For the estimation of α 0 , the BE obtained from the Lindley’s approximation performs better than BE obtained from the Tierney and Kadane technique. Additionally, BE obtained from the Tierney and Kadane technique performs slightly better than MLE. For λ estimation, BE obtained from the Lindley’s approximation performs best, and the BE obtained from the Tierney and Kadane technique and MLE perform similarly. For β estimation, BE obtained from the Lindley’s approximation performs best, followed by BE obtained from Lindley’s approximation and MLE. Overall, BE obtained from the Lindley’s approximation performs well for all three parameters.
(3)
According to the results of the interval estimation in Table 4, Table 5 and Table 6, generally speaking, the asymptotic confidence intervals perform best. The bootstrap percentile confidence intervals and the bootstrap-t confidence intervals perform similarly when the progressive censoring scheme is determined. For all three interval estimation methods, the AIL of λ is longer than those of α 0 and β . The COVP of the intervals of α 0 outperforms those of both λ and β in parametric bootstrap intervals.
(4)
According to Table 7 and Table 8 for the optimal transformed stress level, the values based on A-optimality are smaller than those based on D-optimality. For a similar censoring scheme, the results based on D-optimality become larger as the sample size increases. However, the results based on A-optimality are the opposite. In general, the optimal transformed stress levels do not differ significantly in each of the optimal criteria.

10. Real Data Analysis

A real data set on the breakdown times (in minutes) of a type of insulating fluid from Chapter Three of Nelson [2] is considered to illustrate the methods. The normal stress level is s 0 = 20 kilovolt (kV). In CSALT, the complete failure data under both stress levels are presented in Table 9. To ensure the fit of the data to the Chen distribution, we calculate the Kolmogorov–Smirnov statistics and p-values, also shown in Table 9. Furthermore, the plots of the empirical cumulative distribution function and the fitted cumulative distribution function are shown in Figure 9. It is clear that the data set fits the model well.
Based on the chemical reaction mechanism leading to the failure, the inverse power model is a suitable acceleration model. Therefore, we let ϕ ( s i ) = ln s i , ( i = 1 , 2 ). The value of h is then calculated as h = ( 1 , 1.44966 ) . Under both stress levels s 1 = 30 kV and s 2 = 36 kV, the number of assigned test units n i , i = 1 , 2 , the number of failed test units m i , the progressive Type-II censoring schemes R i and the censored data are shown in Table 10. The point estimates are shown in Table 11 and the intervals estimates are shown in Table 12. The number below each interval in Table 12 is the length of that interval. For Bayesian estimation, the results are calculated with non-informative priors in which the hyperparameters all tend to zero.
From Table 11 and Table 12, the following conclusions can be drawn:
(1)
The point estimates of α 0 and β obtained from the three different point estimation methods are similar, but the point estimates of λ are slightly different. For λ estimation, the MLE is maximal, followed by the BE obtained from the Lindley’s approximation and the BE obtained from the Tierney and Kadane technique;
(2)
In terms of AIL, the results obtained by the three interval estimation methods do not differ significantly. The bootstrap-t confidence intervals outperform slightly the asymptotic confidence intervals and the bootstrap percentile confidence intervals. This is consistent with the pattern that parametric bootstrap intervals outperform asymptotic intervals when the sample size is small.
(3)
For the bounds of the interval estimates of α 0 , the bootstrap percentile confidence intervals and the bootstrap-t confidence intervals are similar. For the interval estimates of λ , the upper and lower bounds of the asymptotic confidence intervals are significantly larger than those of the bootstrap percentile confidence intervals and the bootstrap-t confidence intervals. The bounds of the bootstrap-t confidence intervals are slightly larger than those of the bootstrap percentile confidence intervals. For the bounds on the interval estimates of β , the results obtained by the three methods are similar.

11. Conclusions

This paper considers statistical inference and optimal design of CSALT for the Chen distribution under progressive Type-II censoring. Specifically, maximum likelihood estimates are obtained utilizing the Newton–Raphson algorithm. Bayesian estimation under the squared error loss function is also considered and estimated by the Tierney and Kadane technique and Lindley’s approximation. At the same time, we establish confidence intervals for the lifetime distribution parameters and compare them. We construct the asymptotic intervals based on the observed Fisher information matrix and calculate the bootstrap percentile confidence intervals and the bootstrap-t confidence intervals to address the problem of small sample sizes. In addition, two optimal criteria are applied to determine the optimal transformed stress level that would provide more accurate estimates of the life parameters under the normal stress level. The performance of these estimation methods and the optimal criteria are compared using Monte Carlo simulation methods.
The Chen distribution possessing a bathtub-shaped hazard function is of great significance and practical importance in life experiments. Accelerated life tests and progressive censoring also play important roles. Further research into these models has great potential for estimating lifetime characteristics. This paper can also be extended by considering other distributions or other censoring schemes.

Author Contributions

Investigation, W.Z.; Supervision, W.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Project 202210004146 which was supported by Beijing Training Program of Innovation and Entrepreneurship for Undergraduates.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in [2].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McCool, J.I. Confidence limits for Weibull regression with censored data. IEEE Trans. Reliab. 1980, 29, 145–150. [Google Scholar] [CrossRef]
  2. Nelson, W.B. Accelerated Testing: Statistical Models, Test Plans, and Data Analysis; Wiley: New York, NY, USA, 1990. [Google Scholar]
  3. Miller, R.; Nelson, W. Optimum simple step-stress plans for accelerated life testing. IEEE Trans. Reliab. 1983, 32, 59–65. [Google Scholar] [CrossRef]
  4. Ismail, A.A. Reliability analysis under constant-stress partially accelerated life tests using hybrid censored data from Weibull distribution. Hacet. J. Math. Stat. 2016, 45, 181–193. [Google Scholar] [CrossRef]
  5. El-Din, M.M.M.; Abu-Youssef, S.E.; Ali, N.S.A.; El-Raheem, A.M.A. Optimal plans of constant-stress accelerated life tests for the Lindley distribution. J. Test. Eval. 2017, 45, 1463–1475. [Google Scholar]
  6. El-Din, M.M.M.; Abu-Youssef, S.E.; Ali, N.S.A.; El-Raheem, A.M.A. Parametric inference on step-stress accelerated life testing for the extension of exponential distribution under progressive type-II censoring. Commun. Stat. Appl. Methods 2016, 23, 269–285. [Google Scholar] [CrossRef] [Green Version]
  7. Chandra, N.; Khan, M.A. Analysis and optimum plan for 3-step step-stress accelerated life tests with Lomax model under progressive type-I censoring. Commun. Math. Stat. 2018, 6, 73–90. [Google Scholar]
  8. Abdel-Hamid, A.H.; Abushal, T.A. Inference on progressive-stress model for the exponentiated exponential distribution under type-II progressive hybrid censoring. J. Stat. Comput. Simul. 2015, 85, 1165–1186. [Google Scholar] [CrossRef]
  9. Balakrishnan, N.; Aggarwala, R. Progressive Censoring: Theory, Methods, and Applications; Birkhäuser: Boston, MA, USA, 2000. [Google Scholar]
  10. Balakrishnan, N. Progressive censoring methodology: An appraisal. Test 2007, 16, 211–259. [Google Scholar] [CrossRef]
  11. Jaheen, Z.F.; Moustafa, H.M.; Abd El-Monem, G.H. Bayes inference in constant partially accelerated life tests for the generalized exponential distribution with progressive censoring. Commun. Stat.-Theory Methods 2014, 43, 2973–2988. [Google Scholar] [CrossRef]
  12. Basak, I.; Balakrishnan, N. Prediction of censored exponential lifetimes in a simple step-stress model under progressive Type II censoring. Comput. Stat. 2017, 32, 1665–1687. [Google Scholar] [CrossRef]
  13. Abd El-Raheem, A.M. Inference and optimal design of multiple constant-stress testing for generalized half-normal distribution under type-II progressive censoring. J. Stat. Comput. Simul. 2019, 89, 3075–3104. [Google Scholar] [CrossRef]
  14. Algarni, A.; Almarashi, A.M.; Okasha, H.; Ng, H.K.T. E-bayesian estimation of Chen distribution based on type-I censoring scheme. Entropy 2020, 22, 636. [Google Scholar] [CrossRef] [PubMed]
  15. Chen, Z. A new two-parameter lifetime distribution with bathtub shape or increasing failure rate function. Stat. Probab. Lett. 2000, 49, 155–161. [Google Scholar] [CrossRef]
  16. Hjorth, U. A reliability distribution with increasing, decreasing, constant and bathtub-shaped failure rates. Technometrics 1980, 22, 99–107. [Google Scholar] [CrossRef]
  17. Olkin, I. Life distributions: A brief discussion. Commun. Stat.-Simul. Comput. 2016, 45, 1489–1498. [Google Scholar] [CrossRef]
  18. Sarhan, A.M.; Hamilton, D.C.; Smith, B. Parameter estimation for a two-parameter bathtub-shaped lifetime distribution. Appl. Math. Model. 2012, 36, 5380–5392. [Google Scholar] [CrossRef]
  19. Rastogi, M.K.; Tripathi, Y.M. Estimation using hybrid censored data from a two-parameter distribution with bathtub shape. Comput. Stat. Data Anal. 2013, 67, 268–281. [Google Scholar] [CrossRef]
  20. Ahmed, E.A. Bayesian estimation based on progressive Type-II censoring from two-parameter bathtub-shaped lifetime model: An Markov chain Monte Carlo approach. J. Appl. Stat. 2014, 41, 752–768. [Google Scholar] [CrossRef]
  21. Ahmed, E.A.; Alhussain, Z.A.; Salah, M.M.; Ahmed, H.H.; Eliwa, M.S. Inference of progressively type-II censored competing risks data from Chen distribution with an application. J. Appl. Stat. 2020, 47, 2492–2524. [Google Scholar] [CrossRef]
  22. MacDonald, I.L. Does Newton–Raphson really fail? Stat. Methods Med. Res. 2014, 23, 308–311. [Google Scholar] [CrossRef]
  23. Tierney, L.; Kadane, J.B. Accurate approximations for posterior moments and marginal densities. J. Am. Stat. Assoc. 1986, 81, 82–86. [Google Scholar] [CrossRef]
  24. Lindley, D.V. Approximate Bayesian methods. Trab. Estadística Y Investig. Oper. 1980, 31, 223–245. [Google Scholar] [CrossRef]
  25. Jung, M.; Chung, Y. Bayesian inference of three-parameter bathtub-shaped lifetime distribution. Commun. Stat.-Theory Methods 2018, 47, 4229–4241. [Google Scholar] [CrossRef]
  26. Balakrishnan, N.; Sandhu, R.A. A simple simulational algorithm for generating progressive type-II censored samples. Am. Stat. 1995, 49, 229–230. [Google Scholar]
Figure 1. The pdf when 0 < α < 1 , 0 < β < 1 .
Figure 1. The pdf when 0 < α < 1 , 0 < β < 1 .
Mathematics 10 01609 g001
Figure 2. The pdf when α 1 , 0 < β < 1 .
Figure 2. The pdf when α 1 , 0 < β < 1 .
Mathematics 10 01609 g002
Figure 3. The pdf when β = 1 .
Figure 3. The pdf when β = 1 .
Mathematics 10 01609 g003
Figure 4. The pdf when β > 1 .
Figure 4. The pdf when β > 1 .
Mathematics 10 01609 g004
Figure 5. The hf when 0 < β < 1 .
Figure 5. The hf when 0 < β < 1 .
Mathematics 10 01609 g005
Figure 6. The hf when β 1 .
Figure 6. The hf when β 1 .
Mathematics 10 01609 g006
Figure 7. The influence of α on the mean of the Chen distribution when β is fixed.
Figure 7. The influence of α on the mean of the Chen distribution when β is fixed.
Mathematics 10 01609 g007
Figure 8. For n i assigned units, progressive Type-II censoring scheme R i 1 , R i 2 , , R i m i at stress level s i , i = 1 , 2 , k , in constant-stress accelerated life testing.
Figure 8. For n i assigned units, progressive Type-II censoring scheme R i 1 , R i 2 , , R i m i at stress level s i , i = 1 , 2 , k , in constant-stress accelerated life testing.
Mathematics 10 01609 g008
Figure 9. The empirical cumulative distribution function and the fitted cumulative distribution function for failure data under stress levels s 1 = 30 kV (Left, data 1) and s 2 = 36 kV (Right, data 2).
Figure 9. The empirical cumulative distribution function and the fitted cumulative distribution function for failure data under stress levels s 1 = 30 kV (Left, data 1) and s 2 = 36 kV (Right, data 2).
Mathematics 10 01609 g009
Table 1. The expected values (EV) and mean square errors (MSE) of the maximum likelihood estimates for different sample sizes n i , m i and different censoring schemes R i .
Table 1. The expected values (EV) and mean square errors (MSE) of the maximum likelihood estimates for different sample sizes n i , m i and different censoring schemes R i .
n i , m i R i EV α ^ 0 MSE α ^ 0 EV λ ^ MSE λ ^ EV β ^ MSE β ^
30 , 15 ( 0 7 , 15, 0 7 ) 0.4933 0.0321 2.1014 0.0736 0.7474 0.0142
( 0 5 , 3 5 , 0 5 ) 0.4946 0.0324 2.1050 0.0786 0.7484 0.0151
( 1 15 ) 0.4967 0.0367 2.1238 0.0932 0.7567 0.0182
45 , 25 ( 0 12 , 20, 0 12 ) 0.4728 0.0144 2.0573 0.0330 0.7281 0.0071
( 0 10 , 4 5 , 0 10 ) 0.4729 0.0146 2.0583 0.0368 0.7281 0.0076
( 1 20 , 0 5 ) 0.4703 0.0139 2.0595 0.0355 0.7281 0.0074
60 , 35 ( 0 17 , 25, 0 17 ) 0.4643 0.0096 2.0456 0.0235 0.7217 0.0049
( 0 15 , 5 5 , 0 15 ) 0.4671 0.0100 2.0409 0.0224 0.7205 0.0048
( 1 25 , 0 10 ) 0.4650 0.0098 2.0402 0.0228 0.7199 0.0050
Table 2. The expected values (EV) and mean square errors (MSE) of the Bayes estimates using the Tierney and Kadane technique for different sample sizes n i , m i and different censoring schemes R i .
Table 2. The expected values (EV) and mean square errors (MSE) of the Bayes estimates using the Tierney and Kadane technique for different sample sizes n i , m i and different censoring schemes R i .
n i , m i R i EV α ^ 0 MSE α ^ 0 EV λ ^ MSE λ ^ EV β ^ MSE β ^
30 , 15 ( 0 7 , 15, 0 7 )0.49570.02912.09370.07150.73860.0130
( 0 5 , 3 5 , 0 5 )0.49570.02992.10030.07790.73970.0139
( 1 15 )0.49530.03152.11870.09020.74760.0166
45 , 25 ( 0 12 , 20, 0 12 )0.47440.01472.05410.03300.72290.0067
( 0 10 , 4 5 , 0 10 )0.47470.01502.05630.03500.72360.0069
( 1 20 , 0 5 )0.47500.01422.05560.03450.72440.0071
60 , 35 ( 0 17 , 25, 0 17 )0.46690.00932.03680.02200.71570.0047
( 0 15 , 5 5 , 0 15 )0.46640.00962.03640.02230.71540.0046
( 1 25 , 0 10 )0.46640.00982.04010.02350.71720.0049
Table 3. The expected values (EV) and mean square errors (MSE) of the Bayes estimates using the Lindley’s approximation for different sample sizes n i , m i and different censoring schemes R i .
Table 3. The expected values (EV) and mean square errors (MSE) of the Bayes estimates using the Lindley’s approximation for different sample sizes n i , m i and different censoring schemes R i .
n i , m i R i EV α ^ 0 MSE α ^ 0 EV λ ^ MSE λ ^ EV β ^ MSE β ^
30 , 15 ( 0 7 , 15, 0 7 )0.48690.02702.04240.04910.71170.0098
( 0 5 , 3 5 , 0 5 )0.48250.02652.04250.05140.70990.0105
( 1 15 )0.48260.02582.04620.05930.71190.0125
45 , 25 ( 0 12 , 20, 0 12 )0.46990.01262.02900.02690.70880.0058
( 0 10 , 4 5 , 0 10 )0.47040.01382.02900.02850.70890.0060
( 1 20 , 0 5 )0.47080.01332.02860.03040.70850.0065
60 , 35 ( 0 17 , 25, 0 17 )0.46320.00882.01900.01920.70590.0042
( 0 15 , 5 5 , 0 15 )0.46480.00872.01890.01960.70540.0042
( 1 25 , 0 10 )0.46600.00962.01710.01940.70510.0042
Table 4. The lower bounds (LB), upper bounds (UB), average interval lengths (AIL) and coverage probabilities (COVP) of the 95 % asymptotic confidence intervals for different sample sizes n i , m i and different censoring schemes R i .
Table 4. The lower bounds (LB), upper bounds (UB), average interval lengths (AIL) and coverage probabilities (COVP) of the 95 % asymptotic confidence intervals for different sample sizes n i , m i and different censoring schemes R i .
n i , m i R i ParameterLBUBAILCOVP
30 , 15 ( 0 7 , 15, 0 7 ) α 0 0.20150.78590.58440.9422
λ 1.66092.53890.87790.9556
β 0.55030.94250.39220.9403
( 0 5 , 3 5 , 0 5 ) α 0 0.20200.78760.58560.9451
λ 1.65532.55240.89710.9595
β 0.54700.94950.40250.9436
( 1 15 ) α 0 0.20180.78890.58710.9476
λ 1.62402.62621.00220.9649
β 0.52850.98630.45780.9404
45 , 25 ( 0 12 , 20, 0 12 ) α 0 0.25650.68710.43060.9467
λ 1.73242.38610.65380.9535
β 0.57840.87800.29950.9465
( 0 10 , 4 5 , 0 10 ) α 0 0.25690.68780.43090.9461
λ 1.72962.38830.65870.9546
β 0.57680.87930.30260.9456
( 1 20 , 0 5 ) α 0 0.25530.68400.42870.9446
λ 1.72702.39950.67250.9546
β 0.57420.88470.31060.9470
60 , 35 ( 0 17 , 25, 0 17 ) α 0 0.28600.64480.35880.9461
λ 1.77002.31330.54330.9535
β 0.59400.84540.25140.9474
( 0 15 , 5 5 , 0 15 ) α 0 0.28560.64370.35810.9444
λ 1.76832.31440.54610.9523
β 0.59350.84650.25310.9486
( 1 25 , 0 10 ) α 0 0.28530.64470.35940.9422
λ 1.76752.31440.54680.9562
β 0.59270.84710.25440.9491
Table 5. The lower bounds (LB), upper bounds (UB), average interval lengths (AIL) and coverage probabilities (COVP) of the 95 % bootstrap percentile confidence interval for different sample sizes n i , m i and different censoring schemes R i .
Table 5. The lower bounds (LB), upper bounds (UB), average interval lengths (AIL) and coverage probabilities (COVP) of the 95 % bootstrap percentile confidence interval for different sample sizes n i , m i and different censoring schemes R i .
n i , m i R i ParameterLBUBAILCOVP
30 , 15 ( 0 7 , 15, 0 7 ) α 0 0.27481.00150.72670.9180
λ 1.80412.91661.11250.8520
β 0.60931.06120.45190.8560
( 0 5 , 3 5 , 0 5 ) α 0 0.27621.01050.73430.9320
λ 1.79542.93241.13700.8680
β 0.60501.06850.46350.8540
( 1 15 ) α 0 0.27601.02370.74760.9260
λ 1.79393.08831.29450.8800
β 0.59611.12000.52380.8800
45 , 25 ( 0 12 , 20, 0 12 ) α 0 0.30400.79590.49190.9260
λ 1.82002.56740.74740.8980
β 0.61710.94350.32650.8900
( 0 10 , 4 5 , 0 10 ) α 0 0.30560.80030.49470.9220
λ 1.81252.56070.74820.9040
β 0.61440.94420.32980.9000
( 1 20 , 0 5 ) α 0 0.31340.81930.50590.9120
λ 1.79632.54480.74840.9180
β 0.60680.94140.33460.8960
60 , 35 ( 0 17 , 25, 0 17 ) α 0 0.31880.70940.39060.9460
λ 1.82872.42190.59320.9180
β 0.61880.88540.26660.9060
( 0 15 , 5 5 , 0 15 ) α 0 0.31840.70820.38980.9420
λ 1.83682.44020.60350.9040
β 0.62510.89570.27060.8820
( 1 25 , 0 10 ) α 0 0.31440.70220.38780.9200
λ 1.83672.44000.60320.9020
β 0.62190.89170.26990.8880
Table 6. The lower bounds (LB), upper bounds (UB), average interval lengths (AIL) and coverage probabilities (COVP) of the 95 % bootstrap-t confidence interval for different sample sizes n i , m i and different censoring schemes R i .
Table 6. The lower bounds (LB), upper bounds (UB), average interval lengths (AIL) and coverage probabilities (COVP) of the 95 % bootstrap-t confidence interval for different sample sizes n i , m i and different censoring schemes R i .
n i , m i R i ParameterLBUBAILCOVP
30 , 15 ( 0 7 , 15, 0 7 ) α 0 0.27791.00880.73090.9180
λ 1.80042.89221.09190.8660
β 0.60981.06240.45260.8460
( 0 5 , 3 5 , 0 5 ) α 0 0.28031.02280.74250.9140
λ 1.79862.92781.12930.8640
β 0.60911.07290.46390.8580
( 1 15 ) α 0 0.28621.05850.77230.9260
λ 1.78773.07741.28970.8680
β 0.59811.12470.52660.8540
45 , 25 ( 0 12 , 20, 0 12 ) α 0 0.29910.78430.48510.9300
λ 1.82512.57500.74990.8900
β 0.62030.94610.32570.8760
( 0 10 , 4 5 , 0 10 ) α 0 0.30630.80310.49680.9240
λ 1.81702.56780.75080.9180
β 0.61510.94460.32950.9020
( 1 20 , 0 5 ) α 0 0.31350.81940.50590.9180
λ 1.80882.56590.75710.9100
β 0.61230.94850.33620.8860
60 , 35 ( 0 17 , 25, 0 17 ) α 0 0.32550.72130.39570.9280
λ 1.82542.38760.56230.9060
β 0.61940.86990.25050.9120
( 0 15 , 5 5 , 0 15 ) α 0 0.31820.70990.39170.9460
λ 1.83062.42910.59840.9020
β 0.61990.88860.26870.9040
( 1 25 , 0 10 ) α 0 0.31730.70770.39040.9260
λ 1.83392.43160.59770.9160
β 0.62270.89270.27000.9120
Table 7. The average optimal transformed stress level h k using D-optimality for different sample sizes n i , m i and different censoring schemes R i .
Table 7. The average optimal transformed stress level h k using D-optimality for different sample sizes n i , m i and different censoring schemes R i .
n i , m i R i h k R i h k R i h k
30 , 15 ( 0 7 , 15, 0 7 )8.2396( 0 5 , 3 5 , 0 5 )8.1901( 1 15 )7.7286
45 , 25 ( 0 12 , 20, 0 12 )8.6145( 0 10 , 4 5 , 0 10 )8.5932( 1 20 , 0 5 )8.5267
60 , 35 ( 0 17 , 25, 0 17 )8.7662( 0 15 , 5 5 , 0 15 )8.7874( 1 25 , 0 10 )8.7212
Table 8. The average optimal transformed stress level h k using A-optimality for different sample sizes n i , m i and different censoring schemes R i .
Table 8. The average optimal transformed stress level h k using A-optimality for different sample sizes n i , m i and different censoring schemes R i .
n i , m i R i h k R i h k R i h k
30 , 15 ( 0 7 , 15, 0 7 )6.7736( 0 5 , 3 5 , 0 5 )6.7567( 1 15 )6.7330
45 , 25 ( 0 12 , 20, 0 12 )6.7087( 0 10 , 4 5 , 0 10 )6.7015( 1 20 , 0 5 )6.6857
60 , 35 ( 0 17 , 25, 0 17 )6.6858( 0 15 , 5 5 , 0 15 )6.6738( 1 25 , 0 10 )6.6786
Table 9. The complete failure data and the fitting results.
Table 9. The complete failure data and the fitting results.
Stress LevelData SetK-S Statisticp-Value
30 kV7.74, 17.05, 20.46, 21.02, 22.66, 43.40, 47.30,
139.07, 144.12, 175.88, 194.90
0.22030.5858
36 kV0.35, 0.59, 0.96, 0.99, 1.69, 1.97, 2.07, 2.58,
2.71, 2.90, 3.67, 3.99, 5.35, 13.77, 25.50
0.21790.4154
Table 10. The constant-stress progressive Type-II censored data.
Table 10. The constant-stress progressive Type-II censored data.
Stress Level h i n i , m i R i Censored Data
30 kV1 11 , 10 ( 0 4 , 1, 0 5 )7.74, 17.05, 20.46, 21.02, 22.66, 47.30, 139.07,
144.12, 175.88, 194.90
36 kV1.44966 15 , 14 ( 0 6 , 1, 0 7 )0.35, 0.59, 0.96, 0.99, 1.69, 1.97, 2.07, 2.58,
2.90, 3.67, 3.99, 5.35, 13.77, 25.50
Table 11. The maximum likelihood estimates (MLE) and the Bayes estimates using the Tierney and Kadane technique ( B E T K ) and Lindley’s approximation ( B E L A ).
Table 11. The maximum likelihood estimates (MLE) and the Bayes estimates using the Tierney and Kadane technique ( B E T K ) and Lindley’s approximation ( B E L A ).
ParameterMLE BE TK BE LA
α 0 0.00250.00320.0026
λ 22.806320.431021.3963
β 0.26390.26050.2673
Table 12. The 95% asymptotic confidence intervals (ACI), the 95% bootstrap percentile confidence intervals (Boot-p) and the 95% bootstrap-t confidence intervals (Boot-t).
Table 12. The 95% asymptotic confidence intervals (ACI), the 95% bootstrap percentile confidence intervals (Boot-p) and the 95% bootstrap-t confidence intervals (Boot-t).
ParameterACIBoot-pBoot-t
α 0 (0.0012, 0.0039)
0.0027
(0.0025, 0.0218)
0.0193
(0.0025, 0.0043)
0.0018
λ (9.9884, 35.6241)
25.6356
(4.8530, 28.3640)
23.5111
(6.0894, 28.2857)
22.1964
β (0.2292, 0.2986)
0.0694
(0.2162, 0.2848)
0.0686
(0.2259, 0.2812)
0.0553
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, W.; Gui, W. Statistical Inference and Optimal Design of Accelerated Life Testing for the Chen Distribution under Progressive Type-II Censoring. Mathematics 2022, 10, 1609. https://doi.org/10.3390/math10091609

AMA Style

Zhang W, Gui W. Statistical Inference and Optimal Design of Accelerated Life Testing for the Chen Distribution under Progressive Type-II Censoring. Mathematics. 2022; 10(9):1609. https://doi.org/10.3390/math10091609

Chicago/Turabian Style

Zhang, Wenjie, and Wenhao Gui. 2022. "Statistical Inference and Optimal Design of Accelerated Life Testing for the Chen Distribution under Progressive Type-II Censoring" Mathematics 10, no. 9: 1609. https://doi.org/10.3390/math10091609

APA Style

Zhang, W., & Gui, W. (2022). Statistical Inference and Optimal Design of Accelerated Life Testing for the Chen Distribution under Progressive Type-II Censoring. Mathematics, 10(9), 1609. https://doi.org/10.3390/math10091609

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop