Next Article in Journal
Towards Interval Type-3 Intuitionistic Fuzzy Sets and Systems
Next Article in Special Issue
Modeling Income Data via New Parametric Quantile Regressions: Formulation, Computational Statistics, and Application
Previous Article in Journal
Influenza-like Illness Detection from Arabic Facebook Posts Based on Sentiment Analysis and 1D Convolutional Neural Network
Previous Article in Special Issue
Estimation of the Six Sigma Quality Index
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Process Performance Index for the Weibull Distribution with a Type-I Hybrid Censoring Scheme

1
Department of Statistics, Tamkang University, Tamsui District, New Taipei City 251301, Taiwan
2
Department of Mathematical Sciences, University of South Dakota, Vermillion, SD 57069, USA
3
School of Statistics, Southwestern University of Finance and Economics, Chengdu 611130, China
*
Authors to whom correspondence should be addressed.
Mathematics 2022, 10(21), 4090; https://doi.org/10.3390/math10214090
Submission received: 13 September 2022 / Revised: 28 October 2022 / Accepted: 30 October 2022 / Published: 2 November 2022

Abstract

:
A new life performance index is proposed for evaluating the quality of lifetime products. The maximum likelihood estimation method and the Bayesian approaches using informative and non-informative prior distributions are utilized to infer the parameters of the Weibull distribution and the proposed new life performance index under a Type-I hybrid censoring scheme. Monte Carlo simulation results show that two Bayesian approaches outperform the maximum likelihood estimation method in terms of the measures of relative bias, relative mean square error, and coverage probability for the point and confidence interval estimators, respectively. The Bayesian approach using a non-informative prior distribution is recommended if the knowledge of setting up the hyper-parameters in the informative prior distribution is not available. Two real data sets are provided for illustration.

1. Introduction

Process capability analysis is one of the important issues during the statistical quality control applications. It allows the producer to know whether the product meets the specification and quality requirements. Process capability index (PCI) is a method to quantify process capability. The concept of PCI is to measure the degree of the quality characteristic meeting the requirements of a given specification limit. At present, many commonly used PCIs have been suggested in the literature to improve product quality, such as C p , C p k , C p m , and C p m k (Juran [1], Kane [2], Hsiang and Taguchi [3], Pearn et al. [4] and Montgomery [5]), which are, respectively, defined as follows,
C p = U L 6 σ ,
C p k = min U μ 3 σ , μ L 3 σ ,
C p m = U L 6 σ 2 + ( μ T ) 2 ,
C p m k = min U μ 3 σ 2 + ( μ T ) 2 , μ L 3 σ 2 + ( μ T ) 2 ,
where U, L, T, σ and μ are the upper, lower specification limits, target value, and process standard deviation and mean, respectively. The aforementioned PCIs have two specification limits to quantify the product quality. Practically, producers and customers could only consider the quality of products to meet above the lower specification limit. Kane [2] proposed a PCI, named C p l , which only considers the lower limit of the specification and is defined by
C p l = μ L 3 σ .
The process performance index based on the version of C p l is expressed by
C L = μ L σ .
Readers can refer to the work of Montgomery [5] for more comprehensive discussions. It is worth noting that the premise of using these traditional PCIs is the normality assumption. If lifetime data do not follow a normal distribution, reliable results could not be achieved. For example, measurement data for drilling processes, coating processes, and chemical processes often follow an asymmetric distribution. For improving the quality of traditional PCIs for asymmetrically distributed data, Clements [6] proposed to replace the C p and C p k by C p ( q ) and C p k ( q ) , respectively. The two new PCIs for asymmetric distributions are, respectively, defined by
C p ( q ) = U L U 0.99865 L 0.00135
and
C p k ( q ) = min U M U 0.99865 M , M L M L 0.00135 ,
where L 0.00135 , U 0.99865 , and M are, respectively, the 0.00135 -quantile, 0.99865 -quantile, and median of quality characteristic measurements. Similarly, C L should be modified for asymmetric distributions. A new modification will be proposed in the study for the Weibull distribution.
The manufacturing technology advancing has successfully prolonged the lifetime of products. In order to save the time and cost of life testing, Type-I and Type-II censoring schemes have been popular for implementing life tests. Type-I censoring uses a time censoring scheme, and the life test is terminated at the predetermined time T, then the failure times less than or equal to T are recorded as a type-I censoring sample. Type-II censoring scheme is also called failed number censoring scheme. r failed lifetimes are collected as a Type-II censored sample in which the positive integer, r, is predetermined by engineers before the life test. To make the censoring scheme more flexible, Epstein [7] introduced a Type-I hybrid censoring scheme. The life test can be terminated at the minimal termination time of the Type-I and Type-II censoring schemes. When the life test is terminated at the maximal termination time of the Type-I and Type-II censoring schemes, the censoring scheme is named Type-II hybrid censoring scheme, see Childs et al. [8].
At present, many researchers have made contributions to the life performance index for non-normal distribution. For example, the life performance index based on random sample from Rayleigh distribution by Lee [9], based on type II censored sample from Pareto distribution by Hong et al. [10], based on progressively type II censored sample from exponential distribution by Lee et al. [11], and based on progressively type I interval censored sample from the Weibull distribution by Wu and Lin [12]. Although these studies have insight contributions into the life performance index for asymmetric distributions, it still needs to be improved for reflecting the Weibull distribution data with the traditional life performance index. Hence, a new life performance index is proposed to evaluate the life quality of the Weibull distribution data. Based on our knowledge, no one has studied the inference method for the new proposed life performance index using a Type-I hybrid censoring scheme. In this study, we proposed the maximum likelihood estimation and Bayesian inference methods to obtain the point estimates of the parameters and the new proposed life performance index based on Type-I hybrid censored sample from Weibull distribution. Because the Type-I hybrid censoring scheme is used for life testing to save testing time, the Fisher information matrix is difficult to obtain because the number of failures and stopping time of the life testing are random variables. The facts of the number of failures and stopping time of the life testing are random variables that also have an impact on the quality of the observed Fisher information matrix. The second derivatives of the log-likelihood function based on Type-I hybrid censored samples from the Weibull distribution are complicated, and the Fisher information matrix exists only if some regular conditions are true, see Wang and He [13]. Hence, we suggest using a bootstrap method to find a confidence interval based on the maximum likelihood estimation method. The highest posterior density interval (HPDI) is used to develop a Bayesian credible interval for the new proposed life performance index. Some studies have also investigated the parameter estimation methods based on Type-I hybrid censored samples. For example, Kundu and Pradhan [14] for the generalized exponential distribution paramters; Lin et al. [15] for the Weibull distribution with a progressive hybrid censoring scheme; Cho et al. [16] for the estimate of Weibull distribution entropy with a generalized progressive hybrid censoring scheme; and Okasha and Mustafa [17] for the E-Bayesian estimation of rate parameter based on adaptive Type-I progressive hybrid censored sample from the Weibull distribution assuming the shape parameter as a constant. However, these four works do not deal with the life performance index.
The rest of this paper is organized as follows. In Section 2, a new life performance index is proposed and its properties are discussed. Section 3 introduces the parameter estimation methods based on a Type-I hybrid censored sample from the Weibull distribution by utilizing the maximum likelihood estimation and Bayesian inference methods. Section 4 reports the results of an intensive simulation study for the evaluation of two proposed estimation methods. Two real examples are illustrated in Section 5. Finally, some concluding remarks will be made in Section 6.

