Next Article in Journal
Adaptive Sliding Mode Control Anticipating Proportional Degradation of Actuator Torque in Uncertain Serial Industrial Robots
Previous Article in Journal
Chaotic Time Series Forecasting Approaches Using Machine Learning Techniques: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian Inferential Approaches and Bootstrap for the Reliability and Hazard Rate Functions under Progressive First-Failure Censoring for Coronavirus Data from Asymmetric Model

by
Rashad M. EL-Sagheer
1,
Muqrin A. Almuqrin
2,*,
Mahmoud El-Morshedy
3,4,
Mohamed S. Eliwa
5,6,
Fathy H. Eissa
7,8 and
Doaa A. Abdo
9
1
Mathematics Department, Faculty of Science, Al-Azhar University, Cairo 11884, Egypt
2
Department of Mathematics, College of Science in Zulfi, Majmaah University, Al-Majmaah 11952, Saudi Arabia
3
Department of Mathematics, College of Science and Humanities in Al-Kharj, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
4
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
5
Department of Statistics and Operation Research, College of Science, Qassim University, P.O. Box 6644, Buraydah 51482, Saudi Arabia
6
Department of Statistics and Computer Science, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
7
College of Science and Arts-Rabigh, King Abdulaziz University, Rabigh 21911, Saudi Arabia
8
Department of Mathematics, Faculty of Science, Damanhur University, Damanhur 22511, Egypt
9
Department of Applied Statistics and Insurance, Faculty of Commerce, Mansoura University, Mansoura 35511, Egypt
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(5), 956; https://doi.org/10.3390/sym14050956
Submission received: 15 March 2022 / Revised: 21 April 2022 / Accepted: 5 May 2022 / Published: 7 May 2022
(This article belongs to the Section Mathematics)

Abstract

:
This paper deals with the estimation of the parameters for asymmetric distribution and some lifetime indices such as reliability and hazard rate functions based on progressive first-failure censoring. Maximum likelihood, bootstrap and Bayesian approaches of the distribution parameters and reliability characteristics are investigated. Furthermore, the approximate confidence intervals and highest posterior density credible intervals of the parameters are constructed based on the asymptotic distribution of the maximum likelihood estimators and Markov chain Monte Carlo technique, respectively. In addition, the delta method is implemented to obtain the variances of the reliability and hazard functions. Moreover, we apply two methods of bootstrap to construct the confidence intervals. The Bayes inference based on the squared error and LINEX loss functions is obtained. Extensive simulation studies are conducted to evaluate the behavior of the proposed methods. Finally, a real data set of the COVID-19 mortality rate is analyzed to illustrate the estimation methods developed here.

1. Introduction

In life testing and reliability analysis, some items may fail (called a censored sample) or be removed before failure. Hence, censoring occurs when the exact lifetimes for only some of the items in the test are known. One of the major reasons for the removal of an experimental item is saving the working experimental items for future use, saving the cost and time associated with testing.
In practice, there are usually two random variables: the time and the number of failures of items. This strategy of censoring schemes shows how the examiner imagines the experiment based on a predetermined time. A random number of items is accounted for the first Type-I of a censoring scheme, which means the exact time of stopping an experiment may be assumed, while there is a predetermined number of failure items and a random time in the type-II censoring scheme. In these two types of censoring schemes, items cannot be removed from an experiment until the final stage, or the number of items fail. These schemes allow the detection of some items that are defective after running the experiment. The proposed approach that has more flexibility than type-II censoring in allowing items to be withdrawn from the test at different observed failure times is the progressive type-II censoring.
Balakrishnan and Sandhu [1] generated a progressive type-II censoring, which is useful to achieve the desired aims of the used censoring schemes. The general idea of this type is to remove some items from the test at the observed failure times within the duration of the test until the test is terminated. Schematically: The independent items n are placed on a life test and the observed failure times m ( m n ) , which is called a progressive sample as X 1 : m : n < X 2 : m : n < < X m : m : n , is prefixed. In addition, the scheme R = ( R 1 , R 2 , , R m ) is a prefixed censoring plan. This is achieved as follows: At the first failure occurrence X 1 , the surviving items R 1 are randomly withdrawn from the test. Then, at the second failure occurrence X 2 , the surviving items R 2 are randomly withdrawn from the test and so on until the occurrence of the m-th failure and therefore the remaining surviving items R m = n m R 1 R 2 R m 1 are withdrawn from the test and the test is terminated. Several authors have discussed inference under a progressive type-II censoring scheme with applications; see, for example, Fu et al. [2], Chen et al. [3], Xu et al. [4] and Luo et al. [5]. Recently, El-Sagheer [6] discussed the estimation problem for the Weibull–Fréchet model based on the progressive type-II censoring schemes with applications on gastric cancer survival times.
Wu et al. [7] described a life test in which the experimenter can decide to divide the items under test into several groups and then run all the items simultaneously until the occurrence of the first failure in each group. Such a censoring scheme is called first-failure censoring. However, using this censoring scheme does not able the experimenter to remove experimental groups from the test until the first failure is observed. For this reason, Wu and Kuş [8] introduced life testing, which combines first-failure censoring with progressive type-II censoring, namely as a progressive first-failure censoring (PFFC) scheme. Several authors have discussed inference under a PFFC scheme for different lifetime distributions; see, for example, Haj et al. [9], Abushal [10], Soliman et al. [11,12], Mahmoud et al. [13], Ahmed [14], Xie and Gui [15] and Shi and Shi [16]. This censoring scheme has advantages in terms of reducing test time in which more items are used but only m of n × k items are failures. Note that some censoring rules can be accommodated: the complete sample case for k = 1 , m = n , R i = 0 , and i = 1 , 2 , 3 , , m ; the first-failure censoring sample for k 1 , m = n , R i = 0 , and i = 1 , 2 , 3 , , m ; the progressive type-II censoring sample for k = 1 ; and the usual type-II censored sample case by taking k = 1 , R m = n m , and R i = 0 for i = 1 , 2 , 3 , , m 1 .
Afify et al. [17] introduced the four parameters Weibull–Fréchet distribution, WFr ( α , β , δ , λ ) , which is a generalization of both Weibull and Fréchet distributions. It is a significant continuous lifetime distribution and plays a key role in reliability problems. The probability density function (PDF) of a random variable X W F r ( α , β , δ , λ ) is given by
f x = δ λ β α β x β + 1 e λ α x β 1 e α x β λ + 1 e δ e α x β 1 λ ; x > 0 ,
and cumulative distribution function (CDF) is given by
F x = 1 e δ e α x β 1 λ ; x > 0 ,
where α is a scale parameter and β ,   δ ,   λ are shape parameters. The reliability and hazard rate functions of X are given by the following expressions
r x = e δ e α x β 1 λ ; x > 0 ,
and
h x = δ λ β α β x β + 1 e λ α x β 1 e α x β λ + 1 ; x > 0 .
This distribution is a very exible model that approaches different distributions when its parameters are changed. It contains the following special models:
  • WFr ( α , 1 , δ , λ ) follows the Weibull inverse exponential model.
  • WFr α , β , δ , 1 is the exponential Frechet model.
  • WFr α , 2 , δ , λ refers to the Weibull inverse Rayleigh model.
  • WFr α , 2 , δ , 1 reduces to the exponential inverse Rayleigh model.
  • WFr α , 1 , δ , 1 follows the exponential inverse exponential model.
The main aim of this article is to focus on the designing problem of a progressive first-failure censoring life test with a Weibull–Fréchet failure time distribution. Four methods—maximum likelihood, Bootstrap-p, Bootstrap-t and Bayes—are used to estimate four unknown parameters as well as reliability and hazard rate functions of the WFrD based on a PFFC scheme. The approximate confidence intervals (ACIs) for α ,   β ,   δ and λ are constructed based on the asymptotic normality of the MLEs. In addition, the delta method is implemented to construct the variances of the reliability r t and hazard rate h t functions. In the Bayesian framework, it is difficult to obtain the Bayes estimators in explicit forms. Thus, Markov chain Monte Carlo (MCMC) techniques are applied to compute the Bayes estimators. To this end, the Metropolis–Hastings (M-H) algorithm within the Gibbs sampler is implemented. The Bayes estimates are obtained under both symmetric (SE) and asymmetric (LINEX) loss functions. Moreover, the corresponding credible intervals (CRIs) are constructed under the MCMC technique. Finally, a simulation study is utilized to assay the performance of the proposed methods.
The rest of the paper is arranged as follows. Section 2 deals with the maximum likelihood estimate and asymptotic confidence intervals. Two parametric bootstrap methods are presented in Section 3. Bayesian estimates using the MCMC technique are provided in Section 4. In Section 5, a simulation study is conducted to compare the performance of these estimation methods. A real data set of COVID-19 is presented to illustrate the application of the proposed inference procedures in Section 6. Finally, a brief conclusion is furnished in Section 7.

