Next Article in Journal
An Inventory Model with Advertisement- and Customer-Relationship-Management-Sensitive Demand for a Product’s Life Cycle
Next Article in Special Issue
Neural Networks Simulation of Distributed SEIR System
Previous Article in Journal
Class of Crosscap Two Graphs Arising from Lattices–I
Previous Article in Special Issue
Effect of Macrophages and Latent Reservoirs on the Dynamics of HTLV-I and HIV-1 Coinfection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Design for Progressive Type I Interval Censoring on the Lifetime Performance Index of Chen Lifetime Distribution

Department of Statistics, Tamkang University, Tamsui District, New Taipei City 251037, Taiwan
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(6), 1554; https://doi.org/10.3390/math11061554
Submission received: 27 February 2023 / Revised: 19 March 2023 / Accepted: 21 March 2023 / Published: 22 March 2023

Abstract

:
The lifetime performance index is commonly utilized to assess the lifetime performance of products. Based on the testing procedure for the lifetime of products following Chen distribution, an experimental design for progressive type I interval censoring is determined to achieve the desired power level while minimizing total experimental cost. For fixed inspection interval lengths and an unfixed number of inspection intervals, the required number of inspection intervals and sample sizes to achieve the minimum experimental costs are computed and presented in a table format. For unfixed termination times, the required number of inspection intervals, minimum sample sizes, and equal interval lengths are obtained and presented in a table format, while the minimum experimental costs are achieved. Finally, a practical example is presented to demonstrate the utilization of this experimental design for collecting samples and conducting a testing procedure to evaluate the lifetime performance of products.

1. Introduction

During the pandemic, there was an increase in demand for advanced technological devices, such as laptops, desktops, and mobile phones. The high lifespan of high-tech products can be a key factor in attracting more consumers and enhancing the brand’s market value. The process capability index C L was introduced by Montgomery [1] to assess the quality of larger-the-better characteristics, such as lifespan, the hardness of smart phone cases, battery capacity, and more. In many cases, the experimenters may not have access to complete data, resulting in the need to handle censored data. The two most common types of censoring are type I censoring and type II censoring. Type I censoring occurs when a life test is terminated at a fixed time point, and the number of failure units is random. Type II censoring occurs when a study is terminated when a predetermined fixed number of failure units is observed so that the termination time is random. Progressive censoring possesses the characteristic of permitting the removal of units at certain time points that may not necessarily be the ultimate termination point. More inferences about the progressive censored data can be seen in Balakrishnan and Aggarwala [2], Aggarwala [3], Balakrishnan [4], and Balakrishnan and Cramer [5]. For the progressive type II censored sample, Laumen and Cramer [6] studied the inferences for the lifetime performance index following the gamma distributions. Bdair et al. [7] studied the estimation and prediction for flexible Weibull distribution based on the progressive type II censored sample. Panahi [8] investigated the interval estimation of Kumaraswamy parameters based on progressively type II censored data and record values. EL-Sagheer [9] studied the estimation of parameters of Weibull–Gamma distribution for the progressively censored sample. Lee et al. [10] assessed the lifetime performance index for the exponential distribution model. Wu et al. [11] tested the lifetime performance index based on the Bayesian approach. The advantage of progressive type I interval censoring is that it is very convenient for quality personnel to conduct the life test and collect the censored data in practical situations. Under this type of censoring, Wu and Lin [12] used the maximum likelihood estimator for the lifetime performance index to develop a testing procedure for exponential lifetime distribution. Wu et al. [13] conducted an experimental analysis for the sampling design of Gompertz life time distribution based on progressive type I interval sampled data. For products following the Chen lifetime distribution, Wu [14] developed a testing procedure for the lifetime performance index under progressive type I interval censoring. The research goal of this study is to investigate the experimental plan for the progressive type I interval censoring design for products following the Chen distribution based on the testing procedure proposed in Wu [14]. In Section 2, we introduce and summarize a testing procedure along with the test power to assess whether the lifetime performance of a production process achieves the desired target index for the lifetime of products following the Chen distribution. Section 3 determines the minimum number of inspection intervals required to minimize the total cost under a pre-specified power level and level of significance for either a fixed or unfixed total experimental time, with the aim of achieving the lowest total cost. Additionally, one real-life example is presented to demonstrate the testing procedure. Finally, we conclude the study by summarizing all relevant findings in Section 4.

2. The Introduction of the Testing Procedure for the Lifetime Performance Index in Wu [14]

