# Estimation for Weibull Parameters with Generalized Progressive Hybrid Censored Data

Authors: Xuanjia Zuo; Liang Wang; Huizhong Lin; Sanku Dey; Li Yan
Journal: Journal of Mathematics
 (2021)

---

## Abstract

In this paper, the interest is in estimating the Weibull products when the available data is obtained via generalized progressive hybrid censoring. The testing scheme conducts products of interest under a more flexible way and allows collecting failure data in efficient and adaptable experimental scenarios than traditional lifetime testing. When the latent lifetime of products follows Weibull distribution, classical and Bayesian inferences are considered for unknown parameters. The existence and uniqueness of maximum likelihood estimates are established, and approximate confidence intervals are also constructed via asymptotic theory. Bayes point estimates as well as the credible intervals of the parameters are obtained, and correspondingly, Monte Carlo sampling technique is also provided for complex posterior computation. Extensive numerical analysis is carried out, and the results show that the generalized progressive hybrid censoring is an adaptive procedure in practical lifetime experiment, both proposed classical and Bayesian inferential approaches perform satisfactorily, and the Bayesian results are superior to conventional likelihood estimates.


---

## Body

## 1. Introduction

In lifetime study and reliability, due to many practical reasons such as time constraint and cost reduction, the experiment cannot be finished with all test unit failures; therefore, censoring appears as a common phenomenon in many life test studies. Generally speaking, there is only a portion of the testing units of which the exact failure times are known when censoring is involved under study. In practice, there are several types of schemes (CSs), and the most used ones are Type-I censoring and Type-II censoring, where experiments are terminated after a prefixed inspecting time or a predetermined number of failures. The conventional CSs, however, do not allow engineers to remove survival units from the experiment at points other than the final termination time. For this reason, a progressive censoring scheme was further introduced in practice, where engineers could remove units from testing at various stages. Like conventional CSs, progressive censoring also includes progressive Type-I censoring and progressive Type-II censoring. Compared to Type-I/-II censoring, it is seen that the progressive CSs are more flexible and efficient in lifetime study, and extensive works have been done on these directions. See, for example, some recent contributions of Yan et al. [1], Laumen and Cramer [2], Gunasekera [3], and Maurya and Tripathi [4] as well as the reference therein. For an exhaustive list of references and details on progressive censoring, one can refer to the monographs of Balakrishnan and Cramer [5].With the improvement in manufacturing design and technology, modern products feature high reliability and long span; under such situations, there may not be enough failure samples or the time of the experiment could be very long in both conventional and progressive censoring. To overcome such problems, hybrid censoring schemes were further introduced in practice, which can be viewed as a mixture of Type-I and Type-II censoring in conventional and progressive censoring procedures, respectively. For instance, Prajapati et al. [6] and Algarni et al. [7] discussed inference based on Type-I and Type-II hybrid censoring, and Basu et al. [8], Gorny and Cramer [9], and Singh et al. [10] studied progressive hybrid censored data under different cases. For an exhaustive amount of references and further details on hybrid censoring, one can see the recent review paper [11] and the monograph of Balakrishnan and Cramer [5] within related sections. Recently, Cho et al. [12] proposed a generalized progressive hybrid censoring scheme (GPHCS) where an adequate number of observations are guaranteed within a reasonable experimental duration which in turn may yield improved estimation results, and the detailed procedure of such test can be described as follows: suppose that n independent items are put in a test. Furthermore, assume that k and m are prefixed constants satisfying 1≤k<m≤n, T is also prespecified monitoring time, and ri≥0, i=1,…,m is prefixed censoring schemes with ∑i=1mri+m=n. Following the similar way as progressive censoring, at time of i th failure Xi:m:n, the ri number of live units is withdrawn from the experiment at random. In this case, the test stops at point T∗=maxXk:m:n,minT,Xm:m:n, i.e.,(1)T∗=Xk:m:n,ifT<Xk:m:n<Xm:m:n,T,ifXk:m:n<T<Xm:m:n,Xm:m:n,ifXk:m:n<Xm:m:n<T.Moreover, under GPHCS, one has the following three types of observations:(2)case I:X1:m:n,X2:m:n,…,Xk:m:n,ifT<Xk:m:n<Xm:m:n,case II:X1:m:n,…,Xk:m:n,…,Xd:m:n,ifXk:m:n<T<Xm:m:n,case III:X1:m:n,…,Xk:m:n,…,Xm:m:n,ifXk:m:n<Xm:m:n<T,where for case II, Xd:m:n<T<Xd+1:m:n. A representation of GPHCS is shown in Figure 1, where rk∗=m−k+∑i=kmri and rd+1∗=m−1+∑i=d+1mri. Therefore, from the test procedure of GPHCS, it is noted that if inadequate observations are recorded, then the test duration can be extended beyond T for obtaining more test information. Moreover, it is also observed that some of the well-known censoring methods including progressive and progressive hybrid censoring are specific cases of GPHCS.Figure 1
A sketch of the generalized progressive hybrid censoring scheme.Denotef⋅ and S⋅ are the PDF and survival functions of the latent failure time X, respectively. Based on (2), the joint density function can be written as(3)Ldata=c1∏i=1kfxi:m:n1−Fxi:m:nri,case I,c2∏i=1dfxi:m:n1−Fxi:m:nriSTrd+1∗,case II,c3∏i=1mfxi:m:n1−Fxi:m:nri,case III,where rd+1∗=n−∑i=1dri+1, rk=n−∑i=1k−1ri+1−1 in case I, and(4)c1=∏j=1k∑k=jmrk+1,ifT<Xk:m:n<Xm:m:n,c2=∏j=1d∑k=jmrk+1,ifXk:m:n<T<Xm:m:n,c3=∏j=1m∑k=jmrk+1,ifXk:m:n<Xm:m:n<T.In this paper, we assume that the latent failure time of units follows a Weibull distribution with shape parameterα>0 and scale parameter β>0, of which the probability density function (PDF) and the cumulative distribution function (CDF) can be expressed as(5)fx;α,β=αβxα−1e−βxα,Fx;α,β=1−e−βxα,x>0.Because of its versatility in fitting time-to-failure distributions of a rather extensive variety of complex mechanisms, the Weibull distribution is one of the most popular distributions used in statistical inference of reliability, lifetime study, and survival analysis. Various problems associated with this distribution have been discussed by numerous authors. See, for example, the works of Almetwally and Almongy [13], Chen et al. [14], and Kumar and Ram [15] as well as the reference therein.Suppose latent failure time (2) of units of sample size n follows Weibull distribution (5) with parameters α>0 and β>0; the likelihood function of α and β, namely, Lα,β;data, can be expressed as(6)Lα,β;data=c1∏i=1kαβxi:m:nα−1e−βxi:m:nαe−βxi:m:nαri,ifT<xk:m:n<xm:m:n,c2∏i=1dαβxi:m:nα−1e−βxi:m:nαe−βxi:m:nαrie−βTαrd+1∗,ifxk:m:n<T<xm:m:n,c3∏i=1mαβxi:m:nα−1e−βxi:m:nαe−βxi:m:nαri,ifxk:m:n<xm:m:n<T.Motivated by such reasons as mentioned above and due to the extensive applicability and practical significance of Weibull lifetime distribution in many real-life studies, we discussed the statistical inference for the generalized progressive hybrid censoring when the latent failure data follows the Weibull model. Our objectives in this paper are as follows: First, we estimate the parameters of the Weibull distribution under generalized progressive hybrid censoring using classical estimation methods, namely, maximum likelihood estimation (MLE). In addition, we have obtained approximate confidence intervals (ACIs) based on the asymptotic distribution of MLE. The performance of the associated estimates are demonstrated in terms of their mean squared error and root absolute bias for different sample sizes through a simulation study. The second objective is to obtain Bayes estimates based on likelihood function based on independent gamma priors. Since Bayes estimates cannot be obtained in explicit forms, then Markov Chain Monte Carlo (MCMC) sampling is used to compute the complex posterior functions and in turn to calculate Bayes estimates as well as associated highest posterior density (HPD) credible intervals. Besides, one real-life data set is analyzed to illustrate the applicability of the proposed methods.The rest of this paper is organized as follows: Classic and Bayes estimates are presented in Sections2 and 3, respectively. Section 4 provides numerical illustrations, and some concluding remarks are given in Section 5.


## 2. Frequentist Inference

MLEs are established for unknown parameters in this section, and ACIs are constructed based on asymptotic theory.


### 2.1. Maximum Likelihood Estimation

In order to find MLEs of unknown parameters in a concise way, likelihood function (6) can be rewritten in compact expression as(7)Lα,β;data=c∗αd∗βd∗e−βwα∏i=1d∗xi:m:nα−1,where c∗=∏j=1d∗∑k=jmrk+1,(8)d∗=k,for case I,d,for case II,m,for case III,,wα=∑i=1d∗1+rixi:m:nα,for case I,∑i=1d∗1+rixi:m:nα+rd+1∗Tα,for case II,∑i=1d∗1+rixi:m:nα,for case III.From (7), the associated log-likelihood function, namely, ℓα,β, is given by(9)ℓα,β=lnc∗+d∗lnα+d∗lnβ+α−1∑i=1d∗lnxi:m:n−βwα.Theorem 1.
Suppose that observation (2) comes from Weibull population (5); the MLE of β given α can be written as β˜=d∗/wα.Proof.
See Appendix A.Furthermore, substitutingβ˜=β˜α into (9), one has a profile log-likelihood function of α, namely, ℓα, as(10)ℓα=lnc∗+d∗lnα+α−1∑i=1d∗lnxi:m:n+d∗lnd∗wα−d∗.The following result establishes the MLE for parameterα.Theorem 2.
Suppose that observation (2) comes from Weibull population (5); the MLE of α obtained from (10) uniquely exists which can be derived from the following equation:(11)1α+1d∗∑i=1d∗lnxi:m:n−w′αwα=0,where(12)w′α=∑i=1k1+rixi:m:nαlnxi:m:n,for case I,∑i=1d1+rixi:m:nαlnxi:m:n+rd+1∗TαlnT,for case II,∑i=1m1+rixi:m:nαlnxi:m:n,for case III.Proof.
See Appendix B.Denote the MLE of parameterα obtained from (11) in Theorem 2 as α^; the MLEs of β, namely, β^, can be further obtained from Theorem 1 as β^=d∗/wα^. It is noted that the MLE of α cannot be expressed as closed from (11), and the numerical procedure such as fixed-point iterative algorithm can be used to find the associated estimate.