2. Proposed Life Performance Index

Let X be an asymmetrically distributed random variable that has a finite second moment. In this case, the mean, μ = E [ X ] , in Equation (5) is less suitable to be the center representative for the distribution. A new life performance index is proposed and expressed by
C L M = M L E ( X M ) 2 = δ 1 + δ 2 + C L 1 + δ 2
where M is the median of X and
δ = M μ σ .
The relationship between the proposed C L M and the C L is stated in Theorem 1.
 Theorem 1. 
For any distribution that has the finite second moment,
 (a) 
If M L μ L < 1 + δ 2 , then C L M C L < 0 .
 (b) 
If M L μ L > 1 + δ 2 , then C L M C L > 0 .
 (c) 
If M L μ L = 1 + δ 2 , then C L M = C L .
The proof of Theorem 1 is addressed in Appendix A. Because the condition of
M L μ L = 1 + δ 2
is difficult to hold for an asymmetric distribution, C L M cannot easily be equal to C L . This fact indicates that when the lifetimes follow an asymmetric distribution, C L M can better characterize the quality of lifetime products.
When X has the Weibull distribution with the probability density, cumulative distribution, survival, and quantile functions that are, respectively, expressed by
f ( x ; Θ ) = α λ x λ α 1 exp x λ α , x > 0 ,
F ( x ; Θ ) = 1 exp x λ α , x > 0 ,
S ( x ; Θ ) = 1 F ( x ; Θ ) = exp x λ α , x > 0 ,
and
q ( u ; Θ ) = λ ( log ( 1 u ) ) 1 α , 0 < u < 1 ,
where Θ = ( α , λ ) , α > 0 is the shape parameter and λ > 0 is scale parameter. With the wide range of shape parameter, α > 0 , the Weibull distribution is a flexible and skewed one for lifetime modeling. The Weibull distribution has another form, which can be obtained via using reparameterization of τ = λ α . The probability density function of the Weibull distribution with parameters α and τ was represented by Kundu [18] to be
f ( x ; α , τ ) = α τ x α 1 exp { τ x α } , x > 0 .
It can be shown that the mean, μ , standard deviation, σ , and median, M, of the Weibull distribution are
μ = λ Γ 1 + 1 α ,
σ = λ A α ,
and
M = λ ( log ( 2 ) ) 1 α ,
where
A α = Γ 1 + 2 α Γ 1 + 1 α 2 .
Hence, Equation (5) can be represented as
C L = μ L σ = Γ 1 + 1 α L λ A α
that implies the life performance index C L M of Equation (8) can also be rewritten by
C L M = δ 1 + δ 2 + Γ 1 + 1 α L λ A α 1 + δ 2 ,
and
L λ = Γ 1 + 1 α A α × C L .
Using Equation (13), the product yield can be expressed by
P r = P ( X > L ) = exp L λ α = exp Γ 1 + 1 α A α C L α = exp Γ 1 + 1 α A α C L M 1 + δ 2 δ α .

3. Parameter Estimation Methods Based on Hybrid Censored Data

Let { X i , i = 1 , 2 , , n } be independent random lifetimes from the Weibull distribution with the parameter ( Θ ) . Denote { X ( 1 ) X ( 2 ) X ( n ) } as the order statistics of { X i , i = 1 , 2 , , n } . The life test stops immediately when the rth ordered failed time is observed or the predetermined T is up. The realization of the failure time sample can be presented by { x ( 1 ) x ( 2 ) x ( d ) } , where
d = r , x ( r ) T , m , x ( r ) > T ,
and m is the number of failure times by T. Let
c = x ( r ) , x ( r ) T , T , x ( r ) > T .
Then, the Type-I hybrid censored sample can be expressed by D = { x ( 1 ) , x ( 2 ) , , x ( d ) , c } , which contains n d repetitions of c. We only use one c in D for simplification.

3.1. Maximum Likelihood Estimation

According to the expression of Kundu [19], the log-likelihood function can be expressed by
( Θ ) log L Θ | D = d log ( α ) d log ( λ ) + ( α 1 ) i = 1 d log ( x ( i ) ) d ( α 1 ) log ( λ ) i = 1 d x ( i ) λ α ( n d ) × c λ α .
The partial derivatives of ( Θ ) with respect to λ and α can be derived as follows,
( Θ ) λ = α λ α + 1 i = 1 d x ( i ) α + ( n d ) c α d α λ
and
( Θ ) α = d α d log ( λ ) + i = 1 d log ( x ( i ) ) i = 1 d x ( i ) λ α log ( x ( i ) ) log ( λ ) ( n d ) c λ α ( log ( c ) log ( λ ) ) .
Let
( Θ ) λ = 0
and
( Θ ) α = 0
Under the condition of d > 0 , the maximum likelihood estimator (MLE), Θ ^ = ( α ^ , λ ^ ) , of Θ can be obtained by Equations (20) and (21),
λ ^ = 1 d i = 1 d x ( i ) α ^ + ( n d ) c α ^ 1 / α ^
and
α ^ = d 1 u ( α ^ ) i = 1 d x ( i ) α ^ log ( x ( i ) ) + ( n d ) c α ^ log ( c ) i = 1 d log ( x ( i ) ) ,
where u ( α ^ ) = λ ^ α ^ and
u ( α ) = 1 d i = 1 d x ( i ) d + ( n d ) c α .
Since Equations (20) and (21) indicate the MLEs have no closed form, Θ ^ can only be obtained by using numerical method through R function optim to Equation (17). Replacing Θ with Θ ^ in Equation (8), the MLE of C L M is given by
C ^ L M = δ ^ + C ^ L 1 + δ ^ 2 ,
where
C ^ L = Γ 1 + 1 α ^ L λ ^ A ^ α ,
δ ^ = ( log ( 2 ) ) 1 α Γ 1 + 1 α ^ A ^ α ,
and
A ^ α = Γ 1 + 2 α ^ Γ 1 + 1 α ^ 2 .
The parametric bootstrap method with the following steps is suggested to obtain an approximate confidence interval of C L M :
 Step 1: 
