Next Article in Journal
Risk Analysis: Changing the Story with the Statistical Stochastic Process and VaR
Previous Article in Journal
An Analytic Solution for 2D Heat Conduction Problems with General Dirichlet Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inference for Partially Accelerated Life Test from a Bathtub-Shaped Lifetime Distribution with Progressive Censoring

1
School of Mathematics, Yunnan Normal University, Kunming 650500, China
2
Department of Mathematics, Indian Institute of Technology Patna, Bihta 800013, India
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(5), 417; https://doi.org/10.3390/axioms12050417
Submission received: 22 February 2023 / Revised: 20 April 2023 / Accepted: 22 April 2023 / Published: 24 April 2023

Abstract

:
The analysis of the constant-stress partially accelerated life test was considered under progressive Type-II censoring when the lifetime of the products follows a two-parameter bathtub-shaped distribution. The maximum likelihood estimates of the unknown parameters were established, where the expectation–maximization iterative solution is proposed for the estimation. The approximate confidence intervals were also constructed based on asymptotic theory via the Fisher information matrix. For comparison purposes, the bootstrap (i.e., Studentized-t and percentile) confidence intervals of the unknown parameters were also obtained. Finally, simulation studies and a real-life data example are presented to examine the performance of the different results.

1. Introduction

With the development of manufacturing design technology, modern products have the characteristics of high reliability and a long lifespan. Under these usage conditions, life testing becomes increasingly difficult within a reasonable period due to the time and cost limitations. To conduct life testing in an effective way, the accelerated life test (ALT) has been introduced in practice, which provides a low-cost and quick method to obtain the failure information. In the ALT, in order to obtain the failure data, products are tested either under a harsher environment or under more intensive usage than the usual use conditions, and the common stress factors include the temperature, pressure, and voltage, among others. Data collected in such accelerated conditions are then extrapolated through a physically appropriate statistical model to estimate the lifetime distribution under normal use conditions. Generally, there are mainly three types of acceleration tests, namely the constant-stress ALT, the step-stress ALT, and the progressive-stress ALT. All these types of ALTs have been studied by a number of authors [1,2,3,4,5]; for more details about the ALT, one can refer to the monograph by Nelson [6].
Under accelerating life testing, when the acceleration factor has an unknown value, the partially accelerated life test (PALT) is usually conducted as a reasonable alternative to the life test, where the units are investigated under both accelerated and regular  use conditions. There are mainly two types of PALTs, the constant-stress PALT and the step-stress PALT. In the constant-stress PALT, the units are divided into two groups; one group is assigned to run under the regular use conditions, while the other group is tested under accelerated conditions. For the step-stress PALT, the units are initially run under the regularuse conditions for a pre-specified period of time, and if a test unit does not fail for the specified time, it is then run under accelerated conditions until failure occurs or the observation is censored. PALTs have been discussed by several authors. For instance, Abdel-Hamid [7] considered the parameter estimation of the Burr Type-XII distribution in a constant-stress PALT for progressively Type-II-censored data. Abdel-Hamid and Al-Hussaini [8] studied the estimation problem in the step-stress PALT when the lifetime of the tested units under the regular use conditions followed a finite mixture of a general class of distributions. Abdel-Ghaly et al. [9] investigated the maximum likelihood estimation (MLE) method to estimate the parameters of the Weibull distribution in the step-stress PALT. A similar problem was also studied by Ismail [10] under the progressive hybrid Type-II censoring scheme. Cheng and Wang [11] discussed the MLE of the Burr XII distribution under the constant-stress PALT when multiple censored data were observed.
In practice, the lifetime distributions with bathtub-shaped hazard functions have attracted the interest of many authors and provide appropriate conceptual models for some electronic and mechanical products, as well as the lifetimes of humans. Chen [12] proposed a two-parameter distribution with a bathtub-shaped or increasing failure rate function. The survival function (SF), cumulative distribution function (CDF), probability density function (PDF), and hazard rate function (HRF) of the Chen distribution are given by
S 1 ( t ) = e α ( 1 e t β ) , t > 0 , F 1 ( t ) = 1 e α ( 1 e t β ) , t > 0 , f 1 ( t ) = α β t β 1 exp α ( 1 e t β ) + t β , t > 0 , h 1 ( t ) = α β t β 1 e t β , t > 0 ,
where α > 0 is the scale parameter and β > 0 is the shape parameter, respectively. In practice, the applicability of a model may partially be attributed to the fact that its reliability, hazard rate, and probability density functions all have nice expressions. Due to its flexible structural properties and practical significance, the Chen distribution has found wide applications in life test studies. It has been observed that the Chen distribution has a bathtub-shaped hazard function when 0 < β < 1 , and when β 1 , it features an increasing hazard rate function. The  case α = 1 corresponds to the exponential power distribution. Several authors have discussed the Chen distribution in different cases. For example, based on Type-II-censored samples, Chen [12] constructed exact confidence intervals for the shape parameter and also obtained exact confidence regions for both model parameters. Rastogi and Tripathi [13] discussed the parameter estimation problem of the Chen distribution for hybrid censored data. Wu [14] investigated the MLE method to estimate the parameters of the Chen distribution under progressive Type-II censoring and also derived exact confidence intervals and confidence regions for the related parameters. Ahmed [15] presented the Bayesian approach to estimate the parameters of the Chen distribution for progressively Type-II-censored data. Elshahhat and Rastogi [16] considered the Bayesian life analysis of a generalized Chen’s population under progressive censoring.
In life testing, reliability analysis, and other related fields, censoring is a very common phenomenon, and the experiments are often terminated before all units fail due to the cost and time considerations. In such cases, the exact failure times are known for only a portion of the units under study. The most-common censoring schemes are Type-I and Type-II censoring, which, however, can only remove units at the termination point, lacking flexibility in practical life tests. Therefore, progressive censoring is further proposed to conduct the tests, which allows units to be removed at different testing stages. The strategy of progressive censoring is vital to planning duration experiments in the field of reliability and lifetime analysis and includes progressive Type-I and Type-II censoring as its special population cases. For details about the censoring scheme, the reader can refer to the recent review paper of Balakrishnan [17] and the monograph of Balakrishnan and Aggarwala [18], as well as the references therein. In this paper, we considered the constant-stress PALT applied to units whose lifetime under the use conditions was assumed to be the Chen distribution under a progressive Type-II censoring scheme. It is worth mentioning that, although there are many discussions that have focused on the inference of the Chen distribution such as the previously mentioned ones and others (e.g., Zhang and Gui [19], Soliman et al. [20]), the inference for the partially constant-stress accelerated life test of the Chen distribution has not been discussed in the literature. In addition, except for the traditional likelihood-based inferential approach, another aim of this paper was that the expectation–maximization estimation for the Chen model be proposed under the constant-stress PALT.
The rest of this article is organized as follows. Section 2 provides a brief description of the constant-stress PALT and some basic assumptions. Section 3 deals with the estimation problem of the MLEs of the parameters. The corresponding confidence intervals are proposed in Section 4, and this section also discusses two parametric bootstrap confidence intervals. An illustration example and Monte Carlo simulation studies are presented in Section 5 to investigate the performance the proposed results. Finally, some concluding remarks are addressed in Section 6.

2. Model Description and Basic Assumptions

2.1. Model Description