### 2.2. Approximation Confidence Interval Estimation

Denoteθ=θ1,θ2′ with θ1=α,θ2=β; it is known under mild regularity that the asymptotic distribution of the MLE θ^ is θ^−θ⟶N0,I−1θ, where ‘⟶’ means “distributed as,” and I−1θ is the inverse of the observed information matrix given by(13)I−1θ^=−∂2ℓα,β∂θj∂θk2×2−1θ=θ^=Varα^Covα^,β^Covα^,β^Varβ^,where the second derivatives of ℓα,β can be obtained directly which are omitted here for saving space. Therefore, for arbitrary 0<γ<1, the 1001−γ% ACIs of θj,j=1,2, is given by(14)θ^j−zγ/2Varθ^j,θ^j+zγ/2Varθ^j,j=1,2,where zγ is the upper γ-th quantile of the standard normal distribution.Sometimes, the ACIs obtained by the previous procedure may have negative lower bounds. In order to overcome this drawback, the logarithmic transformation and delta methods can be used to obtain the asymptotic normality distribution oflnθ^j as(15)lnθ^j−lnθjVarlnθ^j⟶N0,1,j=1,2.Therefore, a1001−γ% ACI of θj,j=1,2, derived in this manner is provided by(16)θ^jexpzγ/2Var^lnθ^j,θ^jexpzγ/2Var^lnθ^j,j=1,2,where Var^lnθ^j=Var^θ^j/θ^j,j=1,2.


## 2.1. Maximum Likelihood Estimation

In order to find MLEs of unknown parameters in a concise way, likelihood function (6) can be rewritten in compact expression as(7)Lα,β;data=c∗αd∗βd∗e−βwα∏i=1d∗xi:m:nα−1,where c∗=∏j=1d∗∑k=jmrk+1,(8)d∗=k,for case I,d,for case II,m,for case III,,wα=∑i=1d∗1+rixi:m:nα,for case I,∑i=1d∗1+rixi:m:nα+rd+1∗Tα,for case II,∑i=1d∗1+rixi:m:nα,for case III.From (7), the associated log-likelihood function, namely, ℓα,β, is given by(9)ℓα,β=lnc∗+d∗lnα+d∗lnβ+α−1∑i=1d∗lnxi:m:n−βwα.Theorem 1.
Suppose that observation (2) comes from Weibull population (5); the MLE of β given α can be written as β˜=d∗/wα.Proof.
See Appendix A.Furthermore, substitutingβ˜=β˜α into (9), one has a profile log-likelihood function of α, namely, ℓα, as(10)ℓα=lnc∗+d∗lnα+α−1∑i=1d∗lnxi:m:n+d∗lnd∗wα−d∗.The following result establishes the MLE for parameterα.Theorem 2.
Suppose that observation (2) comes from Weibull population (5); the MLE of α obtained from (10) uniquely exists which can be derived from the following equation:(11)1α+1d∗∑i=1d∗lnxi:m:n−w′αwα=0,where(12)w′α=∑i=1k1+rixi:m:nαlnxi:m:n,for case I,∑i=1d1+rixi:m:nαlnxi:m:n+rd+1∗TαlnT,for case II,∑i=1m1+rixi:m:nαlnxi:m:n,for case III.Proof.
See Appendix B.Denote the MLE of parameterα obtained from (11) in Theorem 2 as α^; the MLEs of β, namely, β^, can be further obtained from Theorem 1 as β^=d∗/wα^. It is noted that the MLE of α cannot be expressed as closed from (11), and the numerical procedure such as fixed-point iterative algorithm can be used to find the associated estimate.


## 2.2. Approximation Confidence Interval Estimation

Denoteθ=θ1,θ2′ with θ1=α,θ2=β; it is known under mild regularity that the asymptotic distribution of the MLE θ^ is θ^−θ⟶N0,I−1θ, where ‘⟶’ means “distributed as,” and I−1θ is the inverse of the observed information matrix given by(13)I−1θ^=−∂2ℓα,β∂θj∂θk2×2−1θ=θ^=Varα^Covα^,β^Covα^,β^Varβ^,where the second derivatives of ℓα,β can be obtained directly which are omitted here for saving space. Therefore, for arbitrary 0<γ<1, the 1001−γ% ACIs of θj,j=1,2, is given by(14)θ^j−zγ/2Varθ^j,θ^j+zγ/2Varθ^j,j=1,2,where zγ is the upper γ-th quantile of the standard normal distribution.Sometimes, the ACIs obtained by the previous procedure may have negative lower bounds. In order to overcome this drawback, the logarithmic transformation and delta methods can be used to obtain the asymptotic normality distribution oflnθ^j as(15)lnθ^j−lnθjVarlnθ^j⟶N0,1,j=1,2.Therefore, a1001−γ% ACI of θj,j=1,2, derived in this manner is provided by(16)θ^jexpzγ/2Var^lnθ^j,θ^jexpzγ/2Var^lnθ^j,j=1,2,where Var^lnθ^j=Var^θ^j/θ^j,j=1,2.


## 3. Bayesian Inference

As a powerful and valid alternative to classical estimation, the Bayesian viewpoint has received much attention in statistical inference. The capability of incorporating prior information in analysis makes it very valuable in reliability, lifetime study, and other associated fields where one of the major challenges is the limited availability of data.


### 3.1. Prior Information

In this section, joint prior forα and β is considered as follows:(17)πα,β∝αa1−1e−b1αβa2−1e−b2β,α>0,β>0,where a1>0,b1>0,a2>0,andb2>0 are hyperparameters. It is seen that the priors of parameters α and β are independently gamma-distributed respectively, and the hyperparameters are chosen to reflect prior knowledge. The information abstracted by the priors may be either objective or subjective based on failure data from a history knowledge, personal experience, and belief among others. Moments-matching and empirical Bayesian methods may be used to fit the prior densities which are not reported here. When aj=bj=0,j=1,2, these priors reduce to noninformative and are also nonproper priors. Furthermore, the associated joint posterior density function of α and β, namely, πα,β|data, can be expressed as(18)πα,β|data=πα,βLα,β;data∫0∞∫0∞πα,βLα,β;datadαdβ.Therefore, under squared error loss, the Bayes estimator of any function ofα and β, say ηα,β, is given by(19)η^α,β=∫0∞∫0∞ηα,βπα,βLα,β;datadαdβ∫0∞∫0∞πα,βLα,β;datadαdβ.Clearly, the corresponding Bayes estimator ofηα,β does not exist in closed form; thus, a Markov Chain Monte Carlo (MCMC) sampling method is used in the following section to compute the associated estimate.


### 3.2. Posterior Analysis and MCMC Sampling

From (7) and (17) and ignoring the additive constant, joint posterior PDF (18) of α and β can be rewritten as(20)πα,β|data=πα|dataπβ|α,data,with(21)πα|data∝αa1+d∗−1b2+wαa2+d∗e−αb1−∑i=1d∗lnxi:m:n,α>0,(22)πβ|α,data=b2+wαa2+d∗Γa2+d∗βa2+d∗−1e−βb2+wα,β>0.One can observe that for a givenα, the posterior density function of β is gamma-distributed. In addition, by direct computation, one has the following result for πα|data of which the proof is omitted for concision.Lemma 1.
The posterior marginal PDFπα|data in (21) is log-concave.In order to generate random samples from the posterior marginal PDFπα|data, the Metropolis–Hastings (M-H) algorithm with the gamma proposal distribution is used in our study (see, e.g., Gelman et al. [16]) for the log-concave posterior density. Of course, some other sampling algorithms like the approach for computing the log-concave density (e.g., Devroye [17]) and the adaptive rejection sampling approach (e.g., Gilks and wild [18]) could also be used to generate the sample α from the density πα|data as well. Therefore, the corresponding Bayes estimates as well as the HPD credible intervals can be constructed, and an algorithm is proposed to generate samples and to compute Bayes estimates (Algorithm 1).Algorithm 1: MCMC algorithm for Bayes estimation.
Step 1 Generate α from πα|data in (21) by using the proper sampling approach.Step 2 For given α, generate β from πβ|α,data in (22).Step 3 Repeat above Step 1 and Step 2 N times, and N samples of α,β are generated as αk,βk,k=1,2,…,N.Step 4 The Bayes estimate of ηα,β with respect to squared error loss, can be constructed asη^α,β=1/N−N0∑k=N0+1Nη^kandη^k=ηαk,βk,whereN0 is burn-in period.Step 5 To construct the HPD credible interval of ηα,β, first arrange all estimates η^k, k=N0+1,N0+2,…,N in an ascend order, as η^1,η^2,…,η^N−N0, then for arbitrary 0<γ<1 the 1001−γ credible interval of ηα,β can be obtained asη^k,η^k+N−N0−γN−N0+1,k=1,2,…,N−N0γwherey denotes the greatest integer less than or equal to y. Therefore, the 1001−γ% HPD credible interval can be constructed as the k∗ th one satisfyingη^k∗+N−N0−γN−N0+1−η^k∗=mink=1N−N0γη^k+N−N0−γN−N0+1−η^k