2. Maximum Likelihood Inference

Suppose that x i = x i : m : n : k R , i = 1 , 2 , , m , are PFFC order statistics from the WFrD with progressive censored scheme R = R 1 , R 2 , , R m . According to Wu and Kuş [8], the log-likelihood function = log L α , β , δ , λ ; x ̲ without normalized constant can be written as
m log k + m log δ + m log λ + m log β + m β log α β + 1 i = 1 m log x i λ i = 1 m α x i β λ + 1 i = 1 m log 1 e α x i β δ i = 1 m k R i + 1 e α x i β 1 λ .
Calculating the first-order partial derivative of (5) with respect to α ,   β ,   δ and λ , respectively, and equating each result to zero, we obtain the likelihood equations
l m β α λ i = 1 m β x i α x i β 1 λ + 1 i = 1 m β x i α x i β 1 e α x i β 1 e α x i β + δ λ i = 1 m k R i + 1 e α x i β 1 λ + 1 β x i α x i β 1 e α x i β = 0 ,
l m β + m log α i = 1 m log x i λ i = 1 m α x i β log α x i λ + 1 i = 1 m α x i β e α x i β log α x i 1 e α x i β + δ λ i = 1 m k R i + 1 e α x i β 1 λ + 1 α x i β e α x i β log α x i = 0 ,
m δ i = 1 m k R i + 1 e α x i β 1 λ = 0 ,
and
m λ i = 1 m α x i β i = 1 m log 1 e α x i β δ i = 1 m k R i + 1 e α x i β 1 λ log e α x i β 1 = 0 .
From (8), we obtain MLE of δ as
δ ^ = m i = 1 m k R i + 1 e α ^ x i β ^ 1 λ ^ 1 .
As such, (6), (7) and (9) do not have closed form solutions. Newton–Raphson iteration method is widely used to obtain the desired MLEs in such situations. For more details, see EL-Sagheer [18]. The algorithm is described as follows:
  • Begin with a guess initial values α 0 , β 0 , δ 0 , λ 0 and set k = 0 .
  • Calculate α , β , δ , λ α k , β k , δ k , λ k and the observed Fisher information matrix I 1 α , β , δ , λ , given in Section 2.1.
  • Update α , β , δ , λ as
    α k + 1 , β k + 1 , δ k + 1 , λ k + 1 = α k , β k , δ k , λ k + α , β , δ , λ α k , β k , δ k , λ k × I 1 α , β , δ , λ .
  • Set k = k + 1 and then go back to Step 1.
  • Continue the iterative steps until α k + 1 , β k + 1 , δ k + 1 , λ k + 1 α k , β k , δ k , λ k is smaller than a threshold value. The final estimates of α ,   β ,   δ and λ are the MLE of the parameters, denoted as α ^ ,   β ^ ,   δ ^ and λ ^ .
Thus, the MLEs of r t and h t for given mission time t can be obtained according to the invariant property of the MLEs as
r ^ t = e δ ^ e α ^ x β ^ 1 λ ^ ,
and
h ^ t = δ ^ λ ^ β ^ α ^ β ^ x β ^ + 1 e λ ^ α ^ x β ^ 1 e α ^ x β ^ λ ^ + 1 .

2.1. Asymptotic Confidence Intervals

In this subsection, we propose the asymptotic confidence intervals (ACIs) for unknown parameters Ω = ( Ω 1 , Ω 2 , Ω 3 , Ω 4 ) = ( α , β , δ , λ ) . This requires to compute the Fisher information matrix (FIM), which is defined by the negative expectation of the second-order partial derivative of the log-likelihood function
I Ω = E 2 Ω i Ω j i , j = 1 , 2 , 3 , 4 = I i j 4 × 4 .
Because MLE has asymptotic normality property under certain regularity conditions, the estimator Ω ^ = ( α ^ , β ^ , δ ^ , λ ^ ) has asymptotic distribution Ω ^ Ω N ( 0 , I 1 Ω ) , and the inverse matrix of I Ω is
I 1 Ω = I 1 ( Ω ^ ) = v a r ( α ^ ) c o v ( α ^ , β ^ ) c o v ( α ^ , δ ^ ) c o v ( α ^ , λ ^ ) c o v ( β ^ , α ^ ) v a r ( β ^ ) c o v ( β ^ , δ ^ ) c o v ( β ^ , λ ^ ) c o v ( δ ^ , α ^ ) c o v ( δ ^ , β ^ ) v a r ( δ ^ ) c o v ( δ ^ , λ ^ ) c o v ( λ ^ , α ^ ) c o v ( λ ^ , β ^ ) c o v ( λ ^ , δ ^ ) v a r ( λ ^ ) .
Thus, the 100 ( 1 γ ) % ACIs of Ω is
Ω ^ i z γ 2 v a r ( Ω ^ i ) , Ω ^ i + z γ 2 v a r ( Ω ^ i ) , i = 1 , 2 , 3 , 4 ,
where v a r ( Ω ^ i ) = I i i 1 , i = 1 , 2 , 3 , 4 . z γ 2 is the percentile of the standard normal distribution with right-tail probability γ 2 . In order to construct the ACIs of r ( t ) and h ( t ) , we need to find the variances of them. Thus, Greene [19] implemented the delta method for this purpose as follows:
Assume that G = ( Λ α , Λ β , Λ δ , Λ λ ) is first-order partial derivative of Λ = ( r ( t ) , h ( t ) ) . Then, the asymptotic variances of Λ ^ = ( r ^ ( t ) , h ^ ( t ) ) can be written as
v a r ( Λ ^ ) ( G T I 1 ( Ω ^ ) G ) ( Ω = Ω ^ ) ,
where G T is the transpose matrix of G. Thus, the ( 1 γ ) 100 % ACIs for Λ ^ = ( r ^ ( t ) , h ^ ( t ) ) can be given by
Λ ^ z γ 2 v a r ( Λ ^ ) , Λ ^ + z γ 2 v a r ( Λ ^ ) .
Furthermore, to prevent the negative lower bound of the asymptotic confidence intervals, the use of the normal approximation for the log-transformed MLE is suggested by Meeker and Escobar [20]. Thus, two-sided ( 1 γ ) 100 % normal ACIs for θ ^ = α ^ , β ^ , δ ^ , λ ^ , r ^ t or h ^ t are given by
θ ^ . exp z γ 2 v a r ( θ ^ ) θ ^ , θ ^ . exp z γ 2 v a r ( θ ^ ) θ ^ .

3. Parametric Bootstrap Methods

It is known that the normal approximation is inappropriate when the sample is small. In this section, we propose a resampling technique, bootstrap procedure, to obtain more widely used confidence intervals. For this, two parametric bootstrap procedures are discussed to construct the bootstrap confidence intervals of θ = α , β , δ , λ ,   r t or h t . The percentile bootstrap (boot-p) and bootstrap-t (boot-t) confidence intervals using the same steps, which were mentioned in DiCiccio and Efron [21] and Hall [22], respectively, are discussed. Several authors have studied these two types of bootstrap, for example, Reiser et al. [23] and Besseris [24].

3.1. Percentile Bootstrap

  • From the original data x ̲ x 1 : m : n : k R ,   x 2 : m : n : k R , , x m : m : n : k R , compute the MLEs of the unknown parameters α ,   β ,   δ , and λ by solving the nonlinear Equations (6)–(10), denoted the estimates as θ ^ = ( α ^ , β ^ , δ ^ , λ ^ , r t , h t ) .
  • Use α ^ , β ^ , δ ^ and λ ^ to generate PFFC sample x ̲ * with the same values of R i , m ; ( i = 1 , 2 , , m ) and compute the bootstrap estimates θ ^ * .
  • Repeat Step 2 N times; then, we have θ ^ 1 * , θ ^ 2 * , , θ ^ N * .
  • Arrange all components in ascending order; the bootstrap estimates are θ ^ 1 * , θ ^ 2 * , , θ ^ N * .
  • Let φ x = P θ ^ * x be the CDF of θ ^ * . Define θ ^ b o o t - p x = φ 1 x for given x. Then, two-side 100 1 γ % percentile confidence intervals of θ given by
    θ ^ b o o t - p γ 2 , θ ^ b o o t - p 1 γ 2 .