Chen [15] presented a new two-parameter lifetime distribution with a bathtub shape or increasing failure rate function called Chen distribution. Let U be the lifetime of products following a Chen distribution with the probability density function (pdf) and the failure rate function defined as:
f U ( u ) = k β u β 1 e u β exp { k ( 1 e u β ) } , 0 u , k > 0 , β > 0
and
r U ( u ) = k β u β 1 e u β .
Chen indicated that this distribution has an increasing failure rate function when β 1 and a bathtub shape failure rate function when β < 1. After the transformation from U to Y by Y = e U β 1 ,   β > 0 , the pdf of the new variable Y is an exponential distribution with failure rate k. The mean and standard deviation of Y are μ = 1 / k and σ = 1/k. If LU is the specified lower specification limit for U, then L = e L U β 1 is the specified lower specification limit for Y. The lifetime performance index proposed by Montgomery [11] is defined as:
C L = μ L σ ,
where μ is denoted as the process mean, σ is regarded as the process standard deviation, and L is the specified lower specification limit. Substituting μ and σ with the mean and standard deviation of Y into Equation (3), we obtain the lifetime performance index CL = 1 − kL.
Subsequently, the calculation of the conforming rate is performed as P r = P ( U L U ) = P ( Y L ) = exp ( k L ) = exp ( C L 1 )         ,         < C L < 1 and the value of Pr increases as CL increases. The relationship between Pr and CL is displayed in Figure 1.
To obtain the sample using the progressive type I interval censoring scheme, the following steps are followed:
Step 1: Put n products on a life test at the starting time 0. Set the termination time as T and the number of inspections as m. Then, we decide the observation time points t1, …, tm, where tm = T.
Step 2: Observe the number of failure units Xi and removed Ri units with the removing rate of pi, where Ri follows a binomial distribution denoted as bin( n j = 1 i X j j = 1 i 1 R j , pi), I = 1, …, m. Wu [14] obtained the maximum likelihood of k as the numerical solution of the following log-likelihood equation
d d k ln L ( k ) = i = 1 m X i ( y i y i 1 )   e k ( y i y i 1 ) 1 e k ( y i y i 1 ) i = 1 m ( R i   y i + X i   y i 1 )
= i = 1 m X i ( e t i β e t i 1 β )   e k ( e t i β e t i 1 β ) 1 e k ( e t i β e t i 1 β ) i = 1 m ( R i   ( e t i β 1 ) + X i   ( e t i 1 β 1 ) ) = 0
Its asymptotic variance is the reciprocal of the Fisher’s information, given by
I ( k ) = E [ d 2 ln L ( k ) d k 2 ] = n k 2 i = 1 m ln 2 ( 1 q i ) q i j = 1 i 1 ( 1 p j ) l = 1 i ( 1 q )
where q i = 1 exp (     k ( e t i β e t i 1 β ) ) .
To facilitate data collection, we considered the case of equal interval lengths t i t i 1 = t and p1 = … = pm−1 = p i = 1, …, m. That is, ti = it, i = 1, , m, and q i = 1 exp (     k ( e ( i t ) β e ( ( i 1 ) t ) β ) ) . Then, Equation (4) is reduced to
d d k ln L ( k ) = i = 1 m X i ln ( 1 q i ) ( 1 q i )   q i i = 1 m ( R i   ( e ( i t ) β 1 ) + X i   ( e ( ( i 1 ) t ) β 1 ) ) 0
The information number in (5) is reduced to
I ( k ) = n k 2 i = 1 m ln 2 ( 1 q i ) q i ( 1 p ) i 1 l = 1 i ( 1 q )
Furthermore, we have k ^ d n N ( k , g ( k ) ) where g(k) = I 1 ( k ) is the asymptotic variance of k ^ .
Due to the property of the invariance of MLE, the MLE of C L can be acquired as
C ^ L = 1 k ^ L
Let c 0   be the desired level of the lifetime performance index to make the process capable. Then, we want to test H 0 : C L c 0 (the process is not capable) vs. H a : C L > c 0 (the process is capable). Under the level of significance α , the MLE of C L found as C ^ L = 1 k ^ L is utilized as the test statistic. The critical region for this right-sided test is { C ^ L | C ^ L > C L 0 } , where the critical value C L 0 is determined as C L 0 = 1 L ( k 0 + Z 1 α g ( k 0 ) ) and Z α represents the 1 − α percentile of the standard normal distribution. Moreover, the power function denoted by h ( c 1 ) of this test at the point of C L = c 1 > c 0 is obtained as
h ( c 1 ) = Φ ( k 0 k 1 + Z 1 α g ( k 0 ) g ( k 1 ) )
where Φ ( ) is the cdf of the standard normal distribution, k 0 = 1 c 0 L and k 1 = 1 c 1 L .

3. Reliability Sampling Design

The objective of this section is to identify the optimal sampling design for progressive type I interval sampling for products’ lifetime following the Chen distribution, given that the parameters of the Chen distribution may have different structures. In Section 3.1, for the fixed experimental time T, we determine the required minimal number of inspection intervals based on the criterion of minimum total cost so that the required sample size can be calculated to reach the specified test power of the level α testing procedure. In Section 3.2, for the unfixed experimental time T, we determine the required minimal number of inspection intervals and the equal length of intervals to minimize the total experimental cost so that the required sample size can be calculated under the specified test power of the level α testing procedure.
Consider the following function w ( k ) = I ( k ) / n = 1 k 2 i = 1 m ln 2 ( 1 q i ) q i ( 1 p ) i 1 l = 1 i ( 1 q ) , which is not a function of sample size n. The power function can be rewritten as h ( c 1 ) = Φ ( k 0 k 1 + Z 1 α g ( k 0 ) g ( k 1 ) ) = Φ ( k 0 k 1 + Z 1 α w 1 ( k 0 ) / n w 1 ( k 1 ) / n ) .
In order to attain the pre-specified power 1- β or the probability of type II error β at c 1 under the level of significance α , assign the above power function to 1- β at c 1 , and then the sample size is determined as
n = c e i l i n g ( Z β w 1 ( k 1 ) + Z α w 1 ( k 0 ) k 1 k 0 ) 2
where ceiling(x) is a ceiling function mapping x to the smallest integer, which is greater than or equal to x.

3.1. The Minimal Required m for Fixed T

