Next Article in Journal
Classification of Literary Works: Fractality and Complexity of the Narrative, Essay, and Research Article
Previous Article in Journal
A Novel Model on Reinforce K-Means Using Location Division Model and Outlier of Initial Value for Lowering Data Cost
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

E-Bayesian Estimation for the Weibull Distribution under Adaptive Type-I Progressive Hybrid Censored Competing Risks Data

1
Department of Statistics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Department of Mathematics, Faculty of Science, Al-Azhar University, Cairo 11884, Egypt
3
Department of Mathematics, Faculty of Science, Islamic University of Madinah, Madinah 42351, Saudi Arabia
4
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
*
Author to whom correspondence should be addressed.
Submission received: 28 June 2020 / Revised: 13 August 2020 / Accepted: 13 August 2020 / Published: 17 August 2020

Abstract

:
This article focuses on using E-Bayesian estimation for the Weibull distribution based on adaptive type-I progressive hybrid censored competing risks (AT-I PHCS). The case of Weibull distribution for the underlying lifetimes is considered assuming a cumulative exposure model. The E-Bayesian estimation is discussed by considering three different prior distributions for the hyper-parameters. The E-Bayesian estimators as well as the corresponding E-mean square errors are obtained by using squared and LINEX loss functions. Some properties of the E-Bayesian estimators are also derived. A simulation study to compare the various estimators and real data application is applied to show the applicability of the different estimators are proposed.

1. Introduction

In life-testing and reliability studies, both type I and type II censoring schemes are widely used. These two types of censoring schemes are described as follows: Consider n identical components are placed in the test, in type I censoring, the experiment continues up to a predetermined time τ . However, in the type II censoring scheme, the experiment is terminated when a predetermined number of failures m < n occurs. For a mixture of type I and type II censoring schemes, which is called the type I hybrid censoring scheme, is introduced by Epstein [1], and the life test experiment is terminated at a random time τ * = min { x m : m : n , τ } . Childs et al. [2] proposed a new hybrid censoring scheme called a type-II hybrid censoring scheme in which the experiment would terminate at the random time τ * = max { x m : m : n , τ } . These schemes do not allow for removing the components from the experiment at any time other than the terminal point. A more general censoring scheme called progressive type II censoring is used to deal with this problem.
In the progressive type II censoring, n components are placed in a lifetime testing experiment, and m is a fixed number of components to be failed. At the time of the first failure x 1 : m : n , R 1 components are randomly removed from the remaining n 1 components. Similarly, at the time of the second failure x 2 : m : n , R 2 components of the remaining n 2 R 1 components are randomly removed, and so on. At the time of the m-th failure x m : m : n , all the remaining n m i = 1 m 1 R i components are removed. The R i , i = 1 , , m are fixed and predetermined prior to the study.
The type I progressive hybrid censoring scheme is considered by Kundu and Joarder [3], which is a mixture of type II progressive and hybrid type I censoring schemes—in which n identical components are put under testing with R i , i = 1 , , m , predetermined progressive censoring scheme and the experiment is terminated at random time τ * = min { x m : m : n , τ } . The experiment stops at the time x m : m : n , if the m-th failure occurs before the time τ . However, if the m-th failure does not occur before time τ and only J failures occur before the time τ ; then, at the time τ , the experiment terminates and all the remaining components are removed. In addition, a type II progressive hybrid censoring scheme is discussed by Childs et al. [4], where the experiment is terminated at time τ * = max { x m : m : n , τ } .
A new censoring scheme, called an adaptive type II progressive hybrid censoring, is introduced by Ng et al. [5], where the number of failures m and the corresponding progressively scheme is given, but no components will be removed when the experimental time passes time τ ; see, for example, Balakrishnan and Kundu [6]. Another adaptive progressive hybrid censoring scheme, called the adaptive type I progressive hybrid censoring scheme (AT-I PHCS), which assures the termination of the lifetime testing experiment at a fixed time τ , and results in a higher efficiency in estimations, is proposed by Lin and Huang [7]. The AT-I PHCS can be described as: suppose n identical components are placed under testing with prefixed R i , i = 1 , , m , 1 m n ; the experiment is terminated at a prefixed time τ , where τ ( 0 , ) . At the time x 1 : m : n , R 1 of the remaining components are randomly removed, at the time x 2 : m : n , R 2 of the remaining components are randomly removed, and so on. Let the number of failures that occur before time τ be J. If the m-th failure x m : m : n occurs before time τ , the process will not stop, but continue to observe failures without any further withdrawals until reach time τ . Then, all remaining components R J * = n J i = 1 J R i are removed at time τ , and the experiment is terminated. The progressive censoring scheme in this case will become R 1 , R 2 , , R m , R m + 1 , , R J , where R m = R m + 1 = = R J = 0 . Otherwise, when x m : m : n > τ , the process will have a progressive censoring scheme such as R 1 , R 2 , , R J . AT-I PHCS is important when the time is the main goal in the experiment, and it is requisite to terminate the experiment at a predetermined time in any case of number of failures.
The point and interval estimation for the exponential distribution are studied by Lin and Huang [7] and investigated Bayesian sampling plans under different progressive censoring schemes. The maximum likelihood and Bayesian estimation for a two-parameter Weibull distribution based on AT-I PHCS are discussed by Lin et al. [8]. They derived the Bayes estimates of the unknown parameters by using the approximated form of Lindley [9] and Tierney and Kadane [10].
On the other hand, the loss function is important in Bayesian methods. In the Bayesian inference, the most commonly used loss function is the squared error loss. It is well known that the use of symmetric loss functions may be inappropriate in many circumstances, particularly when positive and negative errors have different consequences. One of the most commonly used asymmetric loss functions is the LINEX (linear exponential) loss function. It was introduced by Varian [11] became popular due to Zellner [12].
The failure of components in reliability analysis at the same time may be attributable to more than one reason. These reasons are competing for the experimental component for the failure. This is known as the competing risks model. The data in the competing risks analysis are comprised of a failure time and the associated reason for failure. The reasons for failure can be assumed to be independent or dependent.
The latent failure times model in this paper is assumed as suggested by Cox [13]. In addition, the failure times are independently distributed. The failure is due to more than one reason, see Crowder [14]. The competing risks data discussed here under AT-I PHCS and Weibull distributions with common shape parameters are assumed to the failure times. The competing risks data under AT-I PHCS with the assumption of exponential distribution are analyzed by Ashour and Nassar [15]. The exponential distribution has a constant failure rate, so it has serious limitations in modeling lifetime data. The inference for Weibull distribution under adaptive type-I progressive hybrid censored competing risks data are investigated by Ashour and Nassar [16].
The main aim of this paper is studying the competing risk model under AT-I PHCS. The lifetimes under the competing risks have independent Weibull distributions with common shape parameters. This paper can be organized as follows:
The model description and the notation are introduced in Section 1. The maximum likelihood estimation of the unknown parameters is established in Section 3. Bayesian estimation of the parameters under squared error (SE) and LINEX loss functions are discussed in Section 4. The expected Bayesian estimation under squared error and LINEX loss functions are derived in Section 5. Some properties of the E-Bayesian estimators are also derived in Section 6. Finally, two examples of the real data set and numerical simulation results are presented in Section 7.

2. Model Description

Suppose n identical components are put into a lifetime test with prefixed progressive censoring scheme R i , i = 1 , , m , 1 m n and the experiment is terminated at time τ , where τ ( 0 , ) . The lifetime for the components are assumed to be the Weibull distribution. Under the adaptive type-I progressive censoring scheme and in the presence of competing risks data, we have the following observation:
( X 1 : m : n , c 1 , R 1 ) , , ( X m 1 : m : n , c m 1 , R m 1 ) , ( X m : m : n , c m , 0 ) , , ( X J : m : n , c J , 0 ) ,
where the indicator c i is denoting the reason of failure, and J is the number of failures before time τ .
Consider that R J * is the number of remaining components left at the time point τ with R m = R m + 1 = = R J = 0 . Let c i ( 1 , 2 ) , here, c i = k , k = 1 , 2 means that the component i has failed due to reason k. Furthermore, we define
I 1 ( c i = 1 ) = 1 , c i = 1 0 , otherwise and I 2 ( c i = 2 ) = 1 , c i = 2 0 , otherwise
Thus, the random variables D k = i = 1 J I 1 ( c i = k ) describe the number of failures due to the reason k , k = 1 , 2 for failure.
The latent failure times X 1 i and X 2 i are assumed to be independent. X k i , i = 1 , , n , have Weibull distributions with parameters β and λ k , ( k = 1 , 2 ) (same shape and different scale parameter). The corresponding survival function F ¯ k and the hazard rate function h k are given, respectively, by
F ¯ k ( x ; λ k , β ) = e λ k x β , x > 0 , λ k , β > 0 ,
and
h k ( x ; λ k , β ) = λ k β x β 1 , x > 0 , λ k , β > 0 .

3. Maximum Likelihood Estimation

