Next Article in Journal
Secure Data Transmission and Image Encryption Based on a Digital-Redesign Sliding Mode Chaos Synchronization
Next Article in Special Issue
Multiple Comparison Procedures for Exponential Mean Lifetimes Compared with Several Controls
Previous Article in Journal
Connections between Campos-Bolanos and Murofushi–Sugeno Representations of a Fuzzy Measure
Previous Article in Special Issue
A Note on a Meshless Method for Fractional Laplacian at Arbitrary Irregular Meshes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study on the Experimental Design for the Lifetime Performance Index of Rayleigh Lifetime Distribution under Progressive Type I Interval Censoring

1
Department of Statistics, Tamkang University, Tamsui, Taipei 251301, Taiwan
2
Department of Computer Science, University of Taipei, Taipei 100234, Taiwan
*
Author to whom correspondence should be addressed.
Submission received: 29 December 2021 / Revised: 1 February 2022 / Accepted: 3 February 2022 / Published: 5 February 2022
(This article belongs to the Special Issue Applications of Partial Differential Equations in Engineering)

Abstract

:
With the rapid development of technology, improving product life performance has become a very important issue in recent decades. The lifetime performance index is used in this research for the assessment of the lifetime performance of products following the Rayleigh distribution. Based on the hypothesis testing procedure with this index, using the maximum likelihood estimator as a testing statistic, the sampling design is determined and the related values are tabulated for practical use to reach the given power level and minimize the total experimental cost under progressive type I interval censoring. When the inspection interval length is fixed and the number of inspection intervals is not fixed, the required number of inspection intervals and sample size with the minimum total cost are determined and tabulated. When the termination time is not fixed, the required number of inspection intervals, sample size and equal interval length reaching the minimum total cost are determined and tabulated. Lastly, a practical example is given to illustrate the use of this sampling design for the testing procedure to determine whether the process is capable.

1. Introduction

The process capability index C L proposed by Montgomery [1] is frequently used to measure the larger-the-better quality characteristics such as lifetime, mpg, tensile strength, durability, etc. For measuring the lifetime of products, this index is called the lifetime performance index. In many cases, the experimenters can only observe censored data. Two censoring types, namely type I censoring and type II censoring, are frequently considered. Type I censoring occurs if the life test of n subjects stops at a predetermined time and the number of observations is random. Type II censoring occurs if the life test stops when a predetermined number of failure times is observed. Progressive censoring has the property of allowing the removal of units at some time points that may not be the final termination point. More inferences about progressive censored data can be seen in Balakrishnan and Aggarwala [2] and Aggarwala [3]. For a progressive type II right censored sample, Wu et al. [4] assessed the lifetime performance index of products with an exponential distribution. Wu et al. [5] proposed a Bayesian test of the lifetime performance index for exponential products based on a progressive type II right censored sample. Lee et al. [6] proposed a decision procedure for the performance assessment of the lifetime index of products for the Gompertz distribution under progressive type II right censoring. Lee et al. [7] assessed the lifetime performance index of Rayleigh products based on Bayesian estimation using progressive type II right censored samples. Wu et al. [8] implemented a testing procedure for the lifetime performance index of Burr XII products with a progressive type II right censored sample. For progressive type I interval censored data from the Gompertz lifetime distribution, a testing procedure for the lifetime performance index was proposed by Wu and Hsieh [9] based on a progressive type I interval censored sample. Based on this testing procedure, a reliability sampling design was developed by Wu et al. [10] for products following Gompertz distribution. For products following Weibull lifetime distribution, Wu and Lin [11] proposed a hypothesis testing procedure for the lifetime performance index using progressive type I interval censored data. Wu et al. [12] investigated the sampling design for the testing of the lifetime performance index from Weibull distribution. This study focuses on the Rayleigh distribution, which is an asymmetric probability distribution. Wu et al. [13] proposed a testing procedure to test whether the lifetime performance index meets the desired target for products following Rayleigh distribution under progressive type I interval censoring, and the related testing procedure is summarized in Section 2. The progressive type I interval censoring scheme is depicted as follows: there are n products that are subjected to a life test at time 0. Let (t1, …, tm) and (p1, …, pm) be the predetermined inspection times and the pre-specified removal percentages of the remaining survival units at time ti, i = 1, …, m, where pm = 1 and tm = T is the termination time for the experiment. At the ith inspection time interval (ti−1,ti), the number of failure units Xi follows a binomial distribution bin ( n j = 1 i 1 X j j = 1 i 1 R j ,qi), where q i = F U ( t i ) F U ( t i 1 ) 1 F U ( t i 1 ) and FU(x) is the cumulative distribution function for the lifetime variable U. Furthermore, Ri units are randomly removed from the remaining n j = 1 i X j j = 1 i 1 R j survival units and Ri follows a binomial distribution bin ( n j = 1 i X j j = 1 i 1 R j ,pi), i = 1, …, m. Then, the observed number of failure units (X1, …, Xm) for m inspection intervals is the progressive type I interval censored sample under a progressive censoring scheme (R1, …, Rm) with removal probabilities p1, …, pm. There is an increasing relationship between the lifetime performance index and the conforming rate. If the desired conforming rate is specified, the desired lifetime performance index c0 is determined. The null hypothesis H 0 : C L c 0 vs. the alternative hypothesis H a : C L > c 0 are set up. Using the maximum likelihood estimator of the lifetime performance index, the testing procedure for Rayleigh lifetime products under a level of significance α was developed by Wu et al. [13]. The power function has also been obtained. Further details of the methodologies for the testing procedure are given in Section 2. Based on the testing procedure proposed in Wu et al. [13], we conduct a study on the experimental design for the lifetime performance index from Rayleigh distributed products based on progressive type I interval censored data. The optimal sampling design consists of two parts and is given in Section 3. The first part is to determine the minimum number of inspection intervals yielding the minimum total cost under the pre-specified power level and level of significance for a fixed total experimental time, which is described in Section 3.1. The second part considers the case when the interval time of the experiment is not fixed. Under this consideration, the minimum number of inspection intervals and the equal interval length reaching the minimum total cost are determined and tabulated in Section 3.2. In Section 3.3, a practical example is given to illustrate the application of the optimal sampling design for the testing procedure to test whether the process is capable. Finally, the conclusions are presented in Section 4.