Obtain the MLE of Θ ^ = α ^ , λ ^ via using Equations (20) and (21) based on the Type-I hybrid censored data D = { x ( 1 ) , x ( 2 ) , , x ( d ) , c } . Then, obtain the MLE of C L M by using Equation (22).
 Step 2: 
Generate a new Type-I hybrid censored data D * = { x 1 : n * , x 2 : n * , , x d * : n * , c * } from the Weibull distribution with parameter, ( Θ ^ ) , where Θ ^ is MLE of Θ from Step 1.
 Step 3: 
Compute MLEs of Θ and C L M based on Type-I hybrid censored data D * and denote the obtained MLEs by Θ ^ * and C ^ LM * , respectively.
 Step 4: 
Repeat Step 2 to Step 3 B * times to obtain the bootstrap sample { C ^ L M , j * , j = 1 , 2 , , B * } . Denote the empirical distribution based on the obtained bootstrap sample, { C ^ L M , 1 * , C ^ L M , 2 * , , C ^ L M , B * * } , by G ^ * .
 Step 5: 
Given a significance level γ , find the ( γ / 2 ) th and ( 1 γ / 2 ) th empirical quantiles of G ^ * as the lower and upper limits of the confidence interval of C L M , respectively.

3.2. Bayesian Estimation

Let the prior distribution of λ be the inverse Gamma, denoted by π 1 ( λ ) , and the prior distribution of α be a non-informative prior distribution, denoted by π 2 ( α ) . Thus, the joint prior distribution of Θ can be expressed by
π ( Θ ) = π 1 ( λ ) × π 2 ( α ) , α > 0 , λ > 0 ,
where
π 1 ( λ ) = b a Γ ( a ) 1 λ a + 1 exp b λ , λ > 0 , a , b > 0
and
π 2 ( α ) 1 α , α > 0 .
Then, the joint posterior distribution of Θ can be expressed by
π Θ | D π ( Θ ) L Θ | D α λ d i = 1 d x ( i ) λ α 1 exp i = 1 d x ( i ) λ α + ( n d ) × c λ α × b a Γ ( a ) 1 λ a + 1 exp b λ 1 α α d 1 λ d α a 1 i = 1 d x ( i ) α κ D , α , λ ,
where
κ D , α , λ = exp i = 1 d x ( i ) λ α + ( n d ) × c λ α + b λ .
Using Equation (29), we can obtain the conditional posterior distributions of λ given α and of α given λ , which are presented by
π λ | D , α λ d α a 1 κ D , α , λ
and
π α | D , λ α d 1 λ d α i = 1 d x ( i ) α κ D , α , λ ,
respectively. Notably, if we lack knowledge to set up the prior of λ , a non-informative prior distribution π 1 ( λ ) 1 λ can be used for λ . Then, the joint posterior distribution of Θ can be given by
π N Θ | D α λ d i = 1 d x ( i ) λ α 1 exp i = 1 d x ( i ) λ α + ( n d ) × c λ α × 1 α λ α d 1 λ d α 1 i = 1 d x ( i ) α κ 1 D , α , λ ,
where
κ 1 D , α , λ = exp i = 1 d x ( i ) λ α + ( n d ) × c λ α .
Using Equation (32), we can derive the conditional post distributions of λ given α and α given λ by
π N ( λ | D , α ) λ d α 1 κ 1 D , α , λ
and
π N α | D , λ α d 1 λ d α i = 1 d x ( i ) α κ 1 D , α , λ ,
respectively. Because the closed form of the marginal posterior distributions π α | D and π λ | D (or π N α | D and π N λ | D ) cannot be obtained, the Markov chain Monte Carlo (MCMC) approach with mixing Gibbs sampling technique and Metropolis–Hastings algorithm is used to obtain the Bayesian estimators. The implementation steps of the MCMC approach are presented as follows:
 Step 1: 
Give initial values ( λ , α ) = ( λ ( 0 ) , α ( 0 ) ) and propose the transition probability distributions, q 1 λ ( * ) | λ ( t ) and q 1 α ( * ) | α ( t ) , where λ ( * ) and α ( * ) are the updates of λ and α for the next step.
 Step 2: 
Implement Step 3 for i = 0 , 1 , 2 , , B , where B is a huge number.
 Step 3: 
 (a) 
Generate λ ( * ) q 1 · | λ ( i ) and u U ( 0 , 1 ) , where U ( 0 , 1 ) is the uniform distribution over the interval ( 0 , 1 ) . Update λ ( i + 1 ) according to the condition,
λ ( i + 1 ) = λ ( * ) if u min 1 , π λ ( * ) | α ( t 1 ) ; D q 1 λ ( t 1 ) | λ ( * ) π λ ( t 1 ) | α ( t 1 ) ; D q 1 λ ( * ) | λ ( t 1 ) , λ ( i ) , o t h e r w i s e .
 (b) 
Generate α ( * ) q 2 ( · | α ( i ) ) and u U ( 0 , 1 ) . Update α ( i + 1 ) according to the condition:
α ( i + 1 ) = α ( * ) if u min 1 , π α ( * ) | λ ( t ) ; D q 2 α ( t 1 ) | α ( * ) π α ( t 1 ) | λ ( t ) ; D q 2 α ( * ) | α ( t 1 ) , α ( i ) , o t h e r w i s e .
 (c) 
Replacing ( α , λ ) with ( α ( i + 1 ) , λ ( i + 1 ) ) in Equation (8), the updated estimator of C L M can be obtain by
C L M ( i + 1 ) = δ ( i + 1 ) + C L ( i + 1 ) 1 + δ ( i + 1 ) 2 ,
where
C L ( i + 1 ) = Γ 1 + 1 α ( i + 1 ) L λ ( i + 1 ) A α ( i + 1 ) , δ ( i + 1 ) = ( log ( 2 ) ) 1 α ( i + 1 ) Γ 1 + 1 α ( i + 1 ) A α ( i + 1 ) ,
and
A α ( i + 1 ) = Γ 1 + 2 α ( i + 1 ) Γ 1 + 1 α ( i + 1 ) 2 .
 Step 4: 