## 3.1. Prior Information

In this section, joint prior forα and β is considered as follows:(17)πα,β∝αa1−1e−b1αβa2−1e−b2β,α>0,β>0,where a1>0,b1>0,a2>0,andb2>0 are hyperparameters. It is seen that the priors of parameters α and β are independently gamma-distributed respectively, and the hyperparameters are chosen to reflect prior knowledge. The information abstracted by the priors may be either objective or subjective based on failure data from a history knowledge, personal experience, and belief among others. Moments-matching and empirical Bayesian methods may be used to fit the prior densities which are not reported here. When aj=bj=0,j=1,2, these priors reduce to noninformative and are also nonproper priors. Furthermore, the associated joint posterior density function of α and β, namely, πα,β|data, can be expressed as(18)πα,β|data=πα,βLα,β;data∫0∞∫0∞πα,βLα,β;datadαdβ.Therefore, under squared error loss, the Bayes estimator of any function ofα and β, say ηα,β, is given by(19)η^α,β=∫0∞∫0∞ηα,βπα,βLα,β;datadαdβ∫0∞∫0∞πα,βLα,β;datadαdβ.Clearly, the corresponding Bayes estimator ofηα,β does not exist in closed form; thus, a Markov Chain Monte Carlo (MCMC) sampling method is used in the following section to compute the associated estimate.


## 3.2. Posterior Analysis and MCMC Sampling

From (7) and (17) and ignoring the additive constant, joint posterior PDF (18) of α and β can be rewritten as(20)πα,β|data=πα|dataπβ|α,data,with(21)πα|data∝αa1+d∗−1b2+wαa2+d∗e−αb1−∑i=1d∗lnxi:m:n,α>0,(22)πβ|α,data=b2+wαa2+d∗Γa2+d∗βa2+d∗−1e−βb2+wα,β>0.One can observe that for a givenα, the posterior density function of β is gamma-distributed. In addition, by direct computation, one has the following result for πα|data of which the proof is omitted for concision.Lemma 1.
The posterior marginal PDFπα|data in (21) is log-concave.In order to generate random samples from the posterior marginal PDFπα|data, the Metropolis–Hastings (M-H) algorithm with the gamma proposal distribution is used in our study (see, e.g., Gelman et al. [16]) for the log-concave posterior density. Of course, some other sampling algorithms like the approach for computing the log-concave density (e.g., Devroye [17]) and the adaptive rejection sampling approach (e.g., Gilks and wild [18]) could also be used to generate the sample α from the density πα|data as well. Therefore, the corresponding Bayes estimates as well as the HPD credible intervals can be constructed, and an algorithm is proposed to generate samples and to compute Bayes estimates (Algorithm 1).Algorithm 1: MCMC algorithm for Bayes estimation.
Step 1 Generate α from πα|data in (21) by using the proper sampling approach.Step 2 For given α, generate β from πβ|α,data in (22).Step 3 Repeat above Step 1 and Step 2 N times, and N samples of α,β are generated as αk,βk,k=1,2,…,N.Step 4 The Bayes estimate of ηα,β with respect to squared error loss, can be constructed asη^α,β=1/N−N0∑k=N0+1Nη^kandη^k=ηαk,βk,whereN0 is burn-in period.Step 5 To construct the HPD credible interval of ηα,β, first arrange all estimates η^k, k=N0+1,N0+2,…,N in an ascend order, as η^1,η^2,…,η^N−N0, then for arbitrary 0<γ<1 the 1001−γ credible interval of ηα,β can be obtained asη^k,η^k+N−N0−γN−N0+1,k=1,2,…,N−N0γwherey denotes the greatest integer less than or equal to y. Therefore, the 1001−γ% HPD credible interval can be constructed as the k∗ th one satisfyingη^k∗+N−N0−γN−N0+1−η^k∗=mink=1N−N0γη^k+N−N0−γN−N0+1−η^k


## 4. Numerical Illustration

### 4.1. Simulation Studies