2. Methodology for Testing the Lifetime Performance Index

The lifetime U of a product following Rayleigh distribution has a probability density function (pdf) and a cumulative distribution function (cdf) as follows:
f U ( u ) = u λ 2 exp { u 2 2 λ 2   } ,   u 0 , λ > 0
and
F U ( u ) = 1 exp { u 2 2 λ 2   } ,   u 0 , λ > 0 .
The failure rate function is defined as
r U ( u ) = f U ( u ) 1 F U ( u ) = u λ 2 .  
The curves for the pdf and failure rate function under λ = 1, 1.5, 2, 3, 5 are displayed in Figure 1a,b, respectively.
Let Y = U 2 . Then, this new random variable Y has a one-parameter exponential distribution and its pdf and cdf are given as follows:
f Y ( y ) = 1 2 λ 2 exp { y 2 λ 2 } ,   y 0 , λ > 0
and
F Y ( y ) = 1 exp { y 2 λ 2 } ,   y 0 , λ > 0 ,  
where λ 2 is the shape parameter and the failure rate function is defined as
r Y ( y ) = f y ( y ) 1 F y ( y ) = 1 2 λ 2 ,  
where λ 2 is the scale parameter.
From Montgomery [1], the lifetime performance index is defined as:
C L = μ L σ ,
where μ denotes the process mean, σ represents the process standard deviation, and L is the known lower specification limit. Suppose that the lower specification limit for U is LU. Then, the lower specification limit for Y is L = LU2.
The mean and standard deviation of Y are μ = 2 λ 2 and σ = 2 λ 2 . Substituting μ and σ by this information, we then have C L = μ L σ = 2 λ 2 L 2 λ 2 = 1 L 2 λ 2 .
The conforming rate is obtained as
P r = P ( U L U ) = P ( Y L ) = exp ( L 2 λ 2 ) = exp ( C L 1 ) ,         < C L < 1 .
The conforming rate is an increasing function of the lifetime performance index C L . If the experimenter desires the conforming rate to be 0.8607, the desired lifetime performance index should be c0 = 0.85.
For a progressive type I interval censoring sample (X1, …, Xm) under a progressive censoring scheme (R1, …, Rm) with removal probabilities p1, …, pm, the likelihood function is
L ( λ ) = i = 1 m ( e t i 1 2 x i 2 λ 2 e t i 2 x i 2 λ 2 ) e t i 2 R i 2 λ 2 .
Wu et al. [13] obtained the maximum likelihood estimator (MLE) for λ as the numerical solution of the following log-likelihood equation:
d d λ ln L ( λ ) = i = 1 m ( x i exp { t i 2 t i 1 2 2 λ 2 } 1 exp { t i 2 t i 1 2 2 λ 2 } ( t i 2 t i 1 2 λ 2 ) + t i 1 2 x i + t 1 2 R i λ 3 ) = 1 λ 3 i = 1 m ( x i ( t i 2 t i 1 2 ) ( exp { t i 2 t i 1 2 2 λ 2 } 1 ) 2 ( t i 1 2 x i + t 1 2 R i ) ) .
Its asymptotic variance is the reciprocal of the following Fisher’s information:
I ( λ ) = E [ d 2 ln L ( λ ) d λ 2 ] = n λ 2 i = 1 m { 1 q i q i [ 4 ln 2 ( 1 q i ) + 6 q i ln ( 1 q i ) ] + 3 ( t i 1 2 q i + t i 2 p i ( 1 q i ) ) λ 2 } j = 1 i 1 ( 1 p j ) ( 1 q j ) ;
where q i = 1 exp ( t i 2 t i 1 2 2 λ 2 ) .
Then, we have λ ^ d n N ( λ , I 1 ( λ ) ) .
For equal interval lengths, t i t i 1 = t and t i = it, i = 1 , , m   are replaced in Equations (8) and (9).
By the property of the invariance of MLE, the MLE of C L can be obtained as
C ^ L = 1 L 2 λ ^ . C ^ L = 1 L 2 λ ^ 2
Let c0 be the desired level of the lifetime performance index to make the process capable. Then, we wish to test H 0 : C L c 0 (the process is not capable) vs. H a : C L > c 0 (the process is capable) using the MLE of C L , C ^ L = 1 L 2 λ ^ 2 . as the test statistic, and the critical region for this test is { C ^ L | C ^ L > C L 0 } . From Wu et al. [9], the critical value is found to be C L 0 = 1 L 2 ( λ 0 + Z α I 1 ( λ 0 ) ) 2 , under the level of significance α , where λ 0 = L 2 ( 1 c 0 ) and Z α represents the right-tailed α percentile of standard normal distribution. We will conclude that the alternative hypothesis is supported if C ^ L > C L 0 is satisfied.
Moreover, the power h ( c 1 ) of this statistical test at the point of C L = c 1 > c 0 is
h ( c 1 ) = P ( C ^ L > C L 0 | c 1 = 1 L λ 1 2 ) = Φ ( λ 0 λ 1 + Z 1 α I 1 ( λ 0 ) I 1 ( λ 1 ) ) ,  
where Φ ( ) is the cdf for the standard normal distribution, λ 0 = L 2 ( 1 c 0 ) and λ 1 = L 2 ( 1 c 1 ) .