Remove the first N 1 Markov chains for the burn-in operation. Re-coding the Markov chains of { α ( i ) , i = N 1 + 1 , N 1 + 2 , . . . , B } = { α ( i ) , i = 1 , 2 , , B 1 } , { λ ( i ) , i = N 1 + 1 , N 1 + 2 , , B } = { λ ( i ) , i = 1 , 2 , , B 1 } , and { C L M ( i ) , i = N 1 + 1 , N 1 + 2 , , B } = { C L M ( i ) , i = 1 , 2 , , B 1 } . Considering the squared loss function for Bayesian estimation, the Bayes estimates α ˜ , λ ˜ and C ˜ L M can be, respectively, obtained by the sample means of the Markov chains of α ( i ) , λ ( i ) and C LM ( i ) , i = 1 , 2 , , B 1 .
The HPDI of C L M can be obtained via using the following steps proposed by Chen and Shao [20].
 Step 1: 
Sorting { C L M ( j ) , j = 1 , 2 , , B 1 } to obtain an ordered sequence of C L M ( j ) by { C L M ( j ) , j = 1 , 2 , , B 1 } , where C L M ( j ) C L M ( j + 1 ) for j = 1 , 2 , ( B 1 1 ) . The ordered sequence of { C L M ( 1 ) , C L M ( 2 ) , , C L M ( B 1 ) } is used to construct the empirical distribution of the Bayes estimator C ˜ L M .
 Step 2: 
Find all 100 × ( 1 γ ) % intervals of C L M that are labeled by ( C L M ( j ) , C L M ( j + ( 1 γ ) B 1 ) ) , j = 1 , 2 , , B 1 ( 1 γ ) B 1 , where [ y ] is the largest integer smaller or equal to y.
 Step 3: 
Find the interval that has the shortest length among all intervals in Step 2. The obtained interval is the HPDI.
Similarly, if we consider that λ follows a non-informative prior distribution, the Bayes estimators of α , λ and C L M can be obtained through of using the above proposed MCMC approach with replacing π · | · ; D by π N · | · ; D . The HPDI of C L M can also be obtained by using the aforementioned steps proposed by Chen and Shao [20].

4. Simulation Study

Monte Carlo simulations were carried out to examine the performance of the MLEs and Bayes estimators of the model parameters and the proposed life performance index. In practical applications, one often evaluates the quality of lifetime products under a limited sample resource. It is important to evaluate the quality of the proposed maximum likelihood estimation method and Bayesian inference approach under small sample cases. In the simulation study, we consider using r = 20 and 25 for the sample of size n = 30 and using r = 35 and 40 for the sample of size n = 50 . Moreover, we use the parameters of ( α , λ ) = ( 2 , 1 ) and ( 5 , 1 ) for the Weibull distribution. When ( α , λ ) = ( 2 , 1 ) , we consider using T = 1 , 1.5 , ( a , b ) = ( 3 , 2 ) and L = 0.05 . When ( α , λ ) = ( 5 , 1 ) , we consider using T = 1 , 1.15 , ( a , b ) = ( 4 , 5 ) and L = 0.3 . The Type-I hybrid censoring schemes can be obtained based on the aforementioned parameter information used for the simulation study. For each combination of parameters mentioned above, simulation study has been conducted for 5000 runs.
For each simulation run, the numbers of B = 10,000 and N 1 = 1000 are used to implement the MCMC for the Bayesian posterior distribution approximation and the number of B * = 1000 is used for the bootstrap sampling distribution approximation of MLE estimator. Two HPDIs of C L M for both Bayesian credible intervals and the bootstrap approximated confidence interval of C L M for MLE are obtained at the confidence level of 0.95; that is, γ = 0.05. The Baysian point estimate is obtained by using the posterior mean under the squared error loss, and the MLE is obtained via a procedure described in Section 3.
The measures of relative bias (rBias) and relative mean squared error (rMSE) will be used to measure the accuracies of MLE and two Bayesian estimators. Assume the target parameter is θ and its estimator is θ ^ , the rBias and rMSE can be defined, respectively, by
rBias = 1 θ × Bias
and
rMSE = 1 θ 2 × MSE ,
where Bias and MSE can be evaluated based on 5000 estimates, θ ^ j , j = 1 , 2 , , 5000 , by
Bias = 1 5000 j = 1 5000 θ ^ j θ
and
MSE = 1 5000 j = 1 5000 ( θ ^ j θ ) 2 ,
respectively. Please note that the rBias and rMSE are scale-free measures. Meanwhile the coverage probability (coverage) of 95 % confidence or credible interval is evaluated as the percentage of intervals covering the true parameter based on 5000 simulation runs. All obtained simulation results of coverage listed in the most right side columns of Table 1 and Table 2 are for all the interval estimators of C L M under study and the next three left columns contain rMSE and rBias for each point estimator, respectively. The symbols “Non-info.” and “Info.” in the tables denote the Bayesian estimation results through using a non-informative and informative prior distributions, respectively.
In scan of Table 1 and Table 2, we find that the rBias and rMSE for all parameters are decreased when the sample size and the number of failed products are increased. We also find that the values of rBias and rMSE of the MLE of C L M are larger than those of the two Bayes estimators. The Bayesian approach outperforms the maximum likelihood estimation method in terms of rBias and rMSE. In summary, we can also find that the Bayesian approach using an informative prior distribution performs better than the Bayesian inference approach using a non-informative prior distribution. However, the users could not have knowledge to set up the hyper-parameters of the prior distribution. In practical applications, the Bayesian approach with an non-informative prior distribution can be used to infer the life performance index of C L M for the Weibull distribution with a Type-I hybrid censoring scheme if the users do not have knowledge to set up the hyper-parameters of the prior distribution.

5. Examples