3.2. Bootstrap-t

1,2.
The same as the percentile bootstrap.
3.
Compute the t-statistics T θ * = θ ^ * θ ^ v a r ( θ ^ * ) , where v a r ( θ ^ * ) obtained using the FIM for α , β , δ , λ , and delta method for r t and h t .
4.
Repeat Steps 2, 3, 4 N times; then, we have T θ * 1 , T θ * 2 , , T θ * N .
5.
Let ψ x = P T θ * x be the CDF of T θ * . Define θ ^ b o o t - t x = θ ^ + ψ 1 x v a r ( θ ^ * ) for given x. Then, two-side 100 1 γ % bootstrap-t confidence intervals of θ = α , β , δ , λ ,   r t or h t given by
θ ^ b o o t - t γ 2 , θ ^ b o o t - t 1 γ 2 .

4. Bayes Inference

Bayes estimation is quite different from MLE and bootstrap methods because it takes into consideration both the information from observed sample data and the prior information, in addition to symmetric and asymmetric loss functions. It can characterize the problems more rationally and reasonably. In this section, Bayesian inference procedure using MCMC technique are proposed to estimate the parameters α ,   β ,   δ and λ as well as r t or h t under both SE and LINEX loss functions. Moreover, the corresponding CRIs are constructed under MCMC technique. It is known that the family of gamma distributions is flexible enough to cover a large variety of prior beliefs of the experimenter. Therefore, we consider that the unknown parameters α ,   β ,   δ and λ are stochastically independently distributed with conjugate gamma prior as follows
π 1 ( α ) = b 1 a 1 Γ a 1 α a 1 1 e b 1 α , α > 0 , π 2 ( β ) = b 2 a 2 Γ a 1 β a 2 1 e b 2 β , β > 0 , π 3 ( λ ) = b 3 a 3 Γ a 3 δ a 3 1 e b 3 δ , δ > 0 , π 4 ( λ ) = b 4 a 4 Γ a 3 λ a 4 1 e b 4 λ , λ > 0 , ,
where a i and b i , i = 1 , 2 , 3 , 4 are the hyper-parameters which reflect the prior knowledge about the unknown parameters, and they are assumed to be known and non-negative. The joint prior of unknown parameters α ,   β ,   δ and λ is
π ( α , β , δ , λ ) α a 1 1 β a 2 1 δ a 3 1 λ a 4 1 e b 1 α + b 2 β + b 3 δ + b 4 λ .
Then, the joint posterior density of the unknown parameters Ω = ( α , β , δ , λ ) is
π * ( Ω | x ̲ ) = π ( Ω ) L ( Ω | x ̲ ) Ω π ( Ω ) L ( Ω | x ̲ ) d α d β d δ d λ .
Under SE loss function φ η ^ , η = η ^ η 2 , the Bayes estimation of any function W ( Ω ) is
W * ( Ω ) = Ω W ( Ω ) π ( Ω ) L ( Ω | x ̲ ) d α d β d δ d λ Ω L ( Ω | x ̲ ) π ( Ω ) d α d β d δ d λ .
Under asymmetric LINEX loss function ψ η ^ , η = e ε η ^ η ε η ^ η 1 , we have
ζ ^ = 1 ε log E e ε ζ , ε 0 ,
where ζ = α , β , δ , λ , r t or h t . For more accurate and important details, see EL-Sagheer et al. [25]. It is clear that the integrals in (15) and (16) cannot be computed in closed forms. Therefore, the Metropolis–Hastings (MH) within Gibbs sampling is a useful MCMC technique to estimate complex Bayes problems. Thus, the conditional posterior densities of α ,   β ,   δ and λ can be written as
π 1 * ( α | β , δ , λ , x ̲ ) α m β + a 1 1 i = 1 m 1 e α x i β λ + 1 × e b 1 α λ i = 1 m α x i β δ i = 1 m k R i + 1 e α x i β 1 λ ,
π 2 * ( β | α , δ , λ , x ̲ ) β m + a 2 1 α m β i = 1 m x i β 1 e α x i β λ + 1 × e b 2 β λ i = 1 m α x i β δ i = 1 m k R i + 1 e α x i β 1 λ ,
π 3 * ( δ | α , β , λ , x ̲ ) δ m + a 3 1 × e δ b 3 + i = 1 m k R i + 1 e α x i β 1 λ ,
and
π 4 * ( λ | α , β , δ , x ̲ ) λ m + a 4 1 i = 1 m 1 e α x i β λ × e λ b 4 + i = 1 m α x i β δ i = 1 m k R i + 1 e α x i β 1 λ .
It is noticeable that the conditional posterior density of δ given in (19) is Gamma density with shape parameter ( m + a 3 ) and scale parameter b 3 + i = 1 m k R i + 1 e α x i β 1 λ . Thus, by implementing any Gamma-generating routine, samples of δ can be simply generated. In addition, the conditional posterior densities of α ,   β and λ cannot be reduced analytically to well-known distributions. Thus, M-H algorithm within Gibbs sampling is implemented to conduct the MCMC methodology; see Tierney [26]. This procedure is illustrated in the following algorithm steps:
  • Use the MLEs as the initial value, denoted by α ^ 0 , β ^ 0 , δ ^ 0 and λ ^ 0 .
  • Set j = 1 .
  • Generate δ ( j ) from Gamma m + a 3 , b 3 + i = 1 m k R i + 1 e α x i β 1 λ .
  • Using M-H algorithm, generate α j , β j and λ j from π 1 * ( α j 1 | β j 1 , δ j , λ j 1 , x ̲ ) , π 2 * ( β j 1 | α j 1 , δ j , λ j 1 , x ̲ ) and π 4 * ( λ j 1 | α j 1 , β j 1 , δ j , x ̲ ) , respectively, with the normal distributions N α j 1 , v a r α ^ , N β j 1 , v a r β ^ and N λ j 1 , v a r λ ^ .
    (a)
    Generate α * from N α j 1 , v a r α ^ , β * from N β j 1 , v a r β ^ and λ * from N λ j 1 , v a r λ ^ .
    (b)
    Evaluate the acceptance probabilities
    Q α = min 1 , π 2 * ( α * | β j 1 , δ j , λ ( j 1 ) , x ̲ ) π 2 * ( α j 1 | β j 1 , δ j , λ ( j 1 ) , x ̲ ) , Q β = min 1 , π 2 * ( β * | α j , δ j , λ ( j 1 ) , x ̲ ) π 2 * ( β j 1 | α j , δ j , λ ( j 1 ) , x ̲ ) , Q λ = min 1 , π 4 * ( λ * | α ( j ) , β j , δ j , x ̲ ) π 4 * ( λ j 1 | α ( j ) , β j , δ j , x ̲ ) .
    (c)
    Generate a ρ 1 , ρ 2 and ρ 3 from a Uniform 0 , 1 distribution.
    (d)
    If ρ 1 < Q α , accept the proposal and set α * = α j , else set α j = α j 1 .
    (e)
    If ρ 2 < Q β , accept the proposal and set β * = β j , else set β j = β j 1 .
    (f)
    If ρ 3 < Q λ , accept the proposal and set λ * = λ j , else set λ j = λ j 1 .
  • Compute r t and h t as
    r j t = e δ j e α j x β j 1 λ j ,
    h j t = δ j λ j β j α j β j x β j + 1 e λ j α j x β j 1 e α j x β j λ j + 1 .
  • Set j = j + 1 .
  • Reiterate Steps (3)–(6) M times. The first M 0 simulated varieties are ignored to remove the affection of the selection of initial value and to guarantee the convergence. Then, the selected sample is α ( j ) , β ( j ) , δ j ,   λ ( j ) , r ( j ) ( t ) and h ( j ) ( t ) ,   j = M 0 + 1 , , M , for sufficiently large M, forms an approximate posterior sample which can be used to develop the Bayesian inference.
  • Under SE and LINEX loss functions, the approximate Bayes estimate of ζ (where ζ = α ,   β ,   δ ,   λ ,   r t and h t ) can be obtained by
    ζ ^ S E = 1 M M 0 j = M 0 + 1 M ζ j , ζ ^ L E = 1 ε log 1 M M 0 j = M 0 + 1 M e ε ζ j , ε 0 ,
    where M 0 is the burn-in period and ζ j = α j , β j , δ j , λ j , r j t and h j t .
  • To compute the CRIs of ζ , order ζ M 0 + 1 , ζ M 0 + 2 , , ζ M as ζ 1 , ζ 2 , , ζ M M 0 . Then, the 1 γ 100 % CRIs of ζ can be given by
    ζ M M 0 γ 2 , ζ M M 0 1 γ 2 .