3. Sampling Design

In this section, the reliability sampling design is investigated under different setups and considerations. As described in Section 3.1, when the termination time of the experiment T is fixed, the required sample size and minimal number of inspection intervals are determined so that the given power of the level α testing procedure can be reached and the total cost of the experiment can be minimized. In Section 3.2, when the termination time of the experiment T is fixed, the required sample size, minimal number of inspection intervals and the equal lengths of the interval are determined so that the given power of the level α testing procedure can be reached and the total cost of the experiment can be minimized. Let
g ( λ ) = I 1 ( λ ) / n = λ 2 [ i = 1 m { 1 q i q i [ 4 ln 2 ( 1 q i ) + 6 q i ln ( 1 q i ) ] + 3 ( t i 1 2 q i + t i 2 p i ( 1 q i ) ) λ 2 } j = 1 i 1 ( 1 p j ) ( 1 q j ) ] 1
The function g ( λ ) is independent of the sample size n. In terms of g ( λ ) , the power function can be rewritten as h ( c 1 ) = Φ ( ( λ 0 λ 1 ) n + Z 1 α g ( λ 0 ) g ( λ 1 ) ) .
In order to attain the pre-specified power 1 − β or the probability of type II error β at c 1 under the level of significance α, we set the above power function to be 1 − β, and then the sample size is determined as
n = c e i l i n g ( Z β g ( λ 1 ) + Z α g ( λ 0 ) λ 0 λ 1 ) 2   ,  
where ceiling(x) is a ceiling function mapping x to the smallest integer greater than or equal to x.

3.1. The Optimal m and n for Fixed T

The smaller the number of inspection intervals m, the more convenient it is for experimenters to collect the progressive type I interval sample. Let m0 be the upper limit of m specified by the experimenter such that m m 0 . The default value of m0 is 20 if the experimenter does not specify it in advance. The sample size n is a function of m. In this section, we wish to determine the optimal (m,n) to yield the minimum total cost incurred during the progressive type I interval censoring procedure. Similar to Huang and Wu [14], we consider the following costs:
  • Installation cost Ca—the cost of installing all test units;
  • Sample cost Cs—the cost per test unit;
  • Inspection cost CI—the cost of using the inspection equipment;
  • Operation cost Co—the cost consisting of the personnel cost, depreciation of test equipment and so on. It is proportional to the length of the experimental time period.