In this section, two examples from the literature will be provided for the illustrative purpose. Because no goodness-of-fit method is available for model fitting based on a Type-I hybrid censored sample. The data sets in two examples are complete samples. Firstly, we use Kolmogorov–Smirnov (K-S) testing method to confirm that these two data sets follow Weibull distributions. Then, reducing these two data sets into two Type-I hybrid censored samples to illustrate the applications of the proposed methods. The Type-I hybrid censoring scheme is a generalized scheme that includes the Type-I censoring and Type-II censoring schemes as special cases. Many lifetime applications, using a Type-I censoring scheme or Type-II censoring scheme to generate censored sample from a complete random one published, can be found in the literature. For example, Bain and Engelhardt [21] and Joarder et al. [22] presented data analysis based on a Type-I censored sample about 22 remission times of leukemia. The data set was collected through an experiment for 210 days scheme on a complete random sample published. Dodson [23] and Balakrishnan and Kateri [24] presented data analysis based on a Type-II censored sample about grinders testing, in which twenty grinders were tested and twelve failure times were observed. The Type-I hybrid censoring scheme can be a more efficient scheme for life testing to save testing time.
 Example 1. 
The example regarding the endurance of deep groove bearings in the life test is considered for illustrating the applications of the proposed estimation methods. This data set contains millions of revolutions before the failure of 23 ball bearings in the life test and is reported in Lieblein and Zelen [25] and listed in Table 3. Fatigue is a key factor which is highly related to the life of ball bearings. The experience of ball bearing manufacturers had led to relating the fatigue life L to the number of balls, ball diameter, number of rows, constant angle, and bearing load. L is the number of million revolutions that a specified percentage of bearing will fail to survive on account of fatigue causes. Using this complete random sample and K-S test to check the Weibull distribution modeling, the MLEs by utilizing R function optim is α ^ = 2 and λ ^ = 80 , and the K-S test statistic for the data set is 0.17391. The corresponding p-value = 0.8775. The K-S testing results indicate that the Weibull distribution with α = 2 and λ = 80 can well characterize this complete sample. Figure 1 is the quantile-to-quantile plot based on the complete sample of Example 1. In view of Figure 1, we see that the points form a line that is roughly straight. The quantile-to-quantile plot also shows that the Weibull distribution can well model the complete sample of Example 1.
The manufacturer wants to ensure the endurance yield of the deep groove ball bearing is 0.9975 with L = 4 , and the target process performance index is C L M = 1.67 . Assume that the manufacturer would like to save the test time and budget of the life testing, the Type-I hybrid censoring scheme can be a good option. Firstly, we can convert the current complete random sample to a Type-I hybrid censored sample. In this example, let T = 87 and r = 10 , 16 and 23. The hyper-parameters of a = 3 and b = 160 are used to set up the informative-prior Bayesian approaches. We would like to compare the estimation results based on the obtained Type-I hybrid censored sample via our proposed methods and based on the complete sample in Table 3 via the typical maximum likelihood estimation method. The estimation results by using the maximum likelihood estimation method with the complete sample can be a reference for performance comparison. If the likelihood function is complicated, it is known that the Bayesian estimation method outperforms typical maximum likelihood estimation method for small sample cases, the proposed Bayesian estimation method for the case of r = n can be anther reference method because the sample is a Type-I censored sample.
The MLEs, Bayes estimates, and 95% confidence intervals of the model parameters and process performance index were obtained and given in Table 4. The MLEs based on complete sample are α ^ = 2.1015 and λ ^ = 81.8743 , and C ^ L M = 1.7771 . The Bayes estimates using non-informative prior distribution based on the Type-I censored sample are α ˜ = 2.3638 , λ ˜ = 77.6118 , and C ˜ L M = 1.9954 . The Bayes estimates using informative prior distribution based on the Type-I censored sample are α ˜ = 2.2771 , λ ˜ = 79.1921 , and C ˜ L M = 1.9138 . Moreover, it shows that the confidence interval widths obtained using the two proposed Bayesian methods in Table 4 are shorter than that based on the maximum likelihood estimation method. The HPDI with informative prior distribution has the shortest width among three confidence intervals. Overall, we conclude that the endurance quality of the produced deep groove ball bearings could not satisfy the requirement due to the lower bounds in two HPDIs are smaller than 1.67.
 Example 2. 
The data of the waiting time (in minutes) of 100 customers before receiving the service in the bank is used as second example. This data is the complete sample data with sample size n = 100 shown in Table 5, see Ghitany et al. [26]. The K-S test statistic with the MLEs α ^ = 2 and λ ^ = 11 is 0.11. The corresponding p-value = 0.5806. The goodness-of-fit based on the K-S test indicates that the Weibull distribution with α = 2 and λ = 11 can well model this complete sample. The quantile-to-quantile plot based on the complete sample of Example 2 is displayed by Figure 2. In view of Figure 2, we can find that Weibull distribution can well model the complete sample of Example 2.
Let L = 0.5 . For the Weibull distribution with α = 2 and λ = 11 , we can obtain C L M = 1.6877 as the target process performance index. Following the similar data transformation method for Example 1 to convert the complete sample in Table 5 into a Type-I hybrid censored sample with T = 11 , and r = 35 , 70 , and 100. The Type-I hybrid censored sample with r = n = 100 is a Type-I censored sample. The hyper-parameters of a = 11 and b = 110 are considered for the informative prior distribution. The confidence level is set up at 95 % for interval estimation.
The obtained estimation results are reported in Table 6. From Table 6, we can see that the MLEs based on the complete sample in Table 5 are α ^ = 1.4618 , λ ^ = 10.9768 , and C ^ L M = 1.1403 . Because the sample size is large, the maximum likelihood estimation results based on the complete sample in Table 5 can be the benchmark, and we can find that the estimation results of the MLE and two proposed Bayesian methods based on the Type-I censoring sample (the case of r = n = 100 ) are close to their benchmarks. When r = 70 , the estimation results in Table 6 are close to that based on the case of r = 100 . In view of the Type-I censoring scheme in the last row of Table 6, we can find that the point estimates based on two Bayesian methods are more close to their benchmarks than the maximum likelihood estimation results. Moreover, the widths of two HDPIs are shorter than the width of the corresponding confidence interval based on the maximum likelihood method. The Bayesian estimation results outperform the estimation results via using the maximum likelihood estimation method. Hence, the Type-I hybrid censoring scheme can help to save testing time, and the proposed Bayesian estimation methods can provide reliable estimation results.

6. Conclusions and Remarks