5. Simulation Study

In this section, to compare the performance of the proposed Bayes estimators with the MLEs and bootstrap, we perform a Monte Carlo simulation study using different combinations of n, m, k and the different censoring scheme R (different R i values). Using the algorithm introduced by Balakrishnan and Sandhu [1], with distribution function 1 1 F x k , we generate PFFC samples from WFrD with the parameters α , β , δ , λ = 0.5 , 1 , 1.5 , 5 , respectively. Then, the true values of r t and h t at time t = 0.4 are evaluated to be 0.984462 and 0.342945 . The performance of estimators are evaluated in terms of mean square error (MSE) which is computed as M S E = 1 N i = 1 N ( ζ ^ j ( i ) ζ j ) 2 , where N = 10,000, j = 1 , 2 , , 6 ,   ζ 1 = α ,   ζ 2 = β ,   ζ 3 = δ ,   ζ 4 = λ ,   ζ 5 = r ( t ) and ζ 6 = h ( t ) for the point estimates and also average lengths (ALs) and coverage probability (CPs) for interval estimates (asymptotic, bootstrap, HPD). Bayes estimates and the CRIs are computed based on 12,000 MCMC samples and discard the first 2000 values as “burn-in”. In addition, we assume the informative Gamma priors for α ,   β ,   δ and λ , that is, when the hyper-parameters are a i = 1 and b i = 2 , i = 1 , 2 , 3 , 4 . Further, 95 % CRIs were computed for each simulated sample. In our study, we consider two different group sizes k = 2 , 5 and the following censoring schemes:
CS I: 
R 1 = n m ,   R i = 0 for i 1 .
CS II: 
R ( m + 1 ) / 2 = n m ,   R i = 0 for i ( m + 1 ) / 2 if m odd; R m / 2 = n m ,   R i = 0 for i m / 2 if m even.
CS III: 
R m = n m ,   R i = 0 for i m .
All the computations were performed by MATHEMATICA ver. 12. The results of the MSE of the estimates are reported in Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6, while the results of the ALs and CPs of the estimates are shown in Table 7, Table 8 and Table 9.
From the results in Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9, we observe the following:
  • As expected, from Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9, as sample sizes n , m increase, the MSEs and ALs decrease.
  • The Bayes estimates have the smallest MSEs and ALs for the unknown quantities ζ = α , β , δ , λ , r t and h t . Hence, Bayes estimates performed better than the MLEs and bootstrap methods.
  • The Bootstrap methods performed better than the ML method of ζ in terms of the MSEs and ALs. Moreover, boot-t performs better than boot-p in terms of the MSEs and ALs.
  • The Bayes estimate under LINEX with ε = 0.5 provides better estimates for ζ because of having the smallest MSEs.
  • The Bayes estimates of ζ under LINEX for the choice ε = 0.5 performed better than their estimates for the choice ε = 0.5 in the sense of having smaller MSEs.
  • From all of the tables, we observe that as the group size k increases, the MSEs and ALs associated with ζ increase.
  • For fixed sample sizes and observed failures, the first scheme I is the best scheme in the sense of having smaller MSEs and ALs.
  • The MLE, bootstrap and Bayesian methods have very close estimates and their ACIs have quite high CPs (around 0.95 ). Moreover, the Bayesian CRIs have the highest CPs.
In general, we can conclude that the best estimation method is the Bayes method using the asymmetric loss function, especially if the prior information about the problem under study is available. If the prior information about the problem under study is not available, it is preferable to use bootstrap methods that depend primarily on the maximum likelihood estimates.

6. Applications of COVID-19 Data

To clarify the inference methods discussed in the previous sections, we present an application of real-life data for the coronavirus. We consider two sets of real-life data, which were reported in Almongy et al. [27]. The first set represents COVID-19 mortality rate data belonging to Mexico over 108 days, recorded from 4 March to 20 July 2020. The second set represents COVID-19 data belonging to the Netherlands over 30 days, recorded from 31 March to 30 April 2020. The initial density shape is explored using the nonparametric kernel density (KD) estimation approach in Figure 1 and Figure 2, and the density is asymmetric and unimodal (bimodal) functions for data sets I and II, respectively, are noted. The normality condition and hazard rate shapes are tested via the quantile–quantile (QQ) and TTT plots in Figure 1 and Figure 2. The extremes are spotted from the box and violin plots in Figure 1 and Figure 2, and it is showed that some extreme observations were found.
For the purpose of the goodness-of-fit test, the K-S distance between the empirical and the fitted distribution functions have been computed. It is 0.05698 and the associated p-value is 0.8744 for the first set, while it is 0.09230 and the associated p-value is 0.93963 for the second set. Hence, the p-value for K-S has the highest value for the first and second sets. This leads us to conclude that the WFr model is the best fit for the two real data sets. Empirical plots are shown in Figure 3, which make it clear that the WFr model fits the data very well. Figure 4 and Figure 5 show that the ML estimators are unique.
It was found that the log-likelihood function is unimodal-shaped for each parameter.
The first set consisting of 108 survival times for 108 days. The data are randomly divided into 27 groups with k = 4 items within each group. The groups can be divided as follows: {1.041, 2.988, 4.344, 8.826}, {4.424, 5.242, 9.284, 10.855}, {6.105, 7.630, 7.903, 14.604}, {1.402, 6.327, 9.391, 14.962}, {5.143, 7.840, 9.550, 9.935}, {1.800, 6.015, 7.267, 13.220}, {3.286, 3.537, 6.370, 6.968}, {1.815, 4.949, 5.406, 6.182}, {1.867, 2.601, 2.926, 3.218}, {3.228, 7.854, 8.325, 8.551}, {1.923, 3.233, 3.257, 3.632}, {6.697, 6.560, 7.151, 8.108}, {2.077, 4.292, 4.557, 5.442}, {2.326, 2.500, 5.317, 6.535}, {3.298, 4.648, 4.697, 5.459}, {2.438, 6.122, 10.035, 10.685}, {3.440, 5.854, 6.656, 8.813}, {2.352, 5.985, 7.260, 10.043}, {3.215, 4.730, 11.665, 16.498}, {4.661, 4.909, 6.140, 6.625}, {3.499, 3.751, 6.814, 7.486}, {3.395, 6.412, 8.696, 12.042}, {2.070, 3.336, 10.158, 12.878}, {2.506, 2.838, 3.029, 5.392}, {2.065, 3.327, 3.359, 4.089}, {1.205, 3.027, 3.778, 10.383}, {2.058, 3.219, 3.922, 4.120}.
A PFFC sample, say x ( 1 ) , is obtained as follows: from the original data, the items are divided into n = 27 (number of groups) based on k = 4 (number of items within each group) with R = ( 2 , 0 , 2 , 0 , 2 , 0 , 2 , 0 , 2 , 0 , 2 , 0 , 2 ) (13 first-failure times are observed ( m = 13 ) and 14 groups are removed). The result in the obtained PFFC sample is x ( 1 ) = (1.041, 1.402, 1.800, 1.815, 1.867, 1.923, 2.058, 2.065, 2.070, 2.077, 2.326, 2.352, 2.438).
The second set consists of 30 survival times for 30 days. The data are randomly divided into 15 groups with k = 2 items within each group. The groups can be divided as follows: {1.273, 6.027}, {10.656, 12.274}, {1.974, 4.960}, {5.555, 7.584}, {3.883, 4.462}, {4.235, 5.307}, {7.968, 13.211}, {3.611, 3.647}, {6.940, 7.498}, {5.928, 7.099}, {2.254, 5.431}, {10.289, 10.832}, {4.097, 5.048}, {1.416, 2.857}, {3.461, 14.918}.
Suppose that a PFFC scheme is given by R = ( 2 , 1 , 1 , 1 , 1 , 1 , 1 ) ; then, a PFFC sample of size 7 out of 15 groups of data is obtained as x ( 2 ) = (1.273, 1.974, 2.254, 3.461, 3.611, 4.097, 4.235).
Under the previous data x ( 1 ) and x ( 2 ) , the MLEs and bootstrap estimates of α ,   β ,   δ ,   λ , r t and h t are determined to be as in Table 10. Moreover, Table 11 and Table 12 show the 95 % ACIs for α , β , δ , λ , r t and h t . Now, to compute the Bayesian estimates, the prior distributions of the parameters need to be specified. Because we have no prior information, we assume the noninformative Gamma priors for α ,   β ,   δ and λ , that is, when the hyper-parameters are a i = 0.0001 and b i = 0.0001 , i = 1 , 2 , 3 , 4 . Under the MCMC technique, the posterior analysis was performed across the combined M-H algorithm within the Gibbs sampler. To conduct the MCMC algorithm, which was described in Section 4, the initial values for the parameters α ,   β ,   δ and λ were taken to be their MLEs. In addition, 12,000 MCMC samples were generated. To avoid the effect of the initial values (starting point), we expunge the first 2000 samples as “burn-in”. Table 10, Table 11 and Table 12 show the Bayesian estimates as well as the 95 % CRIs for α ,   β ,   δ ,   λ ,   r t and h t .
From the results in Table 10, Table 11 and Table 12, we observe the following:
  • The model under study fits the data very well, and this is quite clear from the figures.
  • The estimation of the four methods is somewhat similar, with very slight differences, which gives a good impression to the reader.
  • Approximate confidence intervals are good, and all point estimates fall within them. In addition, there are slight differences in the lengths of the intervals, as expected.