In this section, extensive simulations are carried out for comparing different proposed methods and the performance of the point estimates (i.e., MLEs and Bayes) is evaluated in terms of mean squared error (MSE) and average bias (AB), whereas the confidence interval estimates (i.e., ACIs and HPDs) are compared by criteria coverage probability (CP) and average width (AW), respectively.The evaluaton criteria quantities for point and interval estimates were examined under different choices ofβ1,β2,α,m,n, CS, r1,r2,…,rm, and T values, and the Bayes sampling algorithm was repeated 10,000 times following a burn-in phase of 5,000 iterations. For generating the generalized progressively hybrid censored Weibull samples, it is conducted as the following procedure. Using the algorithm proposed by Balakrishnan and Aggarwala ([19], p. 32), progressively Type-II censored data is generated from Weibull population with parameters α and β. Then, the generalized progressively hybrid censored sample can be obtained by comparing prefixed T and Xk:m:n with those of the generated progressively Type-II censored data.In the simulation procedure, following CSs are considered:CS-1:n=18,m=15,r1=⋯=r14=0, and r15=3CS-2:n=25,m=15,r1=⋯=r14=0, and r15=10CS-3:n=25,m=20,r1=⋯=r19=0, and r20=5CS-4:n=50,m=40,r1=⋯=r39=0, and r40=10CS-5:n=50,m=40,r2=⋯=r40=0, and r1=10CS-6:n=50,m=40,r2=⋯=r39=0, and r1=r40=5In the simulation procedure, two sets of parameter values are chosen and the Bayes estimates are computed under noninformative prior (NIP) and informative prior (IP) for the Weibull parameters, respectively, where in NIP, the hyperparameters are taken to bea1=b1=a2=b2=0 and the hyperparameters of IP are chosen as a1=1=b1=1,a2=3,andb2=5 for α,β=1,0.6 and a1=3,b1=2,anda2=b2=1 for α,β=1.5,1 satisfying that the prior means are as the same as the original mean. The simulation was repeated 10,000 times, and the significance level for interval estimates was taken to be γ=0.05.Before starting the numerous simulations, an important process is to investigate the convergence of Bayesian estimation. For parametersα,β=1,0.6, we take the CS-1 and CS-5 with k=10and35 under T=0.7 as a group of demos. A random sample of GPHC is further generated once under each censoring scheme for use in Bayes manner. For space-saving, here we only discuss the Bayes estimates based on IP. Based on M-H algorithm with gamma proposal, the corresponding iterative trend plots are drawn in Figure 2. In a similar way, for parameter α,β=1.5,1, we select the CS-2 and CS-6 with k=13and30 under T=0.9 as another set of examples. The associated iterative trend plots are also given in Figure 3. It can be seen from Figures 2 and 3 that the mean value of the posterior sample under different parameters and censoring schemes has become stable after N>5000, which means that the sampling procedure has a good convergence. Hence, it is noted that the convergence of the MCMC sampling in Bayesian computations is guaranteed in our work.Figure 2
The sampling iterative trend plots of parameterα,β=1,0.6 under different censoring schemes.Figure 3
The sampling iterative trend plots of parameterα,β=1.5,1 under different censoring schemes.Furthermore, based on the previous settings and analysis results, the Monte Carlo simulations are carried out and all the numerical results are summarized in Tables1 to 4.Table 1
ABs and MSEs (within brackets) of the MLEs and Bayes estimates under generalized progressively hybrid censored data withα,β=1,0.6.
T=0.7MLENIP BayesIP BayesCSkαβαβαβ1100.319 2[0.212 3]0.205 1[0.086 1]0.285 8[0.163 5]0.197 1[0.079 5]0.237 5[0.104 4]0.131 3[0.028 4]130.249 7[0.123 9]0.164 6[0.050 7]0.230 6[0.102 7]0.158 1[0.047 2]0.201 6[0.074 5]0.119 1[0.023 9]2100.300 9[0.188 0]0.185 1[0.066 9]0.270 7[0.145 6]0.180 4[0.063 1]0.227 4[0.095 5]0.121 6[0.024 5]130.260 3[0.135 2]0.170 3[0.060 8]0.239 9[0.111 5]0.165 8[0.057 9]0.209 8[0.080 7]0.121 0[0.024 9]3150.232 6[0.102 3]0.148 4[0.042 8]0.216 7[0.087 1]0.144 1[0.040 6]0.194 7[0.067 9]0.112 5[0.021 8]180.199 2[0.074 5]0.134 9[0.032 0]0.188 3[0.065 4]0.130 9[0.030 3]0.172 4[0.053 1]0.107 1[0.019 3]4300.145 3[0.037 0]0.097 4[0.016 2]0.140 7[0.034 3]0.096 0[0.015 8]0.135 1[0.031 3]0.085 8[0.012 2]350.132 5[0.029 8]0.091 4[0.013 9]0.128 9[0.027 9]0.090 0[0.013 5]0.123 9[0.025 6]0.081 5[0.010 9]5300.137 1[0.032 0]0.099 6[0.016 5]0.133 3[0.030 8]0.097 8[0.016 0]0.127 4[0.027 2]0.086 9[0.012 4]350.120 5[0.024 7]0.096 1[0.015 0]0.117 8[0.023 4]0.094 5[0.014 5]0.112 4[0.021 2]0.084 5[0.011 5]6300.141 1[0.034 4]0.098 0[0.016 3]0.136 9[0.032 1]0.096 4[0.015 8]0.131 2[0.029 2]0.086 1[0.012 4]350.124 7[0.026 2]0.092 5[0.013 8]0.121 5[0.024 7]0.091 1[0.013 4]0.116 5[0.022 5]0.082 1[0.010 8]T=1.2MLENIP BayesIP BayesCSkαβαβαβ1100.278 0[0.156 8]0.159 1[0.044 4]0.252 5[0.123 7]0.153 5[0.041 6]0.219 9[0.088 0]0.115 3[0.022 3]130.246 5[0.119 3]0.161 3[0.045 8]0.227 4[0.099 1]0.154 7[0.042 6]0.198 6[0.072 1]0.117 9[0.023 2]2100.242 0[0.111 8]0.154 1[0.047 4]0.225 1[0.093 5]0.151 9[0.045 7]0.202 9[0.072 8]0.116 6[0.023 4]130.238 2[0.109 5]0.146 3[0.046 8]0.221 5[0.092 3]0.143 0[0.044 8]0.198 8[0.071 0]0.110 2[0.022 0]3150.221 3[0.094 0]0.134 9[0.029 9]0.207 1[0.080 5]0.131 1[0.028 4]0.188 8[0.064 5]0.106 5[0.018 2]180.196 5[0.071 9]0.131 6[0.029 3]0.185 9[0.063 2]0.127 7[0.027 7]0.170 0[0.051 5]0.105 1[0.018 2]4300.144 9[0.036 5]0.092 5[0.013 8]0.140 2[0.033 8]0.091 2[0.013 4]0.134 8[0.031 0]0.082 0[0.010 8]350.132 6[0.030 6]0.090 8[0.013 5]0.129 0[0.028 7]0.089 4[0.013 2]0.124 1[0.026 4]0.081 0[0.010 7]5300.136 9[0.032 0]0.098 5[0.016 1]0.133 3[0.030 1]0.096 8[0.015 6]0.127 2[0.027 2]0.086 0[0.012 1]350.119 6[0.024 1]0.095 9[0.014 7]0.116 8[0.022 8]0.094 3[0.014 3]0.111 3[0.020 5]0.084 3[0.011 3]6300.143 1[0.035 3]0.097 2[0.015 7]0.138 7[0.032 8]0.095 6[0.015 2]0.132 8[0.029 8]0.085 5[0.012 0]350.126 1[0.027 3]0.093 3[0.014 2]0.122 8[0.025 7]0.091 8[0.013 8]0.117 7[0.023 3]0.082 7[0.011 1]Table 2
AWs and CPs (within brackets) of the ACIs and HPD credible intervals under generalized progressively hybrid censored data withα,β=1,0.6.
T=0.4ACINIP HPDIP HPDCSkαβαβαβ1101.370 7[0.961 2]0.893 1[0.940 2]1.265 2[0.954 5]0.848 9[0.945 2]1.155 3[0.966 2]0.688 2[0.975 4]131.099 4[0.955 4]0.740 5[0.929 3]1.036 1[0.948 5]0.715 7[0.936 5]0.968 8[0.958 0]0.625 3[0.968 0]2101.349 4[0.963 5]0.901 8[0.953 5]1.244 1[0.954 8]0.853 5[0.947 4]1.135 6[0.964 8]0.687 2[0.969 7]131.152 5[0.955 1]0.757 7[0.946 0]1.080 2[0.949 3]0.728 4[0.942 2]1.013 4[0.958 5]0.619 5[0.965 8]3151.039 7[0.959 0]0.673 1[0.944 3]0.982 7[0.952 1]0.652 0[0.943 2]0.933 7[0.960 8]0.577 9[0.967 4]180.898 4[0.951 3]0.616 2[0.932 6]0.856 2[0.946 3]0.600 6[0.937 5]0.819 1[0.954 6]0.545 5[0.961 6]4300.686 8[0.952 5]0.453 4[0.946 0]0.659 9[0.945 1]0.445 2[0.943 1]0.646 9[0.949 6]0.421 5[0.956 3]350.617 2[0.952 2]0.428 7[0.937 5]0.595 0[0.942 0]0.421 9[0.940 6]0.583 2[0.945 9]0.402 0[0.951 3]5300.640 9[0.954 0]0.470 0[0.939 5]0.616 8[0.944 6]0.460 9[0.942 7]0.602 4[0.949 0]0.435 2[0.954 8]350.566 7[0.952 0]0.458 3[0.939 8]0.546 1[0.942 9]0.449 9[0.941 2]0.533 0[0.946 8]0.425 9[0.954 8]6300.667 4[0.955 3]0.458 2[0.937 5]0.641 4[0.948 2]0.449 7[0.939 1]0.627 9[0.949 5]0.425 6[0.951 8]350.593 7[0.953 1]0.440 4[0.941 4]0.571 7[0.944 2]0.432 9[0.944 3]0.559 9[0.950 3]0.411 5[0.956 1]T=1.2ACINIP HPDIP HPDCSkαβαβαβ1101.269 7[0.960 4]0.777 4[0.930 2]1.179 0[0.954 9]0.746 3[0.935 1]1.101 0[0.963 6]0.648 5[0.964 9]131.096 7[0.959 1]0.730 9[0.922 8]1.033 7[0.951 2]0.706 9[0.932 7]0.967 7[0.961 4]0.621 4[0.964 1]2101.125 2[0.953 2]0.702 1[0.955 3]1.054 0[0.946 6]0.676 8[0.942 1]1.001 9[0.956 1]0.596 0[0.970 2]131.092 9[0.958 9]0.691 8[0.944 8]1.028 4[0.950 5]0.668 4[0.940 0]0.976 8[0.958 5]0.588 5[0.965 0]3151.012 6[0.957 6]0.644 1[0.939 9]0.958 3[0.950 7]0.625 2[0.942 4]0.915 6[0.958 4]0.564 6[0.964 4]180.897 2[0.954 8]0.611 1[0.936 5]0.855 3[0.950 6]0.595 7[0.941 6]0.818 2[0.955 5]0.542 4[0.964 0]4300.684 9[0.954 6]0.448 2[0.949 2]0.658 4[0.945 8]0.440 3[0.950 6]0.645 5[0.951 7]0.418 0[0.960 4]350.617 6[0.945 6]0.428 7[0.940 9]0.595 3[0.935 4]0.422 1[0.943 3]0.583 3[0.942 7]0.402 1[0.955 6]5300.640 8[0.952 3]0.469 6[0.942 7]0.616 7[0.943 9]0.460 5[0.942 8]0.602 4[0.949 1]0.434 8[0.956 9]350.565 7[0.952 2]0.459 0[0.940 7]0.545 3[0.944 0]0.450 6[0.942 2]0.531 8[0.948 1]0.426 5[0.956 5]6300.668 8[0.953]00.456 7[0.938 9]0.643 2[0.945 9]0.448 3[0.937 9]0.629 8[0.950 2]0.424 6[0.951 5]350.595 2[0.949 2]0.440 3[0.937 2]0.573 1[0.939 4]0.432 9[0.938 0]0.561 3[0.945 9]0.401 6[0.950 7]Table 3
ABs and MSEs (within brackets) of the MLEs and Bayes estimates under generalized progressively hybrid censored data withα,β=1.5,1.
T=0.4MLENIP BayesIP BayesCSkαβαβαβ1100.496 3[0.512 5]0.563 6[0.631 8]0.414 6[0.335 3]0.528 9[0.573 6]0.286 2[0.138 6]0.296 0[0.170 9]130.383 4[0.286 1]0.307 0[0.285 9]0.337 3[0.212 6]0.292 1[0.260 5]0.267 7[0.123 3]0.241 0[0.114 4]2100.487 5[0.495 7]0.651 7[0.755 4]0.405 6[0.320 2]0.602 3[0.581 7]0.263 9[0.115 8]0.303 1[0.167 3]130.395 7[0.311 4]0.390 0[0.573 5]0.345 9[0.227 1]0.374 2[0.526 8]0.260 8[0.116 4]0.262 1[0.134 1]3150.347 2[0.235 2]0.297 6[0.232 9]0.309 9[0.180 6]0.285 7[0.211 1]0.249 4[0.108 0]0.233 0[0.105 6]180.295 1[0.162 3]0.229 5[0.115 1]0.271 1[0.132 4]0.221 9[0.106 7]0.231 9[0.092 5]0.199 9[0.076 2]4300.219 5[0.083 5]0.179 0[0.064 7]0.208 2[0.073 7]0.175 7[0.061 9]0.188 7[0.059 3]0.163 6[0.050 2]350.197 8[0.067 1]0.151 2[0.040 9]0.189 1[0.060 4]0.148 7[0.039 4]0.175 3[0.051 2]0.142 6[0.035 4]5300.202 8[0.069 6]0.163 9[0.048 8]0.193 8[0.062 6]0.160 8[0.046 7]0.178 9[0.052 6]0.153 6[0.041 5]350.181 4[0.055 5]0.144 0[0.035 2]0.175 1[0.050 9]0.141 5[0.033 9]0.164 3[0.044 4]0.136 8[0.031 5]6300.211 1[0.076 6]0.168 3[0.054 0]0.200 9[0.068 3]0.165 2[0.051 7]0.183 5[0.056 0]0.156 2[0.044 1]350.190 4[0.062 3]0.149 1[0.039 0]0.183 1[0.056 6]0.146 5[0.037 5]0.170 7[0.048 6]0.141 3[0.034 3]T=0.9MLENIP BayesIP BayesCSkαβαβαβ1100.388 4[0.295 6]0.269 8[0.156 3]0.338 0[0.209 9]0.259 1[0.140 6]0.262 6[0.115 8]0.223 4[0.092 1]130.359 0[0.247 4]0.275 3[0.156 0]0.318 1[0.186 1]0.262 2[0.140 6]0.257 0[0.113 3]0.229 3[0.095 4]2100.351 7[0.237 6]0.320 5[0.273 7]0.316 0[0.182 5]0.311 4[0.250 9]0.251 0[0.106 1]0.253 2[0.117 7]130.357 7[0.242 0]0.315 4[0.285 2]0.319 8[0.185 2]0.305 7[0.260 6]0.254 4[0.108 7]0.247 1[0.114 9]3150.312 5[0.181 1]0.224 7[0.096 1]0.284 3[0.143 8]0.217 1[0.089 4]0.238 1[0.096 0]0.195 9[0.068 9]180.288 1[0.151 3]0.224 7[0.097 6]0.263 9[0.123 5]0.217 0[0.090 4]0.227 0[0.087 8]0.198 3[0.070 6]4300.209 9[0.076 2]0.150 4[0.038 0]0.200 2[0.068 0]0.147 9[0.036 7]0.183 6[0.056 3]0.141 3[0.033 1]350.196 2[0.065 4]0.147 2[0.037 4]0.187 7[0.058 9]0.144 6[0.036 0]0.174 1[0.050 1]0.139 0[0.032 8]5300.204 4[0.072 4]0.162 4[0.047 4]0.194 7[0.064 8]0.159 3[0.045 4]0.179 9[0.054 5]0.152 2[0.040 6]350.179 5[0.054 7]0.146 8[0.037 5]0.172 9[0.050 1]0.144 2[0.036 1]0.162 4[0.043 8]0.139 2[0.033 4]6300.209 7[0.075 8]0.163 1[0.045 2]0.200 0[0.067 7]0.160 1[0.043 3]0.183 5[0.056 1]0.152 6[0.038 8]350.189 9[0.060 8]0.146 6[0.037 7]0.182 5[0.055 3]0.144 0[0.036 2]0.170 4[0.047 6]0.138 9[0.033 2]Table 4
AWs and CPs (within brackets) of the ACIs and HPD credible intervals under generalized progressively hybrid censored data withβ,α=1.5,1.
T=0.4ACINIP HPDIP HPDCSkαβαβαβ1102.075 9[0.958 9]2.593 2[0.961 0]1.870 0[0.957 4]2.175 4[0.960 5]1.551 9[0.978 9]1.448 8[0.966 4]131.655 2[0.953 3]1.347 5[0.955 3]1.536 5[0.950 2]1.276 1[0.951 2]1.370 6[0.969 0]1.150 3[0.961 6]2102.119 9[0.964 0]3.294 6[0.968 1]1.898 1[0.960 6]2.796 3[0.966 3]1.521 0[0.981 6]1.613 3[0.971 9]131.741 2[0.954 7]1.768 6[0.968 9]1.605 3[0.952 5]1.627 4[0.966 2]1.392 8[0.972 0]1.297 7[0.970 6]3151.555 4[0.956 9]1.343 0[0.964 2]1.451 2[0.954 0]1.268 4[0.957 1]1.303 1[0.970 5]1.128 5[0.964 7]181.339 6[0.953 6]1.049 6[0.959 0]1.266 2[0.948 1]1.012 3[0.949 4]1.172 7[0.964 4]0.961 1[0.956 3]4301.031 2[0.954 4]0.828 3[0.956 1]0.987 4[0.948 0]0.806 5[0.947 5]0.940 2[0.957 6]0.777 6[0.951 4]350.926 3[0.952 7]0.708 0[0.949 7]0.890 5[0.947 5]0.693 7[0.946 4]0.857 4[0.952 7]0.680 3[0.951 4]5300.959 5[0.952 6]0.763 1[0.951 0]0.920 8[0.947 7]0.746 2[0.949 3]0.884 4[0.954 8]0.730 3[0.952 7]350.847 7[0.949 8]0.691 6[0.946 9]0.816 7[0.944 1]0.679 3[0.946 2]0.790 6[0.950 5]0.668 4[0.951 6]6300.998 7[0.955 4]0.789 6[0.958 9]0.956 4[0.950 4]0.770 8[0.952 7]0.915 5[0.957 1]0.749 7[0.956 5]350.891 6[0.950 6]0.695 8[0.948 9]0.857 4[0.943 9]0.682 9[0.946 6]0.828 7[0.951 8]0.671 1[0.950 3]T=0.9ACINIP HPDIP HPDCSkαβαβαβ1102.075 9[0.958 9]2.593 2[0.961 0]1.870 0[0.957 4]2.175 4[0.960 5]1.551 9[0.978 9]1.448 8[0.966 4]131.655 2[0.953 3]1.347 5[0.955 3]1.536 5[0.950 2]1.276 1[0.951 2]1.370 6[0.969 0]1.150 3[0.961 6]2102.119 9[0.964 0]3.294 6[0.968 1]1.898 1[0.960 6]2.796 3[0.966 3]1.521 0[0.981 6]1.613 3[0.971 9]131.741 2[0.954 7]1.768 6[0.968 9]1.605 3[0.952 5]1.627 4[0.966 2]1.392 8[0.972 0]1.297 7[0.970 6]3151.555 4[0.956 9]1.343 0[0.964 2]1.451 2[0.954 0]1.268 4[0.957 1]1.303 1[0.970 5]1.128 5[0.964 7]181.339 6[0.953 6]1.049 6[0.959 0]1.266 2[0.948 1]1.012 3[0.949 4]1.172 7[0.964 4]0.961 1[0.956 3]4301.031 2[0.954 4]0.828 3[0.956 1]0.987 4[0.948 0]0.806 5[0.947 5]0.940 2[0.957 6]0.777 6[0.951 4]350.926 3[0.952 7]0.708 0[0.949 7]0.890 5[0.947 5]0.693 7[0.946 4]0.857 4[0.952 7]0.680 3[0.951 4]5300.959 5[0.952 6]0.763 1[0.951 0]0.920 8[0.947 7]0.746 2[0.949 3]0.884 4[0.954 8]0.730 3[0.952 7]350.847 7[0.949 8]0.691 6[0.946 9]0.816 7[0.944 1]0.679 3[0.946 2]0.790 6[0.950 5]0.668 4[0.951 6]6300.998 7[0.955 4]0.789 6[0.958 9]0.956 4[0.950 4]0.770 8[0.952 7]0.915 5[0.957 1]0.749 7[0.956 5]350.891 6[0.950 6]0.695 8[0.948 9]0.857 4[0.943 9]0.682 9[0.946 6]0.828 7[0.951 8]0.671 1[0.950 3]From Tables1–4, the following can be seen for point estimates:(i)
With increase in effective sample size (i.e.,n or m or k or their combinations), ABs and MSEs of both MLEs and Bayes estimates for parameters α and β decreases(ii)
For the fixed sample size and CS, ABs and MSEs of MLEs and Bayes estimates decrease as monitor timeT increases(iii)
MLEs and Bayes estimates with respect to NIP feature similar ABs and MSEs, whereas the Bayes estimates under IP have better performance than the MLEs and Bayes estimates with respect to NIP in terms of ABs and MSEsMoreover, one can also observe the following for interval estimates:(i)
AWs (CPs) of ACIs and Bayes intervals decrease (increase) as the effective sample size increases(ii)
When inspecting timeT increases, the AWs (CPs) of ACIs and Bayes credible intervals decrease (increase) generally(iii)
For the fixed sample and CS, AWs of Bayes credible intervals under NIP are slightly shorter than ACIs and the Bayes credible intervals with respect to IP feature smallest AWs among three interval estimates(iv)
Bayes credible intervals from IP have higher CPs than ACIs and Bayes credible intervals under NIP in most cases, where the latter two intervals perform similarly in terms of CPsIn conclusion, from simulation, one can observe that both classic and Bayes estimates have satisfactory behavior and that the performance of Bayes estimates with respect to IP are preferable than MLEs and Bayes results under NIP where the latter two estimates perform similarly.