A new life performance index C L M for the Weibull distribution has been proposed, and the maximum likelihood estimation and Bayesian estimation methods have been applied to develop estimation procedures for C L M based on Type-I hybrid censored samples. Point and interval estimation results have been obtained. The proposed life performance index was developed from the traditional life performance index by replacing the population mean with the population median, so that the proposed life performance index can also be applied to the lifetime data that follow a symmetric or asymmetric distribution.
Because no closed form available for the MLEs of the Weibull distribution parameters can be found under a Type-I hybrid censoring scheme. A numerical iterative procedure must be applied to obtained the MLEs for the model parameters and the proposed life performance index. The convergence of iterative procedure is sensitive to the initial inputs of two parameters. The Fisher information matrix under a Type-I hybrid censoring scheme is difficult to obtain due to the number of failures and the stopping time are random variables with unknown probability distribution. The observed Fisher information matrix is also less reliable. Therefore, a bootstrap procedure is proposed to find confidence interval for the MLE of the new life performance index.
The inverse Gamma distribution is used to obtain a conjugate informative prior for the scale parameter λ of Weibull distribution. However, there is no conjugate prior for the shape parameter α that can be found. Therefore, a non-informative prior distribution is suggested for α in this study. Two types of joint priors for α and λ have been considered to implement Bayesian inference. Due to the analytic form of the posterior distribution of the model parameter is difficult to developed, the MCMC approach via using a hybrid Metropolis–Hastings algorithm and Gibbs sampling is established to draw sample from the posterior distribution to update parameter in the numerical computation procedure.
Three different estimators, the MLE and two Bayes estimators, for the target parameters have been evaluated, and the performance of the maximum likelihood estimation method and two Bayesian approaches has been evaluated through using Monte Carlo simulation method. Simulation results show that two Bayesian approaches outperform the maximum likelihood estimation method in terms of the performance metrics of rBias and rMSE. Moreover, the coverage probability of HPDI based on using two proposed Bayesian approaches outperforms the one obtained based on using the maximum likelihood estimation method and bootstrapping procedure. If users lack knowledge to set up the hyper-parameters in the prior distribution of the parameters, a non-informative prior distribution has been suggested to implement the Bayesian approach.
The Type-I hybrid censoring scheme could not produce any lifetime observation, d = 0 , for highly reliable products. Hence, the maximum likelihood estimation results could not exist when no failure observed during the life testing. The proposed life performance index for the Weibull distribution could also be applied to any other widely used lifetime distributions. The current study focused on using the Type-I hybrid censoring scheme to save testing time and cost. The approximate confidence interval for the proposed life performance index by using Fisher information matrix could be available for other censoring competitive scheme. Meanwhile, the current proposed performance index could be improved by using nonparametric methods. All these topics are interesting and would be worth to investigate in the future.

Author Contributions

Conceptualization, T.-R.T. and Y.-J.H.; methodology, T.-R.T. and Y.L.; software, Y.-J.H. and T.-R.T.; validation, T.-R.T., J.-Y.C. and Y.L.; formal analysis, Y.L. and J.-Y.C.; investigation, T.-R.T., Y.L. and J.-Y.C.; resources, T.-R.T.; data curation, T.-R.T. and Y.-J.H.; writing—original draft preparation, J.-Y.C.; writing—review and editing, T.-R.T. and Y.L.; visualization, J.-Y.C.; supervision, T.-R.T., Y.L. and J.-Y.C.; project administration, T.-R.T.; funding acquisition, T.-R.T. 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 grant number MOST 110-2221-E-032-034-MY2.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data has been included in Section 5 of this article.

Acknowledgments

This study is supported by Ministry of Science and Technology, Taiwan MOST 110-2221-E-032-034-MY2. We thank you for the excellent suggestions from all reviewers to improve the quality of paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Proof of Theorem 1

(a)
Based on the condition of
M L μ L < 1 + δ 2 ,
we can obtain the inequality of
M L < 1 + δ 2 ( μ L ) .
Minus ( μ L ) at the two sides of Inequality (A2), we obtain
M L ( μ L ) < 1 + δ 2 ( μ L ) ( μ L ) .
Inequality (A3) can be represented by
M μ 1 + δ 2 1 < μ L .
Dividing σ at the two sides of Inequality (A4), we can show that
C L > δ 1 + δ 2 1 .
Then we can obtain the condition of
δ + C L C L 1 + δ 2 < 0 .
Thus,
C L M C L = δ 1 + δ 2 + C L 1 + δ 2 C L = 1 1 + δ 2 × δ + C L C L 1 + δ 2 .
Using (A6) for Equation (A7), we can show that C L M C L < 0 .
(b)
Based on similar inference procedures as (a), it is easy to show that
δ + C L C L 1 + δ 2 > 0 .
Thus,
C L M C L = 1 1 + δ 2 × δ + C L C L 1 + δ 2 .
Using (A8) for Equation (A9), we can show that C L M C L > 0 .
(c) Based on similar inference procedures as (a), it is easy to show that
δ + C L C L 1 + δ 2 = 0 .
Thus,
C L M C L = 1 1 + δ 2 × δ + C L C L 1 + δ 2 .
Using (A10) for Equation (A11), we can show that C L M = C L .