Integrating all these costs, we have the total cost of
TC(m,n) = Ca + nCs + m CI +T Co,
where n is given in Equation (12).
The Algorithm 1 using the numeration method to search the optimal (m,n) is given as follows:
Algorithm 1: Search the optimal (m,n)
Step 1:
Give the pre-assigned values of m = 1, c0, c1, α, β, p, T, L and m0 (the default value is 20) and the four costs Ca = aCo, Cs = bCo, CI = cCo, Co.
Step 2:
Compute the sample size n in Equation (12) and then compute the total cost TC(m,n) in Equation (13).
Step 3:
If m < m0, then m = m + 1 and go to Step 2; otherwise, go to Step 4.
Step 4:
For an array of total costs, the optimal solution of m value denoted by m* is the m value with the minimum total cost. Then, the sample size n* in Equation (12) can be computed.
Once the optimal values of m* and n* are determined, the critical value of C L 0 = 1 L 2 ( λ 0 + Z α I 1 ( λ 0 ) ) 2 can be obtained as well. Consider c0 = 0.85, Co = 1 and a = b = c = 1. For testing H 0 : C L 0.85 when β = 0.25, α = 0.05, p = 0.05, c1 = 0.90, m0 = 20, L = 0.05, T = 1.0, we plot the total cost versus m = 1:m0 in Figure 2a. From this figure, we can see that the minimum total cost is incurred when m = 3, with the total cost of 43. For a different setup of parameters, β = 0.15, α = 0.05, p = 0.01, c1 = 0.95, another total cost curve with m = 1:m0 is given in Figure 2b. From this figure, we can see that the curve is a concave upwards curve and the minimum total cost occurs when m = 2, with the total cost of 13.
The minimum suggested inspection intervals m* and the corresponding sample size n* to attain the minimum total cost TC (m*,n*) under the constraint of m < 20 for testing H 0 : C L 0.85 are tabulated in Table 1 and Table 2 for c1 = 0.875, 0.90 and c1 = 0.925, 0.95, respectively, under the condition of α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15, L = 0.05, T = 1.0, p = 0.05, 0.075, 0.1. The corresponding critical values are also tabulated in these two tables.
For example, suppose that the users wish to conduct the level 0.05 hypothesis testing of H 0 : C L 0.85 under the power of 0.85 at c 1 = 0.90, p = 0.05 and m0 = 20. According to Table 1, the minimum required sample size is 53 and the minimum number of inspection intervals is 3. The minimum total cost can also be found as TC = 58 and the critical value is 0.881256.
According to Table 1 and Table 2, the optimal number of inspection intervals m is nonincreasing when c 1 is increasing and the range of m is 1~4. In Figure 3, we plot the minimum total cost TC* vs. c 1 for α = 0.01, 0.05, 0.1 at β = 0.25 and p = 0.05. In Figure 4, we plot the minimum total cost TC* vs. c 1 for 1 − β = 0.75, 0.80, 0.85 at α = 0.1 and p = 0.05. In Figure 5, we plot the minimum total cost TC* vs. c 1 for p = 0.05, 0.075, 0.1 at α = 0.1 and β = 0.25. From Figure 3, we have the finding that the minimum total cost TC* decreases when the level of significance increases. From Figure 4, we find that the minimum total cost TC* increases when the test power increases. From Figure 5, we find that the minimum total cost TC* increases when the removal probability increases. Overall, the minimum total cost TC* becomes smaller when the value of c 1 becomes larger.

3.2. The Optimal m, t and n When the Interval Time of the Experiment Is Not Fixed

In this subsection, the interval time of experiment t is not fixed. The inspection time t for each inspection interval is random. We wish to determine the optimal (m,t,n) to yield the minimum total cost incurred during the type I interval censoring procedure. The total cost is
TC(m,t,n) = Ca + nCs + m CI + mt Co,
where n is given in Equation (12).
The Algorithm 2 using the numeration method to search the optimal (m,t,n) is given as follows:
Algorithm 2: Search the optimal (m,t,n)
Step 1:
Give the pre-assigned values of m = 1, c0, c1, α, β and p, L and m0 (the default value is 20) and the four costs Ca = aCo, Cs = bCo, CI = cCo, Co.
Step 2:
Find the optimal value of t* such that TC(m,t,n) in Equation (14) is minimized (the optimize function in R is used to find t*). Use t = t* to compute the sample size n in Equation (12) and compute the corresponsong total cost TC(m,t*,n) in Equation (14).
Step 3:
If m < m0, then m = m + 1 and go to Step 2; otherwise, go to Step 4.
Step 4:
For an array of total costs, we determine the optimal value of m denoted by m* such that TC(m*,t*,n) = TC* *= min m m 0 TC(m,t*,n) is reached. With m = m* and t = t*, the related sample size n* in Equation (12) can be computed.
Once the optimal values of m*, n* and t* are determined, the critical value of C L 0 = 1 L 2 ( λ 0 + Z α I 1 ( λ 0 ) ) 2 can be determined as well.
Consider c0 = 0.85, Co = 1 and a =b = c = 1. When β = 0.15, α = 0.05, p = 0.1, c1 = 0.9, m0 = 20, L = 0.05, T = 1.0, we plot m = 1:m0 against its corresponding total cost in Figure 6a. We can see that it is a concave upwards curve and the minimum total cost occurs when m = 3, with the total cost of 58.14254. For a different setup of parameters, β = 0.25, α = 0.05, p = 0.1, c1 = 0.875, another total cost curve with m = 1:m0 is given in Figure 6b. We can see that this is a concave upwards curve and the minimum total cost is incurred when m = 4, with the total cost of 198.2394. For other combinations of setups, we also find similar patterns.
The minimum suggested inspection intervals m*, inspection interval time length t* and sample size n* to yield the minimum total cost TC(m*,t*,n*) denoted by TC** under the constraint of m < m0, where m0 = 20 for testing H 0 : C L 0.85 , are tabulated in Table 3 and Table 4 at α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15, L = 0.05, T = 1.0, p = 0.05, 0.075, 0.1 for c1 = 0.875, 0.90 and c1 = 0.925, 0.95, respectively. The corresponding critical values are also tabulated.
For example, suppose that the user wishes to conduct the level 0.05 hypothesis testing of H 0 : C L 0.85 under power of 0.75 at c 1 = 0.875, p = 0.05 and m0 = 20. According to Table 3, the minimum required sample size is 245, the minimum number of inspection intervals is 5 and the inspection interval time length is 0.25. The total cost is calculated as TC = 252.254 and the critical value is 0.8656.
For any other setup of testing procedure, a software program is provided by the authors for users to input L, T, m, c0, c1, α, β to output the minimum required sample size n; or input L, T, c0, c1, α, β to output the minimum suggested number of inspections m; or input L, m, c0, c1, α, β to output the minimum recommended number of inspections m and the equal length of time t for each inspection interval.
According to Table 3 and Table 4, the optimal number of inspection intervals m is nonincreasing when c 1 is increasing and the range of m is 2~6. The optimal length of inspection interval t* is within 0.24 and 0.32 units of time for c 1 = 0.875. The value of t* is within 0.28 and 0.39 units of time for c 1 = 0.90. The value of t* is within 0.33 and 0.55 units of time for c 1 = 0.925. The value of t* is within 0.45 and 0.74 units of time for c 1 = 0.95. In Figure 7, we plot the minimum total cost TC* vs. c 1 for α = 0.01, 0.05, 0.1 at β = 0.25 and p = 0.05. In Figure 8, we plot the minimum total cost TC* vs. c 1 for 1 − β = 0.75, 0.80, 0.85 at α = 0.1 and p = 0.05. In Figure 9, we plot the minimum total cost TC* vs. c 1 for p = 0.05, 0.075, 0.1 at α = 0.1 and β = 0.25. From Figure 7, we find that the minimum total cost TC** decreases when the value of α increases. From Figure 8, we find that the minimum total cost TC** increases when the test power 1 − β increases. From Figure 9, we find that the minimum total cost TC** becomes larger when the removal probability becomes larger. Overall, the minimum total cost TC** becomes smaller when the value of c 1 grows larger.