In the constant-stress PALT, assuming that there are n identical test units, n 1 units are randomly selected from n test units for testing in the regular use conditions and the remaining n 2 = n n 1 units are tested under accelerated stress. Progressive Type-II censoring is applied as follows. For  j = 1 , 2 , when the first failure, say T j 1 , has occurred, R j 1 units are randomly removed from the remaining n j 1 surviving units. When the second failure, say X j 2 , has occurred, R j 2 units from the remaining n j 2 R j 1 units are randomly removed. The test proceeds at different stress levels until the m j th failure, say T j m j , has occurred, at  which time all remaining R j m j = n j m j k = 1 m j 1 R j m j units are withdrawn and the constant-stress PALT is terminated. In our discussion, the values of the censoring numbers R j i , j = 1 , 2 , i = 1 , 2 , , m j were predetermined with m j < n j . It is noted that the complete sample and conventional Type-II censoring are special cases of the progressive Type-II censoring scheme.
In our study, the progressively Type-II-censored samples T j 1 < T j 2 < < T j m j , j = 1 , 2 were from two testing conditions, of  which the CDFs and the PDFs are F j ( t ) and f j ( t ) , respectively, with censoring schemes R j = ( R j 1 , R j 2 , , R j m j ) . Denote t j 1 < t j 2 < < t j m j as the observed values of T j 1 , T j 2 , , T j m j ; the joint PDF based on the two progressively Type-II-censored samples can be expressed as
L ( α , β , λ ; t ) = j = 1 2 C j i = 1 m j f j ( t j i ) ( 1 F j ( t j i ) ) R j i ,
where t = ( t 1 , t 2 ) , t j = ( t j 1 , t j 2 , , t j m j ) , j = 1 , 2 and C j = n j i = 1 m j 1 [ n j k = 1 i ( 1 + R j k ) ] .

2.2. Basic Assumptions

1.
The lifetime of a unit tested under the regular use conditions follows a Chen distribution with the CDF and PDF presented in (1).
2.
The hazard rate of a unit tested under accelerated conditions is given by h 2 ( t ) = λ h 1 ( t ) , where λ is an acceleration factor satisfying λ > 1 . Therefore, the HRF, SF, CDF, and PDF under the accelerated conditions are given by
h 2 ( t ) = α β λ t β 1 e t β , t > 0 , S 2 ( t ) = exp 0 t h 2 ( u ) d u = exp α λ ( 1 e t β ) , t > 0 , F 2 ( t ) = 1 e α λ ( 1 e t β ) , t > 0 , f 2 ( t ) = α β λ t β 1 exp α λ ( 1 e t β ) + t β , t > 0 .
3.
There are n j , j = 1 , 2 units allocated under the regular use conditions with j = 1 and the accelerated conditions with j = 2 , of which the lifetimes T i j , i = 1 , 2 , , n j , j = 1 , 2 are mutually independent.

3. Maximum Likelihood Estimation