### 4.2. Real-Life Illustration

In this part, a set of lifetime data is provided to show the applications of the Weibull model under censored data. From Lawless [20], the number of million revolutions before failure for 25 ball bearings in a life test is used in this illustration and the details are provided in Table 5.Table 5
The number of million revolutions before failure for ball bearings.
0.178 80.289 20.330 00.415 20.421 20.456 00.484 80.518 40.519 60.541 20.555 60.678 00.678 00.678 00.686 40.686 40.688 80.841 20.931 20.986 41.051 21.058 41.279 21.280 41.734 0Before further proceeding, we investigate whether the Weibull distribution can be employed or not to analyze this data. The MLEs ofα and β under complete data are 0.814 3 and 2.173 2, and the Kolmogorov–Smirnov distance and the corresponding p value (within bracket) are 0.14670.6772, respectively. The results indicate that the Weibull distribution provides a reasonable model for this data. Similarly, using Gini statistic, Lee [21] also concluded that the Weibull model can be fitted to this data set. Furthermore, based on the complete failure data from Table 5, censored scheme r1=⋯=r19=0,r20=5 with k=16,T=1.2 and m=20 is used in advance for the testing; then, the following GPHCS data is generated as(23)0.17880.28920.33000.41520.42120.45600.48480.51840.51960.54120.55560.67800.67800.67800.68640.68640.68880.84120.93120.9864.It is noted that the obtained GPHCS data is the common Type-II censored sample which is the most occurred failure observation in practice. Therefore, the corresponding point and interval estimates are calculated, and the results are listed in Table6 with a significance level γ=0.05. To compare the performance of different estimates, the estimated standard errors (ESEs) for point estimates and the interval lengths for interval estimates are shown in squared brackets, respectively. Since there is no information about the unknown parameters, NIP with a1=b1=a2=b2=0 is used in this illustrative example. Based on the results in Table 6, one could observe that although the MLEs and Bayes estimates of the parameters α and β are close to each other, the ESE of Bayes results is smaller than that of MLE. In addition, the Bayes credible intervals under NIP also have better performance than ACIs in terms of interval length. Besides, to illustrate Theorem 2 and Lemma 1, the curves of the profile log-likelihood function (10) and the logarithmic function of posterior marginal density (21) with respect to NIP are plotted in Figure 4 under the censored failure data. It is observed that the profile log-likelihood function of α has a unique maximum point and the posterior marginal density function πα|data is log-concave. Finally, the trace plots and iterative trend plots of Bayes sampling algorithm are given in Figure 5. One could also note that the MCMC sampling in Bayesian computations converges in the real-life illustration.Table 6
Point and interval estimates ofα and β under real application data.
αβMLE2.423 3[0.454 6]1.777 2[0.431 0]Bayes2.372 3[0.440 1]1.751 9[0.428 2]ACI(1.518 7, 3.289 7)[1.781 9](0.932 4, 2.622 1)[1.689 7]HPD(1.512 6, 3.217 7)[1.705 1](0.941 7, 2.558 1)[1.616 4]Figure 4
The profile log-likelihood (a) and the posterior marginal density function (b) ofα under ball-bearings failure data.
(a)(b)Figure 5
The iterative trend plots (a) and the trace plots (b) of sampling algorithm for parameterα and β under ball-bearings failure data.
(a)(b)