3.3. Example

We use the data in Caroni [15] consisting of the failure times (number of cycles in 1000 times) of n = 25 ball bearings in an automatic life test, listed as follows:
0.1788 0.2892 0.3300 0.4152 0.4212 0.4560 0.4848 0.5184 0.5196 0.5412
 0.5556 0.6780 0.6780 0.6780 0.6864 0.6864 0.6888 0.8412 0.9312 0.9864
 1.0512 1.0584 1.2792 1.2804 1.7340                 
In Figure 10, we plot the empirical cumulative distribution function (ecdf) for this data set. The G test based on the Gini statistic (see Gail and Gastwirth [16]) is computed as 0.5052237 with p value 0.9293607 > 0.05. Therefore, we conclude that this data set fits the Rayleigh distribution well. We will now use this example to illustrate the implementation of the previous subsections.
For Section 3.1, we consider the case of α = 0.01, the power 1 – β = 0.85 at c 1 = 0.95 under fixed T = 1.0. From Table 1, we can find that the optimal sampling design is m* = 2, n* = 12, yielding the minimum total cost of 20 units under the cost setup of Co = 1 and a = b = c = 1. The critical value C L 0 = 0.921702 is also found in this table.
We can thus start to perform the testing procedure for H 0 : C L 0.85 as follows:
Step 1:
Take a random sample of size n = 12 from the data set. Observe the progressive type I interval censored sample (X1,X2) = (3,4) at the pre-set times (t1,t2) = (0.5,1.0) with censoring schemes of (R1,R2) = (2,3).
Step 2:
Obtain the MLE of λ as λ ^ = 0.6625991 and then obtain the test statistic C ^ L = 1 0.05 2 ( 0 . 6625991 ) 2 = 0.9430573.
Step 3:
Compared with the critical value, we have C ^ L = 0.9557158 > C L 0 = 0.9042. Thus, we can conclude that the lifetime performance index of the product meets the required level of 0.85.
For Section 3.2, we consider the case of α = 0.05, the power 1 – β = 0.85 at c 1 = 0.925. From Table 4, we can find that the optimal sampling design is m* = 2, n* = 14 and t* = 0.42, yielding the minimum total cost of 17.843 units under the cost setup of Co = 1 and a = b = c = 1. The corresponding critical value C L 0 = 0.9042 is also found in this table.
We can thus start to perform the testing procedure for H 0 : C L 0.85 as follows:
Step 1:
Take a random sample of size n = 14 from the data set. Observe the progressive type I interval censored sample (X1,X2) = (1,5) at the pre-set times (t1,t2) = (0.42,0.84) with censoring schemes of (R1,R2) = (2,6).
Step 2:
Obtain the MLE of λ as λ ^ = 0.7513559 and then obtain the test statistic C ^ L = 1 0 . 05 2 ( 0 . 7513559 ) 2 = 0.9557158.
Step 3:
Compared with the critical value, we have C ^ L = 0.9557158 > C L 0 = 0.9042. Therefore, we reach the same conclusion and reject the null hypothesis.

4. Conclusions