Based on (1)–(3), the log-likelihood function of L ( α , β , λ ; t ) , say , can be expressed as
( m 1 + m 2 ) log ( α β ) + m 2 log ( λ ) + β j = 1 2 i = 1 m j log t j i + j = 1 2 i = 1 m j t j i β   + α j = 1 2 i = 1 m j λ j 1 ( 1 + R j i ) ( 1 e t j i β ) ,
where the notation ∝ means “be proportional to” without additive constant terms.
Taking the derivatives of with respect to α , β , λ and equating them to zero, one can derive the MLEs of the parameters α , β , and λ from the following likelihood equations:
α = 0 , β = 0 , and λ = 0 ,
where
α = α = m 1 + m 2 α + j = 1 2 i = 1 m j λ j 1 ( 1 + R j i ) ( 1 e t j i β ) , β = β = m 1 + m 2 β + j = 1 2 i = 1 m j log t j i + j = 1 2 i = 1 m j t j i β log t j i α j = 1 2 i = 1 m j λ j 1 ( 1 + R j i ) e t j i β t j i β log t j i , λ = λ = m 2 λ + α i = 1 m 2 ( 1 + R 2 i ) ( 1 e t 2 i β ) .
Based on the likelihood equations α = 0 and λ = 0 in (5), one has
α = m 1 i = 1 m 1 ( 1 + R 1 i ) ( e t 1 i β 1 ) , λ = m 2 i = 1 m 1 ( 1 + R 1 i ) ( e t 1 i β 1 ) m 1 i = 1 m 2 ( 1 + R 2 i ) ( e t 2 i β 1 ) ,
which implies that
m 1 + m 2 β + j = 1 2 i = 1 m j ( 1 + t j i β ) log t j i j = 1 2 m j i = 1 m j ( 1 + R j i ) e t j i β t j i β log t j i i = 1 m j ( 1 + R j i ) ( e t j i β 1 ) = 0 .
Therefore, using numerical iterative programs such as the Newton–Raphson algorithm to solve the nonlinear Equation (7), the MLE of β can be calculated. Further, by substituting the MLE of β into (6), the MLEs of α and λ can be obtained. Alternatively, one can also use the expectation–maximization (EM) algorithm to deduce the MLEs of the unknown parameters. The EM algorithm is extensively used to the iterative computation of maximum likelihood estimates and is very useful in a variety of fields such as survival analysis, reliability theory, and other fields. The EM iteration alternates between performing an expectation step or E-step and a maximization step or M-step. The E-step creates a function of the expectation of the log-likelihood of the current estimation evaluation using the parameters, and the M-step calculates the parameters that maximize the expected log-likelihood found in the E-step. For more details, one can refer to the work of Dempster et al. [21], where the authors first introduced the EM algorithm to handle some missing or incomplete data situations, and the monograph by McLachlan and Krishnan [22], as well as the references therein.
In the constant-stress PALT, the progressively Type-II-censored samples can be viewed as an incomplete dataset in each life test stage. Therefore, the EM algorithm will provide a good alternative to the conventional iterative method in the process of numerically computing the MLEs.
For j = 1 , 2 , let X = ( X j 1 , X j 2 , , X j m j ) with X j k = ( X j k 1 , X j k 2 , , X j k R j k ) , k = 1 , 2 , , m j , represent the censored data under the normal use and accelerated conditions, respectively. We treated the censored observations as missing data. Thus, the combination of ( T , X ) forms the complete constant-stress PALT failure dataset, for which the likelihood function can be expressed as
L ( α , β , λ ; t ) = j = 1 2 i = 1 m j f j ( t j i ) · j = 1 2 k = 1 m j s = 1 R j k f j ( x j k s ) .
The log-likelihood function of L ( α , β , λ ; t ) , say , can be expressed as
m 1 + m 2 + j = 1 2 k = 1 m j R j k ln α + m 1 + m 2 + j = 1 2 k = 1 m j R j k ln β + m 2 + k = 1 m 2 R 2 k ln λ + β j = 1 2 i = 1 m j ln t j i + j = 1 2 k = 1 m j s = 1 R j k ln x j k s + j = 1 2 i = 1 m j [ α λ j 1 ( 1 e t j i β ) + t j i β ] + j = 1 2 k = 1 m j s = 1 R j k [ α λ j 1 ( 1 e ( x j k s ) β ) + ( x j k s ) β ] .
The MLEs of the parameters α , β , and λ for the complete failure sample ( T , X ) can be derived by taking the derivatives for the log-likelihood function with respect to α , β , and λ and setting them to zero, which can be expressed as follows:
α = 1 α m 1 + m 2 + j = 1 2 k = 1 m j R j k + j = 1 2 i = 1 m j λ j 1 ( 1 e t j i β ) + j = 1 2 k = 1 m j s = 1 R j k λ j 1 ( 1 e ( x j k s ) β = 0 , β = 1 β m 1 + m 2 + j = 1 2 k = 1 m j R j k + j = 1 2 i = 1 m j ln t j i + j = 1 2 k = 1 m j s = 1 R j k ln x j k s   + j = 1 2 i = 1 m j t j i β ( 1 α λ j 1 e t j i β ) ln t j i + j = 1 2 k = 1 m j s = 1 R j k ( x j k s ) β ( 1 α λ j 1 e ( x j k s ) β ) ln x j k s = 0 , λ = 1 λ m 2 + k = 1 m 2 R 2 k + α i = 1 m 2 ( 1 e t 2 i β ) + k = 1 m 2 s = 1 R 2 k ( 1 e ( x 2 k s ) β ) = 0 .
Furthermore, one has
α α ( β ) = B ( e t 2 i β , e ( x 2 k s ) β ) E ( e t 2 i β , e ( x 2 k s ) β ) A D D C ( e t 1 i β , e ( x 1 k s ) β ) , λ λ ( β ) = C ( e t 1 i β , e ( x 1 k s ) β ) D 2 E ( e t 2 i β , e ( x 2 k s ) β ) [ B ( e t 2 i β , e ( x 2 k s ) β ) E ( e t 2 i β , e ( x 2 k s ) β ) A D ] ,
where
A = m 1 + m 2 + j = 1 2 k = 1 m j R j k , B ( b 2 i , b 2 k s ) = i = 1 m 2 ( 1 b 2 i ) + k = 1 m 2 s = 1 R 2 k ( 1 b 2 k s ) , C ( c 1 i , c 1 k s ) = i = 1 m 1 ( 1 c 1 i ) + k = 1 m 1 s = 1 R 1 k ( 1 c 1 k s ) , D = m 2 + k = 1 m 2 R 2 k , E ( e 2 i , e 2 k s ) = i = 1 m 2 ( 1 e 2 i ) + k = 1 m 2 s = 1 R 2 k ( 1 e 2 k s ) .
and α α ( β ) represents α equivalent to α ( β )
By substituting (10) to β = 0 in (9), one has
β = m 1 + m 2 + j = 1 2 k = 1 m j R j k j = 1 2 i = 1 m j ln t j i + j = 1 2 k = 1 m j s = 1 R j k ln x j k s + j = 1 2 i = 1 m j t j i β ( 1 α ( β ) λ j 1 ( β ) e t j i β ) ln t j i + j = 1 2 k = 1 m j s = 1 R j k ( x j k s ) β ( 1 α ( β ) λ j 1 ( β ) e ( x j k s ) β ) ln x j k s .
Therefore, it can be observed that the pseudo likelihood equations can be rewritten as a nonlinear function of β .
Given T = ( T j 1 , T j 2 , , T j m j ) , j = 1 , 2 , for  k = 1 , 2 , , m j , the conditional distribution of X j k s , s = 1 , 2 , , R j k follows a truncated Chen distribution with scale parameter α λ j 1 and shape parameter β , of which the density can be expressed as
f X | T ( x j k | t ) = f j ( x j k ) 1 F j ( t j k ) = α β λ j 1 x j k β 1 e α λ j 1 ( 1 e x j k β ) + x j k β e α λ j 1 ( 1 e t j k β ) + t j k β , x j k > t j k .
For j = 1 , 2 and k = 1 , 2 , , m j , one has
E j 1 ( α , β , λ ) = E j ( e ( X j k s ) β | X j k > t j k ) = 1 + α λ j 1 e t j k β α λ j 1 , E j 2 ( α , β , λ ) = E j ( ln X j k s | X j k > t j k ) = ln t j k + e α λ j k e t j k β β t j k β t 1 e α λ j 1 e t d t , E j 3 ( α , β , λ ) = E j ( ( X j k s ) β ln X j k s | X j k > t j k ) = α λ j 1 β S j ( t j k ) t j k β t e α λ j 1 ( 1 e t ) + t ln t d t , E j 4 ( α , β , λ ) = E j ( ( X j k s ) β e ( X j k s ) β ln X j k s | X j k > t j k ) = α β 1 λ j 1 e α λ j 1 e t i j β t j k β t e 2 t α λ j 1 e t ln t d t .
Hence, an EM algorithm (Algorithm 1) to calculate the MLEs of the parameters α , β , and λ is proposed as follows.
Algorithm 1 EM iterative algorithm.
Step 1
Let α ( 0 ) , β ( 0 ) , and λ ( 0 ) be the initial guess values of the MLEs for α , β , and λ .
Step 2
In the ( d + 1 ) th iterative:
  • (E-step) Replace any function of X j k s (say h X j k s ) by E ( h ( X j k s ) | X j k s > T j k ) , and the likelihood equations are replaced by
    β = m 1 + m 2 + j = 1 2 k = 1 m j R j k j = 1 2 i = 1 m j ln t j i + j = 1 2 k = 1 m j s = 1 R j k E j 2 ( α ( d ) , β ( d ) , λ ( d ) ) + j = 1 2 i = 1 m j t j i β ( 1 [ α ( β ) ] ( d ) ( [ λ ( β ) ] ( d ) ) j 1 e t j i β ) ln t j i + j = 1 2 k = 1 m j s = 1 R j k E j 3 ( α ( d ) , β ( d ) , γ ( d ) ) j = 1 2 k = 1 m j s = 1 R j k E j 4 ( α ( d ) , β ( d ) , γ ( d ) ) [ α ( β ) ] ( d ) ( [ λ ( β ) ] ( d ) ) j 1 ,
    where
    [ α ( β ) ] ( d ) = B ( e t 2 i β , E 21 ( α ( d ) , β ( d ) , λ ( d ) ) ) E ( e t 2 i β , E 21 ( α ( d ) , β ( d ) , λ ( d ) ) ) A D D C ( e t 1 i β , E 11 ( α ( d ) , β ( d ) , λ ( d ) ) ) , [ λ ( β ) ] ( d ) = C ( e t 1 i β , e ( x 1 k s ) β ( d ) ) D 2 E ( e t 2 i β , E 21 ( α ( d ) , β ( d ) , λ ( d ) ) ) · [ B ( e t 2 i β , E 21 ( α ( d ) , β ( d ) , λ ( d ) ) ) E ( e t 2 i β , E 21 ( α ( d ) , β ( d ) , λ ( d ) ) ) A D ] .
  • (M-step) The estimate of β , namely β ( k + 1 ) , can be derived iteratively by solving the following nonlinear equations:
    β ( d + 1 ) = m 1 + m 2 + j = 1 2 k = 1 m j R j k j = 1 2 i = 1 m j ln t j i + j = 1 2 k = 1 m j R j k E j 2 ( α ( d ) , β ( d ) , λ ( d ) ) + j = 1 2 i = 1 m j t j i β ( d + 1 ) ( 1 [ α ( β ) ] ( d ) ( [ λ ( β ) ] ( d ) ) j 1 e t j i β ( d + 1 ) ) ln t j i + j = 1 2 k = 1 m j R j k E j 3 ( α ( d ) , β ( d ) , γ ( d ) ) j = 1 2 k = 1 m j R j k E j 4 ( α ( d ) , β ( d ) , γ ( d ) ) [ α ( β ) ] ( d ) ( [ λ ( β ) ] ( d ) ) j 1 ,
    and the estimates of α and λ , namely α ( d + 1 ) , λ ( d + 1 ) , can be further derived as
    α ( d + 1 ) = B ( e t 2 i β ( d + 1 ) , E 21 ( α ( d ) , β ( d ) , γ ( d ) ) ) E ( e t 2 i β ( d + 1 ) , E 21 ( α ( d ) , β ( d ) , γ ( d ) ) ) A D D C ( e t 1 i β ( d + 1 ) , E 11 ( α ( d ) , β ( d ) , γ ( d ) ) ) , λ ( d + 1 ) = C ( e t 1 i β ( d + 1 ) , E 11 ( α ( d ) , β ( d ) , γ ( d ) ) ) D 2 E ( e t 2 i β ( d + 1 ) , E 21 ( α ( d ) , β ( d ) , γ ( d ) ) ) · [ B ( e t 2 i β ( d + 1 ) , E 21 ( α ( d ) , β ( d ) , γ ( d ) ) ) E ( e t 2 i β ( d + 1 ) , E 21 ( α ( d ) , β ( d ) , γ ( d ) ) ) A D ] .
Step 4
Stop the iteration, and find the EM-based estimates of α , β , and λ when | α ( d + 1 ) α ( d ) | < ε , | β ( d + 1 ) β ( d ) | < ε and | λ ( d + 1 ) λ ( d ) | < ε for some given tolerance limit ε , for example ε = 0.001 .
Remark 1.
From Wu [14], it is observed that, for samples T j 1 < T j 2 < < T j m j , j = 1 , 2 , one has that
T j ( β ) = 1 n j ( m j 1 ) i = 1 m j ( R j i + 1 ) e t j i β 1 e t j 1 β 1 1 m j 1
follows the F distribution with 2 m j 2 and 2 degrees of freedom and T j ( β ) is a strictly increasing function of β. Hence, for 0 < γ < 1 , one can choose the initial estimate of β from the following 100 ( 1 γ ) % confidence interval:
ρ ( t , F ( 2 ( n 1 ) , 2 ) ( 1 γ 2 ) ) , ρ ( t , F ( 2 ( n 1 ) , 2 ) ( γ 2 ) ) ,
where F ( a , b ) ( p ) is the 100 p % right-tail percentile of the F distribution with a and b degrees of freedom and ρ ( t , y ) is the solution of β for equation T j ( β ) = y . Meanwhile, one can also let T j ( β ) , j = 1 , 2 equal the median of the F ( 2 ( n 1 ) , 2 ) distribution and find the root of β, which can also be utilized as the initial estimate for β.

4. Confidence Interval Estimation

In this section, common large-sample-based asymptotic confidence intervals of the model parameters were constructed via the Fisher information. In addition, the bootstrap sampling technique was also used to obtain the bootstrap confidence intervals for comparison.

4.1. Asymptotic Confidence Intervals

Under some mild regularity (see, e.g., Casella and Berger [23]), the asymptotic confidence intervals of α , β , and λ can be derived from the usually asymptotic normality of the maximum likelihood estimation with empirical variances estimated from the inverse of the observed Fisher information matrix.
By direct calculation, the second partial derivatives of the log-likelihood function in (5) can be expressed as
α α = 2 α 2 = m 1 + m 2 α 2 , α β = 2 α β = j = 1 2 i = 1 m j λ j 1 ( 1 + R j i ) e t j i β t j i β log t j i , α λ = 2 α λ = i = 1 m 2 ( 1 + R 2 i ) ( 1 e t 2 i β ) , β β = 2 β 2 = m 1 + m 2 β 2 + j = 1 2 i = 1 m j t j i β log 2 t j i α j = 1 2 i = 1 m j λ j 1 ( 1 + R j i ) e t j i β t j i β ( 1 + t j i β ) log 2 t j i , β λ = 2 β λ = α i = 1 m 2 ( 1 + R 2 i ) e t 2 i β t 2 i β log t 2 i , λ λ = 2 λ 2 = m 2 λ 2 .
Thus, the observed Fisher information matrix I = I ( α , β , λ ) can be expressed as
I ( α ^ , β ^ , λ ^ ) = α α α β α λ α β β β β λ α λ β λ λ λ | α ^ , β ^ , λ ^ = Var ( α ^ ) Cov ( α ^ , β ^ ) Cov ( α ^ , λ ^ ) Cov ( α ^ , β ^ ) Var ( β ^ ) Cov ( β ^ , λ ^ ) Cov ( α ^ , λ ^ ) Cov ( β ^ , λ ^ ) Var ( λ ^ ) 1 .
Based on the asymptotic theory, the asymptotic distribution of ( α ^ , β ^ , λ ^ ) is normally distributed with mean vector ( α , β , λ ) and variance–covariance matrix I 1 ( α ^ , β ^ , λ ^ ) , i.e.,
α ^ β ^ λ ^ α β λ d N 0 , I 1 ( α ^ , β ^ , λ ^ ) ,
where d denotes the convergence in the distribution and I 1 ( α ^ , β ^ , λ ^ ) is the inverse of the Fisher information matrix I ( α ^ , β ^ , λ ^ ) . For 0 < γ < 1 , the asymptotic 100 ( 1 γ ) % confidence intervals of α , β , and λ are given by
θ ^ i z γ / 2 Var ( θ ^ i ) , θ ^ i + z γ / 2 Var ( θ ^ i ) , i = 1 , 2 , 3 ,
where θ 1 = α , θ 2 = β and θ 3 = λ , and z γ is the upper γ th quantile of the standard normal distribution.
Sometimes, previous asymptotic confidence intervals may have negative lower bounds. In order to overcome this drawback, the logarithmic transformation and delta methods can be employed to obtain the asymptotic normality distribution of ln θ ^ i as ( ln θ ^ i ln θ i ) d N 0 , Var ( θ ^ i ) / θ i 2 , i = 1 , 2 , 3 , which implies that the asymptotic 100 ( 1 γ ) % confidence interval of ln θ i is
( ln θ ^ i ± z λ / 2 Var ( θ ^ i ) / θ ^ i ) ( A i 1 , A i 2 ) , i = 1 , 2 , 3 .
Hence, the asymptotic 100 ( 1 γ ) % confidence interval of θ i , i = 1 , 2 , 3 can be constructed as ( e A i 1 , e A i 2 ) , respectively.

4.2. Bootstrap Confidence Intervals

The bootstrap method [24] is known to enlarge the sample by simulation so as to overcome the shortage of experimental data. Compared with the exact method, the bootstrap estimation is widely used due to its great properties such as its easy calculation. Generally, one can derive the bootstrap confidence interval by both nonparametric and parametric approaches. Following the suggestion of Kundu et al. [25], the parametric bootstrap confidence intervals of the unknown model parameters were constructed in this part, where the Studentized-t bootstrap confidence interval (SBCI) from Hall [26] and the percentile bootstrap confidence interval (PBCI) suggested by Efron [24] were obtained, respectively.
Under the constant-stress PALT, the following procedure of the progressively Type-II-censored samples is provided to generate the associated bootstrap data (e.g., Balakrishnan and Sandhu [27]):
  • Based on the original progressively Type-II-censored sample t j 1 , t j 2 , , t j m j , j = 1 , 2 and censoring scheme R j = ( R j 1 , R j 2 , , R j m j ) , obtain the estimates of α , β , and λ , say α ^ , β ^ , λ ^ .
  • For the given m j and n j , j = 1 , 2 , generate m j independent observations ( W j 1 , W j 2 , , W j m j ) from uniform distribution U ( 0 , 1 ) .
  • Set E j i = 1 / ( i + k = m j i + 1 m j R j k ) , for j = 1 , 2 , i = 1 , 2 , , m j .
  • Set V j i = W j i E j i , for j = 1 , 2 , i = 1 , 2 , , m j .
  • Set U j i = 1 k = m j i + 1 m j V j k , i = 1 , 2 , , m j . Then, U j 1 , U j 2 , , U j m j , j = 1 , 2 is the progressively Type-II-censored data from U ( 0 , 1 ) .
  • Based on α ^ , β ^ , and λ ^ , two random samples ( t j 1 , t j 2 , , t j m j ) , j = 1 , 2 from CDFs F 1 ( t ) and F 2 ( t ) presented in (1) and (3) are generated as follows:
    t j i = ln 1 ln ( 1 U j i ) α ^ λ ^ j 1 1 / β ^ .
  • As in Step 1, based on ( t j 1 , t j 2 , , t j m j ) , j = 1 , 2 , the bootstrap sample estimates of α , β , λ are computed, namely α ^ , β ^ , and λ ^ .
  • Repeat Steps 2–7 N times; we can generate N different bootstrap samples for α , β , and λ , i.e., θ ^ k 1 , θ ^ k 2 , , θ ^ k N , k = 1 , 2 , 3 , where θ 1 = α , θ 2 = β , and θ 3 = λ , respectively.
  • Arrange all α ^ , β ^ , and λ ^ in an ascending order; one can obtain the bootstrap samples for α , β , and λ as θ ^ k [ 1 ] , θ ^ k [ 2 ] , , θ ^ k [ N ] , k = 1 , 2 , 3 .
  • Percentile bootstrap confidence intervals:
Let G ( t ) = P ( θ ^ k t ) be the CDF of θ ^ k and θ ^ k b o o t = G 1 ( t ) for given t. For 0 < γ < 1 , the approximate 100 ( 1 γ ) % PBCI of θ k is given by
θ ^ k b o o t γ 2 , θ ^ k b o o t 1 γ 2 , k = 1 , 2 , 3 .
  • Studentized-t bootstrap confidence intervals:
In order to construct the Studentized-t bootstrap confidence intervals for parameters α , β , and λ , we derive the order statistics:
η k [ 1 ] < η k [ 2 ] < < η k [ N ] , with η k [ d ] = θ ^ k [ d ] θ ^ k Var ( θ ^ k [ d ] ) , k = 1 , 2 , 3 , d = 1 , 2 , , N .
Denote G ( t ) = P ( η k t ) as the CDF of η k and η k b o o t t = θ ^ k + Var ( θ ^ k ) G 1 ( t ) for given t. Thus, for 0 < γ < 1 , the 100 ( 1 γ ) % SBCI of θ k can be expressed as
η ^ k b o o t t γ 2 , η ^ k b o o t t 1 γ 2 , k = 1 , 2 , 3 ,
where Var ( θ ^ k ) can be estimated by the asymptotic variance from the observed Fisher information matrix I ( α ^ , β ^ , λ ^ ) .

5. Numerical Studies

In this section, a real data example and simulation studies are presented for illustrative purposes.

5.1. Data Analysis

The dataset from a light-emitting diode (LED) life test was analyzed, which was also analyzed by Cheng and Wang [11]. Since the original data included 155 samples with censoring level 0.25 , for the sake of simplicity, we just extracted the observed failure samples obtained under the use or stress conditions as the set of complete constant-stress PALT data, for which the details are listed in Table 1.
For the purposes of comparison, we fit the Chen distribution, as well as the normal, the Weibull, gamma, and log-normal to fit this real dataset. Since the exponential distribution is a special case of the Weibull and gamma distributions, the exponential distribution was also considered as a candidate model as well. In order to select the best-fitting model from the candidate distributions, the goodness-of-fit criteria quantities are presented in Table 2, including the Akaike information criterion (AIC), Bayesian information criterion (BIC), the second-order Akaike information criterion (AICc), as well as the values of the log-likelihood. Let k be the number of free parameters in the related model and n be the sample size, the related criteria are defined as
AIC = 2 k 2 ln ( likelihood ) , BIC = k ln n 2 ln ( likelihood ) , AICc = AIC + 2 k ( k + 1 ) n k 1 .
The best model is the one with the least values of the AIC, BIC, and AICc. The results presented in Table 2 indicate that the Chen distribution provided the best fit and can be used to analyze these data.
Based on the complete data given in Table 1, a group of progressively Type-II-censored constant-stress PALT data was generated, and the details are presented in Table 3. Therefore, different point and interval estimates were obtained, and the results are provided in Table 4 with the 90 % significance level for the interval estimates. It can be seen from Table 4, the bootstrap confidence intervals outperformed the ACIs, where the SBCIs of α , β and λ also featured a narrower interval length than those of PBCIs.

5.2. Simulation Studies

In this part, extensive numerical studies were conducted to investigate the performance of our methods. To generate the progressively Type-II-censored samples under the constant-stress PALT, the algorithm developed by Balakrishnan and Aggarwala [18] was implemented. In addition, in order to compare the performance of different estimates, four quantiles were examined as follows:
(a)
The mean-squared errors (MSEs) of the estimates of the parameters α , β , and λ computed as 1 M d = 1 M ( ϕ l ϕ ^ l d ) 2 , repeated M-times, where ϕ 1 = α , ϕ 2 = β , and ϕ 3 = λ denote the selected values of the parameters α , β , and λ and ϕ ^ l d , d = 1 , 2 , , M is the related estimates of ϕ l , l = 1 , 2 , 3 .
(b)
The average estimates (AVEs) of the MLEs for the parameters α , β , and λ , which are defined as ϕ ^ l = 1 M ϕ ^ l d .
(b)
The coverage probabilities (CPs) of the confidence intervals of the parameters α , β , and λ , which are defined as the probability that the interval estimate contains the true parameters.
(c)
The average widths (AWs) of the confidence intervals of the parameters α , β , and λ .
For j = 1 , 2 , three progressive censoring schemes were considered as follows:
Scheme I:
R j 1 = R j 2 = = R j ( m j 1 ) = 0 and R j m j = n j m j ;
Scheme II:
R j 1 = n j m j and R j 2 = = R j m j = 0 ;
Scheme III:
R j 1 = R j 2 = R j ( m j 1 ) = 1 and R j m j = n j 2 m j + 1 .
Under different choices of the parameter values, sample sizes, and censoring schemes, the results of criteria quantities are presented in Table 5, Table 6, Table 7 and Table 8 based on 10,000 repetitions, where the numerical computational work was conducted by using the MATLAB 2017a software; the interval results were obtained under a 90 % significance level, and these same censoring schemes were used under both the normal use and accelerated stress stages. It is noted that the criteria quantities, the MSEs and AVEs, of the parameters α , β , and λ decreased when the effective sample size m , n increased, which indicated that the MLEs were consistent and worked satisfactorily under the simulated scenarios. In addition, the point estimates obtained from Censoring Scheme II outperformed the results from Censoring Schemes I and III in terms of the MSEs and AVEs in general. Furthermore, it was seen that the different AWs of the various interval estimates decreased when the sample size increased, whereas the CPs increased correspondingly with the same trend. Moreover, it was also observed that the two bootstrap confidence intervals were superior to the likelihood-based ACIs according to the quantities’ AWs, whereas the SBCI provided relatively shorter interval widths than those of the PBCI in most cases, which is in agreement with the results of Hall [26]. Similar to the point results and compared to the other two censoring schemes, the interval estimates obtained from Censoring Scheme II also had relatively shorter interval lengths under fixed simulated scenarios.

6. Concluding Remarks

In this paper, statistical inference was discussed for a constant-stress partially accelerated life test. Under progressive Type-II censoring, parameter estimation was conducted, when the lifetime of the units followed the two-parameter bathtub-shaped Chen distribution. The expectation–maximization procedure was implemented for maximum likelihood estimation, and the confidence intervals of the model parameters were derived by using the asymptotic theory and bootstrap techniques. The performance of the point and interval estimates was investigated by simulation studies and a real-life example, and the results indicated that the proposed methods worked satisfactorily. Although the constant-stress partially accelerated life test was considered for the bathtub-shaped Chen distribution in this paper, the results could be extended to other lifetime models with a proper modification. For further works, the testing design and optimal sampling planning for the constant-stress partially accelerated life test seem also interesting and important in practice, which will be discussed in the future.

Author Contributions

Funding acquisition, L.W.; formal analysis, Y.N., L.W., and Y.M.T.; methodology, L.W.; writing—original draft, L.W.; writing—review and editing, Y.M.T.; investigation, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Yunnan Fundamental Research Project (No. 202101AT070103).

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Editor and Reviewers for their valuable comments and suggestions, which improved this paper significantly.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abdel-Hamid, A.H.; Al-Hussaini, E.K. Progressive stress accelerated life tests under finite mixture models. Metrika 2007, 66, 213–231. [Google Scholar] [CrossRef]
  2. Balakrishnan, N.; Kundu, D.; Ng, H.K.T.; Kannan, N. Point and interval estimation for a simple step-stress model with Type-II censoring. J. Qual. Technol. 2007, 39, 35–47. [Google Scholar] [CrossRef]
  3. Kim, C.M.; Bai, D.S. Analysis of accelerated life test data under two failure modes. International Journal of Reliability. Qual. Saf. Eng. 2002, 9, 111–125. [Google Scholar] [CrossRef]
  4. Nelson, W. Residuals and their analysis for accelerated life tests with step and varying stress. IEEE Trans. Reliab. 2008, 57, 360–368. [Google Scholar] [CrossRef]
  5. Watkins, A.J.; John, A.M. On constant stress accelerated life tests terminated by Type-II censoring at one of the stress levels. J. Stat. Plan. Inference 2008, 138, 768–786. [Google Scholar] [CrossRef]
  6. Nelson, W. Accelerated Testing: Statistical Models, Test Plans, and Data Analysis; John Wiley & Sons: New York, NY, USA, 1990. [Google Scholar]
  7. Abdel-Hamid, A.H. Constant-partially accelerated life tests for Burr Type-XII distribution with progressive Type-II censoring. Comput. Stat. Data Anal. 2009, 53, 2511–2523. [Google Scholar] [CrossRef]
  8. Abdel-Hamid, A.H.; Al-Hussaini, E.K. Step partially accelerated life tests under finite mixture models. J. Stat. Comput. Simul. 2008, 78, 911–924. [Google Scholar] [CrossRef]
  9. Abdel-Ghaly, A.A.; Attia, A.F.; bdel-Ghani, M.M. The maximum likelihood estimates in step partially accelerated life tests for the Weibull parameters in censored data. Commun. Stat. Theory Methods 2002, 31, 551–573. [Google Scholar] [CrossRef]
  10. Ismail, A.A. Inference for a step-stress partially accelerated life test model with an adaptive Type-II progressively hybrid censored data from Weibull distribution. J. Comput. Appl. Math. 2014, 260, 533–542. [Google Scholar] [CrossRef]
  11. Cheng, Y.F.; Wang, F.K. Estimating the Burr XII parameters in constant-stress partially accelerated life tests under multiple censored data. Commun. Stat.-Simul. Comput. 2012, 41, 1711–1727. [Google Scholar] [CrossRef]
  12. Chen, Z. A new two parameter lifetime distribution with bathtub shaped or increasing failure rate function. Stat. Probab. Lett. 2000, 49, 155–161. [Google Scholar] [CrossRef]
  13. 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]
  14. Wu, S.J. Estimation of the two parameter bathtub-shaped lifetime distribution with progressive censoring. J. Appl. Stat. 2008, 35, 1139–1150. [Google Scholar] [CrossRef]
  15. 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]
  16. Elshahhat, A.; Rastogi, M.K. Bayesian Life Analysis of Generalized Chen’s Population Under Progressive Censoring. Pak. Stat. Oper. Res. 2022, 18, 675–702. [Google Scholar] [CrossRef]
  17. Balakrishnan, N. Progressive censoring methodology: An appraisal. TEST 2007, 16, 211–259. [Google Scholar] [CrossRef]
  18. Balakrishnan, N.; Aggarwala, R. Progressive Censoring: Theory, Methods and Applications; Birkhause: Boston, MA, USA, 2000. [Google Scholar]
  19. 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. [Google Scholar] [CrossRef]
  20. Soliman, A.A.; Abd-Elmougod, G.A.; Al-Sobhi, M.M. Estimation in step-stress partially accelerated life tests for the Chen distribution using progressive Type-II censoring. Appl. Math. Inf. Sci. 2017, 11, 325–332. [Google Scholar] [CrossRef]
  21. Dempster, A.P.; Laird, N.M.; Rubin, D.B. Maximum likelihood estimation from incomplete data via the EM algorithm. J. R. Stat. Soc. Ser. B 1977, 39, 1–22. [Google Scholar]
  22. McLachlan, G.J.; Krishnan, T. The EM Algorithm and Extensions, 2nd ed.; John Wiley & Sons: New York, NY, USA, 2008. [Google Scholar]
  23. Casella, G.; Berger, R.L. Statistical Inference, 2nd ed.; Duxbury Press: Pacific Grove, CA, USA, 2002. [Google Scholar]
  24. Efron, B. Computers and the theory of statistics: Thinking the unthinkable. SIAM Rev. 1979, 21, 460–480. [Google Scholar] [CrossRef]
  25. Kundu, D.; Kannan, N.; Balakrishnan, N. Analysis of progressively censored competing risks data. In Handbook of Statistics; Balakrishnan, N., Rao, C.R., Eds.; Elsevier: New York, NY, USA, 2004; Volume 23. [Google Scholar]
  26. Hall, P. Theoretical comparison of bootstrap confidence intervals. Ann. Stat. 1988, 16, 927–953. [Google Scholar] [CrossRef]
  27. Balakrishnan, N.; Sandhu, R.A. A simple simulational algorithm for generating progressive Type-II censored samples. Am. Stat. 1995, 49, 229–230. [Google Scholar]