Based on adaptive type-I progressive hybrid censoring with a sample of size m obtained from a life test experiment of n items from the Weibull distribution, the likelihood function of the observed data ( x 1 , c 1 ) , , ( τ , R J * ) , for a given scheme R 1 , R 2 , , R m 1 , 0 , 0 , R J * , can be written as
L = c i = 1 J f 1 ( x i ) F ¯ 2 ( x i ) I ( c i = 1 ) f 2 ( x i ) F ¯ 1 ( x i ) I ( c i = 2 ) F ¯ 1 ( x i ) F ¯ 2 ( x i ) R i F ¯ 1 ( τ ) F ¯ 2 ( τ ) R J * ,
where x i = x i : m : n and c is a constant that doesn’t depend on the parameters.
The likelihood function can be represented by applying the identity f k = h k F ¯ k as follows:
L = c i = 1 J h 1 ( x i ) I ( c i = 1 ) h 2 ( x i ) I ( c i = 2 ) F ¯ 1 ( x i ) F ¯ 2 ( x i ) R i + 1 F ¯ 1 ( τ ) F ¯ 2 ( τ ) R J * ,
where J = D 1 + D 2 and J > 0 .
From Equations (2), (3) and (5), the likelihood function ignoring the normalized constant can be written as follows:
L ( λ 1 , λ 2 | x _ ) = λ 1 D 1 λ 2 D 2 φ ( β ; x _ ) e λ 1 + λ 2 T ,
where x _ = ( x 1 , x 2 , , x r ) , φ ( β ; x _ ) = β J i = 1 J x i β 1 , T = i = 1 J ( 1 + R i ) x i β + R J * τ β , R J * = n J i = 1 J R i , and x i = x i : m : n for simplicity of notation.
Assume that the parameter β is known and takes the natural logarithm of (6); then, differentiating with respect to λ 1 and λ 2 , we can obtain the MLEs of λ 1 and λ 2 as
λ ^ 1 = D 1 T and λ ^ 2 = D 2 T .
Before progressing further, we assume the case of 1 D i m , i = 1 , 2 in all of the next sections. Three ways to get an estimation performance assessment are the expectation, variance, and mean-square error (MSE). It can be shown that T has a Gamma ( D i , λ i , i = 1 , 2 ) distribution. Therefore, the expectations and mean square errors for all estimators mentioned above are obtained respectively as follows: For MLE, λ ^ M L E , of λ ,
E λ i ^ M L E = 0 D i t λ i D i Γ ( D i ) t D i 1 e λ i t d t = D i D i 1 λ i , i = 1 , 2 ,
and
M S E λ i ^ M L E = E λ i λ i ^ M L E 2 = λ i 2 2 λ i E ( λ i ^ M L E ) + E λ i ^ M L E 2 = D i + 2 ( D i 1 ) ( D i 2 ) λ i 2 , i = 1 , 2 .

4. Bayesian Estimation

In this section, we obtain the Bayes estimators of the parameters λ 1 and λ 2 based on SE and LINEX loss functions. For developing the Bayesian estimation, we assume that the parameters λ 1 and λ 2 are independently distributed according to gamma distribution. Let λ j , j = 1 , 2 have a gamma prior with scale parameters c j and shape parameters d j . The joint prior density of λ 1 and λ 2 can be written as follows:
g ( λ 1 , λ 2 ) λ 1 d 1 1 λ 2 d 2 1 e c 1 λ 1 c 2 λ 2 , c j , d j > 0 , j = 1 , 2 .
The joint prior (9), can be derived as a special case from the dependent prior proposed by Kundu and Pradhan [17]. Based on the likelihood function (6) and the joint prior density (9), the joint posterior density of λ 1 and λ 2 given β , and the data can be written as
q ( λ 1 , λ 2 | x _ ) λ 1 D 1 + d 1 1 λ 2 D 2 + d 2 1 e λ 1 ( T + c 1 ) λ 2 ( T + c 2 ) .
From (10), we observe that the posterior density functions of λ 1 and λ 2 are g a m m a ( D 1 + d 1 , T + c 1 ) and g a m m a ( D 2 + d 2 , T + c 2 ) , respectively. Based on the SE loss function, the Bayes estimators of λ 1 and λ 2 can be obtained as the posterior means with the following forms:
λ ^ j B S = D j + d j T + c j , j = 1 , 2 .
Considering d 1 = d 2 = c 1 = c 2 = 0 , the Bayes estimators in (11) coincide with the MLEs in (7).
Now, we obtain the Bayes estimators based on LINEX loss function proposed by Varian [11]. The LINEX loss function with parameter ω is given by
ψ ( μ B L , μ ) = e ω ( μ B L μ ) ω ( μ B L μ ) 1 ,
where ω 0 is a constant, μ is the parameter to be estimated, and μ B L is the Bayes estimator under the LINEX loss function. From (12), the Bayes estimator μ B L is
μ B L = 1 ω ln E ( e ω μ ) , ω 0 .
From (10) and (13), the Bayes estimators of the parameters λ 1 and λ 2 can be obtained as
λ ^ j B L = D j + d j ω ln 1 + ω T + c j , j = 1 , 2 .
Under the SE loss function, the MSE of the parameters, λ j , j = 1 , 2 , can be obtained as
M S E λ ^ j B S ( c j , d j ) | T = E ( λ j λ ^ j B S ( c j , d j ) ) 2 | T = E ( λ j 2 | T ) 2 λ ^ j B S ( c j , d j ) E ( λ j | T ) + ( λ ^ j B S ( c j , d j ) ) 2 = ( D j + d j + 1 ) ( D j + d j ) ( T + c j ) 2 2 D j + d j T + c j D j + d j T + c j + D j + d j T + c j 2 = D j + d j ( T + c j ) 2 , j = 1 , 2 .
In addition, when 1 < w T + c j < 1 , the Bayesian estimator, λ ^ j B L ( c j , d j ) , j = 1 , 2 , has the following properties:
M S E λ ^ j B L ( c j , d j ) | T = E ( λ j λ ^ j B L ( c j , d j ) ) 2 | T = E ( λ j 2 | T ) 2 λ ^ j B L ( c j , d j ) E ( λ j | T ) + ( λ ^ j B L ( c j , d j ) ) 2 = ( D j + d j + 1 ) ( D j + d j ) ( T + c j ) 2 2 ( D j + d j ) 2 w ( T + c j ) ln ( 1 + w T + c j )         + D j + d j w ln ( 1 + w T + c j ) 2 = D j + d j ( T + c j ) 2 + D j + d j T + c j 2 i = 2 w ( i 1 ) ( 1 ) ( i 1 ) i ( T + c j ) ( i 1 ) 2 , j = 1 , 2 .

5. E-Bayesian Estimation

E-Bayesian (Expected Bayesian) estimation was first introduced in literature by Han [18]. He obtained the estimate of the scale parameter of the Weibull distribution based on SE loss function and also derived the properties of the E-Bayesian estimation. E-Bayesian based on three different prior distributions of hyper parameter are used in this section to investigate the influence of different prior distributions on the E-Bayesian of λ j , j = 1 , 2 . For more relevant research about the E-Bayesian estimation, see Han [19], Jaheen and Okasha [20], Azimi et al. [21], Okasha [22], Okasha and Wang [23], and Abdallah and Jumping [24].
Han [18] stated that the prior distribution of d j and c j , j = 1 , 2 should be determined to ensure that the prior distribution g ( λ j ) , j = 1 , 2 is a decreasing function in λ j , j = 1 , 2 . To be sure from this condition, we find the first derivative of g ( λ j ) , j = 1 , 2 with respect to λ j , j = 1 , 2 as
g ( λ j ) λ j λ j d j 2 e c j λ j [ d j 1 c j λ j ] , j = 1 , 2 .
Thus, for 0 < d j < 1 and c j > 0 , j = 1 , 2 , the function g ( λ j ) λ j < 0 , j = 1 , 2 , and therefore g ( λ j ) , j = 1 , 2 is a decreasing function of λ j , j = 1 , 2 . Suppose that d j and c j , j = 1 , 2 are independent with bivariate PDF given by
π ( d j , c j ) = π j ( d j ) π j ( c j ) , j = 1 , 2 .
Then, the E-Bayesian estimate of the parameter λ j , j = 1 , 2 , (expectation of the Bayesian estimate of λ j , j = 1 , 2 ), according to Han [18] can be obtained as follows:
λ j ^ E B = D λ j ^ B ( d j , c j ) π ( d j , c j ) d d j d c j , j = 1 , 2 ,
where λ j ^ B ( d j , c j ) is the Bayesian estimator of λ j given by Equations (11) and (14), λ j ^ E B ( d j , c j ) is the expected Bayes estimate of λ j , j = 1 , 2 under any loss function. The value range of hyper parameter d j and c j , j = 1 , 2 satisfy D = ( a , b ) | 0 < d j < 1 , 0 < c j < s , j = 1 , 2 . Suppose the prior distribution of d j and c j , j = 1 , 2 are beta distribution and uniform distribution in ( 0 , s ) , respectively. For more details about E-Bayesian estimation, see Han [18], Okasha [22] and Okasha, and Wang [23].
The expected mean square error (E-MSE) of the parameter λ j , j = 1 , 2 according to Han [19] can be obtained as follows:
E - MSE ( λ j ^ ) = D M S E ( λ j ^ ) ( d j , c j ) π ( d j , c j ) d d j d c j , j = 1 , 2 ,
where M S E ( λ j ^ ) ( d j , c j ) is the MSE of λ j , j = 1 , 2 under any loss function.

5.1. E-Bayesian Estimation under SE Loss Function

Based on three different prior distributions of the hyper-parameters d j and c j , the E-Bayesian estimates of the parameter λ j , j = 1 , 2 can be obtained. Accordingly, these prior distributions are selected to show the effect of the different prior distributions on the E-Bayesian estimation of the parameter λ j , j = 1 , 2 . The selected priors distributions are given by
π 1 ( d j , c j ) = 1 s B ( a , b ) d j a 1 ( 1 d j ) b 1 0 < d j < 1 , 0 < c j < s π 2 ( d j , c j ) = 2 s 2 B ( a , b ) ( s c j ) d j a 1 ( 1 d j ) b 1 , 0 < d j < 1 , 0 < c j < s π 3 ( d j , c j ) = 2 c j s 2 B ( a , b ) d j a 1 ( 1 d j ) b 1 , 0 < d j < 1 , 0 < c j < s .
These prior distributions are used to guarantee that g ( λ j ) , j = 1 , 2 is a decreasing function in λ j , j = 1 , 2 . Now, under SE loss function, the E-Bayesian estimates of the parameter λ j , j = 1 , 2 can be obtained from (11), (17), and (19). Using the prior distribution, π 1 ( d j , c j ) is given by
λ j ^ E B S 1 = D λ j ^ B S ( d j , c j ) π 1 ( d j , c j ) d c j d d j = 1 s D j + a a + b ln 1 + s T , j = 1 , 2 .
Similarly, under SE loss function, the E-Bayesian estimates of λ j , j = 1 , 2 based on π 2 ( d j , c j ) and π 3 ( d j , c j ) are given, respectively, by
λ j ^ E B S 2 = 2 s D j + a a + b 1 + T s ln 1 + s T 1 , j = 1 , 2 ,
and
λ j ^ E B S 3 = 2 s D j + a a + b 1 T s ln 1 + s T , j = 1 , 2 .
The E-MSE of the parameter λ j , j = 1 , 2 can be obtained from (15), (18), and (19). Using the prior distribution, π 1 ( d j , c j ) is given by
E - MSE ( λ j ^ B S 1 ) = D MSE ( λ j ^ B S ) ( d j , c j ) π 1 ( d j , c j ) d c j d d j = D j + a a + b 1 T ( T + s ) , j = 1 , 2 .
Similarly, under SE loss function, the E-MSE of λ j , j = 1 , 2 based on π 2 ( d j , c j ) and π 3 ( d j , c j ) are given, respectively, by
E - MSE ( λ j ^ B S 2 ) = 2 s 2 D j + a a + b s T ln 1 + s T , j = 1 , 2 ,
and
E - MSE ( λ j ^ B S 3 ) = 2 s 2 D j + a a + b ln 1 + s T s T + s , j = 1 , 2 .