Process capability indices are widely applied by manufacturers to evaluate the capability performance of a specific process when the lifetime of a product follows Rayleigh distribution. For many practical cases, a progressive type I interval censored sample is collected instead of a complete sample. We investigate the required minimum number of inspection intervals when the termination time of the experiment is fixed to reach given power and the minimum total cost for a level α test. When the termination time of the experiment is not fixed, the required minimum sample size, number of inspection intervals and the inspection interval time length are determined in this paper to reach given power and the minimum total cost for a level α test under progressive type I interval censoring. In the future, we plan to propose a testing procedure and study the experimental design for other lifetime distributions such as Lomax distribution. Other than the procreative type I interval sampling scheme, our studies can be applied to type II progressively hybrid censored samples.

Author Contributions

Conceptualization, S.-F.W.; methodology, S.-F.W.; software, S.-F.W., T.-H.L., Y.-H.L. and W.-T.C.; validation, T.-H.L., Y.-H.L. and W.-T.C.; formal analysis, S.-F.W.; investigation, S.-F.W., T.-H.L., Y.-H.L. and W.-T.C.; resources, S.-F.W.; data curation, S.-F.W., T.-H.L., Y.-H.L. and W.-T.C.; writing—original draft preparation, S.-F.W. and W.-T.C.; writing—review and editing, S.-F.W.; visualization, T.-H.L., Y.-H.L. and W.-T.C.; 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 was funded by the Ministry of Science and Technology, Taiwan under MOST 108-2118-M-032-001 and MOST 109-2118-M-032-001-MY2 and the APC was funded by MOST 109-2118-M-032-001-MY2.

Data Availability Statement

Data are available in a publicly accessible repository The data presented in this study are openly available in Caroni [15].

Conflicts of Interest

The authors declare no conflict of interest.

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] [CrossRef]
  4. 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]
  5. 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]
  6. Lee, W.C.; Wu, J.W.; Hong, C.W.; Pan, H.Y.; Hung, W.L. Decision procedure of performance assessment of lifetime index of products for the Gompertz distribution. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2010, 224, 493–499. [Google Scholar] [CrossRef]
  7. Lee, W.C.; Wu, J.W.; Hong, M.L.; Lin, L.S.; Chan, R.L. Assessing the lifetime performance index of Rayleigh products based on the Bayesian estimation under progressive type II right censored samples. J. Comput. Appl. Math. 2011, 235, 1676–1688. [Google Scholar] [CrossRef]
  8. Wu, J.W.; Lee, W.C.; Hong, C.W.; Yeh, S.Y. Implementing Lifetime Performance Index of Burr XII Products with Progressively Type II Right Censored Sample. Int. J. Innov. Comput. Inf. Control. 2014, 10, 671–693. [Google Scholar]
  9. Wu, S.F.; Hsieh, Y.T. The assessment on the lifetime performance index of products with Gompertz distribution based on the progressive type I interval censored sample. J. Comput. Appl. Math. 2019, 351, 66–76. [Google Scholar] [CrossRef]
  10. Wu, S.F.; Xie, Y.J.; Liao, M.F.; Chang, W.T. Reliability sampling design for the lifetime performance index of Gompertz lifetime distribution under progressive type I interval censoring. Mathematics 2021, 9, 2109. [Google Scholar] [CrossRef]
  11. Wu, S.F.; Lin, M.J. Computational testing algorithmic procedure of assessment for lifetime performance index of products with weibull distribution under progressive type I interval censoring. J. Comput. Appl. Math. 2017, 311, 364–374. [Google Scholar] [CrossRef]
  12. Wu, S.F.; Wu, Y.C.; Wu, C.H.; Chang, W.T. Sampling design for the lifetime performance index of Weibull lifetime distribution under progressive type I interval censoring. Symmetry 2021, 13, 1691. [Google Scholar]
  13. Wu, S.F.; Lin, Y.T.; Chang, W.J.; Chang, C.W.; Lin, C. A computational algorithm for the evaluation on the lifetime performance index of products with Rayleigh distribution under progressive type I interval censoring. J. Comput. Appl. Math. 2018, 328, 508–519. [Google Scholar] [CrossRef]
  14. 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]
  15. Caroni, C. The correct “ball bearings” data. Lifetime Data Anal. 2002, 8, 395–399. [Google Scholar] [CrossRef] [PubMed]
  16. 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. (a) The curve of pdf for Rayleigh distribution; (b) the curve of failure rate function for Rayleigh distribution.