In numerous practical situations, the experimenters aim to minimize the number of inspection intervals m so that they do not need to frequently gather data for the progressive type I interval sampling. Suppose that the upper limit of m is m0 for experimenters, the value of m must satisfy m m 0 . If the value for m0 is not predetermined, the default value of 30 is utilized for m. In this subsection, our aim is to find the optimal value of m, denoted as m*, that minimizes the total cost incurred during the progressive type I interval censoring procedure. Similar to Huang and Wu [16], we consider the following costs:
  • Inspection cost CI: the cost for operating a single inspection station;
  • Sample cost Cs: the cost for obtaining one unit of sample;
  • Operation cost Co: the cost incurred for conducting the experiment per unit of time, which encompasses expenses such as the personnel cost, the depreciation of test equipment, and other related costs.
Taking into account all of these expenses, the overall cost of conducting this experiment is
TC ( m )   =   m C I + c e i l i n g ( Z β w 1 ( k 1 ) + Z α w 1 ( k 0 ) k 1 k 0 ) 2 C s   +   T C 0
Here is the Algorithm 1 that utilizes the numeration method to search for the optimal (m, n):
Algorithm 1: Utilize the numeration method to search for the optimal (m, n)
Step 1: Specify the pre-assigned values of m = 1, c0, c1, α, β, p, T, L, and m0 (the default value is 30) and CI = aCo, Cs = bCo, Co.
Step 2: Compute the sample size n in Equation (11) first and then compute the related total cost TC(m) in Equation (12).
Step 3: If m < m0, then m = m + 1 and go to Step 2; otherwise go to Step 4.
Step 4: For a array of total costs TC(1),…, TC(m0), The optimal solution of m* is the minimum m value, such that TC(m*) = TC* = min m m 0 TC(m), and then the related sample size n* in Equation (11) is computed.
Step 5: Calculate the value of k 0 = 1 c 0 L followed by determining the critical value C L 0 = 1 L k 0 + Z 1 ? α g k 0 .
Consider Co = 1, a = 2, and b = 1. For testing H 0 : C L 0.8 with β = 0.15, α = 0.05, p = 0.01, c1 = 0.95, L = 0.1, and T = 0.8, the curve of total cost with m = 1:m0 is displayed in Figure 2a. It can be seen that the minimum total cost occurred at m = 2, with a total cost of 13.8. For a different set up of parameters β = 0.25, α = 0.05, p = 0.05, c1 = 0.90, another curve of total cost with m = 1:m0 is displayed in Figure 2b. You can see that the minimum total cost occurred at m = 2, with a total cost of 25.8.
For testing H 0 : C L 0.80 , the required minimal inspection intervals m* and the related sample size n* to yield the minimum total cost TC(m*) with m < 50 are tabulated in Table A1 and Table A2 at the conditions of α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15, p = 0.01, 0.025, 0.05, 0.15, 0.25 for c1 = 0.825, 0.85 and c1 = 0.875, 0.90 respectively. Table A1 and Table A2 also contain the relevant critical values.
Looking at Table A2, suppose that the experimenter wants to conduct the level 0.01 hypothesis test with c0 = 0.80 and 1 − β = 0.85 at c 1 = 0.90, p = 0.05. We find that the required minimal number of inspection intervals is 2 and the sample size is determined as 14 with the minimum total cost TC* = 18.8 and the critical value of 0.877235.
From Table A1 and Table A2, the optimal number of inspection intervals m is nonincreasing when c 1 is increasing and the range of m is 2~6. In Figure 3, the plot of the minimum total cost TC* vs. c 1 for α = 0.01, 0.05, 0.1 at β = 0.25 and p = 0.05 is displayed. In Figure 4, the plot of the minimum total cost TC* vs. 1 − β = 0.75, 0.80, 0.85 at α = 0.1 and p = 0.05 is displayed. In Figure 5, the plot of the minimum total cost TC* vs. p = 0.05, 0.075, 0.1 is displayed at α = 0.1 and β = 0.25. From Figure 3, it can be observed that, as the level of significance increases, there is a decrease in the minimum total cost TC*. From Figure 4, it can be observed that, as the test power increases, there is an increase in the minimum total cost TC*. From Figure 5, it can be observed that, as the removal rate p increases, there is an increase in the minimum total cost TC*. Furthermore, these three figures show that the minimum total cost TC* is a deceasing function of c 1 . According to Table A1 and Table A2, the required minimal number of inspection intervals is inversely proportional to c 1 .

3.2. The minimal Required m, t, and n When the Interval Time of the Experiment Is Unfixed