5.2. E-Bayesian Estimation under LINEX Loss Function

Under LINEX loss function, the E-Bayesian estimation of λ j , j = 1 , 2 can be obtained by using the different prior distributions of the hyper-parameters given by (19). For the prior distribution, π 1 ( d j , c j ) , j = 1 , 2 and, based on (14), (17), and (19), the E-Bayesian estimate of λ j , j = 1 , 2 is obtained as
λ j ^ E B L 1 = D λ j ^ B L ( d j , c j ) π 1 ( d j , c j ) d c j d d j = 1 s ω D j + a a + b s ln 1 + ω s + T + ( T + ω ) ln 1 + s T + ω         T ln 1 + s T , j = 1 , 2 .
Similarly, under LINEX loss function, the E-Bayesian estimates of λ j , j = 1 , 2 using π 2 ( d j , c j ) and π 3 ( d j , c j ) are given, respectively, by
λ j ^ E B L 2 = D j + a a + b 1 ω ln 1 + ω T ( T + s ) 2 s 2 ω ln 1 + s T + ( T + ω + s ) 2 s 2 ω × ln 1 + s T + ω 1 s , j = 1 , 2 ,
and
λ j ^ E B L 3 = D j + a a + b 1 ω ln 1 + ω s + T + T 2 s 2 ω ln 1 + s T ( T + ω ) 2 s 2 ω × ln 1 + s T + ω + 1 s , j = 1 , 2 .
The E-MSE of the parameter λ j , j = 1 , 2 can be obtained from (16), (18), and (19). Using the prior distribution, π 1 ( d j , c j ) is given by
E - MSE ( λ j ^ B L k ) = D MSE ( λ j ^ B L ) ( d j , c j ) π k ( d j , c j ) d c j d d j , j = 1 , 2 , k = 1 , 2 , 3 ,
where λ j ^ B L , j = 1,2, is the Bayesian estimator of λ j given by Equation (14), M S E ( λ j ^ B L ) is the mean square error of Bayesian estimator of λ j given by Equation (16) and D is the domain of d j and c j for which the prior density is decreasing in λ j .
The E-MSE of the parameter λ j , j = 1 , 2 can be obtained from (14), (18), and (19). Using the prior distribution, π 1 ( d j , c j ) is given by
E - MSE ( λ j ^ B L 1 ) = D MSE ( λ j ^ B L ) ( d j , c j ) π 1 ( d j , c j ) d c j d d j = 1 s { s T ( T + s ) I 1 j + 2 s ω I 2 j Li 2 ( ω T ) Li 2 ( ω s + T ) 1 s ω 2 I 2 j [ ω + 2 ω log ( T ) log [ T + ω ] 2 + 2 ω log ( T + ω ) log T + ω T + T log T + ω T 2 × log ( T + ω ) log T + ω ω 2 ω log ( s + T ) log ( s + T + ω ) + ω log ( s + T + ω ) 2 + 2 ω log ( s + T ) log s + T + ω ω 2 ω log ( s + T + ω ) × log 1 + ω s + T ( s + T ) log 1 + ω s + T 2 2 ω Li 2 ( T ω ) + 2 ω Li 2 ( s + T ω ) ] } ,
where Li n ( z ) = k = 1 z k k n , | z | < 1 is the polylogarithm function,
I 1 j = D j 2 + 3 a + b a + b D j + a ( 2 + 2 a + b ) ( a + b ) ( a + b + 1 ) , j = 1 , 2 ,
and
I 2 j = D j 2 + 2 a a + b D j + a ( a + 1 ) ( a + b ) ( a + b + 1 ) , j = 1 , 2 .
Similarly, under LINEX loss function, the E-MSE of λ j , j = 1 , 2 based on π 2 ( d j , c j ) and π 3 ( d j , c j ) are given, respectively, by
E - MSE ( λ j ^ B L 2 ) = 2 I 1 j s 2 s T log 1 + s T 4 I 2 j s 2 ω ω log T + ω + T log T + ω T ω log s + T + ω ( s + T ) log 1 + ω s + T ( s + T ) Li 2 ( ω T ) + ( s + T ) Li 2 ( ω s + T ) 1 4 ω ω ( 2 ( s + T ) + ω ) log T + ω 2 + 2 T ω log T + ω T T ( 2 s + T ) log T + ω T 2 + 2 ω log T + ω ω ( 2 ( s + T ) + ω ) log ( T ) + log T + ω T + ω 2 ( 2 ( s + T ) + ω ) log ( T ) log T + ω ω + 2 ω + ( 2 ( s + T ) + ω ) 2 log ( s + T ) log ( s + T + ω ) log ( s + T + ω ) 2 ( 2 ( s + T ) + ω ) log ( s + T ) log s + T + ω ω + 2 ω s T + ( 2 ( s + T ) + ω ) × log ( s + T + ω ) log 1 + ω s + T + ( s + T ) 2 log 1 + ω s + T 2 + 2 ω ( 2 ( s + T ) + ω ) Li 2 ( T ω ) Li 2 ( s + T ω ) .
E - MSE ( λ j ^ B L 3 ) = 2 I 1 j s 2 log ( 1 + s T ) s s + T + 4 I 2 j s 2 ω ω log ( T + ω ) + T log T + ω T ω log ( s + T + ω ) ( s + T ) log 1 + ω s + T T Li 2 ( ω T ) + T Li 2 ( ω s + T ) + I 2 j s 2 ω 2 ω ( 2 T + ω ) log ( T + ω ) 2 2 T ω log T + ω T + T 2 log T + ω T 2 + 2 ω log ( T + ω ) ω + ( 2 T + ω ) log ( T ) + log T + ω T + ω 2 ( 2 T + ω ) log ( T ) log T + ω ω + log ( s + T + ω ) 2 ω 2 ( 2 T + ω ) log ( s + T ) + ( 2 T + ω ) log ( s + T + ω ) ) + 2 ( 2 T + ω ) log ( s + T ) log s + T + ω ω + 2 ω ( s + T ( 2 T + ω ) log ( s + T + ω ) ) log 1 + ω s + T + ( s T ) ( s + T ) × log 1 + ω s + T 2 + 2 ω ( 2 T + ω ) Li 2 ( T ω ) + Li 2 ( s + T ω ) .

6. Properties of E-Bayesian Estimation Based on SE Loss Function

Now, the relations among λ j ^ E B S i and E-MSE ( λ j ^ E B S i ) , ( j = 1 , 2 , i = 1 , 2 , 3 ) estimations will be discussed.
I. Relations among λ j ^ E B S i , ( j = 1 , 2 , i = 1 , 2 , 3 )
Proposition 1.
Let 0 < s < T , a > 0 , b > 0 and λ j ^ E B S i , ( j = 1 , 2 , i = 1 , 2 , 3 ) be given by (20)–(22). Then, the following inequalities are
(i) 
λ j ^ E B S 3 < λ j ^ E B S 1 < λ j ^ E B S 2 .
(ii) 
lim T λ j ^ E B S 1 = lim T λ j ^ E B S 2 = lim T λ j ^ E B S 3 .
Proof. 
See Appendix A. □
II. Relations among λ j ^ E B L i , ( j = 1 , 2 , i = 1 , 2 , 3 )
Proposition 2.
Let 0 < s < T , a > 0 , b > 0 and λ j ^ E B L i , ( j = 1 , 2 , i = 1 , 2 , 3 ) be given by (26)–(28). Then, the following inequalities are
(i) 
λ j ^ E B L 3 < λ j ^ E B L 1 < λ j ^ E B L 2 .
(ii) 
lim T λ j ^ E B L 1 = lim T λ j ^ E B L 2 = lim T λ j ^ E B L 3 .
Proof. 
See Appendix A. □
III. Relations among E-MSE ( λ j ^ E B S i ) , ( j = 1 , 2 , i = 1 , 2 , 3 )
Proposition 3.
Let 0 < s < T , a > 0 , b > 0 , and E-MSE ( λ j ^ E B S i ) ( i = 1 , 2 , 3 ) be given by (23)–(25). Then, the following inequalities are
(i) 
E - MSE ( λ j ^ E B S 3 ) < E - MSE ( λ j ^ E B S 1 ) < E - MSE ( λ j ^ E B S 2 ) .
(ii) 
lim T E - MSE ( λ j ^ E B S 1 ) = lim T E - MSE ( λ j ^ E B S 2 ) = lim T E - MSE ( λ j ^ E B S 3 ) .
Proof. 
See Appendix A. □

7. Numerical Results

In this section, the numerical discussion is performed on the hypothesis that the scale parameters and common shape parameter are unknown. MATHCAD 2007 was used to perform all calculations. Numerical clarification is made to evaluate these estimates and to clarify the behavior of the proposed methods. We reanalyze the real data set analyzed by Hoel [25], and reused by Kundu et al. [26]. The simulation study used a comparison of performance of the different estimators values and different schemes.

7.1. Real Data