Table 1. Complete LED constant-stress PALT failure data.
Table 1. Complete LED constant-stress PALT failure data.
Regular Use Condition
0.18, 0.19, 0.19, 0.34, 0.36, 0.40 0.44, 0.44, 0.45, 0.46, 0.47, 0.53, 0.57, 0.57, 0.63, 0.65, 0.70, 0.71,
0.71, 0.75, 0.76, 0.76, 0.79, 0.80, 0.85, 0.98, 1.01, 1.07, 1.12, 1.14, 1.15, 1.17, 1.20, 1.23, 1.24, 1.25,
1.26, 1.32, 1.33, 1.33, 1.39, 1.42, 1.50, 1.55, 1.58, 1.59, 1.62, 1.68, 1.70, 1.79, 2.00, 2.01, 2.04, 2.54,
3.61, 3.76, 4.65, 8.97.
Stress condition
0.13, 0.16, 0.20, 0.20, 0.21, 0.25, 0.26, 0.28, 0.28, 0.30, 0.31, 0.33, 0.35, 0.35, 0.35, 0.39, 0.50, 0.52,
0.58, 0.60, 0.60, 0.62, 0.63, 0.67, 0.71, 0.73, 0.75, 0.75, 0.78, 0.80, 0.80, 0.86, 0.90, 0.91, 0.93, 0.93,
0.94, 0.98, 0.99, 1.01, 1.03, 1.06, 1.06, 1.10, 1.22, 1.22, 1.24, 1.28, 1.39, 1.39, 1.46, 1.48, 1.52, 1.74,
1.95, 2.46, 3.02, 5.16.
Table 2. Summary of the goodness-of-fit test for the fit distributions.
Table 2. Summary of the goodness-of-fit test for the fit distributions.
Results under Data from Regular Use Conditions
Modellog-likelihoodAICBICAICc
Chen−52.9318109.8636113.9845110.0818
Weibull−71.6092147.2185151.3394147.4367
Gamma−69.2207142.4414146.5623142.6596
normal−126.8569257.7138261.8347257.9320
log-normal−93.9304191.8608195.9817192.0790
exponential−74.3596150.7193152.7797150.7907
results under data from stress conditions
Modellog-likelihoodAICBICAICc
Chen−41.127186.254290.375186.4724
Weibull−49.0693102.1386106.2595102.3568
Gamma−47.470298.9404103.061399.1586
normal−97.2140198.4279202.5488198.6461
log-normal−93.9592191.9184196.0393192.1366
exponential−53.4458108.8916110.9520108.9630
Table 3. Progressively Type-II-censored LED constant-stress PALT data.
Table 3. Progressively Type-II-censored LED constant-stress PALT data.
Regular Use Conditions
T 1 = ( 0.18 , 0.40 , 0.47 , 0.65 , 0.79 , 1.07 , 1.20 , 1.32 , 1.50 , 1.68 , 2.04 , 4.65 )
R 1 = ( 4 , 4 , 4 , 6 , 4 , 4 , 4 , 4 , 4 , 4 , 3 , 1 )
Stress conditions
T 2 = ( 0.13 , 0.21 , 0.30 , 0.35 , 0.58 , 0.63 , 0.75 , 0.86 , 0.94 , 1.03 , 1.22 , 1.39 , 1.95 , 5.16
R 2 = ( 3 , 4 , 3 , 4 , 3 , 4 , 3 , 4 , 3 , 4 , 3 , 4 , 2 , 0 )
Table 4. Estimates under LED constant-stress PALT data.
Table 4. Estimates under LED constant-stress PALT data.
MLEACIPBCISBCI
α 0.0833(0.0389,0.2177)(0.0496,0.2224)(0.0244,0.1931)
0.17880.17280.1687
β 0.7286(0.6137,0.8435)(0.6234,0.8460)(0.5461,0.7459)
0.22980.22260.1998
λ 1.6067(0.5604,2.6530)(1.1400,2.5212)(1.2593,2.6011)
2.09261.38121.3418
Table 5. AVEs and MSEs (with brackets) for the parameters with ( α , β , λ ) = ( 1 , 0.5 , 2 ) at n 1 = n 2 = n and m 1 = m 2 = m .
Table 5. AVEs and MSEs (with brackets) for the parameters with ( α , β , λ ) = ( 1 , 0.5 , 2 ) at n 1 = n 2 = n and m 1 = m 2 = m .
nmCS α β λ
3010I1.2217 (0.3254)0.6041 (0.2322)1.8104 (0.2416)
II1.2034 (0.3135)0.5783 (0.2158)1.8329 (0.2337)
III1.2516 (0.3217)0.5926 (0.2411)2.2021 (0.2504)
20I1.1510 (0.3126)0.4124 (0.2258)2.1765 (0.2300)
II1.1379 (0.2989)0.5473 (0.2122)2.1341 (0.2213)
III1.2037 (0.3097)0.5721 (0.2194)1.8647 (0.2321)
5015I0.8153 (0.2205)0.5678 (0.1601)1.8790 (0.1818)
II1.1536 (0.2114)0.5545 (0.1534)1.9444 (0.1726)
III1.1764 (0.2198)0.5581 (0.1643)1.8426 (0.1769)
25I0.8437 (0.2143)0.5454 (0.1532)2.1095 (0.1745)
II0.8924 (0.2009)0.4792 (0.1521)2.0423 (0.1708)
III0.9609 (0.2132)0.5562 (0.1530)2.1111 (0.1723)
10035I1.1111 (0.1421)0.5333 (0.1532)2.0899 (0.1251)
II1.0613 (0.1346)0.5201 (0.1467)2.0351 (0.1139)
III1.1453 (0.1399)0.5400 (0.1499)2.0932 (0.1300)
50I1.0712 (0.1314)0.5204 (0.1406)2.0532 (0.1126)
II1.0221 (0.1289)0.5172 (0.1375)2.0241 (0.1089)
III0.9105 (0.1357)0.5198 (0.1437)1.9613 (0.1231)
Table 6. AWs and CPs (with brackets) for the parameters with ( α , β , λ = ( 1 , 0.5 , 2 ) at n 1 = n 2 = n and m 1 = m 2 = m ( γ = 0.1 ).
Table 6. AWs and CPs (with brackets) for the parameters with ( α , β , λ = ( 1 , 0.5 , 2 ) at n 1 = n 2 = n and m 1 = m 2 = m ( γ = 0.1 ).
nmCS α β λ
ACISBCIPBCIACISBCIPBCIACISBCIPBCI
3010I0.82130.53060.61892.32161.53241.60411.42790.90780.9111
(0.8722)(0.8631)(0.8624)(0.8701)(0.8623)(0.8712)(0.8671)(0.8641)(0.8650)
II0.80110.49460.57172.20751.45611.51021.24360.85680.8389
(0.8693)(0.8556)(0.8579)(0.8643)(0.8656)(0.8697)(0.8632)(0.8637)(0.8725)
III0.93110.50000.60342.24531.56561.53221.29790.87690.8662
(0.8731)(0.8615)(0.8611)(0.8679)(0.8611)(0.8714)(0.8574)(0.8679)(0.8734)
20I0.78910.50080.58691.97851.25621.30081.31530.83590.8533
(0.8760)(0.8698)(0.8658)(0.8697)(0.8630)(0.8620)(0.8602)(0.8635)(0.8639)
II0.72220.46250.56341.80111.10541.21211.11370.78860.8216
(0.8663)(0.8571)(0.8567)(0.8635)(0.8627)(0.8654)(0.8551)(0.8612)(0.8654)
III0.80100.47730.57431.88841.34091.25641.25320.80420.8480
(0.8723)(0.8673)(0.8712)(0.8722)(0.8543)(0.8635)(0.8563)(0.8627)(0.8640)
5015I0.80590.46720.49181.55341.20741.23000.91570.79930.7899
(0.8924)(0.8731)(0.8710)(0.8822)(0.8717)(0.8652)(0.8791)(0.8663)(0.8711)
II0.74740.45510.46351.53071.19351.21150.89930.72690.7425
(0.8812)(0.8678)(0.8689)(0.8769)(0.8642)(0.8657)(0.8790)(0.8659)(0.8661)
III0.77130.47500.50011.60891.26171.30210.93340.76190.7713
(0.8929)(0.8754)(0.8723)(0.8843)(0.8754)(0.8690)(0.8814)(0.8660)(0.8700)
25I0.74510.43000.44191.22300.96110.98870.86560.68900.6998
(0.8809)(0.8754)(0.8665)(0.8669)(0.8713)(0.8743)(0.8691)(0.8658)(0.8671)
II0.70090.42140.43271.21320.87940.92360.81020.65720.6627
(0.8807)(0.8629)(0.8624)(0.8614)(0.8610)(0.8689)(0.8678)(0.8654)(0.8665)
III0.75130.42670.46781.23450.89990.95420.84230.73140.7320
(0.8816)(0.8742)(0.8741)(0.8736)(0.9645)(0.8720)(0.8684)(0.8637)(0.8702)
10035I0.53270.28110.28410.86780.62310.64230.62170.50130.5134
(0.9132)(0.8956)(0.9012)(0.9231)(0.8942)(0.8972)(0.8963)(0.8956)(0.9142)
II0.52020.26790.27000.75910.61240.62220.59360.48970.4999
(0.8965)(0.8817)(0.8835)(0.8922)(0.8879)(0.8938)(0.8856)(0.8901)(0.9009)
III0.52200.27440.27860.77770.65150.62310.61100.51100.5183
(0.9111)(0.8924)(0.8971)(0.9045)(0.9023)(0.8939)(0.8942)(0.9014)(0.9125)
50I0.47120.18490.20130.61350.52140.53270.51210.42170.3868
(0.9121)(0.9121)(0.9178)(0.9127)(0.9011)(0.9126)(0.8971)(0.9131)(0.9063)
II0.46570.17350.19890.56700.51730.51990.47850.35610.3620
(0.8934)(0.8903)(0.8997)(0.8930)(0.8954)(0.8979)(0.8849)(0.8943)(0.8974)
III0.46990.23520.22470.59780.53360.52310.49900.37540.3926
(0.9023)(0.9156)(0.9234)(0.8965)(0.9087)(0.9118)(0.8937)(0.9072)(0.9125)
Table 7. AVE and MSEs (with brackets) for the parameters with ( α , β , λ ) = ( 0.7 , 1 , 1.5 ) at n 1 n 2 and m 1 m 2 .
Table 7. AVE and MSEs (with brackets) for the parameters with ( α , β , λ ) = ( 0.7 , 1 , 1.5 ) at n 1 n 2 and m 1 m 2 .
n 1 m 1 CS α β λ
n 2 m 2
3010I0.9755 (0.3257)1.3221 (0.2214)1.7987 (0.2032)
4015II0.9537 (0.3216)1.3024 (0.2185)1.7611 (0.2011)
III0.9701 (0.3309)1.3176 (0.2210)1.8432 (0.2058)
16I0.9314 (0.3203)1.3010 (0.2179)1.7534 (0.1978)
18II0.8973 (0.3168)1.2877 (0.2124)1.7518 (0.1963)
III0.9102 (0.3200)1.3112 (0.2213)1.7642 (0.1991)
5015I0.8333 (0.2372)1.3124 (0.1843)1.7013 (0.1816)
6020II0.8259 (0.2341)1.2013 (0.1785)1.6752 (0.1807)
III0.8367 (0.2405)1.2657 (0.1841)1.6986 (0.1821)
25I0.7910 (0.2316)1.1211 (0.1742)1.6923 (0.1709)
30II0.7832 (0.2254)1.0799 (0.1699)1.3485 (0.1636)
III0.7932 (0.2332)1.1154 (0.1721)1.6000 (0.1687)
10035I0.7543 (0.1550)1.0799 (0.1174)1.6032 (0.1125)
12040II0.6815 (0.1537)1.0231 (0.1112)1.5649 (0.1062)
III0.8011 (0.1589)1.1014 (0.1217)1.6101 (0.1198)
45I0.7385 (0.1421)1.0473 (0.1069)1.5892 (0.1054)
50II0.7112 (0.1378)1.0117 (0.1052)1.5437 (0.1012)
III0.7361 (0.1436)1.0762 (0.1110)1.6045 (0.1103)
Table 8. AWs and CPs (with brackets) for the parameters with ( α , β , λ = ( 0.7 , 1 , 1.5 ) at n 1 n 2 and m 1 m 2 ( γ = 0.1 ).
Table 8. AWs and CPs (with brackets) for the parameters with ( α , β , λ = ( 0.7 , 1 , 1.5 ) at n 1 n 2 and m 1 m 2 ( γ = 0.1 ).
n 1 m 1 CS α β λ
n 2 m 2 ACISBCIPBCIACISBCIPBCIACISBCIPBCI
3010I0.95280.76000.81792.62231.12531.26381.38791.10221.2135
4015 (0.8864)(0.8872)(0.8912)(0.8804)(0.8815)(0.8872)(0.8726)(0.8810)(0.8821)
II0.84910.73060.75252.57621.10941.21141.35640.94230.9745
(0.8793)(0.8765)(0.8823)(0.8765)(0.8724)(0.8769)(0.8683)(0.8659)(0.8721)
III1.03720.79620.82132.59901.11021.24711.40110.99871.1104
(0.8910)(0.8898)(0.8913)(0.8796)(0.8793)(0.8834)(0.8795)(0.8746)(0.8793)
16I0.92000.69830.72322.10320.93450.98141.11230.89710.9214
18 (0.8845)(0.8849)(0.8842)(0.8799)(0.8829)(0.8823)(0.8728)(0.8753)(0.8812)
II0.87190.68540.71211.99670.89720.94231.10590.88650.8937
(0.8762)(0.8754)(0.8834)(0.8734)(0.8721)(0.8800)(0.8667)(0.8640)(0.8711)
III0.89890.71560.73582.11240.92910.98781.12070.89920.9148
(0.8906)(0.8900)(0.8873)(0.8812)(0.8794)(0.8907)(0.8779)(0.8801)(0.8734)
5015I0.86610.63740.63921.12310.73140.75690.79800.58960.6474
6020 (0.8976)(0.8823)(0.8859)(0.8821)(0.8797)(0.8894)(0.8809)(0.8809)(0.8902)
II0.82190.60110.62381.09890.61110.67770.72370.56190.5938
(0.8810)(0.8748)(0.8795)(0.8769)(0.8730)(0.8825)(0.8710)(0.8726)(0.8788)
III0.84320.64550.62991.12370.73680.74950.75950.60350.6371
(0.8931)(0.8847)(0.8810)(0.8824)(0.8806)(0.8876)(0.8769)(0.8864)(0.8869)
25I0.83240.60110.62130.98710.64570.65340.73510.56000.5990
30 (0.8865)(0.8789)(0.8902)(0.8812)(0.8879)(0.8907)(0.8735)(0.8798)(0.8847)
II0.81070.58460.59970.96580.58920.60040.70320.54780.5612
(0.8812)(0.8736)(0.8769)(0.8753)(0.8734)(0.8798)(0.8694)(0.8731)(0.8762)
III0.84850.61090.61441.02160.62350.62780.74440.55920.6349
(0.8927)(0.8824)(0.8873)(0.8835)(0.8812)(0.8859)(0.8738)(0.8765)(0.8924)
10035I0.51240.32830.36110.72410.46530.52240.55310.34590.3618
12040 (0.8910)(0.8875)(0.8913)(0.8967)(0.8923)(0.9078)(0.8952)(0.8924)(0.9078)
II0.48370.32410.35480.69730.43270.44190.51360.33270.3412
(0.8867)(0.8856)(0.8879)(0.8952)(0.8856)(0.8978)(0.8911)(0.8874)(0.8942)
III0.49650.33600.37160.75360.50610.53720.54290.34460.3547
(0.8881)(0.8907)(0.8942)(0.9011)(0.9108)(0.9116)(0.9035)(0.8979)(0.9011)
45I0.46110.32370.35720.68930.43250.46110.48030.33150.3600
50 (0.8889)(0.8912)(0.8931)(0.9015)(0.8897)(0.9174)(0.8977)(0.8993)(0.9057)
II0.43520.31200.33450.64790.42340.43300.47600.32440.3309
(0.8845)(0.8849)(0.8908)(0.8947)(0.8861)(0.8921)(0.8914)(0.8868)(0.8937)
III0.44580.32960.34780.66420.44960.45370.49930.32780.3469
(0.8876)(0.8924)(0.8920)(0.8969)(0.8974)(0.9123)(0.9123)(0.8895)(0.8986)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Niu, Y.; Wang, L.; Tripathi, Y.M.; Liu, J. Inference for Partially Accelerated Life Test from a Bathtub-Shaped Lifetime Distribution with Progressive Censoring. Axioms 2023, 12, 417. https://doi.org/10.3390/axioms12050417

AMA Style

Niu Y, Wang L, Tripathi YM, Liu J. Inference for Partially Accelerated Life Test from a Bathtub-Shaped Lifetime Distribution with Progressive Censoring. Axioms. 2023; 12(5):417. https://doi.org/10.3390/axioms12050417

Chicago/Turabian Style

Niu, Yingzi, Liang Wang, Yogesh Mani Tripathi, and Jia Liu. 2023. "Inference for Partially Accelerated Life Test from a Bathtub-Shaped Lifetime Distribution with Progressive Censoring" Axioms 12, no. 5: 417. https://doi.org/10.3390/axioms12050417

APA Style

Niu, Y., Wang, L., Tripathi, Y. M., & Liu, J. (2023). Inference for Partially Accelerated Life Test from a Bathtub-Shaped Lifetime Distribution with Progressive Censoring. Axioms, 12(5), 417. https://doi.org/10.3390/axioms12050417

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