When the equal interval time t is not fixed, we would like to determine the optimal (m,t) to yield the minimum total cost incurred for the progressive type I interval censored sampling. The total cost becomes
TC ( m , t ) =   m C I + c e i l i n g ( Z β w 1 ( k 1 ) + Z α w 1 ( k 0 ) k 1 k 0 ) 2 C s   +   mt C o .
We use the numeration method to search for the optimal (m,t), and the steps of the algorithm are as follows:
Step 1: Specify the pre-assigned values of m = 1, c0, c1, α , β and p, L and m0 (the default value is 30), and the costs CI = aCo, Cs = bCo, Co.
Step 2: Determine the optimal solution, denoted as t*, to minimize the total cost TC(m,t), as described in Equation (12). Put m = m and t = t* in Equation (11) so that the sample size n can be computed. Subsequently, the corresponding total cost TC(m,t*) in Equation (12) is computed.
Step 3: If m < m0, then let m = m + 1 and go to Step 2; otherwise go to Step 4.
Step 4: For an array of total costs TC(1,t*), …, TC(m0,t*), the optimal solution of m* is the minimum value of m such that TC(m*,t*) = TC** = min m m 0 TC(m,t*) is achieved. Put m = m* and t = t* in Equation (11), and then the related sample size n* in Equation (11) can be computed.
Step 5: Calculate the value of k 0 = 1 c 0 L followed by determining the critical value of C L 0 = 1 L ( k 0 + Z 1 α g ( k 0 ) ) .
Consider Co = 1, a = 2, and b = 1. For testing H 0 : C L 0.8 when β = 0.25, α = 0.05, p = 0.025, c1 = 0.9, m0 = 50, L = 0.05, T = 0.8, we plot m = 1:m0 against its corresponding total cost in Figure 6a. We find that the curve is a concave upward curve and the minimum total cost occurred at m = 2 with a total cost of 25.563. For another set up of parameters β = 0.15, α = 0.01, p = 0.05, c1 = 0.875, another curve of total cost with m = 1:m0 is given in Figure 6b. It can be seen that it is a concave upward curve and the minimum total cost occurred at m = 4 with a total cost of 81.962. For other combinations of setups, we can also find similar patterns.
For testing H 0 : C L 0.8 , the required minimum inspection intervals m*, the inspection interval time length t*, and sample size n* to yield the minimum total cost TC(m*,t*) are tabulated in Table A3 and Table A4 at α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15, p = 0.01, 0.025, 0.05, 0.15, 0.25 for c1 = 0.825, 0.85, and c1 = 0.875, 0.90, respectively, under the constraint of m < m0, with m0 = 50. Table A3 and Table A4 also contain the relevant critical values.
Looking at Table A4, if the experimenter wants to conduct a level 0.05 hypothesis test under a power of 0.75 at c 1 = 0.90 and p = 0.05, the minimum required sample size is obtained as 21, the minimum number of inspection intervals is obtained as 2, and the optimal inspection interval time length is 0.28. For this case, the minimum total cost is TC** = 25.55 and the relevant critical value is 0.881071.
From Table A3 and Table A4, the optimal required minimal number of inspection intervals is inversely proportional to c 1 and the range of m is 2~10. The optimal length of inspection interval t* is within 0.15 and 0.22 unit of times for c 1 = 0.825. The values of t* are within 0.18 and 0.28 units of time for c 1 = 0.875. The values of t* are within 0.21 and 0.37 units of time for c 1 = 0.875. The values of t* are within 0.23 and 0.37 units of time for c 1 = 0.90. Figure 7 displays a graph showing the relationship between the minimum total cost TC** and c 1 for α = 0.01, 0.05, 0.1 at β = 0.25, and p = 0.05. Figure 8 displays a graph showing the relationship between the minimum total cost TC** and c 1 for 1 − β = 0.75, 0.80, 0.85 at α = 0.1, and p = 0.05. Figure 9 displays a graph showing the relationship between the minimum total cost TC** and c 1 for p = 0.01, 0.025, 0.05, 0.15, 0.25 at α = 0.1 and β = 0.25. From Figure 7 and Figure 8, it can be seen that the minimum total cost TC** is a decreasing function of α or an increasing function of 1 − β . From Figure 9, it can be seen that, as the removal rate p increases, there is an increase in the minimum total cost TC**. Furthermore, the minimum total cost TC** deceases as c 1 increases.

3.3. Example

The dataset utilized in this study, obtained from Xie and Lai [17], comprises the failure times (measured in number of cycles in 100,000 times) of n = 18 electronic devices, which are provided as 0.05, 0.11, 0.21, 0.31, 0.46, 0.75, 0.98, 1.22, 1.45, 1.65, 1.95, 2.24, 2.45, 2.93, 3.21, 3.30, 3.50, and 4.20.
To test the goodness of fit of the Chen distribution, we employ the Gini statistic suggested by Gill and Gastwirth [18]. The p-value of this test is a function of β and the p-value versus the β value from 0 to 1.0 is given in Figure 10. From Figure 10, the value of β = 0.64 is determined with the largest p-value of 0.9788521. The large p-value indicated that the data fitted the Chen distribution very well. We also conducted the Kolmogorov-Smirnov test (ks. test in R) with a p-value of 0.781, which fitted the Chen distribution as well.
Using this example, the implementation of Section 3.1 and Section 3.2 is given as follows: Suppose we want to test H 0 : C L 0.8 . Refer to Section 3.1, the case of α = 0.1, the power 1 − β = 0.75 at c 1 = 0.90, p = 0.05 and T = 0.8 is considered, where the termination time of experiment T is fixed. From Table A2, we can find that the optimal sampling design is m* = 2, n* = 17 with critical value C L 0 = 0.870090 and a minimum total cost of 21.8 units under the cost setup of Co = 1, a = 2, and b = 1.
The procedure for testing is executed in the following manner:
Step 1
Take a random sample of size n = 17 from the data set. Observe the progressive type I interval censored sample (X1,X2) = (4,1) at the pre-set observation time points (t1,t2) = (0.4,0.8) with censoring schemes of (R1,R2) = (1,11).
Step 2
Obtain the MLE of k as k ^   =   0.3155534 , and then we can obtain the test statistic C ^ L = 1 k ^ L = 0.9684447.
Step 3
Compare the test statistic with the critical value. We have C ^ L = 0.9684447 > C L 0 = 0.870090. It can be inferred that the lifetime performance index of product surpasses the required level of 0.80.
Refer to Section 3.2, the case of α = 0.10, the power 1 − β = 0.85 at c 1 = 0.95 is considered. We can find that the optimal sampling design is m* = 2, n* = 16, and t* = 0.34 with critical value C L 0 = 0.9266319 and a minimum total cost of 20.682 units under the cost setup of Co = 1, a = 2, and b = 1 from our software.
The procedure for testing is executed in the following manner:
Step 1
Take a random sample of size n = 17 from the data set. Observe the progressive type I interval censored sample (X1,X2) = (4,1) at the pre-set observation time points (t1,t2) = (0.34,0.68) with censoring schemes of (R1,R2) = (0,12).
Step 2
Obtain the MLE of k as k ^   =   0.3052468 , and then we can obtain the test statistic C ^ L = 1 k ^ L = 0.9694753.
Step 3
Comparing the test statistic with the critical value, we have C ^ L = 0.9694753 > C L 0 = 0.9266319. As a result, we arrived at the same conclusion of rejecting the null hypothesis.