In this subsection, we reanalyze a real data set, which was originally reported in Hoel [25] and later by Kundu et al. [26], Pareek et al. [27], and Cramer and Schmiedt [28]. The data were obtained from a laboratory experiment in which male mice received a radiation dose of 300 roentgens at 35 days to 42 days (5–6 weeks) of age. The cause of death for each mouse was determined by reticulum cell sarcoma as cause 1 and other causes of death as cause 2; there were n = 77 observations that remain in the analysis.
Example 1.
Suppose τ = 600 , the censoring scheme D = 25 , R 1 = R 2 = = R 19 = 2 and R 20 = R 21 = = R D = 0 . The adaptive type I progressive hybrid censored sample from the original data is given as follows: (40,2), (42,2), (51,2), (62,2), (163,2), (179,2), (206,2), (222,2), (228,2), (252,2), (317,1), (324,2), (341,2), (385,2), (399,1), (407,2), (420,2), (517,2), (524,2), (536,1), (554,1), (571,1), (586,1), (586,2), (594,1).
The first component of the data denotes the lifetime and the second indicate the reason of failure. There were D = 25 , D 1 = 7 deaths due to reason 1 and D 2 = 18 deaths due to reason 2. The maximum likelihood (ML), the Bayesian estimates (BE) for Prior 0: d 1 = d 2 = c 1 = c 2 = 0 and Prior 1: d 1 = 0.5 , c 1 = 1.2 , d 2 = 0.6 , c 2 = 1.3 , and the E-Bayesian estimate (E-BE) are computed for the following initial values β = 0.018 , ω = 1.5 , s = 15 , a = 2 , b = 3 , see Table 1.
Example 2.
The same data are used, m and R i ’s are the same as before, while τ = 650 instead of τ = 600 . The adaptive type-I progressive hybrid sample from the original data in this case is given as: (40,2), (42,2), (51,2), (62,2), (163,2), (179,2), (206,2), (222,2), (228,2), (252,2), (317,1), (324,2), (341,2), (385,2), (399,1), (407,2), (420,2), (517,2), (524,2), (536,1), (554,1), (571,1), (586,1), (586,2), (594,1), (619,2), (621,1), (622,2).
Here, D = 28 , D 1 = 8 and D 2 = 20 . The MLE, BEs, and E-BE are computed for the following initial values: β = 0.018 , d 1 = 0.5 , c 1 = 1.2 , d 2 = 0.6 , c 2 = 1.3 , ω = 1.5 , s = 3 , a = 2 , b = 3 , see Table 2.
From two examples of previous real data, we noticed that:
  • Exact time τ plays an important role in estimating unknown parameters.
  • Based on the SE and LINEX loss functions for Prior 0, we note that MLE performance is very close to BEs and E-BEs.
  • The E-BEs of λ 1 and λ 2 under LINEX loss function with Prior III give smaller E-MSEs.
  • The ordering of performance of estimates of λ 1 and λ 2 in terms of minimum MSEs (E-MSEs) (from best to worst) are E-BEs with prior III, E-BEs with Prior I, E-BEs with prior II, BEs with prior 1, with BEs with prior 0 and MLs.
  • Based on prior 0, it can be seen that the BEs are quite close to the MLEs.
  • Based on LINEX loss function and minimum MSE’s (E-MSE’s), the E-BEs and BEs are better than those under SE loss function; this due to the use of SE which is symmetric loss function and based on the assumption that the loss is the same in any direction.

7.2. Simulation Study

In this subsection, we performed a simulation study to evaluate the performance of different methods presented above. The simulation was conducted according to the following steps:
  • For given values, the sample size n = 30 , 50 and different effective number of failures m = 5 , 10 and τ = 0.8 , 1.2 .
  • For given values, the parameters ( λ 1 , λ 2 , β ) = ( 0.4 , 0.6 , 1.5 ) for each case.
  • Determine three different sampling schemes:
    • Scheme 1: R 1 = = R m 1 = 0 and R m = n m ,
    • Scheme 2: R 1 = = R m 1 = 1 and R m = n 2 m + 1 , and
    • Scheme 3: R 1 = = R m 1 = R m = ( n m ) / m .
  • Determine n , m , R i s , T and the value of the parameters λ j , j = 1 , 2 .
  • If x m : m : n < τ , the progressive censoring scheme will become R 1 , R 2 , , R m , R m + 1 , , R D , where R m = R m + 1 = = R D = 0 .
  • All Bayesian estimates are calculated by two types of priests:
    • Prior 0 described the case of hyper-parameter values, d 1 = d 2 = c 1 = c 2 = 0 .
    • Prior 1 described the case of hyper-parameter values, d 1 = 0.8 , c 1 = 2 , d 2 = 0.6 , c 2 = 1 , see Table 3, Table 4, Table 5 and Table 6.
  • Generate the conventional progressive type-I censored sample from the Weibull model according to the method proposed by Kemp and Kemp [29], by using X = 1 λ k log 1 1 U 1 / β , k = 1 , 2 , where U is uniform ( 0 , 1 ) .
  • The assumption that the number of failures due to each reason of failures is at least one.
  • Under the SE loss function, the estimates λ j ^ B S , M S E ( λ j ^ B S ) , λ j ^ E B S i and E - MSE ( λ j ^ B S i ) , j = 1 , 2 , and i = 1 , 2 , 3 are computed from (11), (15), (20)–(25).
  • Using the LINEX loss function by choosing ( ω = 1.5 ) in all the cases, the estimates λ j ^ B L , M S E ( λ j ^ B L ) , λ j ^ E B L i and E - MSE ( λ j ^ B L i ) , j = 1 , 2 , i = 1 , 2 , 3 are computed from (14), (16), and (26)–(32).
  • Repeat the above steps 10,000 times. The average of all 10,000 estimated values from Steps 9–10 are respectively calculated and summarized.
  • The computational results are displayed in Table 3, Table 4, Table 5 and Table 6.
From Table 3, Table 4, Table 5 and Table 6, we have the following observations:
  • The MSE and E-MSE of the different estimators decrease as n increases.
  • For fixed n, the MSE of the different estimators decreases as m increases.
  • For fixed n and m, the MSE of the different estimators of λ 1 decrease as T increases.
  • For fixed n and m, the MSE of the different estimators of λ 2 increases as T increases.
  • As m increases for fixed n, the MSE and E-MSE decrease.
  • For fixed n and m, the MSE and E-MSE of λ 1 decrease as T increases.
  • For fixed n and m, the MSE and E-MSE of λ 2 increase as T increases.
  • The Bayesian and E-Bayesian estimates of λ 1 and λ 2 perform better than MLEs in terms of minimum MSE.
  • The E-Bayesian estimates of λ 1 and λ 2 have the minimum MSE among all other estimates.
  • The E-Bayesian estimates of λ 1 and λ 2 using prior distribution 3 perform better than other estimates in terms of minimum MSE.
  • The E-Bayesian estimates of λ 1 and λ 2 based on SE loss function under prior distribution 3 have the minimum MSE comparing with all other estimates.
  • The E-posterior risk of E-Bayesian estimation of λ 1 and λ 2 using prior distribution 3 under LINEX loss function have the minimum values among all other prior distributions.
Combining all the above results, we recommend using the E-Bayesian procedure to estimate the parameters λ 1 and λ 2 for the Adaptive Type-I Progressive Hybrid Censoring scheme based on the prior distribution 3 which performs better than other estimates in terms of minimum MSE and E-MSE.

8. Conclusions

In this article, we have investigated the E-Bayesian estimation of the parameter and the reliability function of the Weibull distribution based on A-I PHCS. The E-Bayesian estimation is considered by using three different prior distributions under two loss functions, namely the SE and LINEX loss functions. The properties of the E-Bayesian estimation as well as the E-posterior risk are also derived. We compare the performance of the E-Bayesian estimation with the maximum likelihood and Bayesian estimators via an extensive simulation study. The simulation results reveal that the E-Bayesian estimation performs better than the maximum likelihood and Bayesian estimators in terms of minimum biases and MSEs. Moreover, we analyze two real data sets for illustration purposes, and the results coincide with those in the simulation part.

Author Contributions

All authors conceived of the study, participated in its design and coordination, drafted the manuscript, participated in the sequence alignment. All authors have read and agreed to the published version of the manuscript.

Funding

This project was funded by Deanship Scientific Research (DSR), King Abdulaziz Universuty, Jeddah grant number D-334-130-1441.

Acknowledgments

The authors would like to thank the anonymous reviewers for their helpful and constructive comments that greatly contributed to improving the final version of the paper. The authors would like to extend their sincere appreciation to the Deanship of Scientific Research King Abdulaziz university for funding this work through the Research Group No. (D-334-130-1441).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of Proposition 1.
 
(i) From (22)–(24), we have
λ j ^ E B S 2 λ j ^ E B S 1 = λ j ^ E B S 1 λ j ^ E B S 3 = 1 s D j + a a + b s + 2 T s ln T + s T 2 , j = 1 , 2 .
For 1 < x < 1 , we have: ln ( 1 + x ) = x x 2 2 + x 3 3 x 4 4 + = k = 1 ( 1 ) k 1 x k k . Let x = s T , when 0 < s < T , 0 < s T < 1 , we get:
s + 2 T s ln T + s T 2 = s + 2 T s s T 1 2 s T 2 + 1 3 s T 3 1 4 s T 4 + 1 5 s T 5 2 = s T 1 2 s T 2 + 1 3 s T 3 1 4 s T 4 + 1 5 s T 5 2 + 2 s T + 2 3 s T 2 2 4 s T 3 + 2 5 s T 4 = s 2 6 T 2 s 3 6 T 3 + 3 c 4 6 T 4 2 c 5 15 T 5 + = s 2 6 T 2 1 s T + s 4 60 T 4 9 8 s T + > 0 .
According to (A1) and (A2), we have
λ j ^ E B S 2 λ j ^ E B S 1 = λ j ^ E B S 1 λ j ^ E B S 3 > 0 ,
that is,
λ j ^ E B S 3 < λ j ^ E B S 1 < λ j ^ E B S 2 , j = 1 , 2 .
(ii) From (A1) and (A2), we get
lim T λ j ^ E B S 2 λ j ^ E B S 1 = lim T λ j ^ E B S 1 λ j ^ E B S 3 = 1 s D j + a a + b lim T s 2 6 T 2 1 s T + s 4 60 T 4 9 8 c T + = 0 .
That is, lim T λ j ^ E B S 1 = lim T λ j ^ E B S 2 = lim T λ j ^ E B S 3 , j = 1 , 2 .
Thus, the proof is complete. □
Proof of Proposition 2.
 