Figure 1. (a) The curve of pdf for Rayleigh distribution; (b) the curve of failure rate function for Rayleigh distribution.
Mathematics 10 00517 g001
Figure 2. (a): Total cost curve at m = 1:m0 for the first case; (b) total cost curve at m = 1:m0. for the second case.
Figure 2. (a): Total cost curve at m = 1:m0 for the first case; (b) total cost curve at m = 1:m0. for the second case.
Mathematics 10 00517 g002
Figure 3. Minimum total cost curve for α = 0.01, 0.05, 0.1 at β = 0.25 and p = 0.05.
Figure 3. Minimum total cost curve for α = 0.01, 0.05, 0.1 at β = 0.25 and p = 0.05.
Mathematics 10 00517 g003
Figure 4. Minimum total cost curve for 1 − β = 0.75, 0.80, 0.85 at α = 0.1 and p = 0.05.
Figure 4. Minimum total cost curve for 1 − β = 0.75, 0.80, 0.85 at α = 0.1 and p = 0.05.
Mathematics 10 00517 g004
Figure 5. Minimum total cost curve for p = 0.05, 0.075, 0.1 at α = 0.1 and β = 0.25.
Figure 5. Minimum total cost curve for p = 0.05, 0.075, 0.1 at α = 0.1 and β = 0.25.
Mathematics 10 00517 g005
Figure 6. (a) Total cost curve at m = 1:m0 for the first case; (b) Total cost curve at m = 1:m0. for the second case.
Figure 6. (a) Total cost curve at m = 1:m0 for the first case; (b) Total cost curve at m = 1:m0. for the second case.
Mathematics 10 00517 g006
Figure 7. Minimum total cost curve for α = 0.01, 0.05, 0.1 at β = 0.25 and p = 0.05.
Figure 7. Minimum total cost curve for α = 0.01, 0.05, 0.1 at β = 0.25 and p = 0.05.
Mathematics 10 00517 g007
Figure 8. Minimum total cost curve for 1 − β = 0.75, 0.80, 0.85 at α = 0.1 and p = 0.05.
Figure 8. Minimum total cost curve for 1 − β = 0.75, 0.80, 0.85 at α = 0.1 and p = 0.05.
Mathematics 10 00517 g008
Figure 9. Minimum total cost curve for p = 0.05, 0.075, 0.1 at α = 0.1 and β = 0.25.
Figure 9. Minimum total cost curve for p = 0.05, 0.075, 0.1 at α = 0.1 and β = 0.25.
Mathematics 10 00517 g009
Figure 10. The ecdf for the data set in Caroni [15].
Figure 10. The ecdf for the data set in Caroni [15].
Mathematics 10 00517 g010
Table 1. The optimal (m*,n*), total cost TC and critical value for c1 = 0.875, 0.90, α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15 and p = 0.05, 0.075, 0.1 under m0 = 20, L = 0.05 and c0 = 0.85.
Table 1. The optimal (m*,n*), total cost TC and critical value for c1 = 0.875, 0.90, α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15 and p = 0.05, 0.075, 0.1 under m0 = 20, L = 0.05 and c0 = 0.85.
c10.8750.90
αβpm* m * n* n * T C T C C L 0   C L 0 m* m * n* n * TC T C C L 0   C L 0
0.010.250.05043093150.869497361660.889173
0.07533193240.869534362670.889269
0.10033273320.869515364690.889127
0.200.05043473530.868502369740.887289
0.07533593640.868519371760.887197
0.10033683730.868502372770.887327
0.150.05043954010.867445379840.885299
0.07534084130.867474381860.885262
0.10034184230.867462383880.885237
0.050.250.05041871930.867883338430.885858
0.07531941990.867877240440.886332
0.10031982030.867895241450.886105
0.200.05042172230.86671344490.883769
0.07532252300.866709345500.883769
0.10032302350.866714346510.883778
0.150.05042552610.865517353580.881256
0.07532642690.865527354590.881305
0.10032702750.865528355600.881362
0.100.250.05031361410.866571229330.883662
0.07531391440.866575328330.883426
0.10031431480.866532230340.883447
0.200.05031621670.865291333380.880926
0.07531661710.865277235390.881249
0.10031701750.865269236400.881003
0.150.05041931990.864012341460.878209
0.07532002050.864015243470.878647
0.10032052100.864002244480.878479
Table 2. The optimal (m*,n*), total cost TC and critical value for c1 = 0.925, 0.95, α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15 and p = 0.05, 0.075, 0.1 under m0 = 20, L = 0.05 and c0 = 0.85.
Table 2. The optimal (m*,n*), total cost TC and critical value for c1 = 0.925, 0.95, α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15 and p = 0.05, 0.075, 0.1 under m0 = 20, L = 0.05 and c0 = 0.85.
c10.9250.95
αβpm* n* TC C L 0 m* n* TC C L 0
0.010.250.050221250.90972529130.927812
0.075222260.90897229130.928015
0.100222260.90915829130.92822
0.200.050323280.906423210140.925489
0.075225290.906429210140.92569
0.100225290.906611210140.925894
0.150.050327320.903317212160.921503
0.075229330.90355212160.921702
0.100229330.903728212160.921903
0.050.250.050214180.90404426100.921499
0.075214180.90422126100.921699
0.100214180.90439926100.9219
0.200.050216200.90151427110.918176
0.075216200.90168627110.918373
0.100217210.90072927110.918571
0.150.050220240.89744729130.912881
0.075220240.8976129130.913071
0.100220240.89777529130.913263
0.100.250.050210140.9009711690.921343
0.075210140.9011411690.921343
0.100211150.8995561690.921343
0.200.050213170.89624317100.918022
0.075213170.89640417100.918022
0.100213170.89656617100.918022
0.150.050215190.89378927110.907841
0.075216200.89286327110.908024
0.100216200.89301727110.908208
Table 3. The optimal (m*,t*,n*), total cost TC and critical value for c1 = 0.875, 0.90, α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15 and p = 0.05, 0.075, 0.1 under m0 = 20, L = 0.05 and c0 = 0.85.
Table 3. The optimal (m*,t*,n*), total cost TC and critical value for c1 = 0.875, 0.90, α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15 and p = 0.05, 0.075, 0.1 under m0 = 20, L = 0.05 and c0 = 0.85.
c10.8750.90
αβpm* t*n*TC** C L 0 m* t*n*TC** C L 0
0.010.150.05050.26297304.2850.869640.285864.1270.8895
0.07550.27307314.3520.869630.356166.0520.8894
0.10040.31317323.2560.869630.366267.0930.8895
0.200.05060.24332340.4640.868640.286672.1180.8875
0.07550.27345352.3430.868640.326773.2930.8876
0.10040.32356362.2630.868630.397075.1630.8876
0.250.05060.25377385.4840.867540.37581.2190.8856
0.07550.28391398.4160.867540.327783.2710.8856
0.10040.32404410.2880.867540.337985.3170.8856
0.050.150.05050.25180187.250.867930.343742.0180.8862
0.07540.3187193.20.867930.343843.0060.8861
0.10040.31192198.2390.867930.383843.1420.8864
0.200.05050.25209216.2360.866830.364348.0840.8840
0.07550.28215222.4150.866830.364449.0790.8840
0.10040.31223229.2230.866830.364550.080.8840
0.250.05050.25245252.2540.865630.375156.0980.8817
0.07540.31254260.2440.865630.375257.120.8817
0.10040.32261267.2640.865630.385358.1430.8817
0.100.150.05050.26129136.3140.866630.342732.0190.8836
0.07540.32134140.2820.866630.332832.9780.8834
0.10040.31138144.240.866630.362833.0670.8836
0.200.05060.24153161.4210.865330.383237.1390.8813
0.07540.3160166.1950.865330.363338.0690.8811
0.10040.31164170.260.865330.353439.0420.8810
0.250.05050.25185192.2750.864130.383944.1320.8788
0.07540.31192198.2450.864130.374045.1040.8787
0.10040.32197203.2980.864130.364146.0930.8786
Table 4. The optimal (m*,t*,n*), total cost TC and critical value for c1 = 0.925, 0.95, α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15 and p = 0.05, 0.075, 0.1 under m0 = 20, L = 0.05 and c0 = 0.85.
Table 4. The optimal (m*,t*,n*), total cost TC and critical value for c1 = 0.925, 0.95, α = 0.01, 0.05, 0.1, β = 0.25, 0.20, 0.15 and p = 0.05, 0.075, 0.1 under m0 = 20, L = 0.05 and c0 = 0.85.
c10.9250.95
αβpm* t*n*TC** C L 0 m* t*n*TC** C L 0
0.010.150.05020.482124.9580.909720.52812.0470.9306
0.07530.352025.0430.909520.55812.1080.9310
0.10030.392025.1660.909920.61812.220.9317
0.200.05030.332327.9880.906420.461013.9190.9254
0.07530.362328.0650.906720.471013.9430.9256
0.10030.42328.1920.907220.481013.970.9259
0.250.05030.372631.1110.904020.611115.2130.9245
0.07530.352732.0540.903620.481215.9540.9216
0.10030.382732.1510.904020.491215.9780.9219
0.050.150.05020.421417.8430.904210.7479.7360.9220
0.07520.431417.8660.904410.7479.7360.9220
0.10020.451417.8920.904410.7479.7360.9220
0.200.05020.471619.9470.901420.45710.9060.9181
0.07520.491619.9870.901720.46710.9260.9183
0.10020.521620.0420.901920.47710.9470.9185
0.250.05020.511923.0140.898420.53812.0670.9155
0.07520.541923.0830.898820.55812.10.9159
0.10030.421823.2560.898720.57812.1370.9162
0.100.150.05020.471013.9350.900920.5448.0850.9198
0.07520.491013.9710.901120.5648.1260.9201
0.10020.511014.0170.901320.5948.1730.9205
0.200.05020.511216.0170.897720.5659.1130.9152
0.07520.541216.0840.898020.5859.1520.9156
0.10020.441316.8850.896720.659.1970.9159
0.250.05020.491518.990.893820.47710.9380.9078
0.07520.521519.0320.894020.48710.9570.9080
0.10020.551519.090.894320.49710.9770.9082
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://0-doi-org.brum.beds.ac.uk/10.3390/math10030517

AMA Style

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(3):517. https://0-doi-org.brum.beds.ac.uk/10.3390/math10030517

Chicago/Turabian Style

Wu, Shu-Fei, Tzu-Hsuan Liu, Yu-Hua Lai, and Wei-Tsung Chang. 2022. "A Study on the Experimental Design for the Lifetime Performance Index of Rayleigh Lifetime Distribution under Progressive Type I Interval Censoring" Mathematics 10, no. 3: 517. https://0-doi-org.brum.beds.ac.uk/10.3390/math10030517

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