4. Conclusions

The evaluation of the lifetime performance index for products is a crucial subject in various manufacturing industries, particularly when the product’s lifetime follows a Chen distribution. To facilitate the collection of data, a sample was collected using the progressive type I interval censoring scheme. Our investigation aimed to determine the minimum number of inspection intervals required to achieve the given test power with a minimum total cost for a level α test when the total experimental time was fixed. When the total experimental time was not fixed, the required minimum sample size, number of inspection intervals, and the equal inspection interval time length were determined to achieve the given test power with a minimum total cost for a level α test under progressive type I interval censoring. The influences of various structures of level α , the power, and p on the minimum total cost were analyzed for the given c1 value. Nine figures for the total cost vs. c1 value in the alternative hypothesis were displayed and analyzed. We also observed that, in all cases, the minimum total cost decreased as c1 increased.

Author Contributions

Conceptualization, S.-F.W.; methodology, S.-F.W.; software, S.-F.W. and M.-Z.S.; validation, M.-Z.S.; formal analysis, S.-F.W.; investigation, S.-F.W. and M.-Z.S.; resources, S.-F.W.; data curation, S.-F.W. and M.-Z.S.; writing—original draft preparation, S.-F.W. and M.-Z.S.; writing—review and editing, S.-F.W.; visualization, M.-Z.S.; supervision, S.-F.W.; project administration, S.-F.W.; funding acquisition, S.-F.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC are funded by [National Science and Technology Council, Taiwan] NSTC 111-2118-M-032-003-MY2.

Data Availability Statement