7. Conclusions

The main aim of this paper was to develop different methods to estimate the unknown quantities ζ = α , β , δ , λ , r t and h t of the WFrD based on a PFFC scheme, which was introduced by Wu and Kuş [8]. The MLEs as well as the ACIs using asymptotic distributions were obtained. Furthermore, to obtain the CIs of the reliability and hazard functions, we used the delta method. Two parametric bootstrap procedures (boot-p and boot-t) were discussed to obtain more widely used confidence intervals. It is clear that, after studying the Bayesian estimates, the posterior distribution equations of the unknown quantities are complicated and so hard to reduce analytically to well-known forms. For this reason, we applied the MCMC method to compute the Bayes estimators. The Bayes estimates were computed under both the SE and LINEX loss functions. For illustrative purposes, two real data sets of COVID-19 mortality rates were considered. To check and compare the performance of the proposed methods, a simulation study was implemented with different sample sizes n , m , k and different CSs ( I , I I , I I I ) .

Author Contributions

Conceptualization, R.M.E.-S., M.A.A. and D.A.A.; methodology, M.E.-M. and M.S.E.; software, R.M.E.-S., M.E.-M. and M.S.E.; validation, F.H.E. and D.A.A.; formal analysis, R.M.E.-S. and M.S.E.; investigation, M.E.-M. and D.A.A.; resources, M.A.A., F.H.E. and D.A.A.; data curation, R.M.E.-S., M.E.-M. and M.S.E.; writing—original draft preparation, R.M.E.-S. and D.A.A.; writing—review and editing, R.M.E.-S. and M.S.E.; supervision, M.E.-M. and F.H.E.; funding acquisition, M.A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Majmaah University, grant number R-2022-110. The APC was funded by Muqrin A. Almuqrin.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Deanship of Scientific Research at Majmaah University for supporting this work under Project R-2022-110.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Balakrishnan, N.; Sandhu, R.A. A simple simulation algorithm for generating progressively type-II censored samples. Am. Stat. 1995, 49, 229–230. [Google Scholar]
  2. Fu, J.; Xu, A.; Tang, Y. Objective Bayesian analysis of Pareto distribution under progressive Type-II censoring. Stat. Probab. Lett. 2012, 82, 1829–1836. [Google Scholar] [CrossRef]
  3. Chen, P.; Xu, A.; Ye, Z. Generalized fiducial inference for accelerated life tests with Weibull distribution and progressively type-II censoring. IEEE Trans. Reliab. 2016, 65, 1737–1744. [Google Scholar] [CrossRef]
  4. Xu, A.; Zhou, S.; Tang, Y. A unified model for system reliability evaluation under dynamic operating conditions. IEEE Trans. Reliab. 2021, 70, 65–72. [Google Scholar] [CrossRef]
  5. Luo, C.; Shen, L.; Xu, A. Modelling and estimation of system reliability under dynamic operating environments and lifetime ordering constraints. Reliab. Eng. Syst. Saf. 2022, 218, 108136. [Google Scholar] [CrossRef]
  6. EL-Sagheer, R.M.; Shokr, E.M.; Mahmoud, M.A.W.; El-Desouky, B.S. Inferences for Weibull Fréchet distribution using a Bayesian and Non-Bayesian methods on gastric cancer survival times. Comput. Math. Methods Med. 2021, 9965856. [Google Scholar] [CrossRef] [PubMed]
  7. Wu, J.-W.; Hung, W.-L.; Tsai, C.-H. Estimation of the parameters of the Gompertz distribution under the first failure-censored sampling plan. Statistics 2003, 37, 517–525. [Google Scholar] [CrossRef]
  8. Wu, S.J.; Kuş, C. On estimation based on progressive first-failure-censored sampling. Comput. Stat. Data Anal. 2009, 10, 3659–3670. [Google Scholar] [CrossRef]
  9. Haj Ahmad, H.; Salah, M.M.; Eliwa, M.S.; Ali Alhussain, Z.; Almetwally, E.M.; Ahmed, E.A. Bayesian and non-Bayesian inference under adaptive type-II progressive censored sample with exponentiated power Lindley distribution. J. Appl. Stat. 2021. [Google Scholar] [CrossRef]
  10. Abushal, T.A. Estimation of the unknown parameters for the compound Rayleigh distribution based on progressive first-failure-censored sampling. Open J. Stat. 2011, 1, 161–171. [Google Scholar] [CrossRef] [Green Version]
  11. Soliman, A.A.; Abd-Ellah, A.H.; Abou-Elheggag, N.A.; Modhesh, A.A. Estimation of the coefficient of variation for non-normal model using progressive first-failure-censoring data. J. Appl. Stat. 2012, 12, 2741–2758. [Google Scholar] [CrossRef]
  12. Soliman, A.A.; Abd-Ellah, A.H.; Abou-Elheggag, N.A.; EL-Sagheer, R.M. Estimation Based on Progressive First-Failure Censored Sampling with Binomial Removals. Intell. Inf. Manag. 2013, 5, 117–125. [Google Scholar] [CrossRef] [Green Version]
  13. Mahmoud, M.A.W.; Soliman, A.A.; AAAbd-Ellah, A.H.; EL-Sagheer, R.M. Bayesian Inference and Prediction using Progressive First-Failure Censored from Generalized Pareto Distribution. J. Stat. Appl. Probab. 2013, 3, 269–279. [Google Scholar] [CrossRef]
  14. Ahmed, E.A.; Ali Alhussain, Z.; Salah, M.M.; Haj Ahmed, H.; Eliwa, M.S. Inference of progressively type-II censored competing risks data from Chen distribution with an application. J. Appl. Stat. 2020, 47, 2492–2524. [Google Scholar] [CrossRef]
  15. Xie, Y.; Gui, W. Statistical inference of the lifetime performance index with the Log-Logistic distribution based on progressive first-failure-censored data. Symmetry 2020, 12, 937. [Google Scholar] [CrossRef]
  16. Shi, X.; Shi, Y. Inference for Inverse Power Lomax distribution with progressive first-failure censoring. Entropy 2021, 23, 1099. [Google Scholar] [CrossRef]
  17. Afify, A.Z.; Yousof, H.M.; Cordeiro, G.M.; Ortega, E.M.; Nofal, Z.M. The Weibull Frechet distribution and its applications. J. Appl. Stat. 2016, 43, 2608–2626. [Google Scholar] [CrossRef]
  18. EL-Sagheer, R.M. Estimation of parameters of Weibull–Gamma distribution based on progressively censored data. Stat. Pap. 2018, 59, 725–757. [Google Scholar] [CrossRef] [Green Version]
  19. Greene, W.H. Econometric Analysis, 4th ed.; Prentice-Hall: New York, NY, USA, 2000. [Google Scholar]
  20. Meeker, W.Q.; Escobar, L.A. Statistical Methods for Reliability Data; Wiley: New York, NY, USA, 1998. [Google Scholar]
  21. DiCiccio, T.J.; Efron, B. Bootstrap confidence intervals. Stat. Sci. 1996, 11, 189–212. [Google Scholar] [CrossRef]
  22. Hall, P. Theoretical comparison of bootstrap confidence intervals. Ann. Stat. 1988, 16, 927–953. [Google Scholar] [CrossRef]
  23. Reiser, M.; Yao, L.; Wang, X.; Wilcox, J.; Gray, S. A Comparison of Bootstrap Confidence Intervals for Multi-Level Longitudinal Data Using Monte-Carlo Simulation, in ‘Monte-Carlo Simulation-Based Statistical Modeling’; Springer: Berlin/Heidelberg, Germany, 2017; pp. 367–403. [Google Scholar]
  24. Besseris, G.J. Evaluation of robust scale estimators for modified Weibull process capability indices and their bootstrap confidence intervals. Comput. Ind. Eng. 2019, 128, 135–149. [Google Scholar] [CrossRef]
  25. EL-Sagheer, R.M.; Eliwa, M.S.; Alqahtani, K.M.; EL-Morshedy, M. Asymmetric randomly censored mortality distribution: Bayesian framework and parametric bootstrap with application to COVID-19 data. J. Math. 2022. [Google Scholar] [CrossRef]
  26. Tierney, L. Markov chains for exploring posterior distributions (with discussion). Ann. Stat. 1994, 22, 1701–1722. [Google Scholar]
  27. Almongy, H.M.; Almetwally, E.M.; Aljohani, H.M.; Alghamdi, A.S.; Hafez, E.H. A new extended rayleigh distribution with applications of COVID-19 data. Results Phys. 2021, 23, 104012. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The KD, box, TTT, QQ and violin plots for data set I.