References

  1. Juran, J.M. Juran’s Quality Control Handbook, 3rd ed.; McGraw-Hill: New York, NY, USA, 1974. [Google Scholar]
  2. Kane, V.E. Process capability indices. J. Qual. Technol. 1986, 18, 41–52. [Google Scholar] [CrossRef]
  3. Hsiang, T.C.; Taguchi, G. A tutorial on quality control and assurance-the Taguchi mthods. In Joint Meetings of the American Statistical Association; InASA Annual Meeting LA: Las Vegas, NV, USA, 1985; p. 188. [Google Scholar]
  4. Pearn, W.L.; Kotz, S.; Johnson, N.L. Distributional and inferential properties of process capability indices. J. Qual. Technol. 1992, 24, 216–231. [Google Scholar] [CrossRef]
  5. Montgomery, D.C. Introduction to Statistical Quality Control, 6th ed.; John Wiley and Sons: New York, NY, USA, 2013; pp. 374–375. [Google Scholar]
  6. Clements, J.A. Process capability calculations for non-normal distributions. Qual. Prog. 1989, 22, 95–100. [Google Scholar]
  7. Epstein, B. Truncated life tests in the exponential case. Ann. Math. Stat. 1954, 25, 555–564. [Google Scholar] [CrossRef]
  8. Childs, A.; Chandrasekar, B.; Balakrishnan, N.; Kundu, D. Exact likelihood inference based on Type-I and Type-II hybrid censored samples from the exponential distribution. Ann. Inst. Stat. Math. 2003, 55, 319–330. [Google Scholar] [CrossRef]
  9. Lee, W.C. Statistical testing for assessing lifetime performance index of the Rayleigh lifetime products. J. Chin. Inst. Ind. Eng. 2008, 25, 433–445. [Google Scholar] [CrossRef]
  10. Hong, C.W.; Wu, J.W.; Cheng, C.H. Computational procedure of performance assessment of lifetime index of businesses for the Pareto lifetime model with the right type II censored sample. Appl. Math. Comput. 2007, 184, 336–350. [Google Scholar] [CrossRef]
  11. 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]
  12. 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]
  13. Wang, Y.; He, S. Fisher information in censored data. Stat. Probab. Lett. 2005, 73, 199–206. [Google Scholar] [CrossRef]
  14. Kundu, D.; Pradhan, B. Estimating the parameters of the generalized exponential distribution in presence of hybrid censoring. Commun. Stat. Theory Methods 2009, 38, 2030–2041. [Google Scholar] [CrossRef]
  15. Lin, C.T.; Chou, C.C.; Huang, Y.L. Inference for the Weibull distribution with progressive hybrid censoring. Comput. Stat. Data Anal. 2012, 56, 451–467. [Google Scholar] [CrossRef]
  16. Cho, Y.; Sun, H.; Lee, K. Estimating the entropy of a Weibull distribution under generalized progressive hybrid censoring. Entropy 2015, 17, 102–122. [Google Scholar] [CrossRef] [Green Version]
  17. Okasha, H.; Mustafa, A. E-Bayesian estimation for the Weibull distribution under adaptive type-I progressive hybrid censored competing risks data. Entropy 2020, 22, 903. [Google Scholar] [CrossRef]
  18. Kundu, D. Bayesian inference and life testing plan for the Weibull distribution in presence of progressive censoring. Technometrics 2008, 50, 144–154. [Google Scholar] [CrossRef]
  19. Kundu, D. On hybrid censored Weibull distribution. J. Stat. Plan. Inference 2007, 137, 2127–2142. [Google Scholar] [CrossRef] [Green Version]
  20. Chen, M.H.; Shao, Q.M. Monte Carlo estimation of Bayesian credible and HPD intervals. J. Comput. Graph. Stat. 1999, 8, 69–92. [Google Scholar]
  21. Bain, L.J.; Engelhardt, M. Statistical Analysis of Reliability and Life-Testing Models; Marcel Dekker: New York, NY, USA, 1991. [Google Scholar]
  22. Joarder, A.; Krishna, H.; Kundu, D. Infrences on Weibull parameters with conventional type-I censoring. Comput. Stat. Data Anal. 2011, 55, 1–11. [Google Scholar] [CrossRef]
  23. Dodson, B. The Weibull Analysis Handbook, 2nd ed.; ASQ Quality Press: Milwaukee, WI, USA, 2006. [Google Scholar]
  24. Balakrishnan, N.; Kateri, M. On the maximum likelihood estimation of parameters of Weibull distribution based on complete and censored data. Stat. Probab. Lett. 2008, 55, 2971–2975. [Google Scholar] [CrossRef]
  25. Lieblein, J.; Zelen, M. Statistical investigation of the fatigue life of deep-groove ball bearings. J. Res. Natl. Bur. Stand. 1956, 57, 273–316. [Google Scholar] [CrossRef]
  26. Ghitany, M.E.; Atieh, B.; Nadarajah, S. Lindley distribution and its application. Math. Comput. Simul. 2008, 78, 493–506. [Google Scholar] [CrossRef]