Data available in a publicly accessible repository The data presented in this study are openly available in Xie and Lai [17].

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The optimal m*, n*, related total cost TC* and the critical value for c1 = 0.825, 0.85 and p = 0.01, 0.025, 0.05, 0.15, 0.25 under m0 = 30, L = 0.1 and c0 = 0.80.
Table A1. The optimal m*, n*, related total cost TC* and the critical value for c1 = 0.825, 0.85 and p = 0.01, 0.025, 0.05, 0.15, 0.25 under m0 = 30, L = 0.1 and c0 = 0.80.
c1 0.825 0.85
αβp m * n * TC* C L 0 m * n * TC* C L 0
0.010.250.0106745757.80.8179274177185.80.837169
0.0255760770.80.8179264179187.80.837135
0.0505777787.80.8179184181189.80.837218
0.1503828834.80.8179483193199.80.837176
0.2503863869.80.8179283201207.80.837149
0.200.0106668680.80.8189324160168.80.839094
0.0255681691.80.8189374162170.80.839035
0.0505696706.80.8189324165173.80.838981
0.1503743749.80.8189473175181.80.839041
0.2503774780.80.8189313182188.80.839040
0.150.0106605617.80.8198934147155.80.840786
0.0255617627.80.8198954148156.80.840840
0.0505631641.80.8198834151159.80.840748
0.1503673679.80.8199083160166.80.840830
0.2503701707.80.8198923167173.80.840755
0.050.250.0106465477.80.8160444108116.80.833644
0.0255474484.80.8160494109117.80.833648
0.0504487495.80.8160433113119.80.833679
0.1503517523.80.8160603118124.80.833617
0.2503538544.80.8160553123129.80.833577
0.200.0105407417.80.817209495103.80.835873
0.0255413423.80.817194496104.80.835853
0.0504423431.80.8172143100106.80.835801
0.1503450456.80.8172143104110.80.835808
0.2503468474.80.8172143108114.80.835833
0.150.0106356368.80.81833638793.80.838078
0.0255363373.80.81834038894.80.837975
0.0504373381.80.81833138995.80.837949
0.1503396402.80.81835039399.80.837866
0.2503412418.80.818346299103.80.838114
0.100.250.0105345355.80.81456338187.80.830747
0.0255350360.80.81455238187.80.830839
0.0504358366.80.81457838288.80.830804
0.1503381387.80.81457638692.80.830680
0.2503397403.80.81456229195.80.830973
0.200.0105293303.80.81580336975.80.833314
0.0255297307.80.81579737076.80.833174
0.0504305313.80.81579437177.80.833104
0.1503323329.80.81583137480.80.833074
0.2503337343.80.81580527882.80.833455
0.150.0105252262.80.81704036167.80.835431
0.0255255265.80.81704936167.80.835537
0.0504262270.80.81704136268.80.835425
0.1503278284.80.81706436470.80.835564
0.2503290296.80.81703726973.80.835570
Table A2. The optimal m*, n*, related total cost TC* and the critical value for c1 = 0.875, 0.90 and p = 0.01, 0.025, 0.05, 0.15, 0.25 under m0 = 30, L = 0.1 and c0 = 0.85.
Table A2. The optimal m*, n*, related total cost TC* and the critical value for c1 = 0.875, 0.90 and p = 0.01, 0.025, 0.05, 0.15, 0.25 under m0 = 30, L = 0.1 and c0 = 0.85.
c1 0.875 0.90
αβp m * n * TC* C L 0 m * n * TC* C L 0
0.010.250.01037581.80.85800433945.80.880437
0.02537682.80.85779334046.80.879663
0.05037682.80.85808234046.80.880060
0.15038086.80.85774324448.80.880105
0.25028589.80.85817524549.80.879954
0.200.01036975.80.86047323943.80.884001
0.02536975.80.86065433743.80.882829
0.05037076.80.86052033743.80.883242
0.15037379.80.86044824145.80.882984
0.25027882.80.86073024145.80.883764
0.150.01036470.80.86279133440.80.886148
0.02536470.80.86297923741.80.886359
0.05036571.80.86280423741.80.886555
0.15036874.80.86263123842.80.886197
0.25027276.80.86321023943.80.885885
0.050.250.01034551.80.85294622529.80.874182
0.02534551.80.85310422529.80.874283
0.05034652.80.85278622529.80.874452
0.15025054.80.85313122630.80.873680
0.25025155.80.85310322630.80.874373
0.200.01034046.80.85615822327.80.877341
0.02524347.80.85664022327.80.877446
0.05034147.80.85591222327.80.877622
0.15024549.80.85600522327.80.878338
0.25024650.80.85591422428.80.877410
0.150.01033642.80.85919522125.80.880940
0.02523943.80.85947422125.80.881050
0.05033743.80.85885722125.80.881234
0.15024044.80.85940322125.80.881984
0.25024145.80.85922622226.80.880852
0.100.250.01023438.80.84956121721.80.870090
0.02533238.80.84906521721.80.870185
0.05023539.80.84902521822.80.868363
0.15023539.80.84947821822.80.868994
0.25023640.80.84924521822.80.869642
0.200.01023034.80.85276221620.80.872247
0.02523034.80.85283321620.80.872345
0.05023135.80.85209221620.80.872510
0.15023135.80.85257321620.80.873179
0.25023236.80.85223221620.80.873867
0.150.01022731.80.85561621418.80.877235
0.02522731.80.85569121418.80.877340
0.05022731.80.85581821418.80.877516
0.15022832.80.85531821519.80.875579
0.25022832.80.85583821519.80.876289
Table A3. The optimal (m*,t*), n*, total cost TC** and the critical value for c1 = 0.825,0.85 and p = 0.01, 0.025, 0.05, 0.15, 0.25 under m0 = 30, L = 0.1 and c0 = 0.80.
Table A3. The optimal (m*,t*), n*, total cost TC** and the critical value for c1 = 0.825,0.85 and p = 0.01, 0.025, 0.05, 0.15, 0.25 under m0 = 30, L = 0.1 and c0 = 0.80.
c1 0.825 0.85
αβp m * t * n * T C * * C L 0 m * t * n * T C * * C L 0
0.010.150.010100.14702723.360.81797860.18167180.080.837372
0.02580.17719736.380.81798550.21171182.070.837391
0.05070.20739754.370.81798250.22174185.090.837367
0.15050.27796807.340.81798440.27186195.090.837368
0.25050.29836847.440.81797340.31194203.240.837387
0.200.01090.15632651.340.81898460.20151164.180.839282
0.02590.16643662.430.81898350.22155166.100.839262
0.05070.20663678.380.81898550.21158169.070.839238
0.15060.25712725.500.81898040.27169178.060.839203
0.25050.29750761.440.81897540.31176185.260.839252
0.250.01090.15573592.350.81993750.20141152.010.840983
0.02580.17585602.380.81993850.22142153.120.841020
0.05070.18602617.300.81993850.21145156.040.840959
0.15050.26648659.280.81993040.26155164.040.840953
0.25050.29680691.430.81992830.31164170.920.840931
0.050.150.01080.16441458.300.81611150.22103114.090.833883
0.02580.17448465.360.81610950.23104115.140.833897
0.05060.21463476.280.81611140.25108117.000.833837
0.15050.26496507.320.81610730.30116122.910.833820
0.25040.30523532.190.81609830.33120126.980.833832
0.200.01080.16384401.300.81726550.2191102.060.836050
0.02580.17390407.390.81726640.2494102.980.836036
0.05070.19401416.350.81726240.2396104.900.835931
0.15050.26432443.300.81725940.28100109.120.836031
0.25040.31455464.230.81726130.32106112.950.835992
0.250.01070.17341356.190.81839040.218492.850.838071
0.02570.17346361.220.81839340.248492.950.838121
0.05070.20353368.410.81839640.258594.000.838141
0.15050.25381392.240.81838130.279298.820.838054
0.25040.30401410.190.81838430.3095101.890.838044
0.100.150.01070.17327342.200.81463240.227785.890.830936
0.02570.17332347.210.81462940.257785.990.831018
0.05060.20341354.220.81463040.267887.020.831024
0.15050.26365376.300.81462930.308490.900.830965
0.25040.29385394.180.81462130.328793.950.830953
0.200.01070.17278293.200.81586940.246674.950.833370
0.02560.20284297.220.81587430.296975.870.833391
0.05060.20290303.210.81586530.277076.820.833330
0.15040.28313322.120.81586830.287379.830.833249
0.25040.30327336.200.81586230.327581.950.833338
0.250.01070.18239254.230.81711230.286066.830.835706
0.02560.19245258.110.81709640.245866.970.835744
0.05050.22252263.100.81710530.286167.830.835685
0.15050.26267278.280.81710430.326369.950.835767
0.25040.28282291.130.81709330.296672.860.835591
Table A4. The optimal (m*,t*), n*, total cost TC** and the critical value for c1 = 0.875, 0.90 and p = 0.01, 0.025, 0.05, 0.15, 0.25 under m0 = 30, L = 0.1 and c0 = 0.80.
Table A4. The optimal (m*,t*), n*, total cost TC** and the critical value for c1 = 0.875, 0.90 and p = 0.01, 0.025, 0.05, 0.15, 0.25 under m0 = 30, L = 0.1 and c0 = 0.80.
c1 0.875 0.90
αβp m * t * n * T C * * C L 0 m * t * n * T C * * C L 0
0.010.150.01040.217280.860.85815830.263945.780.880480
0.02540.247280.950.85823530.283945.830.880627
0.05040.247381.960.85823030.254046.740.880226
0.15030.287985.830.85801930.294147.860.880481
0.25030.338187.980.85823920.314549.630.880063
0.200.01040.226674.860.86065630.283642.830.883677
0.02540.246674.960.86082530.233743.700.883314
0.05040.246775.960.86078130.253743.760.883415
0.15030.297278.880.86073230.293844.870.883598
0.25030.297581.870.86060720.374145.740.883592
0.250.01040.226169.880.86309330.253440.750.886296
0.02530.266470.770.86301630.273440.800.886387
0.05040.246270.970.86318520.323741.650.886264
0.15030.286773.830.86300120.313842.630.886147
0.25030.316975.930.86310320.313943.620.886001
0.050.150.01030.254551.750.85303720.302529.610.874037
0.02530.264551.780.85313620.312529.620.874068
0.05030.294551.870.85331820.322529.640.874201
0.15030.294753.860.85314820.302630.590.873748
0.25020.345155.680.85301120.342630.690.874244
0.200.01030.264046.770.85618820.282327.550.877466
0.02530.274046.820.85631320.282327.560.877600
0.05030.254147.740.85602820.292327.580.877662
0.15030.294248.860.85622320.342327.680.878102
0.25020.314650.630.85599020.302428.610.877649
0.250.01030.273642.800.85918420.282125.550.881071
0.02520.343943.680.85920520.282125.560.881211
0.05030.253743.740.85897920.292125.580.881276
0.15020.354044.710.85921320.352125.700.881722
0.25020.354145.690.85910020.302226.590.881101
0.100.150.01020.373438.750.84938620.371721.740.869842
0.02530.263238.790.84909420.271822.540.868511
0.05030.293238.860.84926220.281822.560.868540
0.15020.393539.780.84942120.311822.630.868954
0.25020.363640.720.84913520.371822.740.869499
0.200.01020.343034.670.85251320.261620.520.872750
0.02520.353034.700.85259820.271620.530.872667
0.05030.292834.870.85266320.271620.550.872881
0.15020.343135.670.85241520.311620.620.873136
0.25020.323236.640.85223120.371620.730.873715
0.250.01020.292731.590.85559520.291418.580.877207
0.02520.302731.600.85559920.301418.590.877212
0.05020.312731.630.85568220.311418.620.877327
0.15020.292832.580.85547310.431719.430.877249
0.25020.342832.680.85574110.431719.430.877249