Figure 1. The KD, box, TTT, QQ and violin plots for data set I.
Symmetry 14 00956 g001
Figure 2. The KD, box, TTT, QQ and violin plots for data set II.
Figure 2. The KD, box, TTT, QQ and violin plots for data set II.
Symmetry 14 00956 g002
Figure 3. Empirical plots of WFr model for COVID-19 data in Mexico (panel (A)) and Netherlands (panel (B)).
Figure 3. Empirical plots of WFr model for COVID-19 data in Mexico (panel (A)) and Netherlands (panel (B)).
Symmetry 14 00956 g003
Figure 4. The log-likelihood profiles for data set I.
Figure 4. The log-likelihood profiles for data set I.
Symmetry 14 00956 g004
Figure 5. The log-likelihood profiles for data set II.
Figure 5. The log-likelihood profiles for data set II.
Symmetry 14 00956 g005
Table 1. MSE of estimates for the parameter α .
Table 1. MSE of estimates for the parameter α .
(k, n, m)SCMLEBootstrap Bayesian
Boot-pBoot-tSELINEX
ε = 0 . 5 ε = 0 . 5
( 2 , 30 , 15 ) I0.004130.003840.003520.003310.003130.00289
II0.004920.004930.004740.004520.004250.00372
III0.005260.005300.005130.004820.004670.00443
( 2 , 30 , 20 ) I0.003840.003710.003320.003240.002740.00248
II0.004350.004290.003910.003780.003460.00326
III0.004670.004620.004490.003990.003780.00352
( 2 , 50 , 35 ) I0.002380.002360.002190.002080.002010.00189
II0.002570.002490.002390.002320.002280.00215
III0.002950.002880.002730.002650.002590.00239
5 , 30 , 15 I0.005220.005230.004580.004250.004120.00399
II0.005930.005840.005290.004970.004690.00448
III0.006190.006180.005860.005530.005130.00494
5 , 30 , 20 I0.004590.004610.003850.003670.003510.00336
II0.004720.004710.004520.004360.004040.00384
III0.005150.005080.004770.004540.004410.00419
5 , 50 , 35 I0.003120.002990.002870.002660.002590.00228
II0.003750.003680.003470.003340.003250.00249
III0.004280.004210.003840.003710.003640.00291
Table 2. MSE of estimates for the parameter β .
Table 2. MSE of estimates for the parameter β .
(k, n, m)SCMLEBootstrap Bayesian
Boot-pBoot-tSELINEX
ε = 0 . 5 ε = 0 . 5
( 2 , 30 , 15 ) I0.007990.007850.007710.007530.007480.00713
II0.008250.008230.007940.007750.007620.00746
III0.008860.008810.008540.008340.008170.00792
( 2 , 30 , 20 ) I0.007220.007180.006990.006750.006640.00625
II0.007630.007550.007120.006990.006870.00652
III0.007960.007850.007640.007330.007290.00706
( 2 , 50 , 35 ) I0.005610.005570.005130.004980.004860.00445
II0.006150.006080.005790.005640.005580.00498
III0.006690.006610.006170.005990.005840.00537
5 , 30 , 15 I 0.008540.008470.008120.007980.007820.00746
II0.008840.008780.008350.008240.008230.00785
III0.009220.009180.008960.008810.008670.00842
5 , 30 , 20 I 0.008130.008110.007850.007780.007610.00713
II0.008420.008420.008140.008050.007930.00762
III0.008840.008790.008560.008470.008390.00786
5 , 50 , 35 I 0.007510.007490.007240.007050.006920.00665
II0.007920.007830.007640.007510.007260.00689
III0.008330.008280.008140.007830.007750.00744
Table 3. MSE of estimates for the parameter δ .
Table 3. MSE of estimates for the parameter δ .
(k, n, m)SCMLEBootstrap Bayesian
Boot-pBoot-tSELINEX
ε = 0 . 5 ε = 0 . 5
( 2 , 30 , 15 ) I0.009540.009470.009120.009110.008830.00856
II0.009820.009780.009540.009460.009320.00897
III0.011020.010870.009860.009780.009690.00942
( 2 , 30 , 20 ) I0.009120.009080.008830.008790.008540.00817
II0.009550.009510.009420.009350.008940.00876
III0.009840.009820.009610.009540.009250.00913
( 2 , 50 , 35 ) I0.008230.008170.008060.007830.007750.00737
II0.008650.008550.008390.008100.008080.00777
III0.008990.008920.008750.008630.008410.00813
5 , 30 , 15 I 0.009940.009820.009540.009490.009450.00916
II0.012430.012390.009940.009750.009640.00943
III0.013540.013510.011620.009930.009840.00971
5 , 30 , 20 I 0.009470.009450.009250.009180.009140.00873
II0.009840.009850.009610.009540.009430.00915
III0.011950.011790.011010.009940.009820.00960
5 , 50 , 35 I 0.009110.009070.008830.008790.008690.00826
II0.009450.009410.009260.009110.009030.00874
III0.009890.009870.009620.009550.009470.00914
Table 4. MSE of estimates for the parameter λ .
Table 4. MSE of estimates for the parameter λ .
(k, n, m)SCMLEBootstrap Bayesian
Boot-pBoot-tSELINEX
ε = 0 . 5 ε = 0 . 5
( 2 , 30 , 15 ) I0.041280.040640.039440.037840.036940.03351
II0.044560.044410.042910.041230.040980.03756
III0.049220.049130.046210.045230.044280.04154
( 2 , 30 , 20 ) I0.037880.037810.035990.033870.033170.02999
II0.042450.042340.040450.039110.038940.03541
III0.046230.045120.044590.042300.041550.03864
( 2 , 50 , 35 ) I0.025410.025110.023170.022740.021990.01884
II0.028450.027560.025470.024890.023970.02103
III0.032540.031840.029870.027540.026990.02480
5 , 30 , 15 I0.052330.051780.048650.047960.047850.04211
II0.055360.055240.052610.050120.050010.04699
III0.058660.057360.055270.053240.053110.04986
5 , 30 , 20 I0.048640.047990.045620.044670.044130.03950
II0.051550.051320.048350.047930.046810.04356
III0.056870.056460.053450.051890.050990.04786
5 , 50 , 35 I0.029450.028740.026140.024160.023540.02147
II0.032550.031870.030020.029940.027580.02540
III0.037610.036970.033540.033000.032140.02893
Table 5. MSE of estimates for the parameter r t .
Table 5. MSE of estimates for the parameter r t .
(k, n, m)SCMLEBootstrap Bayesian
Boot-pBoot-tSELINEX
ε = 0 . 5 ε = 0 . 5
( 2 , 30 , 15 ) I0.002240.002210.002100.001990.001950.00189
II0.002450.002410.002330.002110.002050.00195
III0.002710.002750.002670.002550.002390.00210
( 2 , 30 , 20 ) I0.002010.001990.001930.001910.001890.00171
II0.002350.002320.002050.002030.001980.00183
III0.002530.002540.002350.002240.002180.00197
( 2 , 50 , 35 ) I0.001590.001490.001350.001310.001290.00119
II0.001760.001710.001520.001490.001380.00127
III0.002020.001980.001760.001690.001550.00139
5 , 30 , 15 I0.003310.003320.003140.002980.002980.00271
II0.003610.003610.003540.003140.003140.00305
III0.004010.004110.003820.003670.003650.00342
5 , 30 , 20 I0.002550.002490.002350.002270.002270.00199
II0.002690.002680.002540.002480.002470.00212
III0.002930.002870.002780.002650.002640.00240
5 , 50 , 35 I0.001980.001910.001850.001850.001810.00164
II0.002210.002130.002040.001900.001900.00171
III0.002510.002490.002380.002380.002370.00219
Table 6. MSE of estimates for the parameter h t .
Table 6. MSE of estimates for the parameter h t .
(k, n, m)SCMLEBootstrap Bayesian
Boot-pBoot-tSELINEX
ε = 0 . 5 ε = 0 . 5
( 2 , 30 , 15 ) I0.001120.001100.001050.001000.000990.00078
II0.001290.001250.001150.001090.001060.00091
III0.001390.001380.001250.001160.001110.0099
( 2 , 30 , 20 ) I0.001010.001010.000970.000910.000910.00071
II0.001200.001220.001100.001010.001010.00082
III0.001330.001310.001270.001250.001150.00093
( 2 , 50 , 35 ) I0.000890.000880.000830.000810.000820.00067
II0.000950.000940.000890.000890.000900.00078
III0.001010.00990.00960.000950.000950.00086
5 , 30 , 15 I0.001250.001240.001210.001210.001200.00191
II0.001460.001440.001360.001320.001320.00119
III0.001750.001680.001540.001420.001410.00128
5 , 30 , 20 I0.001150.001140.001060.001060.001050.00092
II0.001230.001230.001180.001170.001170.00112
III0.001450.001440.001380.001370.001350.00127
5 , 50 , 35 I0.001010.001010.000910.000890.000890.00075
II0.001170.001160.000990.000970.000950.00086
III0.001350.001330.001100.001090.001080.00093
Table 7. ALs (first row) and CPs (second row) of 95 % ACIs for α and β .
Table 7. ALs (first row) and CPs (second row) of 95 % ACIs for α and β .
α β
(k, n, m)SCMLEBootstrapBayesMLEBootstrapBayes
ACIBoot-pBoot-tCRIACIBoot-pBoot-tCRI
( 2 , 30 , 15 ) I1.12451.08991.07450.99872.52632.18911.94661.6469
0.92540.93410.92940.94660.92410.92330.94560.9541
II1.29541.11231.11181.08942.64512.36542.14571.9546
0.93550.94550.93990.94710.93450.92560.93880.9491
III1.37811.24571.20181.11012.74912.53312.24782.1564
0.92890.95110.93420.95220.92870.92930.92380.9487
( 2 , 30 , 20 ) I1.11461.06471.05320.91452.41361.98451.79411.5342
0.94260.94520.93660.96880.92990.93740.94610.9588
II1.23461.09141.06981.02452.54612.21451.99631.7645
0.92340.93450.92840.95440.93240.92870.93560.9455
III1.29351.11411.11261.09632.63142.46652.13281.9746
0.91940.92340.91990.97630.92970.93660.93880.9541
( 2 , 50 , 35 ) I0.99840.98320.88430.75362.21461.74781.53991.3457
0.91480.92480.93620.96120.93540.92470.94550.9593
II1.08741.02280.98740.83452.34981.89991.67941.5617
0.95410.93120.94330.95420.94110.93990.95410.9641
III1.20121.10991.10960.91442.56662.24561.89741.7224
0.92580.92350.93470.93660.93200.93880.94210.9510
5 , 30 , 15 I1.34421.24571.11571.10992.89412.74622.34711.9946
0.94220.94610.95410.96550.92940.93230.94100.9500
II1.44571.35641.24341.11872.93612.90022.59412.2340
0.91990.92780.94710.95470.93420.93260.94780.9544
III1.49631.41001.32101.29142.98332.95412.65742.4562
0.91970.92360.94730.96310.94110.95140.93990.9481
5 , 30 , 20 I1.23541.11351.10861.08942.61242.58412.41892.139
0.94110.93450.91540.95940.93990.94320.94770.9521
II1.35471.29361.113651.10942.79422.71552.60842.3667
0.93440.92470.95360.94960.94120.93240.94870.9527
III1.46121.37451.2364 1.11022.83542.79412.62632.4110
0.94180.95100.93550.95430.91470.92140.92480.9392
5 , 50 , 35 I1.20311.19781.13110.99652.34512.221472.01451.9746
0.93990.93440.92550.94680.93980.93220.94750.9532
II1.31241.23471.2004 1.11392.46522.36542.16572.0197
0.94380.95120.94780.95870.94120.95390.95110.9746
III1.39241.29631.25941.20772.56312.47192.29482.1345
0.92140.93980.94550.96410.92810.92780.93560.9523
Table 8. ALs (first row) and CPs (second row) of 95 % ACIs for δ and λ .
Table 8. ALs (first row) and CPs (second row) of 95 % ACIs for δ and λ .
δ λ
(k, n, m)SCMLEBootstrapBayesMLEBootstrapBayes
ACIBoot-pBoot-tCRIACIBoot-pBoot-tCRI
( 2 , 30 , 15 ) I3.24783.21473.11482.24628.64437.54687.26316.1946
0.92480.93870.93470.95440.92980.92480.93890.9536
II3.52473.34563.2467 2.32338.85617.73147.56156.3210
0.93540.92880.94120.96110.94510.93810.92560.9498
III3.73563.65643.4326 2.48548.92327.87657.66946.5328
0.91990.93660.92550.94980.91930.92880.94810.9611
( 2 , 30 , 20 ) I3.01242.99912.89741.99988.53267.26576.83645.7698
0.92880.91490.92580.93990.92880.92660.93870.9642
II3.23213.21103.1342 2.21108.63457.44587.03265.9461
0.94870.95110.93470.95550.94750.93540.94920.9782
III3.42983.34573.26452.32148.74567.54727.28846.1632
0.94110.94820.95410.96140.92650.91990.93250.9611
( 2 , 50 , 35 ) I2.63122.48712.43571.72458.14576.86356.14735.2654
0.93880.92870.95710.96320.94730.93550.94560.952
II2.75642.56342.5166 1.98468.24656.99846.31455.4231
0.92770.92370.93450.95440.94980.95140.92560.9558
III2.85532.76512.7033 2.11448.36587.19846.53275.5321
0.93980.94120.94250.94990.93540.93550.92810.9521
(k, n, m)SCMLEBootstrapBayesMLEBootstrapBayes
ACIBoot-pBoot-tCRIACIBoot-pBoot-tCRI
5 , 30 , 15 I3.61783.51483.31492.64608.94437.84697.76326.5942
0.93970.93980.94750.95770.95000.93440.93870.9499
II3.71233.64783.5124 2.73419.04568.01137.96356.7251
0.92780.93850.94220.95240.93260.94870.92860.9621
III3.82543.76943.6324 2.82119.12368.26548.03216.8654
0.95410.94970.95140.95780.91560.92360.92440.9762
5 , 30 , 20 I3.32453.25683.14562.27548.66247.35486.93655.9699
0.91990.92870.93450.94870.93450.92880.93560.9566
II3.43583.34573.2147 2.36458.76357.45537.12476.2354
0.93150.92880.94870.95140.92430.92650.91990.9499
III3.56843.49353.39782.43918.83657.56477.35466.4321
0.94840.93470.39810.95770.92550.93660.94780.9516
5 , 50 , 35 I3.14573.11472.82241.99748.34587.06366.44725.5651
0.95000.94870.93900.95350.93660.92640.92870.9578
II3.25863.21143.0051 2.24578.46527.18796.55415.6234
0.92940.94110.93880.95470.93920.94810.93560.9614
III3.34513.33473.19632.35518.53267.33656.76545.7653
0.93790.92990.93780.94990.94210.93870.94560.9526
Table 9. ALs (first row) and CPs (second row) of 95 % ACIs for r t and h t .
Table 9. ALs (first row) and CPs (second row) of 95 % ACIs for r t and h t .
r(t)h(t)
(k, n, m)SCMLEBootstrapBayesMLEBootstrapBayes
ACIBoot-pBoot-tCRIACIBoot-pBoot-tCRI
( 2 , 30 , 15 ) I0.65320.63330.61540.56420.35650.33540.31250.2547
0.93990.94210.94850.96510.94100.93890.93480.9587
II0.66980.64750.6235 0.57230.36540.34620.32470.2645
0.91890.92450.93660.95460.92430.92460.94230.9641
III0.67540.65870.6378 0.58360.37650.35640.33120.2765
0.92450.92890.93550.94870.91230.91880.92650.9499
( 2 , 30 , 20 ) I0.63240.60540.57320.52360.31480.29870.26340.2198
0.94120.94930.95680.97120.924780.92510.92380.9514
II0.6436 0.61350.58460.53660.32540.31580.27940.2365
0.92560.92450.93560.96140.92450.93450.93560.9651
III0.65360.62450.59640.54760.33650.32980.28970.2476
0.91360.92540.92890.95660.92980.93560.93240.9558
( 2 , 50 , 35 ) I0.58640.56680.54630.49770.27360.26990.22980.1899
0.94560.93450.93620.95840.93140.92160.92450.9499
II0.59320.57360.5559 0.51460.28360.27960.23450.1997
0.92460.92580.93560.95830.92560.93610.93840.9544
III0.61230.58640.5678 0.53460.29360.28740.24550.2195
0.91560.91840.93650.95110.91990.92360.93260.9562
5 , 30 , 15 I0.68350.67360.64550.59470.39640.38570.35260.2945
0.93650.93870.93140.96140.92650.93660.92140.9745
II0.69980.68350.6545 0.60140.41100.39840.36760.3009
0.93220.94230.94110.94980.93470.93480.93650.9633
III0.71360.69360.6634 0.62010.42650.41050.37650.3187
0.91870.92650.91990.95470.93210.93540.93990.9584
5 , 30 , 20 I0.66550.63550.61340.56330.36230.35640.32540.2551
0.92560.92640.93450.95820.92650.93240.94780.9614
II0.67740.64680.6247 0.57630.37650.36580.33040.2669
0.92880.92350.94510.95470.93510.94150.93870.9526
III0.68360.65230.63480.58630.38640.37840.34660.2863
0.92840.91230.93540.96450.92470.94810.93690.9485
5 , 50 , 35 I0.61680.59670.57630.52750.32650.31980.27960.2187
0.94710.93950.94550.96100.95120.94980.93990.9587
II0.62580.61450.5889 0.53020.33590.32570.29000.2341
0.93450.92890.93560.95840.92450.93550.93870.9612
III0.63210.62350.59940.54670.34620.33640.31660.2499
0.92630.92650.93710.94980.91870.93250.92890.9499
Table 10. Point estimates of α ,   β ,   δ ,   λ ,   r t and h t .
Table 10. Point estimates of α ,   β ,   δ ,   λ ,   r t and h t .
Data(k, n, m)ParameterMLEBootstrapBayesian
Boot-pBoot-tSELINEX
ε = 0 . 5 ε = 0 . 5
x ( 1 ) ( 4 , 27 , 13 ) α 2.85652.67842.58992.71122.70472.4999
β 1.54811.49771.47231.51891.45661.4119
δ 0.98510.99630.96640.95470.96870.9489
λ 0.35320.37510.34650.33980.34560.3297
r 0.4 0.99970.97670.95670.97670.97440.9653
h 0.4 0.00410.00450.00390.00420.00430.0038
x ( 2 ) ( 2 , 15 , 7 ) α 1.50341.51491.49571.47681.47221.4701
β 5.04845.15564.99744.96434.91074.8777
δ 0.05040.05060.05010.04970.04990.04875
λ 0.57790.61440.59970.56680.57570.5638
r 0.4 0.99870.98870.97670.98560.97770.9686
h 0.4 0.00310.00270.00250.00240.00230.0021
Table 11. The 95 % ACIs and CRIs of α ,   β ,   δ ,   λ ,   r t and h t based on x ( 1 ) .
Table 11. The 95 % ACIs and CRIs of α ,   β ,   δ ,   λ ,   r t and h t based on x ( 1 ) .
Method α β δ
IntervalLengthIntervalLengthIntervalLength
A C I [1.2485,4.4645]3.2160[0.0523,8.1486]8.0963[0.1205,3.0907]2.9702
b o o t - p [1.3234,4.7217]3.3983[0.2357,7.6541]7.4184[0.2456,3.2641]3.0185
b o o t - t [0.9984,4.1568]3.1584[0.0584,6.7542]6.6958[0.0874,2.8456]2.7582
C R I [1.0749,4.2559]3.1810[0.1425,5.3622]5.2197[0.1005,2.5312]2.4307
λ r t h t
IntervalLengthIntervalLengthIntervalLength
A C I [0.0485,0.6579]0.6094[0.9918,1.0076]0.0158[0.0008,1.0076]1.0068
b o o t - p [0.0524,0.6149]0.5625[0.9484,1.0099]0.0615[0.0009,0.9456]0.9447
b o o t - t [0.0147,0.6037]0.5890[0.9473,1.0087]0.0614[0.0007,0.9249]0.9242
C R I [0.0099,0.5831]0.5732[0.9402,1.0035]0.0633[0.0004,0.9189]0.9185
Table 12. The 95 % ACIs and CRIs of α ,   β ,   δ ,   λ ,   r t and h t based on x ( 2 ) .
Table 12. The 95 % ACIs and CRIs of α ,   β ,   δ ,   λ ,   r t and h t based on x ( 2 ) .
Method α β δ
IntervalLengthIntervalLengthIntervalLength
A C I [0.6567,2.3502]1.6935[1.2560,9.3527]8.0967[0.0017,0.2025]0.2008
b o o t - p [0.6241,2.2633]1.6392[1.5476,9.8475]8.2999[0.0014,0.2361]0.2347
b o o t - t [0.5147,2.1685]1.6538[1.8947,9.7566]7.8619[0.0011,0.2287]0.2276
C R I [0.4566,2.0472]1.5906[1.9965,8.1253]6.1288[0.0009,0.2154]0.2145
λ r t h t
IntervalLengthIntervalLengthIntervalLength
A C I [0.0826,1.9817]1.8991[0.9517,1.0023]0.0506[0.0005,1.0047]1.0042
b o o t - p [0.1145,2.1475]2.0330[0.9454,1.0025]0.0571[0.0003,0.9599]0.9596
b o o t - t [0.0967,1.9987]1.9020[0.9321,0.9997]0.0676[0.0003,0.9467]0.9464
C R I [0.0845,1.9632]1.8787[0.9254,0.9873]0.0619[0.0002,0.9354]0.9352
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