## 4.1. Simulation Studies

In this section, extensive simulations are carried out for comparing different proposed methods and the performance of the point estimates (i.e., MLEs and Bayes) is evaluated in terms of mean squared error (MSE) and average bias (AB), whereas the confidence interval estimates (i.e., ACIs and HPDs) are compared by criteria coverage probability (CP) and average width (AW), respectively.The evaluaton criteria quantities for point and interval estimates were examined under different choices ofβ1,β2,α,m,n, CS, r1,r2,…,rm, and T values, and the Bayes sampling algorithm was repeated 10,000 times following a burn-in phase of 5,000 iterations. For generating the generalized progressively hybrid censored Weibull samples, it is conducted as the following procedure. Using the algorithm proposed by Balakrishnan and Aggarwala ([19], p. 32), progressively Type-II censored data is generated from Weibull population with parameters α and β. Then, the generalized progressively hybrid censored sample can be obtained by comparing prefixed T and Xk:m:n with those of the generated progressively Type-II censored data.In the simulation procedure, following CSs are considered:CS-1:n=18,m=15,r1=⋯=r14=0, and r15=3CS-2:n=25,m=15,r1=⋯=r14=0, and r15=10CS-3:n=25,m=20,r1=⋯=r19=0, and r20=5CS-4:n=50,m=40,r1=⋯=r39=0, and r40=10CS-5:n=50,m=40,r2=⋯=r40=0, and r1=10CS-6:n=50,m=40,r2=⋯=r39=0, and r1=r40=5In the simulation procedure, two sets of parameter values are chosen and the Bayes estimates are computed under noninformative prior (NIP) and informative prior (IP) for the Weibull parameters, respectively, where in NIP, the hyperparameters are taken to bea1=b1=a2=b2=0 and the hyperparameters of IP are chosen as a1=1=b1=1,a2=3,andb2=5 for α,β=1,0.6 and a1=3,b1=2,anda2=b2=1 for α,β=1.5,1 satisfying that the prior means are as the same as the original mean. The simulation was repeated 10,000 times, and the significance level for interval estimates was taken to be γ=0.05.Before starting the numerous simulations, an important process is to investigate the convergence of Bayesian estimation. For parametersα,β=1,0.6, we take the CS-1 and CS-5 with k=10and35 under T=0.7 as a group of demos. A random sample of GPHC is further generated once under each censoring scheme for use in Bayes manner. For space-saving, here we only discuss the Bayes estimates based on IP. Based on M-H algorithm with gamma proposal, the corresponding iterative trend plots are drawn in Figure 2. In a similar way, for parameter α,β=1.5,1, we select the CS-2 and CS-6 with k=13and30 under T=0.9 as another set of examples. The associated iterative trend plots are also given in Figure 3. It can be seen from Figures 2 and 3 that the mean value of the posterior sample under different parameters and censoring schemes has become stable after N>5000, which means that the sampling procedure has a good convergence. Hence, it is noted that the convergence of the MCMC sampling in Bayesian computations is guaranteed in our work.Figure 2
The sampling iterative trend plots of parameterα,β=1,0.6 under different censoring schemes.Figure 3
The sampling iterative trend plots of parameterα,β=1.5,1 under different censoring schemes.Furthermore, based on the previous settings and analysis results, the Monte Carlo simulations are carried out and all the numerical results are summarized in Tables1 to 4.Table 1
ABs and MSEs (within brackets) of the MLEs and Bayes estimates under generalized progressively hybrid censored data withα,β=1,0.6.
T=0.7MLENIP BayesIP BayesCSkαβαβαβ1100.319 2[0.212 3]0.205 1[0.086 1]0.285 8[0.163 5]0.197 1[0.079 5]0.237 5[0.104 4]0.131 3[0.028 4]130.249 7[0.123 9]0.164 6[0.050 7]0.230 6[0.102 7]0.158 1[0.047 2]0.201 6[0.074 5]0.119 1[0.023 9]2100.300 9[0.188 0]0.185 1[0.066 9]0.270 7[0.145 6]0.180 4[0.063 1]0.227 4[0.095 5]0.121 6[0.024 5]130.260 3[0.135 2]0.170 3[0.060 8]0.239 9[0.111 5]0.165 8[0.057 9]0.209 8[0.080 7]0.121 0[0.024 9]3150.232 6[0.102 3]0.148 4[0.042 8]0.216 7[0.087 1]0.144 1[0.040 6]0.194 7[0.067 9]0.112 5[0.021 8]180.199 2[0.074 5]0.134 9[0.032 0]0.188 3[0.065 4]0.130 9[0.030 3]0.172 4[0.053 1]0.107 1[0.019 3]4300.145 3[0.037 0]0.097 4[0.016 2]0.140 7[0.034 3]0.096 0[0.015 8]0.135 1[0.031 3]0.085 8[0.012 2]350.132 5[0.029 8]0.091 4[0.013 9]0.128 9[0.027 9]0.090 0[0.013 5]0.123 9[0.025 6]0.081 5[0.010 9]5300.137 1[0.032 0]0.099 6[0.016 5]0.133 3[0.030 8]0.097 8[0.016 0]0.127 4[0.027 2]0.086 9[0.012 4]350.120 5[0.024 7]0.096 1[0.015 0]0.117 8[0.023 4]0.094 5[0.014 5]0.112 4[0.021 2]0.084 5[0.011 5]6300.141 1[0.034 4]0.098 0[0.016 3]0.136 9[0.032 1]0.096 4[0.015 8]0.131 2[0.029 2]0.086 1[0.012 4]350.124 7[0.026 2]0.092 5[0.013 8]0.121 5[0.024 7]0.091 1[0.013 4]0.116 5[0.022 5]0.082 1[0.010 8]T=1.2MLENIP BayesIP BayesCSkαβαβαβ1100.278 0[0.156 8]0.159 1[0.044 4]0.252 5[0.123 7]0.153 5[0.041 6]0.219 9[0.088 0]0.115 3[0.022 3]130.246 5[0.119 3]0.161 3[0.045 8]0.227 4[0.099 1]0.154 7[0.042 6]0.198 6[0.072 1]0.117 9[0.023 2]2100.242 0[0.111 8]0.154 1[0.047 4]0.225 1[0.093 5]0.151 9[0.045 7]0.202 9[0.072 8]0.116 6[0.023 4]130.238 2[0.109 5]0.146 3[0.046 8]0.221 5[0.092 3]0.143 0[0.044 8]0.198 8[0.071 0]0.110 2[0.022 0]3150.221 3[0.094 0]0.134 9[0.029 9]0.207 1[0.080 5]0.131 1[0.028 4]0.188 8[0.064 5]0.106 5[0.018 2]180.196 5[0.071 9]0.131 6[0.029 3]0.185 9[0.063 2]0.127 7[0.027 7]0.170 0[0.051 5]0.105 1[0.018 2]4300.144 9[0.036 5]0.092 5[0.013 8]0.140 2[0.033 8]0.091 2[0.013 4]0.134 8[0.031 0]0.082 0[0.010 8]350.132 6[0.030 6]0.090 8[0.013 5]0.129 0[0.028 7]0.089 4[0.013 2]0.124 1[0.026 4]0.081 0[0.010 7]5300.136 9[0.032 0]0.098 5[0.016 1]0.133 3[0.030 1]0.096 8[0.015 6]0.127 2[0.027 2]0.086 0[0.012 1]350.119 6[0.024 1]0.095 9[0.014 7]0.116 8[0.022 8]0.094 3[0.014 3]0.111 3[0.020 5]0.084 3[0.011 3]6300.143 1[0.035 3]0.097 2[0.015 7]0.138 7[0.032 8]0.095 6[0.015 2]0.132 8[0.029 8]0.085 5[0.012 0]350.126 1[0.027 3]0.093 3[0.014 2]0.122 8[0.025 7]0.091 8[0.013 8]0.117 7[0.023 3]0.082 7[0.011 1]Table 2
AWs and CPs (within brackets) of the ACIs and HPD credible intervals under generalized progressively hybrid censored data withα,β=1,0.6.
T=0.4ACINIP HPDIP HPDCSkαβαβαβ1101.370 7[0.961 2]0.893 1[0.940 2]1.265 2[0.954 5]0.848 9[0.945 2]1.155 3[0.966 2]0.688 2[0.975 4]131.099 4[0.955 4]0.740 5[0.929 3]1.036 1[0.948 5]0.715 7[0.936 5]0.968 8[0.958 0]0.625 3[0.968 0]2101.349 4[0.963 5]0.901 8[0.953 5]1.244 1[0.954 8]0.853 5[0.947 4]1.135 6[0.964 8]0.687 2[0.969 7]131.152 5[0.955 1]0.757 7[0.946 0]1.080 2[0.949 3]0.728 4[0.942 2]1.013 4[0.958 5]0.619 5[0.965 8]3151.039 7[0.959 0]0.673 1[0.944 3]0.982 7[0.952 1]0.652 0[0.943 2]0.933 7[0.960 8]0.577 9[0.967 4]180.898 4[0.951 3]0.616 2[0.932 6]0.856 2[0.946 3]0.600 6[0.937 5]0.819 1[0.954 6]0.545 5[0.961 6]4300.686 8[0.952 5]0.453 4[0.946 0]0.659 9[0.945 1]0.445 2[0.943 1]0.646 9[0.949 6]0.421 5[0.956 3]350.617 2[0.952 2]0.428 7[0.937 5]0.595 0[0.942 0]0.421 9[0.940 6]0.583 2[0.945 9]0.402 0[0.951 3]5300.640 9[0.954 0]0.470 0[0.939 5]0.616 8[0.944 6]0.460 9[0.942 7]0.602 4[0.949 0]0.435 2[0.954 8]350.566 7[0.952 0]0.458 3[0.939 8]0.546 1[0.942 9]0.449 9[0.941 2]0.533 0[0.946 8]0.425 9[0.954 8]6300.667 4[0.955 3]0.458 2[0.937 5]0.641 4[0.948 2]0.449 7[0.939 1]0.627 9[0.949 5]0.425 6[0.951 8]350.593 7[0.953 1]0.440 4[0.941 4]0.571 7[0.944 2]0.432 9[0.944 3]0.559 9[0.950 3]0.411 5[0.956 1]T=1.2ACINIP HPDIP HPDCSkαβαβαβ1101.269 7[0.960 4]0.777 4[0.930 2]1.179 0[0.954 9]0.746 3[0.935 1]1.101 0[0.963 6]0.648 5[0.964 9]131.096 7[0.959 1]0.730 9[0.922 8]1.033 7[0.951 2]0.706 9[0.932 7]0.967 7[0.961 4]0.621 4[0.964 1]2101.125 2[0.953 2]0.702 1[0.955 3]1.054 0[0.946 6]0.676 8[0.942 1]1.001 9[0.956 1]0.596 0[0.970 2]131.092 9[0.958 9]0.691 8[0.944 8]1.028 4[0.950 5]0.668 4[0.940 0]0.976 8[0.958 5]0.588 5[0.965 0]3151.012 6[0.957 6]0.644 1[0.939 9]0.958 3[0.950 7]0.625 2[0.942 4]0.915 6[0.958 4]0.564 6[0.964 4]180.897 2[0.954 8]0.611 1[0.936 5]0.855 3[0.950 6]0.595 7[0.941 6]0.818 2[0.955 5]0.542 4[0.964 0]4300.684 9[0.954 6]0.448 2[0.949 2]0.658 4[0.945 8]0.440 3[0.950 6]0.645 5[0.951 7]0.418 0[0.960 4]350.617 6[0.945 6]0.428 7[0.940 9]0.595 3[0.935 4]0.422 1[0.943 3]0.583 3[0.942 7]0.402 1[0.955 6]5300.640 8[0.952 3]0.469 6[0.942 7]0.616 7[0.943 9]0.460 5[0.942 8]0.602 4[0.949 1]0.434 8[0.956 9]350.565 7[0.952 2]0.459 0[0.940 7]0.545 3[0.944 0]0.450 6[0.942 2]0.531 8[0.948 1]0.426 5[0.956 5]6300.668 8[0.953]00.456 7[0.938 9]0.643 2[0.945 9]0.448 3[0.937 9]0.629 8[0.950 2]0.424 6[0.951 5]350.595 2[0.949 2]0.440 3[0.937 2]0.573 1[0.939 4]0.432 9[0.938 0]0.561 3[0.945 9]0.401 6[0.950 7]Table 3
ABs and MSEs (within brackets) of the MLEs and Bayes estimates under generalized progressively hybrid censored data withα,β=1.5,1.
T=0.4MLENIP BayesIP BayesCSkαβαβαβ1100.496 3[0.512 5]0.563 6[0.631 8]0.414 6[0.335 3]0.528 9[0.573 6]0.286 2[0.138 6]0.296 0[0.170 9]130.383 4[0.286 1]0.307 0[0.285 9]0.337 3[0.212 6]0.292 1[0.260 5]0.267 7[0.123 3]0.241 0[0.114 4]2100.487 5[0.495 7]0.651 7[0.755 4]0.405 6[0.320 2]0.602 3[0.581 7]0.263 9[0.115 8]0.303 1[0.167 3]130.395 7[0.311 4]0.390 0[0.573 5]0.345 9[0.227 1]0.374 2[0.526 8]0.260 8[0.116 4]0.262 1[0.134 1]3150.347 2[0.235 2]0.297 6[0.232 9]0.309 9[0.180 6]0.285 7[0.211 1]0.249 4[0.108 0]0.233 0[0.105 6]180.295 1[0.162 3]0.229 5[0.115 1]0.271 1[0.132 4]0.221 9[0.106 7]0.231 9[0.092 5]0.199 9[0.076 2]4300.219 5[0.083 5]0.179 0[0.064 7]0.208 2[0.073 7]0.175 7[0.061 9]0.188 7[0.059 3]0.163 6[0.050 2]350.197 8[0.067 1]0.151 2[0.040 9]0.189 1[0.060 4]0.148 7[0.039 4]0.175 3[0.051 2]0.142 6[0.035 4]5300.202 8[0.069 6]0.163 9[0.048 8]0.193 8[0.062 6]0.160 8[0.046 7]0.178 9[0.052 6]0.153 6[0.041 5]350.181 4[0.055 5]0.144 0[0.035 2]0.175 1[0.050 9]0.141 5[0.033 9]0.164 3[0.044 4]0.136 8[0.031 5]6300.211 1[0.076 6]0.168 3[0.054 0]0.200 9[0.068 3]0.165 2[0.051 7]0.183 5[0.056 0]0.156 2[0.044 1]350.190 4[0.062 3]0.149 1[0.039 0]0.183 1[0.056 6]0.146 5[0.037 5]0.170 7[0.048 6]0.141 3[0.034 3]T=0.9MLENIP BayesIP BayesCSkαβαβαβ1100.388 4[0.295 6]0.269 8[0.156 3]0.338 0[0.209 9]0.259 1[0.140 6]0.262 6[0.115 8]0.223 4[0.092 1]130.359 0[0.247 4]0.275 3[0.156 0]0.318 1[0.186 1]0.262 2[0.140 6]0.257 0[0.113 3]0.229 3[0.095 4]2100.351 7[0.237 6]0.320 5[0.273 7]0.316 0[0.182 5]0.311 4[0.250 9]0.251 0[0.106 1]0.253 2[0.117 7]130.357 7[0.242 0]0.315 4[0.285 2]0.319 8[0.185 2]0.305 7[0.260 6]0.254 4[0.108 7]0.247 1[0.114 9]3150.312 5[0.181 1]0.224 7[0.096 1]0.284 3[0.143 8]0.217 1[0.089 4]0.238 1[0.096 0]0.195 9[0.068 9]180.288 1[0.151 3]0.224 7[0.097 6]0.263 9[0.123 5]0.217 0[0.090 4]0.227 0[0.087 8]0.198 3[0.070 6]4300.209 9[0.076 2]0.150 4[0.038 0]0.200 2[0.068 0]0.147 9[0.036 7]0.183 6[0.056 3]0.141 3[0.033 1]350.196 2[0.065 4]0.147 2[0.037 4]0.187 7[0.058 9]0.144 6[0.036 0]0.174 1[0.050 1]0.139 0[0.032 8]5300.204 4[0.072 4]0.162 4[0.047 4]0.194 7[0.064 8]0.159 3[0.045 4]0.179 9[0.054 5]0.152 2[0.040 6]350.179 5[0.054 7]0.146 8[0.037 5]0.172 9[0.050 1]0.144 2[0.036 1]0.162 4[0.043 8]0.139 2[0.033 4]6300.209 7[0.075 8]0.163 1[0.045 2]0.200 0[0.067 7]0.160 1[0.043 3]0.183 5[0.056 1]0.152 6[0.038 8]350.189 9[0.060 8]0.146 6[0.037 7]0.182 5[0.055 3]0.144 0[0.036 2]0.170 4[0.047 6]0.138 9[0.033 2]Table 4
AWs and CPs (within brackets) of the ACIs and HPD credible intervals under generalized progressively hybrid censored data withβ,α=1.5,1.
T=0.4ACINIP HPDIP HPDCSkαβαβαβ1102.075 9[0.958 9]2.593 2[0.961 0]1.870 0[0.957 4]2.175 4[0.960 5]1.551 9[0.978 9]1.448 8[0.966 4]131.655 2[0.953 3]1.347 5[0.955 3]1.536 5[0.950 2]1.276 1[0.951 2]1.370 6[0.969 0]1.150 3[0.961 6]2102.119 9[0.964 0]3.294 6[0.968 1]1.898 1[0.960 6]2.796 3[0.966 3]1.521 0[0.981 6]1.613 3[0.971 9]131.741 2[0.954 7]1.768 6[0.968 9]1.605 3[0.952 5]1.627 4[0.966 2]1.392 8[0.972 0]1.297 7[0.970 6]3151.555 4[0.956 9]1.343 0[0.964 2]1.451 2[0.954 0]1.268 4[0.957 1]1.303 1[0.970 5]1.128 5[0.964 7]181.339 6[0.953 6]1.049 6[0.959 0]1.266 2[0.948 1]1.012 3[0.949 4]1.172 7[0.964 4]0.961 1[0.956 3]4301.031 2[0.954 4]0.828 3[0.956 1]0.987 4[0.948 0]0.806 5[0.947 5]0.940 2[0.957 6]0.777 6[0.951 4]350.926 3[0.952 7]0.708 0[0.949 7]0.890 5[0.947 5]0.693 7[0.946 4]0.857 4[0.952 7]0.680 3[0.951 4]5300.959 5[0.952 6]0.763 1[0.951 0]0.920 8[0.947 7]0.746 2[0.949 3]0.884 4[0.954 8]0.730 3[0.952 7]350.847 7[0.949 8]0.691 6[0.946 9]0.816 7[0.944 1]0.679 3[0.946 2]0.790 6[0.950 5]0.668 4[0.951 6]6300.998 7[0.955 4]0.789 6[0.958 9]0.956 4[0.950 4]0.770 8[0.952 7]0.915 5[0.957 1]0.749 7[0.956 5]350.891 6[0.950 6]0.695 8[0.948 9]0.857 4[0.943 9]0.682 9[0.946 6]0.828 7[0.951 8]0.671 1[0.950 3]T=0.9ACINIP HPDIP HPDCSkαβαβαβ1102.075 9[0.958 9]2.593 2[0.961 0]1.870 0[0.957 4]2.175 4[0.960 5]1.551 9[0.978 9]1.448 8[0.966 4]131.655 2[0.953 3]1.347 5[0.955 3]1.536 5[0.950 2]1.276 1[0.951 2]1.370 6[0.969 0]1.150 3[0.961 6]2102.119 9[0.964 0]3.294 6[0.968 1]1.898 1[0.960 6]2.796 3[0.966 3]1.521 0[0.981 6]1.613 3[0.971 9]131.741 2[0.954 7]1.768 6[0.968 9]1.605 3[0.952 5]1.627 4[0.966 2]1.392 8[0.972 0]1.297 7[0.970 6]3151.555 4[0.956 9]1.343 0[0.964 2]1.451 2[0.954 0]1.268 4[0.957 1]1.303 1[0.970 5]1.128 5[0.964 7]181.339 6[0.953 6]1.049 6[0.959 0]1.266 2[0.948 1]1.012 3[0.949 4]1.172 7[0.964 4]0.961 1[0.956 3]4301.031 2[0.954 4]0.828 3[0.956 1]0.987 4[0.948 0]0.806 5[0.947 5]0.940 2[0.957 6]0.777 6[0.951 4]350.926 3[0.952 7]0.708 0[0.949 7]0.890 5[0.947 5]0.693 7[0.946 4]0.857 4[0.952 7]0.680 3[0.951 4]5300.959 5[0.952 6]0.763 1[0.951 0]0.920 8[0.947 7]0.746 2[0.949 3]0.884 4[0.954 8]0.730 3[0.952 7]350.847 7[0.949 8]0.691 6[0.946 9]0.816 7[0.944 1]0.679 3[0.946 2]0.790 6[0.950 5]0.668 4[0.951 6]6300.998 7[0.955 4]0.789 6[0.958 9]0.956 4[0.950 4]0.770 8[0.952 7]0.915 5[0.957 1]0.749 7[0.956 5]350.891 6[0.950 6]0.695 8[0.948 9]0.857 4[0.943 9]0.682 9[0.946 6]0.828 7[0.951 8]0.671 1[0.950 3]From Tables1–4, the following can be seen for point estimates:(i)
With increase in effective sample size (i.e.,n or m or k or their combinations), ABs and MSEs of both MLEs and Bayes estimates for parameters α and β decreases(ii)
For the fixed sample size and CS, ABs and MSEs of MLEs and Bayes estimates decrease as monitor timeT increases(iii)
MLEs and Bayes estimates with respect to NIP feature similar ABs and MSEs, whereas the Bayes estimates under IP have better performance than the MLEs and Bayes estimates with respect to NIP in terms of ABs and MSEsMoreover, one can also observe the following for interval estimates:(i)
AWs (CPs) of ACIs and Bayes intervals decrease (increase) as the effective sample size increases(ii)
When inspecting timeT increases, the AWs (CPs) of ACIs and Bayes credible intervals decrease (increase) generally(iii)
For the fixed sample and CS, AWs of Bayes credible intervals under NIP are slightly shorter than ACIs and the Bayes credible intervals with respect to IP feature smallest AWs among three interval estimates(iv)
Bayes credible intervals from IP have higher CPs than ACIs and Bayes credible intervals under NIP in most cases, where the latter two intervals perform similarly in terms of CPsIn conclusion, from simulation, one can observe that both classic and Bayes estimates have satisfactory behavior and that the performance of Bayes estimates with respect to IP are preferable than MLEs and Bayes results under NIP where the latter two estimates perform similarly.