References

  1. Montgomery, D.C. Introduction to Statistical Quality Control; John Wiley and Sons Inc.: New York, NY, USA, 1985. [Google Scholar]
  2. Balakrishnan, N.; Aggarwala, R. Progressive Censoring: Theory, Methods and Applications; Birkhäuser: Boston, MA, USA, 2000. [Google Scholar]
  3. Aggarwala, R. Progressive interval censoring: Some mathematical results with applications to inference. Commun. Stat.-Theory Methods 2001, 30, 1921–1935. [Google Scholar]
  4. Balakrishnan, N. Progressive Censoring Methodology: An Appraisal (with Discussions). Test 2007, 16, 211–296. [Google Scholar] [CrossRef]
  5. Balakrishnan, N.; Cramer, E. The Art of Progressive Censoring; Birkhäuser: Boston, MA, USA, 2014. [Google Scholar]
  6. Laumen, B.; Cramer, E. Inference for the lifetime performance index with progressively Type-II censored samples from gamma distributions. Econ. Qual. Control 2015, 30, 59–73. [Google Scholar] [CrossRef]
  7. Bdair, O.M.; Abu Awwad, R.R.; Abufoudeh, G.K.; Naser, M.F.M. Estimation and Prediction for Flexible Weibull Distribution Based on Progressive Type II Censored Data. Commun. Math. Stat. 2020, 8, 255–277. [Google Scholar] [CrossRef]
  8. Panahi, H. Interval estimation of Kumaraswamy parameters based on progressively type II censored sample and record values. Miskolc Math. Notes 2020, 21, 319. [Google Scholar] [CrossRef]
  9. EL-Sagheer, R.M. Estimation of parameters of Weibull–Gamma distribution based on progressively censored data. Stat. Pap. 2018, 59, 725–757. [Google Scholar] [CrossRef] [Green Version]
  10. Lee, W.C.; Wu, J.W.; Hong, C.W. Assessing the lifetime performance index of products with the exponential distribution under progressively type II right censored samples. J. Comput. Appl. Math. 2009, 231, 648–656. [Google Scholar] [CrossRef] [Green Version]
  11. Wu, J.W.; Lee, W.C.; Lin, L.S.; Hong, M.L. Bayesian test of lifetime performance index for exponential products based on the progressively type II right censored sample. J. Quant. Manag. 2011, 8, 57–77. [Google Scholar]
  12. Wu, S.F.; Lin, Y.P. Computational testing algorithmic procedure of assessment for lifetime performance index of products with one-parameter exponential distribution under progressive type I interval censoring. Math. Comput. Simul. 2016, 120, 79–90. [Google Scholar] [CrossRef]
  13. Wu, S.F.; Liu, T.H.; Lai, Y.H.; Chang, W.T. A study on the experimental design for the lifetime performance index of Rayleigh lifetime distribution under progressive type I interval censoring. Mathematics 2022, 10, 517. [Google Scholar] [CrossRef]
  14. Wu, S.F. The performance assessment on the lifetime performance index of products following Chen lifetime distribution based on the progressive type I interval censored sample. J. Comput. Appl. Math. 2018, 334, 27–38. [Google Scholar] [CrossRef]
  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. Huang, S.R.; Wu, S.J. Reliability sampling plans under progressive type-I interval censoring using cost functions. IEEE Trans. Reliab. 2008, 57, 445–451. [Google Scholar] [CrossRef]
  17. Xie, M.; Lai, C.D. Reliability analysis using an additiveWeibull model with bathtub-shaped failure rate function. Reliab. Eng. Syst. Saf. 1995, 52, 87–93. [Google Scholar] [CrossRef]
  18. Gail, M.H.; Gastwirth, J.L. A scale-free goodness of fit test for the exponential distribution based on the Gini Statistic. J. R. Stat. Soc. B 1978, 40, 350–357. [Google Scholar] [CrossRef]