EL-Sagheer, R.M.; Almuqrin, M.A.; El-Morshedy, M.; Eliwa, M.S.; Eissa, F.H.; Abdo, D.A. Bayesian Inferential Approaches and Bootstrap for the Reliability and Hazard Rate Functions under Progressive First-Failure Censoring for Coronavirus Data from Asymmetric Model. Symmetry 2022, 14, 956. https://doi.org/10.3390/sym14050956

AMA Style

EL-Sagheer RM, Almuqrin MA, El-Morshedy M, Eliwa MS, Eissa FH, Abdo DA. Bayesian Inferential Approaches and Bootstrap for the Reliability and Hazard Rate Functions under Progressive First-Failure Censoring for Coronavirus Data from Asymmetric Model. Symmetry. 2022; 14(5):956. https://doi.org/10.3390/sym14050956

Chicago/Turabian Style

EL-Sagheer, Rashad M., Muqrin A. Almuqrin, Mahmoud El-Morshedy, Mohamed S. Eliwa, Fathy H. Eissa, and Doaa A. Abdo. 2022. "Bayesian Inferential Approaches and Bootstrap for the Reliability and Hazard Rate Functions under Progressive First-Failure Censoring for Coronavirus Data from Asymmetric Model" Symmetry 14, no. 5: 956. https://doi.org/10.3390/sym14050956

APA Style

EL-Sagheer, R. M., Almuqrin, M. A., El-Morshedy, M., Eliwa, M. S., Eissa, F. H., & Abdo, D. A. (2022). Bayesian Inferential Approaches and Bootstrap for the Reliability and Hazard Rate Functions under Progressive First-Failure Censoring for Coronavirus Data from Asymmetric Model. Symmetry, 14(5), 956. https://doi.org/10.3390/sym14050956

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