(i) From (26)–(28), we have
λ j ^ E B L 2 λ j ^ E B L 1 = λ j ^ E B L 1 λ j ^ E B L 3 = D j + a a + b ( T + ω ) ( T + ω + s ) s 2 ω ln 1 + s T + ω T ( T + s ) s 2 ω ln 1 + s T 1 s .
For 1 < x < 1 , we have: ln ( 1 + x ) = x x 2 2 + x 3 3 x 4 4 + = k = 1 ( 1 ) k 1 x k k . Let x = s T , when 0 < s < T , 0 < s T < 1 , we get:
( T + ω ) ( T + ω + s ) s 2 ω ln 1 + s T + ω T ( T + s ) s 2 ω ln 1 + s T 1 s = ( T + ω ) 2 s 2 ω 1 + s T + ω ln 1 + s T + ω T 2 s 2 ω 1 + s T ln 1 + s T + ω 1 s = ( T + ω ) 2 s 2 ω 1 + s T + ω s T + ω 1 2 s T + w 2 + 1 3 s T + ω 3 1 4 s T + ω 4 + T 2 s 2 ω 1 + s T s T 1 2 s T 2 + 1 3 s T 3 1 4 s T 4 + 1 s = ( T + ω ) 2 s 2 ω { 1 2 s T + w 2 + 1 3 s T + ω 3 1 4 s T + ω 4 + 1 5 s T + ω 5 + 1 2 s T + w 3 + 1 3 s T + ω 4 1 4 s T + ω 5 + 1 5 s T + ω 6 } T 2 s 2 ω { 1 2 s T 2 + 1 3 s T 3 1 4 s T 4 + 1 5 s T 5 + 1 2 s T 3 + 1 3 s T 4 1 4 s T 5 + } = ( T + ω ) 2 s 2 ω 1 2 s T + w 2 1 6 s T + ω 3 + 1 12 s T + ω 4 T 2 s 2 ω 1 2 s T 2 1 6 s T 3 + 1 12 s T 4 = 1 6 1 T + ω s T 1 12 2 T + ω ( T + ω ) 2 s T 2 + > 0 .
According to (A3) and (A4), we have
λ j ^ E B L 2 λ j ^ E B L 1 = λ j ^ E B L 1 λ j ^ E B L 3 > 0 ,
that is,
λ j ^ E B L 3 < λ j ^ E B L 1 < λ j ^ E B L 2 , j = 1 , 2 .
(ii) From (A3) and (A4), we get
lim T λ j ^ E B L 2 λ j ^ E B L 1 = lim T λ j ^ E B L 1 λ j ^ E B L 3 = D j + a a + b lim T 1 6 1 T + ω s T 1 12 2 T + ω ( T + ω ) 2 s T 2 + = 0 .
That is, lim T λ j ^ E B L 1 = lim T λ j ^ E B L 2 = lim T λ j ^ E B L 3 , j = 1 , 2 .
Thus, the proof is complete. □
Proof of Proposition 3.
 