Figure 1. The relationship between Pr and CL.
Figure 1. The relationship between Pr and CL.
Mathematics 11 01554 g001
Figure 2. (a): Total cost curve at m = 1:m0. (b) Total cost curve at m = 1:m0.
Figure 2. (a): Total cost curve at m = 1:m0. (b) Total cost curve at m = 1:m0.
Mathematics 11 01554 g002
Figure 3. The minimum total cost curve for α = 0.01, 0.05, 0.1.
Figure 3. The minimum total cost curve for α = 0.01, 0.05, 0.1.
Mathematics 11 01554 g003
Figure 4. The minimum total cost for 1 − β = 0.75, 0.80, 0.85.
Figure 4. The minimum total cost for 1 − β = 0.75, 0.80, 0.85.
Mathematics 11 01554 g004
Figure 5. The minimum total cost for p = 0.01, 0.025, 0.05.
Figure 5. The minimum total cost for p = 0.01, 0.025, 0.05.
Mathematics 11 01554 g005
Figure 6. (a) The minimum total cost at m = 1:m0. (b) The minimum total cost at m = 1:m0.
Figure 6. (a) The minimum total cost at m = 1:m0. (b) The minimum total cost at m = 1:m0.
Mathematics 11 01554 g006
Figure 7. The minimum total cost curve or α = 0.01, 0.05, 0.1.
Figure 7. The minimum total cost curve or α = 0.01, 0.05, 0.1.
Mathematics 11 01554 g007
Figure 8. The minimum total cost curve or 1 − β = 0.75, 0.80, 0.85.
Figure 8. The minimum total cost curve or 1 − β = 0.75, 0.80, 0.85.
Mathematics 11 01554 g008
Figure 9. The minimum total cost curve for p = 0.01, 0.025, 0.05.
Figure 9. The minimum total cost curve for p = 0.01, 0.025, 0.05.
Mathematics 11 01554 g009
Figure 10. The p-value vs. the β values.
Figure 10. The p-value vs. the β values.
Mathematics 11 01554 g010
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

Wu, S.-F.; Song, M.-Z. Experimental Design for Progressive Type I Interval Censoring on the Lifetime Performance Index of Chen Lifetime Distribution. Mathematics 2023, 11, 1554. https://doi.org/10.3390/math11061554

AMA Style

Wu S-F, Song M-Z. Experimental Design for Progressive Type I Interval Censoring on the Lifetime Performance Index of Chen Lifetime Distribution. Mathematics. 2023; 11(6):1554. https://doi.org/10.3390/math11061554

Chicago/Turabian Style

Wu, Shu-Fei, and Meng-Zong Song. 2023. "Experimental Design for Progressive Type I Interval Censoring on the Lifetime Performance Index of Chen Lifetime Distribution" Mathematics 11, no. 6: 1554. https://doi.org/10.3390/math11061554

APA Style

Wu, S. -F., & Song, M. -Z. (2023). Experimental Design for Progressive Type I Interval Censoring on the Lifetime Performance Index of Chen Lifetime Distribution. Mathematics, 11(6), 1554. https://doi.org/10.3390/math11061554

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