## 4.2. Real-Life Illustration

In this part, a set of lifetime data is provided to show the applications of the Weibull model under censored data. From Lawless [20], the number of million revolutions before failure for 25 ball bearings in a life test is used in this illustration and the details are provided in Table 5.Table 5
The number of million revolutions before failure for ball bearings.
0.178 80.289 20.330 00.415 20.421 20.456 00.484 80.518 40.519 60.541 20.555 60.678 00.678 00.678 00.686 40.686 40.688 80.841 20.931 20.986 41.051 21.058 41.279 21.280 41.734 0Before further proceeding, we investigate whether the Weibull distribution can be employed or not to analyze this data. The MLEs ofα and β under complete data are 0.814 3 and 2.173 2, and the Kolmogorov–Smirnov distance and the corresponding p value (within bracket) are 0.14670.6772, respectively. The results indicate that the Weibull distribution provides a reasonable model for this data. Similarly, using Gini statistic, Lee [21] also concluded that the Weibull model can be fitted to this data set. Furthermore, based on the complete failure data from Table 5, censored scheme r1=⋯=r19=0,r20=5 with k=16,T=1.2 and m=20 is used in advance for the testing; then, the following GPHCS data is generated as(23)0.17880.28920.33000.41520.42120.45600.48480.51840.51960.54120.55560.67800.67800.67800.68640.68640.68880.84120.93120.9864.It is noted that the obtained GPHCS data is the common Type-II censored sample which is the most occurred failure observation in practice. Therefore, the corresponding point and interval estimates are calculated, and the results are listed in Table6 with a significance level γ=0.05. To compare the performance of different estimates, the estimated standard errors (ESEs) for point estimates and the interval lengths for interval estimates are shown in squared brackets, respectively. Since there is no information about the unknown parameters, NIP with a1=b1=a2=b2=0 is used in this illustrative example. Based on the results in Table 6, one could observe that although the MLEs and Bayes estimates of the parameters α and β are close to each other, the ESE of Bayes results is smaller than that of MLE. In addition, the Bayes credible intervals under NIP also have better performance than ACIs in terms of interval length. Besides, to illustrate Theorem 2 and Lemma 1, the curves of the profile log-likelihood function (10) and the logarithmic function of posterior marginal density (21) with respect to NIP are plotted in Figure 4 under the censored failure data. It is observed that the profile log-likelihood function of α has a unique maximum point and the posterior marginal density function πα|data is log-concave. Finally, the trace plots and iterative trend plots of Bayes sampling algorithm are given in Figure 5. One could also note that the MCMC sampling in Bayesian computations converges in the real-life illustration.Table 6
Point and interval estimates ofα and β under real application data.
αβMLE2.423 3[0.454 6]1.777 2[0.431 0]Bayes2.372 3[0.440 1]1.751 9[0.428 2]ACI(1.518 7, 3.289 7)[1.781 9](0.932 4, 2.622 1)[1.689 7]HPD(1.512 6, 3.217 7)[1.705 1](0.941 7, 2.558 1)[1.616 4]Figure 4
The profile log-likelihood (a) and the posterior marginal density function (b) ofα under ball-bearings failure data.
(a)(b)Figure 5
The iterative trend plots (a) and the trace plots (b) of sampling algorithm for parameterα and β under ball-bearings failure data.
(a)(b)


## 5. Concluding Remarks

In this paper, classical and Bayesian inferences are considered for generalized progressive hybrid censored data. When the latent failure times are Weibull distributed, point and interval estimations are discussed for the unknown parameters, respectively. The existence and uniqueness of MLEs are established, and Bayes estimates are obtained by using the Monte Carlo sampling method. Numerical illustration indicates that classical and Bayes estimates have similar performance when noninformative prior is adopted and that the Bayes estimators with respect to informative prior perform better than MLEs and Bayes estimators under noninformative prior. Furthermore, as a future work, optimal design and sampling plan of generalized progressive hybrid censoring under Weibull population seem also of interest which will be studied in the future.



---

*Source: 5561191-2021-12-03.xml*