Figure 1. The quantile-to-quantile plot based on the complete sample of Example 1.
Figure 1. The quantile-to-quantile plot based on the complete sample of Example 1.
Mathematics 10 04090 g001
Figure 2. The quantile-to-quantile plot based on the complete sample of Example 2.
Figure 2. The quantile-to-quantile plot based on the complete sample of Example 2.
Mathematics 10 04090 g002
Table 1. The rBias and rMSE (in parentheses) for ( α , λ ) = ( 2 , 1 ) and L = 0.05 .
Table 1. The rBias and rMSE (in parentheses) for ( α , λ ) = ( 2 , 1 ) and L = 0.05 .
( n , T ) rEstimation α λ C LM Coverage
(30, 1)20MLE0.117 (0.013)−0.022 (0.001)0.118 (0.014)0.864
Non-info.0.063 (0.004)0.025 (0.001)0.051 (0.002)0.945
Info.0.074 (0.005)<0.001 (<0.001)0.062 (0.003)0.945
25MLE0.108 (0.011)−0.023 (0.001)0.108 (0.011)0.862
Non-info.0.056 (0.003)0.022 (0.001)0.043 (0.001)0.948
Info.0.066 (0.004)−0.002 (<0.001)0.054 (0.002)0.949
(30, 1.5)20MLE0.084 (0.007)−0.011 (<0.001)0.084 (0.007)0.895
Non-info.0.037 (0.001)0.022 (0.001)0.025 (0.002)0.950
Info.0.043 (0.001)0.004 (<0.001)0.031 (0.001)0.950
25MLE0.070 (0.005)−0.010 (<0.001)0.072 (0.005)0.899
Non-info.0.035 (0.001)0.006 (<0.001)0.027 (<0.001)0.948
Info.0.034 (0.001)−0.005 (<0.001)0.025 (0.001)0.944
(50, 1)35MLE0.065(0.004)−0.015 (<0.001)0.066 (0.004)0.893
Non-info.0.035 (0.001)0.008 (<0.001)0.028 (0.001)0.949
Info.0.040 (0.001)−0.002 (<0.001)0.033 (0.001)0.949
40MLE0.064 (0.004)−0.013 (<0.001)0.065 (0.004)0.892
Non-info.0.035 (0.001)0.010 (<0.001)0.028 (0.001)0.948
Info.0.040 (0.002)<−0.001 (<0.001)0.033 (0.001)0.950
(50, 1.5)35MLE0.047 (0.002)−0.006 (<0.001)0.048 (0.002)0.905
Non-info.0.022 (0.001)0.009 (<0.001)0.015 (<0.001)0.947
Info.0.024 (0.001)0.001 (<0.001)0.017 (<0.001)0.948
40MLE0.041 (0.001)−0.003 (<0.001)0.042 (0.001)0.913
Non-info.0.020 (<0.001)0.007 (<0.001)0.015 (<0.001)0.949
Info.0.020 (<0.001)<0.001 (<0.001)0.014 (<0.001)0.945
Table 2. The rBias and rMSE (in parentheses) for ( α , λ ) = ( 5 , 1 ) and L = 0.3 .
Table 2. The rBias and rMSE (in parentheses) for ( α , λ ) = ( 5 , 1 ) and L = 0.3 .
( n , T ) rEstimation α λ C LM Coverage
(30, 1)20MLE0.113 (0.012)−0.012 (<0.001)0.097 (0.009)0.874
Non-info.0.059 (0.003)0.002 (<0.001)0.047 (0.002)0.951
Info.0.057 (0.003)0.005 (<0.001)0.046 (0.002)0.953
25MLE0.113 (0.012)−0.011 (<0.001)0.097 (0.009)0.864
Non-info.0.060 (0.003)0.003 (<0.001)0.049 (0.002)0.953
Info.0.058 (0.003)0.006 (<0.001)0.048 (0.002)0.953
(30, 1.15)20MLE0.090 (0.008)−0.006 (<0.001)0.078 (0.006)0.883
Non-info.0.043 (0.001)0.004 (<0.001)0.035 (0.001)0.946
Info.0.041 (0.001)0.007 (<0.001)0.034 (0.001)0.943
25MLE0.073 (0.005)−0.004 (<0.001)0.065 (0.004)0.896
Non-info.0.037 (0.001)0.001 (<0.001)0.031 (0.001)0.947
Info.0.038 (0.001)0.002 (<0.001)0.033 (0.001)0.950
(50, 1)35MLE0.066 (0.004)−0.005 (<0.001)0.057 (0.003)0.900
Non-info.0.036 (0.001)0.002 (<0.001)0.030 (0.001)0.951
Bayes20.035 (0.001)0.003 (<0.001)0.029 (0.001)0.949
40MLE0.064 (0.004)−0.006 (<0.001)0.055 (0.003)0.896
Non-info.0.035 (0.001)0.001 (<0.001)0.028 (0.001)0.946
Info.0.034 (0.001)0.002 (<0.001)0.027 (0.001)0.947
(50, 1.15)35MLE0.050 (0.003)−0.003 (<0.001)0.044 (0.002)0.910
Non-info.0.025 (0.001)0.003 (<0.001)0.020 (<0.001)0.946
Info.0.024 (0.001)0.003 (<0.001)0.020 (<0.001)0.948
40MLE0.039 (0.002)−0.003 (<0.001)0.034 (0.001)0.919
Non-info.0.017 (<0.001)0.001 (<0.001)0.014 (<0.001)0.948
Info.0.017 (<0.001)0.002 (<0.001)0.014 (<0.001)0.951
Table 3. The endurance life of deep groove ball bearings (in million revolutions) will fail to survive.
Table 3. The endurance life of deep groove ball bearings (in million revolutions) will fail to survive.
17.8828.9233.0041.5242.1245.60
48.4851.8451.9654.1255.5667.80
68.6468.6468.8884.1293.1298.64
105.12105.84127.92128.04173.40
Table 4. Estimation results based on Example 1 ( n = 23 ball bearings and T = 87 in million revolutions).
Table 4. Estimation results based on Example 1 ( n = 23 ball bearings and T = 87 in million revolutions).
rEstimation α λ C LM Confidence Interval or HPDI
MLE (Complete data)2.101581.87431.7771(1.1403, 3.2407)
10MLE3.609363.70123.0302(1.9612, 6.5553)
Non-Info.3.298966.83052.7601(1.0169, 4.8253)
Info.3.263868.31202.7290(1.1303, 4.4750)
16MLE2.469276.70072.1059(1.4421, 3.6492)
Non-Info.2.335179.32181.9704(0.8488, 3.1089)
Info.2.394877.17402.0229(0.9891, 3.0498)
23MLE2.105976.70072.4692(1.4518, 3.7428)
Non-Info.2.363877.61181.9954(0.9942, 2.9651)
Info.2.277179.19211.9138(0.9167, 2.9240)
Table 5. Waiting time data of customers before receiving service in the bank (in minutes).
Table 5. Waiting time data of customers before receiving service in the bank (in minutes).
0.80.81.31.51.81.91.92.12.62.7
2.93.13.23.33.53.64.04.14.24.2
4.34.34.44.44.64.74.74.84.94.9
5.05.35.55.75.76.16.26.26.26.3
6.76.97.17.17.17.17.47.67.78.0
8.28.68.68.68.88.88.98.99.59.6
9.79.810.710.911.011.011.111.211.211.5
11.912.412.512.913.013.113.313.613.713.9
14.115.415.417.317.318.118.218.418.919.0
19.920.621.321.421.923.027.031.633.138.5
Table 6. The estimation results of Example 2 (n = 100 waiting times and T = 11 in minutes).
Table 6. The estimation results of Example 2 (n = 100 waiting times and T = 11 in minutes).
rEstimation α λ C LM Confidence Interval or HPDI
MLE (Complete data)1.461810.97681.1403(0.8396, 1.5746)
35MLE1.97368.73901.6373(1.2207, 2.3970)
Non-Info.1.93639.04631.5946(1.0437, 2.1372)
Info.1.90629.18631.5673(1.0496, 2.1682)
70MLE1.671710.21631.3551(0.9573, 2.0706)
Non-Info.1.653910.32011.3338(0.9873, 1.6867)
Info.1.663210.27371.3432(1.0167, 1.6850)
100MLE1.642610.35741.3263(1.0244, 1.7854)
Non-Info.1.627410.42791.3081(1.0028, 1.6236)
Info.1.623210.37461.3031(0.9690, 1.6093)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tsai, T.-R.; Lio, Y.; Chiang, J.-Y.; Huang, Y.-J. A New Process Performance Index for the Weibull Distribution with a Type-I Hybrid Censoring Scheme. Mathematics 2022, 10, 4090. https://doi.org/10.3390/math10214090

AMA Style

Tsai T-R, Lio Y, Chiang J-Y, Huang Y-J. A New Process Performance Index for the Weibull Distribution with a Type-I Hybrid Censoring Scheme. Mathematics. 2022; 10(21):4090. https://doi.org/10.3390/math10214090

Chicago/Turabian Style

Tsai, Tzong-Ru, Yuhlong Lio, Jyun-You Chiang, and Yi-Jia Huang. 2022. "A New Process Performance Index for the Weibull Distribution with a Type-I Hybrid Censoring Scheme" Mathematics 10, no. 21: 4090. https://doi.org/10.3390/math10214090

APA Style

Tsai, T. -R., Lio, Y., Chiang, J. -Y., & Huang, Y. -J. (2022). A New Process Performance Index for the Weibull Distribution with a Type-I Hybrid Censoring Scheme. Mathematics, 10(21), 4090. https://doi.org/10.3390/math10214090

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