(i) From (23)–(25), we have
E - MSE ( λ j ^ E B S 3 ) E - MSE ( λ j ^ E B S 1 ) = E - MSE ( λ j ^ E B S 1 ) E - MSE ( λ j ^ E B S 2 ) = D j + a a + b s + 2 T s T ( T + s ) 2 s 2 ln 1 + s T , j = 1 , 2 .
Let x = s T , when 0 < s < T , 0 < s T < 1 , we get:
s + 2 T s T ( T + s ) 2 s 2 ln 1 + s T = s + 2 T s T ( T + s ) 2 s 2 s T 1 2 s T 2 + 1 3 s T 3 1 4 s T 4 + 1 5 s T 5 = s T 2 ( T + s ) 2 s 2 1 3 s T 3 1 4 s T 4 + 1 5 s T 5 1 6 s T 6 + 1 7 s T 7 = s 1 2 s T 3 T 2 ( T + s ) + 2 s T 3 1 4 1 5 s T + 1 6 s T 2 1 7 s T 3 + = s 1 2 s T 3 T 2 ( T + s ) + 2 s T 3 1 20 5 4 s T + 1 42 s T 2 7 6 s T + 1 72 s T 4 9 8 s T + > 0 ,
implying
E - MSE ( λ j ^ E B S 3 ) < E - MSE ( λ j ^ E B S 1 ) < E - MSE ( λ j ^ E B S 2 ) .
Thus, the proof is complete. ⋄
(ii) From (A5) and (A6), we get
lim T E - MSE ( λ j ^ E B S 2 lim T E - MSE ( λ j ^ E B S 1 ) = lim T E - MSE ( λ j ^ E B S 1 ) lim T E - MSE ( λ j ^ E B S 3 ) = D j + a a + b lim T { s 1 2 s T 3 T 2 ( T + s ) + 2 s T 3 { 1 20 5 4 s T + 1 42 s T 2 × 7 6 s T + 1 72 s T 4 9 8 s T + } } = 0 .
That is, lim T E - MSE ( λ j ^ E B S 1 ) = lim T E - MSE ( λ j ^ E B S 2 ) = lim T E - MSE ( λ j ^ E B S 3 ) .
Thus, the proof is complete. □

References

  1. Epstein, B. Truncated life-test in exponential case. Ann. Math. Stat. 1954, 25, 555–564. [Google Scholar] [CrossRef]
  2. Childs, A.; Chandrasekar, B.; Balakrishnan, N.; Kundu, D. Exact likelihood inference based on Type-I and Type-II hybrid censored samples from the exponential distribution. Ann. Inst. Stat. Math. 2003, 55, 319–330. [Google Scholar] [CrossRef]
  3. Kundu, D.; Joarder, A. Analysis of type-II progressively hybrid censored data. Comput. Stat. Data Anal. 2006, 50, 2509–2528. [Google Scholar] [CrossRef]
  4. Childs, A.; Chandrasekar, B.; Balakrishnan, N. Exact likelihood inference for an exponential parameter under progressive hybrid censoring schemes. In Statistical Models and Methods for Biomedical and Technical Systems; Vonta, F., Nikulin, M., Limnios, N., Huber-Carol, C., Eds.; Birkhäuser: Boston, MA, USA, 2008; pp. 323–334. [Google Scholar]
  5. Ng, H.K.T.; Kundu, D.; Chan, P.S. Statistical analysis of exponential lifetimes under an adaptive Type-II progressively censoring scheme. Nav. Res. Logist. 2009, 56, 687–698. [Google Scholar] [CrossRef] [Green Version]
  6. Balakrishnan, N.; Kundu, D. Hybrid censoring: Models—inferential results and applications. Comput. Stat. Data Anal. 2013, 57, 166–209. [Google Scholar] [CrossRef]
  7. Lin, C.T.; Huang, Y.L. On progressive hybrid censored exponential distribution. J. Stat. Comput. Simul. 2012, 82, 689–709. [Google Scholar] [CrossRef]
  8. Lin, C.T.; Chou, C.C.; Huang, Y.L. Inference for the Weiull distribution with progressive hybrid censoring. Comput. Stat. Data Anal. 2012, 56, 451–467. [Google Scholar] [CrossRef]
  9. Lindley, D.V. Approximate Bayes methods. Trabajos de Estadistica y de Investigacion Operativa. 1980, 31, 223–245. [Google Scholar] [CrossRef]
  10. Tierney, L.; Kadane, J.B. Accurate approximations for posterior moments and marginal densities. J. Am. Stat. Assoc. 1986, 81, 82–86. [Google Scholar] [CrossRef]
  11. Varian, H.R. A Bayesian approach to real state assessment. In Studies in Bayesian Econometrics and Statistics in Honor of Leonard J. Savage; Stephen, E.F., Zellner, A., Eds.; North-Holland: Amsterdam, The Netherlands, 1975; pp. 195–208. [Google Scholar]
  12. Zellner, A. Bayesian estimation and prediction using asymmetric loss functions. J. Am. Stat. Assoc. 1986, 81, 446–451. [Google Scholar] [CrossRef]
  13. Cox, D.R. The analysis of exponentially distributed lifetimes with two types of failure. J. R. Stat. Soc. B 1959, 21, 411–421. [Google Scholar]
  14. Crowder, M.J. Classical Competing Risks; Chapman & Hall: Boca Raton, FL, USA, 2001. [Google Scholar]
  15. Ashour, S.K.; Nassar, M.M.A. Analysis of exponential distribution under adaptive type-I progressive hybrid censored competing risks data. PJSOR 2014, 10, 229–245. [Google Scholar] [CrossRef]
  16. Ashour, S.K.; Nassar, M.M.A. Inference for Weibull distribution under adaptive type-I progressive hybrid censored competing risk data. Commun Stat. Theory Methods 2017, 46, 4756–4773. [Google Scholar] [CrossRef]
  17. Kundu, D.; Pradhan, B. Bayesian analysis of progressively censored competing risks data. Sankhya B 2011, 73, 276–296. [Google Scholar] [CrossRef]
  18. Han, M. E-Bayesian estimation and hierarchical Bayesian estimation of failure rate. Appl. Math. Model. 2009, 33, 1915–1922. [Google Scholar] [CrossRef]
  19. Han, M. E-Bayesian estimation and its E-mean-square error of the exponential distribution parameter based on complete and type I censored samples. Commun Stat. Theory Methods 2019. [Google Scholar] [CrossRef]
  20. Jaheen, Z.R.; Okasha, H.M. E-Bayesian estimation for the Burr type XII model based on type-2 censoring. Appl. Math. Model. 2011, 35, 4730–4737. [Google Scholar] [CrossRef]
  21. Azimi, R.; Yaghmaei, F.; Fasihi, B. E-Bayesian estimation based on generalized half Logistic progressive type-II censored data. IJAMS 2013, 1, 56–63. [Google Scholar] [CrossRef]
  22. Okasha, H.M. E-Bayesian estimation for the Lomax distribution based on type-II censored data. JOEMS 2014, 22, 489–495. [Google Scholar] [CrossRef] [Green Version]
  23. Okasha, H.M.; Wang, J. E-Bayesian estimation for the geometric model based on record statistics. Appl. Math. Model. 2016, 40, 658–670. [Google Scholar] [CrossRef]
  24. Abdalla, R.; Junping, L. E-Bayesian estimation for Burr-X distribution based on generalized type-I hybrid censoring scheme. AJMMS 2019, 39, 41–55. [Google Scholar] [CrossRef] [Green Version]
  25. Hoel, D.G. A representation of mortality data by competing risks. Biometrics 1972, 28, 475–488. [Google Scholar] [CrossRef] [PubMed]
  26. Kundu, D.; Kannan, N.; Balakrishnan, N. Analysis of progressively censored competing risks data. In Handbook of Statistics; Balakrishnan, N., Rao, C.R., Eds.; Elsevier: New York, NY, USA, 2004; Volume 23. [Google Scholar]
  27. Pareek, B.; Kundu, D.; Kumar, S. On progressively censored competing risks data for Weibull distribution. Comput. Stat. Data Anal. 2009, 53, 4083–4094. [Google Scholar] [CrossRef]
  28. Cramer, E.; Schmiedt, A.B. Progressively type-II censored competing risks data from Lomax distribution. Comput. Stat. Data Anal. 2011, 55, 1285–1303. [Google Scholar] [CrossRef]
  29. Kemp, C.D.; Kemp, W. Rapid generation of frequency tables. J. Appl. Stat. 1987, 36, 277–282. [Google Scholar] [CrossRef]
Table 1. The MLE, BE, E-BE, and the MSE(E-MSE).
Table 1. The MLE, BE, E-BE, and the MSE(E-MSE).
BEE-BE
λ ^ j MLELFPrior 0Prior 1Prior IPrior IIPrior III
λ ^ 1 SE0.082174560.086821120.080015690.082177760.07785363
(MSE)0.08217456 (9.8767× 10 4 )(9.8669× 10 4 )(8.6710× 10 4 )(9.1394× 10 4 )(8.2026× 10 4 )
(2.0258× 10 3 )LINEX0.081459450.086075940.079372340.081499840.07724485
(9.8718× 10 4 )(9.8623× 10 4 )(8.6652× 10 4 )(9.1340× 10 4 )(8.1064× 10 4 )
λ ^ 2 SE0.211306020.214158760.198957940.204333880.193582
(MSE)0.21130602 (2.5806× 10 3 )(2.4868× 10 3 )(2.1560× 10 3 )(2.2725× 10 3 )(2.0396× 10 3 )
(3.2831× 10 3 )LINEX0.209467150.212320660.197358260.202648250.19206827
(2.5739× 10 3 )(2.4802× 10 3 )(2.1486× 10 3 )(2.2654× 10 3 )(2.0319× 10 3 )
Table 2. The MLE, BE, E-BE, and the MSE(E-MSE).
Table 2. The MLE, BE, E-BE, and the MSE(E-MSE).
BEE-BS
λ ^ j MLELFPrior 0Prior 1Prior IPrior IIPrior III
λ ^ 1 SE0.093891910.098374670.0908090.0932620.088356
(MSE)0.09389191 (1.1219× 10 3 )(1.1177× 10 3 )(9.8386× 10 4 )(1.0369× 10 3 )(9.3072× 10 4 )
(2.0889× 10 3 )LINEX0.093075010.097530520.0900790.0924930.087665
(1.1216× 10 3 )(1.1174× 10 3 )(9.8339× 10 4 )(1.0366× 10 3 )(9.3020× 10 4 )
λ ^ 2 SE0.234729780.237256550.2205360.2264940.214579
(MSE)0.23472978 (2.7549× 10 3 )(2.7529× 10 3 )(2.3894× 10 3 )(2.5184× 10 3 )(2.2603× 10 3 )
(3.5343× 10 3 )LINEX0.232687530.235220660.2187640.2246260.212901
(2.7491× 10 3 )(2.7471× 10 3 )(2.3825× 10 3 )(2.5119× 10 3 )(2.2532× 10 3 )
Table 3. The average values of the MLE, BE, E-BE, and the MSE(E-MSE) in parentheses, under different censoring schemes and different τ ’s, for n = 30 and m = 5 .
Table 3. The average values of the MLE, BE, E-BE, and the MSE(E-MSE) in parentheses, under different censoring schemes and different τ ’s, for n = 30 and m = 5 .
BEE-BE
S λ ^ j MLELFPrior 0Prior 1Prior IPrior IIPrior III
τ = 0.8
SE0.071130.091680.064270.07020.05923
λ ^ 1 0.07113 (2.3712× 10 3 )(2.3652× 10 3 )(1.6893× 10 3 )(1.9745× 10 3 (1.4041× 10 3 )
(MSE)(2.7518× 10 3 )LINEX0.069410.08960.063480.068760.0582
1 (2.3546× 10 3 )(2.3498× 10 3 )(1.67115× 10 3 )(1.9569× 10 3 )(1.3854× 10 3 )
SE0.095520.11180.083410.090480.07634
λ ^ 2 0.09552 (3.1837× 10 3 )(3.6059× 10 3 )(2.1769× 10 3 )(2.5444× 10 3 )(1.8095× 10 3 )
(MSE)(NA)LINEX0.093210.109180.081820.088630.07502
(3.1694× 10 3 )(3.5931× 10 3 )(2.1599× 10 3 )(2.5283× 10 3 )(1.7915× 10 3 )
SE0.07140.091970.064920.070430.05941
λ ^ 1 0.0714 (2.3846× 10 3 )(2.3792× 10 3 )(1.69710× 10 3 )(1.9840× 10 3 )(1.4102× 10 3 )
(MSE)(4.1442× 10 3 )LINEX0.069670.089880.063680.068980.05838
2 (2.3681× 10 3 )(2.3639× 10 3 )(1.6789× 10 3 )(1.9665× 10 3 )(1.3915× 10 3 )
SE0.095560.111860.083420.090510.07634
λ ^ 2 0.09556 (3.1908× 10 3 )(3.1144× 10 3 )(2.1806× 10 3 )(2.549× 10 3 )(1.8119× 10 3 )
(MSE)(NA)LINEX0.093240.109230.081830.088650.07502
(3.1766× 10 3 )(3.1017× 10 3 )(2.1635× 10 3 )(2.5331× 10 3 )(1.7939× 10 3 )
SE0.071730.092430.065140.070700.05958
λ ^ 1 0.07173 (2.4134× 10 3 )(2.3135× 10 3 )(1.7133× 10 3 )(2.0045× 10 3 )(1.4219× 10 3 )
(MSE)(4.1634× 10 3 )LINEX0.069980.090310.063890.069240.05854
3 (2.3969× 10 3 )(2.3084× 10 3 )(1.6952× 10 3 )(1.9871× 10 3 )(1.4033× 10 3 )
SE0.096480.112870.084070.091240.07689
λ ^ 2 0.09648 (3.2457× 10 3 )(3.1735× 10 3 )(2.2109× 10 3 )(2.5867× 10 3 )(1.8349× 10 3 )
(MSE)(NA)LINEX0.094120.110200.082450.089360.07555
(3.2317× 10 3 )(3.1609× 10 3 )(2.1939× 10 3 )(2.5708× 10 3 )(1.8170× 10 3 )
τ = 1.2
SE0.071080.091640.064680.070170.0592
λ ^ 1 0.07108 (2.3696× 10 3 )(2.3638× 10 3 )(1.6884× 10 3 )(1.9734× 10 3 )(1.4033× 10 3 )
(MSE)(NA)LINEX0.069360.089560.063450.068730.05817
1 (2.3531× 10 3 )(2.3485× 10 3 )(1.6702× 10 3 )(1.9559× 10 3 )(1.3846× 10 3 )
SE0.095570.111840.083450.090520.07638
λ ^ 2 0.09557 (3.1853× 10 3 )(3.1674× 10 3 )(2.1779× 10 3 )(2.5455× 10 3 )(1.8103× 10 3 )
(MSE)(6.7849× 10 2 )LINEX0.093260.109220.081860.088670.07505
(3.1711× 10 3 )(3.1547× 10 3 )(2.1609× 10 3 )(2.5295× 10 3 )(1.7923× 10 3 )
SE0.071170.091760.064740.070240.05925
λ ^ 1 0.07117 (2.3770× 10 3 )(2.3725× 10 3 )(1.6925× 10 3 )(1.9787× 10 3 )(1.4064× 10 3 )
(MSE)(NA)LINEX0.069450.089670.063510.06880.05822
2 (2.3604× 10 3 )(2.3572× 10 3 )(1.6744× 10 3 )(1.9611× 10 3 )(1.3877× 10 3 )
SE0.095790.112080.08360.09070.0765
λ ^ 2 0.09579 (3.1986× 10 3 )(3.1627× 10 3 )(2.1852× 10 3 )(2.5546× 10 3 )(1.8158× 10 3 )
(MSE)(5.0019× 10 3 )LINEX0.093470.109450.08200.088830.07517
(3.1844× 10 3 )(3.1609× 10 3 )(2.1682× 10 3 )(2.5386× 10 3 )(1.7979× 10 3 )
SE0.071940.092620.06530.070880.05973
λ ^ 1 0.07194 (2.4206× 10 3 )(2.4199× 10 3 )(1.7176× 10 3 )(2.0096× 10 3 )(1.4255× 10 3 )
(MSE)(4.1995× 10 3 )LINEX0.070180.09050.064050.069410.05868
3 (2.4042× 10 3 )(2.4041× 10 3 )(1.6995× 10 3 )(1.9922× 10 3 )(1.4068× 10 3 )
SE0.096280.112680.083920.091080.07675
λ ^ 2 0.09628 (3.2392× 10 3 )(3.1674× 10 3 )(2.2069× 10 3 )(2.5822× 10 3 )(1.8317× 10 3 )
(MSE)(NA)LINEX0.093930.110010.08230.08920.07541
(3.2252× 10 3 )(3.1549× 10 3 )(2.2189× 10 3 )(2.5662× 10 3 )(1.8138× 10 3 )
Table 4. The average values of the MLE, BE, E-BE, and the MSE(E-MSE) in parentheses, under different censoring schemes and different τ ’s, for n = 30 and m = 10 .
Table 4. The average values of the MLE, BE, E-BE, and the MSE(E-MSE) in parentheses, under different censoring schemes and different τ ’s, for n = 30 and m = 10 .
BEE-BE
S λ ^ j MLELFPrior 0Prior 1Prior IPrior IIPrior III
τ = 0.8
SE0.142860.158950.119660.129820.10951
λ ^ 1 0.14286 (4.7681× 10 3 )(4.6730× 10 3 )(3.1266× 10 3 )(3.6549× 10 3 )(2.5983× 10 3 )
(MSE)(7.3873× 10 3 )LINEX0.13940.155340.117380.127150.10761
1 (4.7609× 10 3 )(4.6667× 10 3 )(3.1126× 10 3 )(3.6429× 10 3 )(2.5824× 10 3 )
SE0.190860.204070.156430.169710.14316
λ ^ 2 0.19086 (6.3693× 10 3 )(6.2941× 10 3 )(4.0869× 10 3 )(4.7775× 10 3 )(3.3965× 10 3 )
(MSE)(7.6945× 10 3 )LINEX0.186240.199280.153450.166220.14067
(6.3715× 10 3 )(6.2940× 10 3 )(4.0771× 10 3 )(4.7708× 10 3 )(3.3833× 10 3 )
SE0.143650.159750.12020.130430.10997
λ ^ 1 0.14365 (4.8142× 10 3 )(4.0175× 10 3 )(3.1512× 10 3 )(3.6854× 10 3 )(2.6171× 10 3 )
(MSE)(5.1985× 10 3 )LINEX0.140160.15610.11790.127750.10805
2 (4.8073× 10 3 )(4.0115× 10 3 )(3.1374× 10 3 )(3.6735× 10 3 )(2.6012× 10 3 )
SE0.191460.204710.15680.170140.14345
λ ^ 2 0.19146 (6.4159× 10 3 )(6.3375× 10 3 )(4.1103× 10 3 )(4.8070× 10 3 )(3.4136× 10 3 )
(MSE)(2.1338× 10 3 )LINEX0.18680.199890.153790.166640.14095
(6.4185× 10 3 )(6.3215× 10 3 )(4.1005× 10 3 )(4.8005× 10 3 )(3.4005× 10 3 )
SE0.143980.160120.120380.130660.1101
λ ^ 1 0.14398 (4.8463× 10 3 )(4.0498× 10 3 )(3.1672× 10 3 )(3.7059× 10 3 )(2.6286× 10 3 )
(MSE)(5.6554× 10 3 )LINEX0.140460.156450.118070.127960.10817
3 (4.8395× 10 3 )(4.0439× 10 3 )(3.1535× 10 3 )(3.6942× 10 3 )(2.6128× 10 3 )
SE0.192610.205870.157570.171020.14411
λ ^ 2 0.19261 (6.4829× 10 3 )(6.4238× 10 3 )(4.1455× 10 3 )(4.8505× 10 3 )(3.4404× 10 3 )
(MSE)(6.6906× 10 3 )LINEX0.18790.20100.154540.167490.14159
(6.4859× 10 3 )(6.4183× 10 3 )(4.1358× 10 3 )(4.8442× 10 3 )(3.4274× 10 3 )
τ = 1.2
SE0.143440.159490.120110.13030.10992
λ ^ 1 0.14344 (4.7873× 10 3 )(4.6898× 10 3 )(3.1381× 10 3 )(3.6683× 10 3 )(2.6079× 10 3 )
(MSE)(3.3292× 10 3 )LINEX0.139970.155870.117820.127620.10801
1 (4.7802× 10 3 )(4.6837× 10 3 )(3.1242× 10 3 )(3.6564× 10 3 )(2.5919× 10 3 )
SE0.190280.203510.155990.169220.14275
λ ^ 2 0.19028 (6.3498× 10 3 )(6.2721× 10 3 )(4.0753× 10 3 )(4.7639× 10 3 )(3.3868× 10 3 )
(MSE)(7.3309× 10 3 )LINEX0.185670.198730.153010.165750.14027
(6.3519× 10 3 )(6.2656× 10 3 )(4.0653× 10 3 )(4.7571× 10 3 )(3.3735× 10 3 )
SE0.144030.160110.120490.130750.11023
λ ^ 1 0.14403 (4.8272× 10 3 )(4.0289× 10 3 )(3.1589× 10 3 )(3.6945× 10 3 )(2.6235× 10 3 )
(MSE)(5.0914× 10 3 )LINEX0.140530.156450.118180.128050.10831
2 (4.8203× 10 3 )(4.0229× 10 3 )(3.1452× 10 3 )(3.6827× 10 3 )(2.6077× 10 3 )
SE0.19110.204360.156520.169840.14319
λ ^ 2 0.1911 (6.4044× 10 3 )(6.5267× 10 3 )(4.1033× 10 3 )(4.7989× 10 3 )(3.4078× 10 3 )
(MSE)(6.3308× 10 3 )LINEX0.186450.199550.153520.166340.1407
(6.4069× 10 3 )(6.5306× 10 3 )(4.0935× 10 3 )(4.7923× 10 3 )(3.3946× 10 3 )
SE0.143930.160080.120340.130620.11006
λ ^ 1 0.14393 (4.8448× 10 3 )(4.0485× 10 3 )(3.1663× 10 3 )(3.7048× 10 3 )(2.6278× 10 3 )
(MSE)(3.5418× 10 3 )LINEX0.140410.156410.118030.127920.10814
3 (4.8380× 10 3 )(4.0427× 10 3 )(3.1525× 10 3 )(3.6931× 10 3 )(2.6119× 10 3 )
SE0.192670.205930.157610.171070.14415
λ ^ 2 0.19267 (6.4850× 10 3 )(6.4747× 10 3 )(4.1467× 10 3 )(4.8519× 10 3 )(3.4414× 10 3 )
(MSE)(6.8266× 10 3 )LINEX0.187960.201060.154580.167530.14163
(6.4880× 10 3 )(6.4713× 10 3 )(4.1370× 10 3 )(4.8457× 10 3 )(3.4284× 10 3 )
Table 5. The average values of the MLE, BE, E-BE, and the MSE(E-MSE) in parentheses, under different censoring scheme and different τ ’s, for n = 50 and m = 5 .
Table 5. The average values of the MLE, BE, E-BE, and the MSE(E-MSE) in parentheses, under different censoring scheme and different τ ’s, for n = 50 and m = 5 .
BEE-BE
S λ ^ j MLELFPrior 0Prior 1Prior IPrior IIPrior III
τ = 0.8
SE0.042730.056460.042670.045060.04029
λ ^ 1 0.04273 (8.5401× 10 4 )(8.3850× 10 4 )(7.2425× 10 4 )(8.0513× 10 4 )(6.4337× 10 4 )
(MSE)(1.4851× 10 3 )LINEX0.04210.055660.042140.044470.03981
1 (8.5247× 10 4 )(8.6570× 10 4 )(7.2259× 10 4 )(8.0354× 10 4 )(6.4164× 10 4 )
SE0.057190.067820.054830.05790.05177
λ ^ 2 0.05719 (1.4270× 10 3 )(1.3287× 10 3 )(7.3054× 10 4 )(8.0344× 10 4 )(7.2664× 10 4 )
(MSE)(NA)LINEX0.056350.066840.054150.057140.05116
(1.1235× 10 3 )(1.0968× 10 3 )(7.29066× 10 4 )(8.0151× 10 4 )(7.2505× 10 4 )
SE0.042750.05650.042690.04500.0403
λ ^ 1 0.04275 (8.5542× 10 4 )(8.456× 10 4 )(7.2528× 10 4 )(8.0635× 10 4 )(6.4422× 10 4 )
(MSE)(NA)LINEX0.042130.05570.042160.044490.03983
2 (8.5388× 10 4 )(8.4756× 10 4 )(7.2362× 10 4 )(8.0476× 10 4 )(6.4248× 10 4 )
SE0.057270.067910.05490.057980.05183
λ ^ 2 0.05727 (1.4565× 10 3 )(1.3320× 10 3 )(9.3260× 10 4 )(1.0368× 10 3 )(8.2838× 10 4 )
(MSE)(2.4449× 10 3 )LINEX0.056430.066930.054220.057210.05122
(2.1264× 10 3 )(1.3129× 10 3 )(9.3113× 10 4 )(1.0175× 10 3 )(8.2679× 10 4 )
SE0.043160.057010.043040.045470.04061
λ ^ 1 0.04316 (1.7145× 10 3 )(1.1063× 10 3 )(7.3691× 10 4 )(8.1991× 10 4 )(6.5391× 10 4 )
(MSE)(NA)LINEX0.042520.05620.04250.044860.04013
3 (1.7172× 10 3 )(1.0871× 10 3 )(7.3526× 10 4 )(8.1833× 10 4 )(6.5218× 10 4 )
SE0.057770.06850.055310.058430.05219
λ ^ 2 0.05777 (1.6633× 10 3 )(1.3556× 10 3 )(9.4698× 10 4 )(1.0536× 10 3 )(8.4033× 10 4 )
(MSE)(2.4824× 10 3 )LINEX0.056920.067510.054610.057660.05157
(1.1471× 10 3 )(1.3366× 10 3 )(9.4553× 10 4 )(1.0343× 10 3 )(8.3876× 10 4 )
τ = 1.2
SE0.042530.056270.042510.044890.04013
λ ^ 1 0.04253 (8.5009× 10 3 )(1.0814× 10 3 )(7.2145× 10 4 )(8.0201× 10 4 )(6.4088× 10 4 )
(MSE)(9.0935× 10 3 )LINEX0.041910.055480.041980.044290.03966
1 (8.4853× 10 3 )(1.0621× 10 3 )(7.1978× 10 4 )(8.0041× 10 4 )(6.3914× 10 4 )
SE0.057380.068010.05500.058080.05192
λ ^ 2 0.05738 (1.1466× 10 3 )(1.1325× 10 3 )(9.3335× 10 4 )(1.0376× 10 3 )(8.2913× 10 4 )
(MSE)(NA)LINEX0.056540.067030.054310.057310.05131
(1.1274× 10 3 )(1.1135× 10 3 )(9.3188× 10 4 )(1.0182× 10 3 )(8.2755× 10 4 )
SE0.042780.056520.042710.04510.04032
λ ^ 1 0.04278 (8.5592× 10 3 )(1.0873× 10 3 )(7.2564× 10 4 )(8.0674× 10 4 )(6.4453× 10 4 )
(MSE)(1.4880× 10 3 )LINEX0.042150.055720.042180.044510.03985
2 (8.5437× 10 3 )(1.0680× 10 3 )(7.2397× 10 4 )(8.0515× 10 4 )(6.4279× 10 4 )
SE0.057250.067890.054890.057960.05181
λ ^ 2 0.05725 (1.4526× 10 3 )(1.3316× 10 3 )(9.3232× 10 4 )(1.0365× 10 3 )(8.2813× 10 4 )
(MSE)(NA)LINEX0.056410.066920.05420.05720.0512
(1.3603× 10 3 )(1.3126× 10 3 )(9.3085× 10 4 )(1.0172× 10 3 )(8.2654× 10 4 )
SE0.043140.056990.043020.045450.04059
λ ^ 1 0.04314 (8.7087× 10 3 )(1.1058× 10 3 )(7.3649× 10 4 )(8.1944× 10 4 )(6.5355× 10 4 )
(MSE)(NA)LINEX0.04250.056170.042480.044840.04011
3 (8.6934× 10 3 )(1.0965× 10 3 )(7.3484× 10 4 )(8.1786× 10 4 )(6.5182× 10 4 )
SE0.057790.068520.055330.058450.0522
λ ^ 2 0.05779 (1.1665× 10 3 )(1.1557× 10 3 )(9.4711× 10 4 )(1.0538× 10 3 )(8.4046× 10 4 )
(MSE)(2.4932× 10 3 )LINEX0.056930.067520.054630.057670.05158
(1.1473× 10 3 )(1.1368× 10 3 )(9.4566× 10 4 )(1.0344× 10 3 )(8.3888× 10 4 )
Table 6. The average values of the MLE, BE, E-BE, and the MSE(E-MSE) in parentheses, under different censoring scheme and different τ ’s, for n = 50 and m = 10 .
Table 6. The average values of the MLE, BE, E-BE, and the MSE(E-MSE) in parentheses, under different censoring scheme and different τ ’s, for n = 50 and m = 10 .
BEE-BE
S λ ^ j MLELFPrior 0Prior 1Prior IPrior IIPrior III
τ =0.8
SE0.085910.097990.07900.083420.07458
λ ^ 1 0.08591 (1.7174× 10 3 )(1.6835× 10 3 )(1.3411× 10 3 )(1.4909× 10 3 )(1.1913× 10 3 )
(MSE)(1.7436× 10 3 )LINEX0.084650.09660.078010.082320.0737
1 (1.6991× 10 3 )(1.6655× 10 3 )(1.3222× 10 3 )(1.4722× 10 3 )(1.1721× 10 3 )
SE0.113980.123510.102620.108360.09687
λ ^ 2 0.11398 (2.27844× 10 3 )(2.2205× 10 3 )(1.7419× 10 3 )(1.9364× 10 3 )(1.5473× 10 3 )
(MSE)(4.7343× 10 3 )LINEX0.112310.121730.101330.106930.09573
(2.2614× 10 3 )(2.2037× 10 3 )(1.7237× 10 3 )(1.9186× 10 3 )(1.5287× 10 3 )
SE0.086130.098220.079170.083610.07473
λ ^ 1 0.08613 (1.7261× 10 3 )(1.6926× 10 3 )(1.3469× 10 3 )(1.4976× 10 3 )(1.1691× 10 3 )
(MSE)(1.7476× 10 3 )LINEX0.084860.096830.078170.082500.07385
2 (1.7078× 10 3 )(1.6747× 10 3 )(1.3279× 10 3 )(1.4789× 10 3 )(1.1769× 10 3 )
SE0.114270.123820.102840.10860.09707
λ ^ 2 0.11427 (2.290× 10 3 )(2.2325× 10 3 )(1.7495× 10 3 )(1.9453× 10 3 )(1.5537× 10 3 )
(MSE)(2.7584× 10 3 )LINEX0.112590.122030.101550.107170.09592
(2.2729× 10 3 )(2.2158× 10 3 )(1.7313× 10 3 )(1.9275× 10 3 )(1.5351× 10 3 )
SE0.086570.098740.07950.083990.07501
λ ^ 1 0.08657 (1.7484× 10 3 )(1.7167× 10 3 )(1.3615× 10 3 )(1.5149× 10 3 )(1.2082× 10 3 )
(MSE)(2.2551× 10 3 )LINEX0.085290.097330.07850.082870.07412
3 (1.7301× 10 3 )(1.6988× 10 3 )(1.3427× 10 3 )(1.4963× 10 3 )(1.1890× 10 3 )
SE0.115380.124970.103690.109550.09784
λ ^ 2 0.11538 (2.3300× 10 3 )(2.2738× 10 3 )(1.7759× 10 3 )(1.9759× 10 3 )(1.5758× 10 3 )
(MSE)(2.7765× 10 3 )LINEX0.113660.123150.102390.108090.09668
(2.3130× 10 3 )(2.2572× 10 3 )(1.7577× 10 3 )(1.9582× 10 3 )(1.5573× 10 3 )
τ = 1.2
SE0.085610.097690.078740.083150.07434
λ ^ 1 0.08561 (1.7113× 10 3 )(1.6778× 10 3 )(1.3367× 10 3 )(1.4859× 10 3 )(1.1874× 10 3 )
(MSE)(NA)LINEX0.084350.096310.077760.082050.07346
1 (1.6929× 10 3 )(1.6598× 10 3 )(1.3178× 10 3 )(1.4673× 10 3 )(1.1682× 10 3 )
SE0.114290.12380.102870.108620.09711
λ ^ 2 0.11429 (2.2844× 10 3 )(2.2262× 10 3 )(1.7462× 10 3 )(1.9412× 10 3 )(1.5512× 10 3 )
(MSE)(2.7105× 10 3 )LINEX0.112610.122020.101580.107190.09597
(2.2673× 10 3 )(2.2095× 10 3 )(1.7279× 10 3 )(1.9235× 10 3 )(1.5325× 10 3 )
SE0.086130.098220.079170.083610.07473
λ ^ 1 0.08613 (1.7261× 10 3 )(1.6926× 10 3 )(1.3469× 10 3 )(1.4976× 10 3 )(1.1961× 10 3 )
(MSE)(2.4756× 10 3 )LINEX0.084860.096830.078170.08250.07385
2 (1.7078× 10 3 )(1.6747× 10 3 )(1.3279× 10 3 )(1.4789× 10 3 )(1.1769× 10 3 )
SE0.114270.123820.102840.10860.09707
λ ^ 2 0.11427 (2.2900× 10 3 )(2.2325× 10 3 )(1.7495× 10 3 )(1.9453× 10 3 )(1.5537× 10 3 )
(MSE)(2.7584× 10 3 )LINEX0.112590.122030.10150.107170.09592
(2.2729× 10 3 )(2.2158× 10 3 )(1.7313× 10 3 )(1.9275× 10 3 )(1.5351× 10 3 )
SE0.086420.098590.079370.083850.07489
λ ^ 1 0.08642 (1.7454× 10 3 )(1.7139× 10 3 )(1.3594× 10 3 )(1.5125× 10 3 )(1.2062× 10 3 )
(MSE)(1.9812× 10 3 )LINEX0.085130.097180.078370.082740.0740
3 (1.7271× 10 3 )(1.659× 10 3 )(1.3405× 10 3 )(1.4939× 10 3 )(1.1871× 10 3 )
SE0.115540.125130.103830.109690.09797
λ ^ 2 0.11534 (2.3334× 10 3 )(2.2771× 10 3 )(1.7783× 10 3 )(1.9786× 10 3 )(1.5779× 10 3 )
(MSE)(2.6643× 10 3 )LINEX0.113820.123310.102520.108230.09681
(2.31647× 10 3 )(2.2605× 10 3 )(1.7601× 10 3 )(1.9609× 10 3 )(1.5594× 10 3 )

Share and Cite

MDPI and ACS Style

Okasha, H.; Mustafa, A. E-Bayesian Estimation for the Weibull Distribution under Adaptive Type-I Progressive Hybrid Censored Competing Risks Data. Entropy 2020, 22, 903. https://0-doi-org.brum.beds.ac.uk/10.3390/e22080903

AMA Style

Okasha H, Mustafa A. E-Bayesian Estimation for the Weibull Distribution under Adaptive Type-I Progressive Hybrid Censored Competing Risks Data. Entropy. 2020; 22(8):903. https://0-doi-org.brum.beds.ac.uk/10.3390/e22080903

Chicago/Turabian Style

Okasha, Hassan, and Abdelfattah Mustafa. 2020. "E-Bayesian Estimation for the Weibull Distribution under Adaptive Type-I Progressive Hybrid Censored Competing Risks Data" Entropy 22, no. 8: 903. https://0-doi-org.brum.beds.ac.uk/10.3390/e22080903

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