Next Article in Journal
Operator Methods of the Maximum Principle in Problems of Optimization of Quantum Systems
Previous Article in Journal
Factors Influencing Physicians Migration—A Case Study from Romania
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Estimating the Parameters of the Beta Inverted Exponential Distribution under Type-II Censored Samples

Department of Statistics, College of Science, University of Jeddah, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Submission received: 13 November 2021 / Revised: 1 January 2022 / Accepted: 27 January 2022 / Published: 5 February 2022
(This article belongs to the Section Probability and Statistics)

Abstract

:
This article aims to consider estimating the unknown parameters, survival, and hazard functions of the beta inverted exponential distribution. Two methods of estimation were used based on type-II censored samples: maximum likelihood and Bayes estimators. The Bayes estimators were derived using an informative gamma prior distribution under three loss functions: squared error, linear exponential, and general entropy. Furthermore, a Monte Carlo simulation study was carried out to compare the performance of different methods. The potentiality of this distribution is illustrated using two real-life datasets from difference fields. Further, a comparison between this model and some other models was conducted via information criteria. Our model performs the best fit for the real data.

1. Introduction

In many life-testing and reliability studies for engineering or medical sciences, the information on failure times for all experimental units may not be obtained ultimately by the experimenter. Due to this, there are many situations in which it is pre-planned to remove units before failure, and these obtained data are called censored data. The most common censoring schemes in life-testing experiments are type-I and type-II censoring schemes. The type-II censoring scheme is used often in toxicology experiments and life-testing applications, where it has proven to save time and money. Many authors have addressed Bayesian and non-Bayesian estimations based on type-II censored samples or different types of samples, including [1], who derived maximum likelihood estimation (MLE) and Bayes estimation under different types of loss functions for exponentiated Weibull distribution based on type-II censored samples. Prakash [2] discussed the properties of the Bayes estimator and the minimax estimator of the parameter of the inverted exponential distribution. The moments of the lower record value and the estimation of the parameter were presented, based on a series of observed record values by the maximum likelihood (ML). Furthermore, Dey and Dey [3] derived the MLE of the generalized inverted exponential distribution parameters in the case of the progressive type-II censoring scheme with binomial removals. Singh and Goel [4] studied a three-parameter beta inverted exponential distribution (BIED). They derived the non-central moments, inverse moments, moment-generating function, inverse-moment-generating function, and mode. Furthermore, they examined the distributional properties of order statistics. Moreover, a statistical inference about the distribution parameters based on a complete sample was investigated. Garg et al. [5] studied the MLE of the parameters and the expected Fisher information under a random censoring model of the generalized inverted exponential distribution. Bakoban and Abu-Zinadah [6] considered the four-parameter beta generalized inverted exponential distribution for complete samples. In their research, the MLE, the Fisher information matrix, and the confidence interval were found. Besides that, the Monte Carlo simulation was discussed to illustrate the theoretical results of the estimation. Finally, applications on real datasets were provided. Aldahlan [7] applied the ML method to estimate the inverse Weibull inverse exponential distribution parameters. One real dataset about time between failures for repairable items was applied. This article focuses on estimation methods based on the type-II censoring scheme. Two estimation methods were used to estimate the unknown parameters for the beta inverted exponential distribution (BIED): MLE and Bayes estimation. The proposed distribution has three parameters (scale parameter λ and shape parameters α and β ). The cumulative distribution function (CDF) and the probability density function (PDF) of BIED, respectively, are:
F ( x ) = 1 B ( α , β ) 0 e λ x ω α 1 ( 1 ω ) β 1 d ω , x > 0 , α , β , λ > 0 ,
and:
f ( x ) = λ x 2 B ( α , β ) e α λ x [ 1 e λ x ] β 1 , x > 0 , α , β , λ > 0 ,
where B ( α , β ) = 0 1 w α 1 ( 1 w ) β 1 d w is the beta function. Equation (1) could also be written as a regularized incomplete beta function:
F ( x ) = I e λ x ( α , β ) = B ( e λ x , α , β ) B ( α , β ) , x > 0 , α , β , λ > 0 ,
where B ( y , α , β ) is the incomplete beta function, such that:
B ( y , α , β ) = 0 y ω α 1 ( 1 ω ) β 1 d ω , 0 y 1 , α β > 0 .
The inverse of CDF is called the quantile function and is given by:
x = Q ( u ) = λ log [ I u 1 ( α , β ) ] , 0 < u < 1 .
The survival and hazard function of the BIED, respectively, are given by:
S ( x ) = 1 I e λ x ( α , β ) , x > 0 , α , β , λ > 0 ,
and:
h ( x ) = λ x 2 B ( α , β ) e α λ x [ 1 e λ x ] β 1 1 I e λ x ( α , β ) , x > 0 , α , β , λ > 0 .
The layout of this article is as follows: In Section 2, the estimation of the unknown parameters for the BIED under type-II censored samples is introduced. A simulation study is discussed in Section 3. In Section 4, an application with real data is provided. Finally, the conclusion is given in Section 5.

2. Method of Estimation

In this section, we derive the ML and Bayesian estimators for the unknown parameters of the BIED based on type-II censored samples.

2.1. Maximum Likelihood Estimation

Assume that X 1 , X 2 , , X n are a random sample from the BIED and the order statistics of this sample are X 1 : n < X 2 : n < < X n : n . The rth observations were chosen in type-II censored sample ( r < n ) . The likelihood function of the type-II censored sample is given by (see [8]):
L ( Θ ̲ | x ̲ ) = n ! ( n r ) ! i = 1 r f ( x i ) [ 1 F ( x r ) ] n r ,
By substituting Equations (2) and (3) in (8), the likelihood function for the vector Θ ̲ = ( α , β , λ ) is given by:
L ( Θ ̲ | x ̲ ) = n ! ( n r ) ! λ r [ B ( α , β ) ] r i = 1 r 1 x i 2 e i = 1 r α λ x i i = 1 r [ 1 e λ x i ] β 1 1 I e λ x r ( α , β ) ( n r ) .
Then, the log-likelihood function can be expressed as follows:
= log n ! ( n r ) ! + r log λ r log [ B ( α , β ) ] 2 i = 1 r log x i i = 1 r α λ x i + ( β 1 ) i = 1 r log [ 1 e λ x i ] + ( n r ) log 1 I e λ x r ( α , β ) ,
The partial derivatives of the log-likelihood function with respect to α , β , and λ are given as:
α = r B ( α , β ) α B ( α , β ) i = 1 r λ x i + ( n r ) 1 I e λ x r ( α , β ) 1 I e λ x r ( α , β ) α ,
where:
α B ( α , β ) = Γ ( α ) Γ ( β ) Γ ( α + β ) Γ ( α ) Γ ( β ) α Γ ( α + β ) [ Γ ( α + β ) ] 2 = B ( α , β ) [ ψ ( α ) ψ ( α + β ) ] ,
According to Abramowitz and Stegun [9], ψ ( z ) = d [ log Γ ( z ) ] d z = Γ ( z ) Γ ( z ) is called the Psi (or digamma) function. Then:
α = r [ ψ ( α ) ψ ( α + β ) ] i = 1 r λ x i + ( n r ) 1 I e λ x r ( α , β ) 1 I e λ x r ( α , β ) α ,
By using the Leibniz integral rule to find 1 I e λ x r ( α , β ) α (see [9]), then:
α = r [ ψ ( α ) ψ ( α + β ) ] i = 1 r λ x i + ( n r ) 1 I e λ x r ( α , β ) { 1 B ( α , β ) × 0 e λ x r w α 1 ( 1 w ) β 1 log ( w ) d w I e λ x r ( α , β ) B ( α , β ) [ ψ ( α ) ψ ( α + β ) ] } .
Let U = ( 1 W ) . Then:
α = r [ ψ ( α ) ψ ( α + β ) ] i = 1 r λ x i + ( n r ) 1 I e λ x r ( α , β ) { 1 B ( α , β ) × 1 e λ x r 1 ( 1 u ) α 1 u β 1 log ( 1 u ) d u I e λ x r ( α , β ) B ( α , β ) [ ψ ( α ) ψ ( α + β ) ] } .
By using  log ( 1 u ) = k = 1 u k k and taking Y = 1 U , we obtain the partial derivatives of α :
α = r [ ψ ( α ) ψ ( α + β ) ] i = 1 r λ x i + ( n r ) 1 I e λ x r ( α , β ) { 1 B ( α , β ) × k = 1 1 k B ( e λ x r ; α , β + k ) + I e λ x r ( α , β ) [ ψ ( α ) ψ ( α + β ) ] } .
Equation (12) can be rewritten as:
α = r [ ψ ( α ) ψ ( α + β ) ] i = 1 r λ x i + ( n r ) 1 I e λ x r ( α , β ) ( Γ ( α ) β ( α ) ( e λ x r ) α × 3 F 2 ( α , α , 1 β ; α + 1 , α + 1 ; e λ x r ) [ log ( e λ x r ) ψ ( α ) + ψ ( α + β ) ] I e λ x r ( α , β ) ) ,
where p F q ( a 1 , , a p ; b 1 , , b q ; z ) = p F q ( a 1 , , a p ; b 1 , , b q ; z ) Γ ( b 1 ) Γ ( b q ) is called the regularized hypergeometric function and p F q ( a 1 , , a p ; b 1 , , b q ; z ) = k = 0 ( a 1 ) k , , ( a p ) k ( b 1 ) k , , ( b q ) k z k k ! is the generalized hypergeometric function, and ( a ) n = a ( a + 1 ) , , ( a + n 1 ) = Γ ( a + n ) Γ ( a ) denotes the ascending factorial (Pochhammer symbol) (see [10]).
Next,
β = r B ( α , β ) β B ( α , β ) + i = 1 r log [ 1 e λ x i ] + ( n r ) 1 I e λ x r ( α , β ) 1 I e λ x r ( α , β ) β ,
where:
β B ( α , β ) = Γ ( α ) Γ ( β ) Γ ( α + β ) Γ ( α ) Γ ( β ) β Γ ( α + β ) [ Γ ( α + β ) ] 2 = B ( α , β ) [ ψ ( β ) ψ ( α + β ) ] ,
then:
β = r [ ψ ( β ) ψ ( α + β ) ] + i = 1 r log [ 1 e λ x i ] + ( n r ) 1 I e λ x r ( α , β ) 1 I e λ x r ( α , β ) β , = r [ ψ ( β ) ψ ( α + β ) ] + i = 1 r log [ 1 e λ x i ] + ( n r ) 1 I e λ x r ( α , β ) { 1 B ( α , β ) × [ 0 e λ x r w α 1 ( 1 w ) β 1 log ( 1 w ) d w I e λ x r ( α , β ) B ( α , β ) × [ ψ ( β ) ψ ( α + β ) ] ] } .
By using  log ( 1 u ) = k = 1 u k k , we obtain:
β = r [ ψ ( β ) ψ ( α + β ) ] + i = 1 r log [ 1 e λ x i ] + ( n r ) 1 I e λ x r ( α , β ) × 1 B ( α , β ) k = 1 1 k B ( e λ x r ; α + k , β ) + I e λ x r ( α , β ) [ ψ ( β ) ψ ( α + β ) ] .
Equation (15) can be rewritten as:
β = r [ ψ ( β ) ψ ( α + β ) ] + i = 1 r log [ 1 e λ x i ] + ( n r ) 1 I e λ x r ( α , β ) { ( 1 e λ x r ) β × Γ ( β ) α ( β ) 3 F 2 ( β , β , 1 α ; 1 + β , 1 + β ; 1 e λ x r ) I ( 1 e λ x r ) ( β , α ) × [ log ( 1 e λ x r ) + ψ ( β ) ψ ( α + β ) ] } .
and:
λ = r λ i = 1 r α x i + ( β 1 ) i = 1 r e λ x i [ 1 e λ x i ] x i + ( n r ) 1 I e λ x r ( α , β ) 1 I e λ x r ( α , β ) λ ,
λ = r λ i = 1 r α x i + ( β 1 ) i = 1 r x i 1 ( e λ x i 1 ) 1 + ( n r ) 1 I e λ x r ( α , β ) e α λ x r [ 1 e λ x r ] β 1 x r B ( α , β ) .
After equating Equations (13), (16) and (18) with zero and solving them, simultaneously, the MLE of α , β , and λ could be found using the Newton–Raphson method via Mathematica 11. Furthermore, the invariance property of the ML is used to estimate S ( x 0 ) and h ( x 0 ) .

2.2. Bayes Estimation

Bayes estimators for the BIED are obtained based on type-II censored samples in this subsection. Singh and Goel [4] derived the Bayes estimators for the BIED based on complete samples under the SE loss function. They considered the gamma prior distribution for the unknown BIED parameters. The prior distribution is denoted by π ( θ ) , which tells us what is known about θ without observing the data. Bayes theorem is based on the posterior distribution, which is defined as π * ( θ | x ̲ ) and given by (see [11]):
π * ( θ | x ̲ ) = L ( x ̲ | θ ) π ( θ ) L ( x ̲ | θ ) π ( θ ) d θ ,
where θ is continuous and L ( x ̲ | θ ) is the likelihood function. Furthermore, Equation (19) could be written as:
π * ( θ | x ̲ ) = k L ( x ̲ | θ ) π ( θ ) ,
where k is called the normalizing constant, necessary to ensure that the posterior distribution π * ( θ | x ̲ ) integrates or sums to one.
Here, we derive the Bayes estimates for α , β , and λ under three types of loss functions: squared error (SE), linear exponential (LINEX), and general entropy (GE). Moreover, four cases are considered first when α is unknown, while β and λ are known. A second case is when β is unknown, while α , λ are known. A third case is when the scale parameter λ is unknown. Finally, a fourth case is when both β and λ are unknown. Two techniques are used to compute the estimates: the standard Bayes and importance sampling techniques for the first three cases. The last case is computed via the importance sampling technique.

2.2.1. Case 1: Bayes Estimators When α Is Unknown

Assume α is unknown and has the following prior distribution α G a m m a ( a , b ) ; thus, the prior for α is given by:
π ( α ) = b a Γ ( a ) α a 1 e b α , α > 0 , a , b > 0 .
By combining (9) and (21), the posterior distribution of the unknown parameter α is given by:
π * ( α | x ̲ ) = k 1 α a 1 Γ ( α + β ) Γ ( α ) r e i = 1 r λ x i + b α e ( n r ) log 1 I e λ x r ( α , β ) ,
where k 1 is the normalizing constant, defined as:
k 1 1 = 0 α a 1 Γ ( α + β ) Γ ( α ) r e i = 1 r λ x i + b α e ( n r ) log 1 I e λ x r ( α , β ) d α .
Therefore, the Bayes estimator of α , denoted by φ ( α ) , is obtained under three types of loss functions and two techniques as follows.
  • SE Loss Function
The symmetric loss function SE is defined as:
φ ^ S E ( θ ) = E ( φ ( θ ) | x ̲ ) = φ ( θ ) π * ( θ | x ̲ ) d θ
Then, the Bayes estimator of φ ( α ) under the SE loss function, denoted by φ ^ S S E C ( α ) , and can be found using Equations (24) and (22).
φ ^ S S E C ( α ) = k 1 0 φ ( α ) α a 1 Γ ( α + β ) Γ ( α ) r e i = 1 r λ x i + b α e ( n r ) log 1 I e λ x r ( α , β ) d α .
where k 1 1 is defined in Equation (23).
ii.
LINEX Loss Function
Varian [12] proposed the LINEX loss function as follows:
φ ^ L E ( θ ) = 1 τ log E θ ( e τ φ ( θ ) | x ̲ ) = 1 τ log e τ φ ( θ ) π * ( θ | x ̲ ) d θ ,
The Bayes estimator of φ ( α ) under the LINEX loss function, denoted by φ ^ S L E C ( α ) , can be found by using Equations (26) and (22).
φ ^ S L E C ( α ) = 1 τ log [ k 1 0 e τ φ ( α ) α a 1 Γ ( α + β ) Γ ( α ) r e i = 1 r λ x i + b α × e ( n r ) log 1 I e λ x r ( α , β ) d α ] ,
where k 1 1 is defined in (23).
iii.
GE Loss Function
According to Calabria and Pulcini [13], the GE loss function of φ ( θ ) can be defined as:
φ ^ G E ( θ ) = E θ ( φ ( θ ) q | x ̲ ) 1 q = φ ( θ ) q π ( θ | x ̲ ) d θ 1 q ,
The Bayes estimator of φ ( α ) under the GE loss function, denoted by φ ^ S G E C ( α ) , can be found by using Equations (22) and (28).
φ ^ S G E C ( α ) = [ k 1 0 [ φ ( α ) ] q α a 1 Γ ( α + β ) Γ ( α ) r e i = 1 r λ x i + b α × e ( n r ) log 1 I e λ x r ( α , β ) d α ] 1 q ,
where k 1 1 is defined in (23).
The Bayes estimator of α using the importance sampling technique can be derived by rewriting the posterior density function in Equation (22); thus:
π * ( α | x ̲ ) ( i = 1 r λ x i + b ) a Γ ( a ) α a 1 e i = 1 r λ x i + b α Γ ( α + β ) Γ ( α ) r e ( n r ) log 1 I e λ x r ( α , β ) .
The posterior density function of α can be considered as:
π * ( α | x ̲ ) G a m m a ( a , i = 1 r λ x i + b ) g 1 ( α | x ̲ ) ,
where:
g 1 ( α | x ̲ ) = Γ ( α + β ) Γ ( α ) r e ( n r ) log 1 I e λ x r ( α , β ) .
The Bayes estimators of φ ( α ) under the SE, LINEX, and GE loss functions based on the importance sampling technique, denoted by φ ^ I S E C ( α ) , φ ^ I L E C ( α ) , and φ ^ I G E C ( α ) , respectively, could be found using the following Algorithm 1.
Algorithm 1 Importance sampling technique when α is unknown based on type-II censored samples.
  • Generate α i Gamma (a, i = 1 r λ x i + b ) .
  • Repeat Step 1 to obtain α 1 , α 2 , , α N .
  • Calculate the values.
    φ ^ I S E C ( α ) = j = 1 N φ ( α j ) g 1 ( α j | x ̲ ) j = 1 N g 1 ( α j | x ̲ )
    φ ^ I L E C ( α ) = 1 τ log j = 1 N e τ φ ( α j ) g 1 ( α j | x ̲ ) j = 1 N g 1 ( α j | x ̲ )
    φ ^ I G E C ( α ) = j = 1 N [ φ ( α j ) ] q g 1 ( α j | x ̲ ) j = 1 N g 1 ( α j | x ̲ ) 1 q
where
g 1 ( α j | x ̲ ) = Γ ( α j + β ) Γ ( α j ) r e ( n r ) log 1 I e λ x r ( α j , β )
The Bayes estimators of α are found numerically under these three loss functions by the NIntegrate function via Mathematica 11.

2.2.2. Case 2: Bayes Estimators When β Is Unknown

Suppose β is unknown and has the following prior distribution β G a m m a ( c , d ) , given by:
π ( β ) = d c Γ ( c ) β c 1 e d β , β > 0 , c , d > 0 .
By combining (9) and (35), the posterior distribution of the unknown parameter β is given by:
π * ( β | x ̲ ) = k 2 β c 1 Γ ( α + β ) Γ ( β ) r e d i = 1 r log [ 1 e λ x i ] β e ( n r ) log 1 I e λ x r ( α , β ) ,
where k 2 is the normalizing constant and can be written as:
k 2 1 = 0 β c 1 Γ ( α + β ) Γ ( β ) r e d i = 1 r log [ 1 e λ x i ] β e ( n r ) log 1 I e λ x r ( α , β ) d β .
Therefore, the Bayes estimator of β , denoted by φ ( β ) , is obtained under three types of loss function and two techniques as follows.
  • SE Loss Function
The Bayes estimator of φ ( β ) under the SE loss function, denoted by φ ^ S S E C ( β ) , can be found by using Equations (24) and (36).
φ ^ S S E C ( β ) = k 2 0 φ ( β ) β c 1 Γ ( α + β ) Γ ( β ) r e d i = 1 r log [ 1 e λ x i ] β × e ( n r ) log 1 I e λ x r ( α , β ) d β ,
where k 2 1 is defined in Equation (37).
ii.
LINEX Loss Function
The Bayes estimator of φ ( β ) under the LINEX loss function, denoted by φ ^ S L E C ( β ) , can be found by using Equations (26) and (36).
φ ^ S L E ( β ) = 1 τ log [ k 2 0 e τ φ ( β ) β c 1 Γ ( α + β ) Γ ( β ) r × e d i = 1 r log [ 1 e λ x i ] β e ( n r ) log 1 I e λ x r ( α , β ) d β ] ,
where k 2 1 is defined in (37).
iii.
GE Loss Function
The Bayes estimator of φ ( β ) under the GE loss function, denoted by φ ^ S G E C ( β ) , can be found by using Equations (28) and (36).
φ ^ S G E C ( β ) = [ k 2 0 [ φ ( β ) ] q β c 1 Γ ( α + β ) Γ ( β ) r × e d i = 1 r log [ 1 e λ x i ] β e ( n r ) log 1 I e λ x r ( α , β ) d β ] 1 q ,
where k 2 1 is defined in (37).
The Bayes estimator of β using the importance sampling technique can be derived by rewriting the posterior density function in Equation (36); thus:
π * ( β | x ̲ ) ( d i = 1 r log [ 1 e λ x i ] ) c Γ ( c ) β c 1 Γ ( α + β ) Γ ( β ) r × e d i = 1 r log [ 1 e λ x i ] β e ( n r ) log 1 I e λ x r ( α , β ) .
The posterior density function of β can be considered as:
π * ( β | x ̲ ) G a m m a ( c , d i = 1 r log [ 1 e λ x i ] ) g 2 ( β | x ̲ ) ,
where:
g 2 ( β | x ̲ ) = Γ ( α + β ) Γ ( β ) r e ( n r ) log 1 I e λ x r ( α , β ) .
The Bayes estimators of φ ( β ) under the SE, LINEX, and GE loss functions based on the importance sampling technique, denoted by φ ^ I S E C ( β ) , φ ^ I L E C ( β ) , and φ ^ I G E C ( β ) , respectively, could be found using the following Algorithm 2.
Algorithm 2 Importance sampling technique when β is unknown based on type-II censored samples.
  • Generate β i Gamma(c, d- i = 1 r log [ 1 e λ x i ] ) .
  • Repeat Step 1 to obtain β 1 , β 2 , , β N .
  • Calculate the values.
    φ ^ I S E C ( β ) = j = 1 N φ ( β j ) g 2 ( β j | x ̲ ) j = 1 N g 2 ( β j | x ̲ )
    φ ^ I L E C ( β ) = 1 τ log j = 1 N e τ φ ( β j ) g 2 ( β j | x ̲ ) j = 1 N g 2 ( β j | x ̲ )
    φ ^ I G E C ( β ) = j = 1 N [ φ ( β j ) ] q g 2 ( β j | x ̲ ) j = 1 N g 2 ( β j | x ̲ ) 1 q
where
g 2 ( β j | x ̲ ) = Γ ( α + β j ) Γ ( β j ) r e ( n r ) log 1 I e λ x r ( α , β j )
The Bayes estimators of β under the three loss functions cannot be computed analytically through the two techniques used. They can be found numerically using the NIntegrate function via Mathematica 11.

2.2.3. Case 3: Bayes Estimators When λ Is Unknown

Suppose that the scale parameter λ is unknown and has the following prior distribution λ G a m m a ( f , ν ) , given by:
π ( λ ) = ν f Γ ( f ) λ f 1 e ν λ , λ > 0 , f , ν > 0 .
By combining (9) and (46), the posterior distribution of the unknown parameter λ is given by:
π * ( λ | x ̲ ) = k 3 λ f + r 1 e i = 1 r α x i + ν λ e ( β 1 ) i = 1 r log [ 1 e λ x i ] e ( n r ) log 1 I e λ x r ( α , β ) ,
where k 3 is the normalizing constant and can be written as:
k 3 1 = 0 λ f + r 1 e i = 1 r α x i + ν λ e ( β 1 ) i = 1 r log [ 1 e λ x i ] e ( n r ) log 1 I e λ x r ( α , β ) d λ .
The Bayes estimator of λ is indicated by φ ( λ ) , and this estimator is obtained under three types of loss function and two techniques as follows.
  • SE Loss Function
The Bayes estimator of φ ( λ ) under the SE loss function, denoted by φ ^ S S E C ( λ ) , can be found by using Equations (24) and (47).
φ ^ S S E C ( λ ) = k 3 0 φ ( λ ) λ f + r 1 e i = 1 r α x i + ν λ e ( β 1 ) i = 1 r log [ 1 e λ x i ] × e ( n r ) log 1 I e λ x r ( α , β ) d λ ,
where k 3 1 is defined in Equation (48).
ii.
LINEX Loss Function
The Bayes estimator of φ ( λ ) under the LINEX loss function, denoted by φ ^ S L E C ( λ ) , can be found by using Equations (26) and (47).
φ ^ S L E C ( λ ) = 1 τ log [ k 3 0 e τ φ ( λ ) λ f + r 1 e i = 1 r α x i + ν λ × e ( β 1 ) i = 1 r log [ 1 e λ x i ] e ( n r ) log 1 I e λ x r ( α , β ) d λ ] ,
where k 3 1 is defined in (48).
iii.
GE Loss Function
The Bayes estimator of φ ( λ ) under the GE loss function, denoted by φ ^ S G E C ( λ ) , can be found by using Equations (28) and (47).
φ ^ S G E C ( λ ) = [ k 3 0 [ φ ( λ ) ] q λ f + r 1 e i = 1 r α x i + ν λ × e ( β 1 ) i = 1 r log [ 1 e λ x i ] e ( n r ) log 1 I e λ x r ( α , β ) d λ ] 1 q .
where k 3 1 is defined in (48).
The Bayes estimator of λ using the importance sampling technique can be obtained by rewriting the posterior density function in Equation (47); thus:
π * ( λ | x ̲ ) i = 1 r α x i + ν f + r Γ ( f + r ) λ f + r 1 e i = 1 r α x i + ν λ × e ( β 1 ) i = 1 r log [ 1 e λ x i ] e ( n r ) log 1 I e λ x r ( α , β ) .
The posterior density function of λ can be considered as:
π * ( λ | x ̲ ) G a m m a ( f + r , i = 1 r α x i + ν ) g 3 ( λ | x ̲ ) ,
where:
g 3 ( λ | x ̲ ) = e ( β 1 ) i = 1 r log [ 1 e λ x i ] e ( n r ) log 1 I e λ x r ( α , β ) .
The Bayes estimators of φ ( λ ) under the SE, LINEX, and GE loss functions based on the importance sampling technique, denoted by φ ^ I S E C ( λ ) , φ ^ I L E C ( λ ) , and φ ^ I G E C ( λ ) , respectively, could be found using the following Algorithm 3.
Algorithm 3 Importance sampling technique when λ is unknown based on type-II censored samples.
  • Generate λ i Gamma(f + r, i = 1 r α x i + ν ) .
  • Repeat Step 1 to obtain λ 1 , λ 2 , , λ N .
  • Calculate the values.
    φ ^ I S E C ( λ ) = j = 1 N φ ( λ j ) g 3 ( λ j | x ̲ ) j = 1 N g 3 ( λ j | x ̲ )
    φ ^ I L E C ( λ ) = 1 τ log j = 1 N e τ φ ( λ j ) g 3 ( λ j | x ̲ ) j = 1 N g 3 ( λ j | x ̲ )
    φ ^ I G E C ( λ ) = j = 1 N [ φ ( λ j ) ] q g 3 ( λ j | x ̲ ) j = 1 N g 3 ( λ j | x ̲ ) 1 q
where
g 3 ( λ j | x ̲ ) = e ( β 1 ) i = 1 r log [ 1 e λ j x i ] e ( n r ) log 1 I e λ j x r ( α , β ) .
The Bayes estimators of λ under the three loss functions cannot be computed analytically through the two techniques used. It can be found numerically using the NIntegrate function via Mathematica 11.

2.2.4. Case 4: Bayes Estimators When λ and β Are Unknown

Consider that both parameters λ and β are unknown. Suppose β G a m m a ( c , d ) and λ G a m m a ( f , ν ) . Therefore, the joint prior distribution is given by:
π ( λ , β ) = ν f Γ ( f ) λ f 1 e ν λ d c Γ ( c ) β c 1 e d β , λ , β > 0 , f , ν , c , d > 0 .
By combining (9) and (57), the joint posterior distribution of the unknown parameters λ and β is given by:
π * ( λ , β | x ̲ ) = k 4 λ f + r 1 e i = 1 r α x i + ν λ β c 1 Γ ( α + β ) Γ ( β ) r e d i = 1 r log [ 1 e λ x i ] β × e i = 1 r log [ 1 e λ x i ] e ( n r ) log 1 I e λ x r ( α , β ) ,
where k 4 is the normalizing constant and can be written as:
k 4 1 = 0 0 λ f + r 1 e i = 1 r α x i + ν λ β c 1 Γ ( α + β ) Γ ( β ) r e d i = 1 r log [ 1 e λ x i ] β × e i = 1 r log [ 1 e λ x i ] e ( n r ) log 1 I e λ x r ( α , β ) d λ d β .
The Bayes estimators of λ and β are derived under three types of loss function as follows.
  • SE Loss Function
The Bayes estimators of φ ( λ , β ) under the SE loss function, denoted by φ ^ S S E C ( λ , β ) , can be found by using Equations (24) and (58).
φ ^ S S E C ( λ , β ) = k 4 0 0 φ ( λ , β ) λ f + r 1 e i = 1 r α x i + ν λ β c 1 Γ ( α + β ) Γ ( β ) r × e d i = 1 r log [ 1 e λ x i ] β e i = 1 r log [ 1 e λ x i ] × e ( n r ) log 1 I e λ x r ( α , β ) d λ d β ,
where k 4 1 is defined in Equation (59).
Different forms of the Bayes estimator for φ ( λ , β ) are obtained from Equation (60):
  • When φ ( λ , β ) = λ , we obtain the Bayes estimator for λ , denoted by λ ^ S S E C , as:
    λ ^ S S E C = k 4 0 0 λ f + r e i = 1 r α x i + ν λ β c 1 Γ ( α + β ) Γ ( β ) r × e d i = 1 r log [ 1 e λ x i ] β e i = 1 r log [ 1 e λ x i ] × e ( n r ) log 1 I e λ x r ( α , β ) d λ d β ,
    where k 4 1 is defined in Equation (59);
  • When φ ( λ , β ) = β , we obtain the Bayes estimator for β , denoted by β ^ S S E C , as:
    β ^ S S E C = k 4 0 0 λ f + r 1 e i = 1 r α x i + ν λ β c Γ ( α + β ) Γ ( β ) r e d i = 1 r log [ 1 e λ x i ] β e i = 1 r log [ 1 e λ x i ] × e ( n r ) log 1 I e λ x r ( α , β ) d λ d β ,
    where k 4 1 is defined in Equation (59).
ii.
LINEX Loss Function
The Bayes estimator of φ ( λ , β ) under the LINEX loss function, denoted by φ ^ S L E C ( λ , β ) , can be found by using Equations (26) and (58).
φ ^ S L E C ( λ , β ) = 1 τ log [ k 4 0 0 e τ φ ( λ , β ) λ f + r 1 e i = 1 r α x i + ν λ β c 1 × Γ ( α + β ) Γ ( β ) r e d i = 1 r log [ 1 e λ x i ] β e i = 1 r log [ 1 e λ x i ] × e ( n r ) log 1 I e λ x r ( α , β ) d λ d β ] ,
where k 4 1 is defined in (59).
Different forms of the Bayes estimator for φ ( λ , β ) are obtained from Equation (63):
  • When φ ( λ , β ) = λ , we obtain the Bayes estimator for λ , denoted by λ ^ S L E C , as:
    λ ^ S L E C = 1 τ log [ k 4 0 0 λ f + r 1 e i = 1 r α x i + ν + τ λ β c 1 × Γ ( α + β ) Γ ( β ) r e d i = 1 r log [ 1 e λ x i ] β e i = 1 r log [ 1 e λ x i ] × e ( n r ) log 1 I e λ x r ( α , β ) d λ d β ] ,
    where k 4 1 is defined in Equation (59);
  • When φ ( λ , β ) = β , we obtain the Bayes estimator for β , denoted by β ^ S L E C , as:
    β ^ S L E C = 1 τ log [ k 4 0 0 λ f + r 1 e i = 1 r α x i + ν λ β c 1 Γ ( α + β ) Γ ( β ) r × e d i = 1 r log [ 1 e λ x i ] + τ β e i = 1 r log [ 1 e λ x i ] × e ( n r ) log 1 I e λ x r ( α , β ) d λ d β ] ,
    where k 4 1 is defined in Equation (59).
iii.
GE Loss Function
The Bayes estimator of φ ( λ , β ) under the GE loss function, denoted by φ ^ S G E C ( λ , β ) , can be found by using Equations (28) and (58).
φ ^ S G E C ( λ , β ) = [ k 4 0 0 [ φ ( λ , β ) ] q λ f + r 1 e i = 1 r α x i + ν λ β c 1 × Γ ( α + β ) Γ ( β ) r e d i = 1 r log [ 1 e λ x i ] β e i = 1 r log [ 1 e λ x i ] × e ( n r ) log 1 I e λ x r ( α , β ) d λ d β ] 1 q ,
where k 4 1 is defined in (59).
Different forms of the Bayes estimator for φ ( λ , β ) are obtained from Equation (66):
  • When φ ( λ , β ) = λ , we obtain the Bayes estimator for λ , denoted by λ ^ S G E C , as:
    λ ^ S G E C = [ k 4 0 0 λ f + r q 1 e i = 1 r α x i + ν λ β c 1 Γ ( α + β ) Γ ( β ) r × e d i = 1 r log [ 1 e λ x i ] β e i = 1 r log [ 1 e λ x i ] × e ( n r ) log 1 I e λ x r ( α , β ) d λ d β ] 1 q ,
    where k 4 1 is defined in Equation (59);
  • When φ ( λ , β ) = β , we obtain the Bayes estimator for β , denoted by β ^ S G E C , as:
    β ^ S G E C = [ k 4 0 0 λ f + r 1 e i = 1 r α x i + ν λ β c q 1 × Γ ( α + β ) Γ ( β ) r e d i = 1 r log [ 1 e λ x i ] β e i = 1 r log [ 1 e λ x i ] × e ( n r ) log 1 I e λ x r ( α , β ) d λ d β ] 1 q ,
    where k 4 1 is defined in Equation (59).
The Bayes estimators of β and λ can be obtained using the importance sampling technique by rewriting the joint posterior density function of λ and β in Equation (58) as follows:
π * ( λ , β | x ̲ ) i = 1 r α x i + ν f + r Γ ( f + r ) λ f + r 1 e i = 1 r α x i + ν λ × ( d i = 1 r log [ 1 e λ x i ] ) c Γ ( c ) β c 1 e d i = 1 r log [ 1 e λ x i ] β × Γ ( α + β ) Γ ( β ) r e i = 1 r log [ 1 e λ x i ] e ( n r ) log 1 I e λ x r ( α , β ) ( d i = 1 r log [ 1 e λ x i ] ) c .
The joint posterior density function of β and λ can be considered as:
π * ( λ , β | x ̲ ) G a m m a ( f + r , i = 1 r α x i + ν ) × G a m m a ( c , d i = 1 r log [ 1 e λ x i ] ) × g 4 ( λ , β | x ̲ ) ,
where:
g 4 ( λ , β | x ̲ ) = Γ ( α + β ) Γ ( β ) r e i = 1 r log [ 1 e λ x i ] e ( n r ) log 1 I e λ x r ( α , β ) ( d i = 1 r log [ 1 e λ x i ] ) c
The Bayes estimators of φ ( λ , β ) under the SE, LINEX, and GE loss functions based on the importance sampling technique, denoted by φ ^ I S E C ( λ , β ) , φ ^ I L E C ( λ , β ) , and φ ^ I G E C ( λ , β ) , respectively, could be found using the following Algorithm 4.
Algorithm 4 Importance sampling technique when λ and β are unknown based on type-II censored samples.
  • Generate λ i Gamma(f + r, i = 1 r α x i + ν ) and β i Gamma(c, d- i = 1 r log [ 1 e λ x i ] ) .
  • Repeat Step 1 to obtain ( λ 1 , β 1 ) , ( λ 2 , β 2 ) , , ( λ N , β N ) .
  • Calculate the values.
    φ ^ I S E C ( λ , β ) = j = 1 N φ ( λ j , β j ) g 4 ( λ j , β j | x ̲ ) j = 1 N g 4 ( λ j , β j | x ̲ )
    φ ^ I L E C ( λ , β ) = 1 τ log j = 1 N e τ φ ( λ j , β j ) g 4 ( λ j , β j | x ̲ ) j = 1 N g 4 ( λ j , β j | x ̲ )
    φ ^ I G E C ( λ , β ) = j = 1 N [ φ ( λ j , β j ) ] q g 4 ( λ j , β j | x ̲ ) j = 1 N g 4 ( λ j , β j | x ̲ ) 1 q
where
g 4 ( λ j , β j | x ̲ ) = Γ ( α + β j ) Γ ( β j ) r e i = 1 r log [ 1 e λ j x i ] e ( n r ) log 1 I e λ j x r ( α , β j ) ( d i = 1 r log [ 1 e λ j x i ] ) c .
The Bayes estimators of λ and β under the three loss functions can be found numerically using the NIntegrate function via Mathematica 11.

3. Simulation Study

Simulation studies were conducted using Mathematica 11 to clarify the performance of the proposed estimators. Simulation results are given for the ML and Bayesian methods based on type-II censored samples. Furthermore, biases and mean-squared errors (MSE) were considered to illustrate the performance of the different estimators, defined as:
B i a s ( θ ^ ) = E [ θ ^ ] θ M S E ( θ ^ ) = E ( θ ^ θ ) 2
The ML estimates of the parameters α , β λ S ( x 0 ) and h ( x 0 ) could be found using the following Algorithm 5.
Algorithm 5 ML method of the parameters α , β λ S ( x 0 ) and h ( x 0 ) based on type-II censored samples.
  • For given true values selected as ( α , β , λ ) , generate a random sample of size n from Equation (5).
  • Arrange Step 1 in ascending order to obtain X 1 : n < X 2 : n < X 3 : n < < X n : n .
  • Obtain the censored sample according to the censoring percentage.
  • Estimate the ML of the parameters α , β , and λ using the Newton–Raphson method to solve the equations given in (13), (16) and (18), simultaneously.
  • Compute the estimators of S ( x 0 ) and h ( x 0 ) from (6) and (7), respectively, using the estimates in Step 4.
  • Repeat Steps 1–5 1000 times.
  • Calculate the mean, bias, and MSE for each estimate.
Three different sets of true parameter values were selected to perform the simulation. Furthermore, different sample sizes n = 30 , 50 , and 100 and two censoring percentage of 80 % and 90 % were selected. Table 1 shows the ML estimates for α , β , λ , S ( x 0 ) , h ( x 0 ) , the bias, and the MSE for the true values ( α = 0.8 , β = 4 , λ = 3 ) and S ( x 0 ) = 0.72904 , h ( x 0 ) = 0.81786 at x 0 = 1 . The ML estimates for α , β , λ , S ( x 0 ) , h ( x 0 ) , the bias, and the MSE for the true values ( α = 3 , β = 0.8 , λ = 3 ) and S ( x 0 ) = 0.87604 , h ( x 0 ) = 0.05352 at x 0 = 5 are summarized in Table 2. Moreover, for the true values ( α = 3 , β = 8 , λ = 2 ) and S ( x 0 ) = 0.22471 , h ( x 0 ) = 1.60828 at x 0 = 2 , the ML estimates for α , β , λ , S ( x 0 ) , h ( x 0 ) , the bias, and the MSE are listed in Table 3.
It is clear from Table 1, Table 2 and Table 3 that the MSEs of the estimates decrease as the sample size increases. Based on biases, the parameter λ was underestimated for all sets of true values, except in Table 2, the parameter λ was overestimated. Likewise, the parameter α was overestimated. Besides, we noticed that when the percentage of censoring was 90 % , the MSE of the estimates in most cases was better.
Bayes estimates of the parameters α , β λ S ( x 0 ) , and h ( x 0 ) could be found using the following Algorithm 6.
Algorithm 6 Bayesian method of the parameters α , β λ S ( x 0 ) , and h ( x 0 ) based on type-II censored samples.
  • For given true parameter values selected as ( α , β , λ ) , generate a random sample of size n from Equation (5).
  • Arrange Step 1 in ascending order to obtain X 1 : n < X 2 : n < X 3 : n < X n : n .
  • Obtain the censored sample according to the censoring percentage.
  • For given values of the hyperparameters parameters a , b , c , d , f , ν , compute the Bayes estimates of the parameters α , β , λ via the standard Bayes technique for the first three cases using the NIntegrate function under the SE, LINEX, and GE loss functions as shown in Table 4.
  • For given values of the hyperparameters parameters a , b , c , d , f , ν , compute the Bayes estimates of the parameters α , β , λ via the importance sampling technique for all cases under the SE, LINEX, and GE loss functions, as shown in Table 4.
  • Compute the Bayes estimates of S ( x 0 ) and h ( x 0 ) for the four cases from (6) and (7), respectively, using the estimates in the previous steps.
  • Repeat Steps 1–6 1000 times.
  • Calculate the mean, bias, and MSE for each estimate.
The simulation study was carried out with the true value of parameter ( α = 0.8 , β = 4 , λ = 3 ) and different sample sizes n = 30 , 50 , and 100 for all four cases. Furthermore, we considered three different values of the LINEX shape parameter ( τ = 0.001 , τ = 2 , τ = 5 ) and three values of the GE shape parameter ( q = 1 , q = 3 , q = 3 ). The Bayes estimates for first three cases were derived based on the standard Bayes and importance sampling techniques. For the last case, this was obtained via the importance sampling technique. The values of the hyperparameters for the standard Bayes technique are ( a = 2 , b = 4 , c = 5 , d = 2 , f = 2 , ν = 5 ) . Further, the values of the hyperparameters for the importance sampling technique are ( a = 80 , b = 0.1 , c = 14 , d = 2 , f = 30 , ν = 0.001 ) with a sample size of N = 1000 . For given S ( x 0 ) = 0.72904 and h ( x 0 ) = 0.81786 at x 0 = 1 , the ML and Bayes estimates of the S ( x 0 ) and h ( x 0 ) were computed. The four cases were calculated as follows.
Table 5, Table 6, Table 7, Table 8, Table 9, Table 10 and Table 11 show the Bayes and ML estimates of the parameters.
The results of the two estimation methods from Table 5, Table 6, Table 7, Table 8, Table 9, Table 10 and Table 11 are summarized as follows:
  • From Table 5, Table 6, Table 7, Table 8, Table 9, Table 10 and Table 11, the MSEs of the ML estimates and Bayes estimates of α , β , λ , S ( x 0 ) , and h ( x 0 ) decrease as the sample size increases;
  • From Table 5, Table 6, Table 7, Table 8, Table 9, Table 10 and Table 11, the mean, bias, and MSE values for the Bayes estimates under the SE loss function, the LINEX loss function with ( τ = 0.001 ), and the GE loss function with ( q = 1 ) are very similar;
  • For Case 1, the Bayes estimates via the standard Bayes technique of α , S ( x 0 ) , and h ( x 0 ) perform the estimates better than the ML estimates under the different loss functions, as shown in Table 5. According to Table 6, we note that the ML estimates give better values than the Bayes estimates via the importance sampling technique;
  • From Table 5, the Bayes estimates under the GE loss function ( q = 3 ) are considered the best estimates of S ( x 0 ) ;
  • From Table 7, the Bayes estimates of β , S ( x 0 ) , and h ( x 0 ) via the standard Bayes technique perform the best based on MSEs and biases at n = 30 . Furthermore, when n = 50 , 100 , the ML estimates of β , S ( x 0 ) , and h ( x 0 ) perform the estimates better than the Bayes estimates;
  • When β is unknown, the Bayes estimates of β via the importance sampling technique perform the best at n = 30 , under the LINEX loss function ( τ = 5 ) . For n = 50 , 100 , the ML estimates of β give the best estimates. Furthermore, based on the MSEs and biases, the ML estimates of S ( x 0 ) and h ( x 0 ) give the best estimates (see Table 8);
  • Based on the MSEs, the ML estimates of λ , S ( x 0 ) , and h ( x 0 ) perform the estimates better than the Bayes estimates based on the two techniques (see Table 9 and Table 10);
  • From Table 11, the ML estimates of λ , S ( x 0 ) , and h ( x 0 ) perform the best based on the smallest MSEs. Besides, the Bayes estimates of β via the importance sampling technique perform the best under the LINEX loss function ( τ = 5 ).

4. Application

The performance of the BIED based on type-II censored samples is illustrated through two real datasets. The BIED model was compared with other lifetime models, such as the inverse exponential distribution (IED) as a special case of the BIED, introduced by Lin et al. [14], the inverse Weibull distribution (IWD) introduced by Keller and Kamath [15], the Weibull inverted exponential distribution (WIED) defined by [16], the Weibull exponential distribution (WED) (see [17]), and the odd Fréchet inverse exponential distribution (OFIED) introduced by Alrajhi [18]. Moreover, the model selection criteria were considered, which included the Akaike information criterion (AIC), log-likelihood ( ) , Bayesian information criterion (BIC), consistent Akaike information criterion (CAIC), and Hannan–Quinn information criterion (HQIC). The smallest values of the AIC, BIC, CAIC, and HQIC, and the highest value determine the best-fit model for the data. For more details about these criteria and their uses, see [19,20].
A I C = 2 ( Θ ^ ) + 2 p
B I C = 2 ( Θ ^ ) + p log n
C A I C = A I C + 2 p ( p + 1 ) n p 1
H Q I C = 2 ( Θ ^ ) + 2 p log ( log n )
where ( Θ ^ ) denotes the log-likelihood function, p is the number of parameters, and n is the sample size:
  • Aluminum coupons’ cut:
The following data consisting of 102 observations were used by Birnbaum and Saunders [21] and correspond to the fatigue life of 6061-T6 aluminum coupons in cycles ( × 10 3 ) with the maximum pressure of 26,000 psi. These coupons were cut parallel to the direction of rolling and oscillated at 18 cycles per second.
  • 233, 258, 268, 276, 290, 310, 312, 315, 318, 321, 321, 329, 335, 336,
    338, 338, 342, 342, 342, 344, 349, 350, 350, 351, 351, 352, 352, 356,
    358, 358, 360, 362, 363, 366, 367, 370, 370, 372, 372, 374, 375, 376,
    379, 379, 380, 382, 389, 389, 395, 396, 400, 400, 400, 403, 404, 406,
    408, 408, 410, 412, 414, 416, 416, 416, 420, 422, 423, 426, 428, 432,
    432, 433, 433, 437, 438, 439, 439, 443, 445, 445, 452, 456, 456, 460,
    464, 466, 468, 470, 470, 473, 474, 476, 476, 486, 488, 489, 490, 491,
    503, 517, 540, 560.
Descriptive statistics of these data are presented in Table 12.
Based on the descriptive statistics in Table 12, we observed that the skewness value was close to zero; thus, the distribution of the fatigue life of the 6061-T6 aluminum coupons’ cut dataset was approximately normal, while the variance was 3884.30, which indicate high variability in the dataset.
According to Figure 1, we can note that there were no outliers in the fatigue life of the 6061-T6 aluminum coupons’ cut data. The ML estimates and Bayesian estimates via the standard Bayes technique of the BIED parameters for the fatigue life of the 6061-T6 aluminum coupons’ cut data are presented in Table 13 at two censoring percentages of 80 % and 90 % .
The ML estimates of the model parameters and the performance of the BIED against other models for the fatigue life of the 6061-T6 aluminum coupons’ cut data are presented in Table 14 at two censoring percentages of 80 % and 90 % .
In Table 14, the values of the AIC, BIC, CAIC, and HQIC show that the BIED was the best model for analyzing the fatigue life of the 6061-T6 aluminum coupons’ cut data. Furthermore, we can consider that the WED is a good alternative model for these data. The estimated PDF and estimated CDF of the models for the fatigue life of the 6061-T6 aluminum coupons’ cut data at two censoring percentages of 80 % and 90 % are shown in Figure 2 and Figure 3.
  • Patients suffering from acute myelogenous leukemia:
The following data consisting of 33 observations were studied by Feigl and Zelen [22] and represent the survival times (in weeks) of patients suffering from acute myelogenous leukemia.
  • 1, 1, 2, 3, 3, 3, 4, 4, 4, 4, 5, 7, 8, 16, 16, 17, 22, 22, 26, 30,
    39, 43, 56, 56, 65, 65, 65, 100, 108, 121, 134, 143, 156
Descriptive statistics of these data are presented in Table 15.
According to the descriptive statistics in Table 15, we observed that the distribution of the survival times (in weeks) of patients suffering from acute myelogenous leukemia data was positively skewed, while the variance was 2181.17, which indicates high variability in the dataset.
According to Figure 4, we can note that there were no outliers in the survival times (in weeks) of patients suffering from acute myelogenous leukemia data. Additionally, the ML estimates and Bayesian estimates via the standard Bayes technique of the BIED for this dataset are presented in Table 16 at two censoring percentages of 80 % and 90 % .
The ML estimates of the model parameters and the performance of the BIED against other models for the survival times (in weeks) of patients suffering from acute myelogenous leukemia data are shown in Table 14 at two censoring percentages of 80 % and 90 % .
Table 17, based on the values of the AIC, BIC, CAIC, and HQIC, shows that the BIED was the best model for fitting the survival times (in weeks) of patients suffering from acute myelogenous leukemia data. Further, the estimated PDF and estimated CDF of the models for this dataset at two censoring percentages of 80 % and 90 % are shown in Figure 5 and Figure 6.

5. Conclusions

In this article, the maximum likelihood and Bayes estimators of the BIED were derived based on type-II censored samples. The invariance property was used to estimate the survival and hazard functions. Furthermore, in the Bayesian estimation, three loss functions were used with two techniques. The gamma distribution was assumed as a prior distribution for the shape and scale parameters. Besides, it can be concluded that the MSEs of the ML estimates and Bayes estimates for the unknown parameters decreased as the sample size increased. Furthermore, when α was unknown, the Bayes estimates gave better estimates via the standard Bayes technique. The ML estimates gave better estimates than the Bayes estimates using the two techniques when λ was unknown. For the third case, when β was unknown, the ML estimates gave better results than the Bayes estimates as the sample size increased. Likewise, when β and λ were unknown, the Bayes estimates via the importance sampling technique under the LINEX loss function ( τ = 5 ) gave better results than the ML estimates as the sample size increased. Two real datasets were applied, and the BIED provided a better fit than the other compared distributions.

Author Contributions

Software, L.S.A.; supervision and writing—review, R.A.B. and M.A.A.; methodology and writing—original draft, L.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Singh, U.; Gupta, P.K.; Upadhyay, S.K. Estimation of parameters for exponentiated Wibull family under type-II censoring scheme. Comput. Stat. Data Anal. 2005, 48, 509–523. [Google Scholar] [CrossRef]
  2. Prakash, G. Some estimation procedures for the inverted exponential distribution. S. Pac. J. Nat. Appl. Sci. 2009, 27, 71–78. [Google Scholar] [CrossRef] [Green Version]
  3. Dey, S.; Dey, T. On progressively censored generalized inverted exponential distribution. J. Appl. Stat. 2014, 41, 2557–2576. [Google Scholar] [CrossRef]
  4. Singh, B.; Goel, R. The beta inverted exponential distribution: Properties and applications. Int. J. Appl. Sci. Math. 2015, 2, 132–141. [Google Scholar]
  5. Garg, R.; Dube, M.; Kumar, K.; Krishna, H. On randomly censored generalized inverted exponential distribution. Am. J. Math. Manag. Sci. 2016, 35, 361–379. [Google Scholar] [CrossRef]
  6. Bakoban, R.; Abu-Zinadah, H.H. The beta generalized inverted exponential distribution with real data applications. REVSTAT-Stat. J. 2017, 15, 65–88. [Google Scholar]
  7. Aldahlan, M.A. The inverse Weibull inverse exponential distribution with application. Int. J. Contemp. Math. Sci. 2019, 14, 17–30. [Google Scholar] [CrossRef] [Green Version]
  8. Arnold, B.C.; Balakrishnan, N.; Nagaraja, H.N. A First Course in Order Statistics; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2008. [Google Scholar]
  9. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 2nd ed.; Dover Publication: New York, NY, USA, 1972. [Google Scholar]
  10. Olver, F.W.; Lozier, D.W.; Boisvert, R.F.; Clark, C.W. NIST Handbook of Mathematical Functions, 1st ed.; Cambridge University Press: New York, NY, USA, 2010. [Google Scholar]
  11. Box, G.E.; Tiao, G.C. Bayesian Inference in Statistical Analysis, 1st ed.; John Wiley & Sons: Hoboken, NJ, USA, 1992. [Google Scholar]
  12. Varian, H.R. A Bayesian approach to real estate assessment. Stud. Bayesian Econom. Stat. Honor. Leonard J. Savage 1975, 5, 195–208. [Google Scholar]
  13. Calabria, R.; Pulcini, G. An engineering approach to Bayes estimation for the Weibull distribution. Microelectron. Reliab. 1994, 34, 789–802. [Google Scholar] [CrossRef]
  14. Lin, C.; Duran, B.; Lewis, T. Inverted gamma as a life distribution. Microelectron. Reliab. 1989, 29, 619–626. [Google Scholar] [CrossRef]
  15. Keller, A.Z.; Kamath, A.R.R. Alternate reliability models for mechanical systems. In Proceedings of the 3rd International Conference on Reliability and Maintainability, Toulouse, France, 18–21 October 1982; pp. 411–415. [Google Scholar]
  16. Pelumi, E.O.; Adebowale, O.A.; Enahoro, A.O. The Weibull inverted exponential distribution: A generalization of the inverse exponential distribution. In Proceedings of the World Congress on Engineering, London, UK, 5–7 July 2017; Volume 1, pp. 16–19. [Google Scholar]
  17. Oguntunde, P.; Balogun, O.; Okagbue, H.; Bishop, S. The Weibull exponential distribution: Its properties and applications. J. Appl. Sci. 2015, 15, 1305–1311. [Google Scholar] [CrossRef]
  18. Alrajhi, S. The odd Fréchet inverse exponential distribution with application. J. Nonlinear Sci. Appl. 2019, 12, 535–542. [Google Scholar] [CrossRef] [Green Version]
  19. Christie, M.; Cliffe, A.; Dawid, P.; Senn, S.S. Simplicity, Complexity and Modelling; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  20. Hannan, E.J.; Quinn, B.G. The determination of the order of an autoregression. J. R. Stat. Soc. 1979, 41, 190–195. [Google Scholar] [CrossRef]
  21. Birnbaum, Z.W.; Saunders, S.C. Estimation for a family of life distributions with applications to fatigue. J. Appl. Probab. 1969, 6, 328–347. [Google Scholar] [CrossRef]
  22. Feigl, P.; Zelen, M. Estimation of exponential survival probabilities with concomitant information. Biometrics 1965, 21, 826–838. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Boxplot for the fatigue life of the 6061-T6 aluminum coupons’ cut data.
Figure 1. Boxplot for the fatigue life of the 6061-T6 aluminum coupons’ cut data.
Mathematics 10 00506 g001
Figure 2. Plots of the estimated PDF and estimated CDF of the models for the fatigue life of the 6061-T6 aluminum coupons’ cut data for r = 82 .
Figure 2. Plots of the estimated PDF and estimated CDF of the models for the fatigue life of the 6061-T6 aluminum coupons’ cut data for r = 82 .
Mathematics 10 00506 g002
Figure 3. Plots of the estimated PDF and estimated CDF of the models for the fatigue life of the 6061-T6 aluminum coupons’ cut data for r = 92 .
Figure 3. Plots of the estimated PDF and estimated CDF of the models for the fatigue life of the 6061-T6 aluminum coupons’ cut data for r = 92 .
Mathematics 10 00506 g003
Figure 4. Boxplot for the survival times (in weeks) of patients suffering from acute myelogenous leukemia data.
Figure 4. Boxplot for the survival times (in weeks) of patients suffering from acute myelogenous leukemia data.
Mathematics 10 00506 g004
Figure 5. Plots of the estimated PDF and estimated CDF of the models for the patients suffering from acute myelogenous leukemia data for r = 26 .
Figure 5. Plots of the estimated PDF and estimated CDF of the models for the patients suffering from acute myelogenous leukemia data for r = 26 .
Mathematics 10 00506 g005
Figure 6. Plots of the estimated PDF and estimated CDF of the models for the patients suffering from acute myelogenous leukemia data for r = 30 .
Figure 6. Plots of the estimated PDF and estimated CDF of the models for the patients suffering from acute myelogenous leukemia data for r = 30 .
Mathematics 10 00506 g006
Table 1. The ML estimates, bias, and MSE of the parameters α , β λ S ( x 0 ) , and h ( x 0 ) for true parameter values ( α = 0.8 , β = 4 , λ = 3 ).
Table 1. The ML estimates, bias, and MSE of the parameters α , β λ S ( x 0 ) , and h ( x 0 ) for true parameter values ( α = 0.8 , β = 4 , λ = 3 ).
nr α ^ β ^ λ ^ S ^ ( x 0 ) h ^ ( x 0 )
3024MLE1.277914.056042.451340.735190.83076
Bias0.477910.05604−0.548660.006150.01291
MSE0.650261.076370.822830.004020.02662
27MLE0.899244.011432.859280.725550.82323
Bias0.099240.01143−0.14072−0.003490.00537
MSE0.076090.926760.320390.003960.02497
5040MLE0.881023.970672.904040.731320.81202
Bias0.08102−0.02934−0.095960.00227−0.00584
MSE0.056730.286390.207240.002500.01299
45MLE0.874033.959372.902660.732910.80861
Bias0.07403−0.04063−0.097340.00386−0.00925
MSE0.039690.275140.199960.002580.01454
10080MLE0.863523.991162.925670.733800.80855
Bias0.06352−0.00884−0.074330.00475−0.00931
MSE0.035420.220320.138270.001310.00643
90MLE0.841683.983372.963890.731890.81021
Bias0.04168−0.01664−0.036110.00285−0.00765
MSE0.027520.200840.156400.001340.00786
Table 2. The ML estimates, bias, and MSE of the parameters α , β λ S ( x 0 ) and h ( x 0 ) for true parameter values ( α = 3 , β = 0.8 , λ = 3 ) .
Table 2. The ML estimates, bias, and MSE of the parameters α , β λ S ( x 0 ) and h ( x 0 ) for true parameter values ( α = 3 , β = 0.8 , λ = 3 ) .
nr α ^ β ^ λ ^ S ^ ( x 0 ) h ^ ( x 0 )
3024MLE3.286950.882413.087630.877650.05352
Bias0.286950.082410.087630.00161−0.00000
MSE0.962470.048960.188240.002140.00020
27MLE3.202030.879133.070170.872300.05497
Bias0.202030.079130.07017−0.003740.00144
MSE0.869520.046670.157280.002030.00017
5040MLE3.176750.861903.026400.873150.05506
Bias0.176750.061910.02631−0.002880.00154
MSE0.521610.029950.108010.001320.00012
45MLE3.093160.843953.041720.873720.05457
Bias0.093160.043950.04172−0.002320.00104
MSE0.332680.023470.088880.001020.00009
10080MLE3.095050.834873.021980.874450.05438
Bias0.095040.034870.02198−0.001590.00086
MSE0.256490.013830.074760.000770.00006
90MLE3.098430.828403.021110.878050.05333
Bias0.098430.028390.021110.00201−0.00020
MSE0.149810.012120.055920.000460.00003
Table 3. The ML estimates, bias, and MSE of the parameters α , β λ S ( x 0 ) and h ( x 0 ) for true parameter values ( α = 3 , β = 8 , λ = 2 ) .
Table 3. The ML estimates, bias, and MSE of the parameters α , β λ S ( x 0 ) and h ( x 0 ) for true parameter values ( α = 3 , β = 8 , λ = 2 ) .
nr α ^ β ^ λ ^ S ^ ( x 0 ) h ^ ( x 0 )
3024MLE3.762318.041031.830680.228961.62124
Bias0.762310.04103−0.169320.004250.01296
MSE2.047440.266250.173990.002790.02465
27MLE3.719197.997491.842190.230801.61052
Bias0.71919−0.00251−0.157820.006090.00224
MSE2.046870.293520.168680.002890.02392
5040MLE3.260327.985091.937540.226861.60902
Bias0.26032−0.01491−0.062460.002150.00074
MSE0.594330.218770.072950.001420.01337
45MLE3.275497.959981.937060.229601.59955
Bias0.27549−0.04002−0.062940.00489−0.00873
MSE0.613350.200420.074510.001320.01260
10080MLE3.172477.975331.970100.230021.59641
Bias0.17247−0.02467−0.029900.00531−0.01187
MSE0.416010.171260.053370.000720.00807
90MLE3.189837.983821.957480.227951.60322
Bias0.18983−0.01618−0.042530.00324−0.00506
MSE0.394860.165810.049620.000680.00795
Table 4. The four cases for calculating the Bayes estimators via the standard Bayes and importance sampling techniques.
Table 4. The four cases for calculating the Bayes estimators via the standard Bayes and importance sampling techniques.
Case 1 When α
Is Unknown
Case 2 When β
Is Unknown
standard
Bayes
technique
The Bayes estimates of φ ( α ) under SE,
LINEX, and GE are obtained by computing
Equations (25), (27) and (29).
The Bayes estimates of φ ( β ) are obtained
numerically under the SE, LINEX, and GE loss functions
by evaluating (38), (39) and (40), respectively.
importance
sampling
technique
According to Algorithm 1, the Bayes estimates
of φ ( α ) are obtained under the SE, LINEX, and
GE loss functions by computing
Equations (31), (32) and (33), respectively.
Based on Algorithm 2, the Bayes estimates of
φ ( β ) are obtained numerically by computing
Equations (42)–(44) under three loss functions.
Case 3 When λ
Is Unknown
Case 4 When β
and λ Are Unknown
standard
Bayes
technique
The Bayes estimates of φ ( λ ) are obtained
under the SE, LINEX, and GE loss functions
by evaluating (49), (50) and (51), respectively.
importance
sampling
technique
Based on Algorithm 3, the Bayes estimates of
φ ( λ ) are obtained numerically by computing
Equations (53)–(55) under three loss functions.
The Bayes estimates of φ ( β , λ ) are obtained
numerically according to Algorithm 4
under the SE, LINEX, and GE loss functions by
computing Equations (70), (71) and (72), respectively.
Table 5. The ML and Bayesian estimates of the unknown parameters α , S ( x 0 ) , and h ( x 0 ) via the standard Bayes technique.
Table 5. The ML and Bayesian estimates of the unknown parameters α , S ( x 0 ) , and h ( x 0 ) via the standard Bayes technique.
nrParameters ML
Estimates
Bayes Estimates
SELINEXGE
τ = 0.001 τ = 2 τ = 5 q = 1 q = 3 q = 3
3024 α Mean0.818650.788550.788550.778280.766920.788550.754480.82343
Bias0.01865−0.01148−0.01148−0.02172−0.03308−0.01148−0.045520.02343
MSE0.018870.013530.013530.013620.013430.013530.015160.01556
S ( x 0 ) Mean0.731070.711390.711390.709880.709220.711390.701860.72454
Bias0.00203−0.01766−0.01766−0.01916−0.01982−0.01766−0.02718−0.00450
MSE0.004260.003620.003620.003840.004010.003620.004430.00340
h ( x 0 ) Mean0.807310.843510.843510.822780.788460.843510.796880.84694
Bias−0.010550.025660.025660.00492−0.029400.02566−0.020980.02908
MSE0.018480.014400.014400.014430.015700.014400.016650.01468
27 α Mean0.819170.792160.792160.780440.767010.792160.756750.82325
Bias0.01917−0.00784−0.00784−0.01956−0.03299−0.00784−0.043250.02325
MSE0.018340.013590.013590.013850.013590.013590.015280.01574
S ( x 0 ) Mean0.731550.713450.713450.710920.709210.713450.702960.72447
Bias0.00250−0.01559−0.01559−0.01813−0.01983−0.01559−0.02608−0.00457
MSE0.004150.003550.003550.003860.003980.003550.004430.00336
h ( x 0 ) Mean0.806520.839540.839540.820740.788780.839540.794850.84697
Bias−0.011330.021680.021680.00288−0.029080.02168−0.023000.02911
MSE0.017970.014250.014250.014690.015620.014250.017040.01467
5040 α Mean0.814290.800570.800570.779730.781160.800570.765340.81616
Bias0.014290.000570.00057−0.02028−0.018840.00057−0.034660.01616
MSE0.010370.008530.008530.009170.008780.008530.009940.00976
S ( x 0 ) Mean0.732150.722020.722020.713400.717590.722020.708690.72696
Bias0.00311−0.00702−0.00702−0.01564−0.01146−0.00702−0.02035−0.00209
MSE0.002410.002180.002180.002580.002460.002180.002820.00223
h ( x 0 ) Mean0.807890.825770.825770.828300.797680.825770.813160.83430
Bias−0.009970.007910.007910.01045−0.020170.00791−0.004690.01644
MSE0.010430.008970.008970.009740.010060.008970.010350.00950
45 α Mean0.812240.800570.800570.789360.781580.800570.775030.81643
Bias0.012240.000570.00057−0.01064−0.018420.00057−0.024970.01643
MSE0.010460.008530.008530.009190.008600.008530.009690.00958
S ( x 0 ) Mean0.731050.722020.722020.718350.717880.722020.713750.72721
Bias0.00201−0.00702−0.00702−0.01069−0.01116−0.00702−0.01529−0.00184
MSE0.002490.002180.002180.002430.002380.002180.002630.00216
h ( x 0 ) Mean0.810070.825770.825770.818290.797340.825770.802870.83380
Bias−0.007790.007910.007910.00044−0.020520.00791−0.014980.01594
MSE0.010660.008970.008970.009720.009820.008970.010660.00925
10080 α Mean0.804130.799390.799390.795280.788550.799390.788030.80636
Bias0.00413−0.00061−0.00061−0.00472−0.01145−0.00061−0.011970.00636
MSE0.005060.004680.004680.004540.004760.004680.004660.00496
S ( x 0 ) Mean0.729070.724860.724860.724000.722200.724860.721760.72698
Bias0.00003−0.00419−0.00419−0.00504−0.00684−0.00419−0.00728−0.00207
MSE0.001260.001190.001190.001180.001280.001190.001230.00121
h ( x 0 ) Mean0.815970.822980.822980.817160.809360.822980.809480.82824
Bias−0.001890.005120.00512−0.00079−0.008490.00512−0.008370.01038
MSE0.005320.004970.004970.004850.005230.004970.005090.00514
90 α Mean0.809250.799390.799390.789860.791950.799390.782640.80978
Bias0.00925−0.00061−0.00061−0.01014−0.00805−0.00061−0.017360.00981
MSE0.005580.004680.004680.004430.004280.004680.004630.00458
S ( x 0 ) Mean0.731480.724860.724860.721290.724140.724860.719030.72887
Bias0.00244−0.00419−0.00419−0.00775−0.00490−0.00419−0.01001−0.00017
MSE0.001340.001190.001190.001170.001150.001190.001230.00110
h ( x 0 ) Mean0.810880.822980.822980.822790.805650.822980.815250.82449
Bias−0.006980.005120.005120.00494−0.012210.00512−0.002610.00663
MSE0.005740.004970.004970.004680.004860.004970.004830.00464
Table 6. The ML and Bayesian estimates of the unknown parameters α , S ( x 0 ) , and h ( x 0 ) via the importance sampling technique.
Table 6. The ML and Bayesian estimates of the unknown parameters α , S ( x 0 ) , and h ( x 0 ) via the importance sampling technique.
nrParameters ML
Estimates
Bayes Estimates
SELINEXGE
τ = 0.001 τ = 2 τ = 5 q = 1 q = 3 q = 3
3024 α Mean0.818651.676611.676611.670451.676541.676611.670181.68018
Bias0.018650.876610.876610.870450.876540.876610.870180.88018
MSE0.018870.805990.805990.796220.808720.805990.795800.81512
S ( x 0 ) Mean0.731070.950510.950510.949890.950500.950510.949880.95056
Bias0.002030.221470.221470.220850.221460.221470.220830.22152
MSE0.004260.049440.049440.049170.049460.049440.049170.04948
h ( x 0 ) Mean0.807310.236780.236780.238750.235950.236780.237780.23714
Bias−0.01055−0.58108−0.58108−0.57910−0.58190−0.58108−0.58008−0.58072
MSE0.018480.342840.342840.340640.344110.342840.341750.34278
27 α Mean0.819171.592181.592181.582241.584601.592181.581981.58762
Bias0.019170.792180.792180.782240.784600.792180.781980.78726
MSE0.018340.660920.660920.645750.6522210.660920.645390.65639
S ( x 0 ) Mean0.731550.941330.941330.940170.940190.941330.940150.94024
Bias0.002500.212280.212280.211120.211140.212280.211100.21120
MSE0.004150.045530.045530.045050.045140.045530.045040.04516
h ( x 0 ) Mean0.806520.269850.269850.273690.272730.269850.272900.27374
Bias−0.01133−0.54800−0.54800−0.54417−0.54512−0.54800−0.54496−0.54412
MSE0.017970.305970.305970.301830.303620.305970.302680.30258
5040 α Mean0.814291.028881.028881.021991.028421.028881.021891.02878
Bias0.014290.228880.228880.221990.228420.228880.221890.22879
MSE0.010370.061840.061840.058130.061640.061840.058090.06180
S ( x 0 ) Mean0.732150.822550.822550.820340.822410.822550.820320.82245
Bias0.003110.093500.093500.091300.093370.093500.091280.09342
MSE0.002410.009820.009820.009380.009840.009820.009370.00985
h ( x 0 ) Mean0.807890.608540.608540.613870.608540.608540.613730.60884
Bias−0.00997−0.20932−0.20932−0.20399−0.20932−0.20932−0.20413−0.20902
MSE0.010430.050110.050110.047610.050250.050110.047670.05014
45 α Mean0.812241.017551.017551.017601.017951.017551.017501.01824
Bias0.012240.217550.217550.217600.217950.217550.217500.21824
MSE0.010460.056010.056010.055800.056970.056010.055770.05709
S ( x 0 ) Mean0.731050.818800.818800.818910.818740.818800.818890.81878
Bias0.002010.089760.089760.089860.089700.089760.089840.08974
MSE0.002490.009090.009090.009090.009200.009090.009090.00920
h ( x 0 ) Mean0.810070.617680.617680.617410.617360.617680.617280.61761
Bias−0.00779−0.20017−0.20017−0.20045−0.20050−0.20017−0.20058−0.20025
MSE0.010660.046040.046040.045990.046750.046040.046030.04666
10080 α Mean0.804130.646560.646560.647570.646850.646560.647560.64688
Bias0.00413−0.15344−0.15344−0.15243−0.15315−0.15344−0.15245−0.15312
MSE0.005060.025460.025460.025260.025310.025460.025270.02530
S ( x 0 ) Mean0.729070.640000.640000.640600.640220.640000.640590.64024
Bias0.00003−0.08904−0.08904−0.08844−0.08882−0.08904−0.08845−0.08880
MSE0.001260.008720.008720.008650.008650.008720.008650.00865
h ( x 0 ) Mean0.815970.989450.989450.988290.989030.989450.988280.98906
Bias−0.001890.171590.171590.170430.171170.171590.170420.17120
MSE0.005320.032130.032130.031870.031890.032130.031870.03190
90 α Mean0.809250.651270.651270.651300.649580.651270.651280.64961
Bias0.00925−0.14873−0.14873−0.14871−0.15042−0.14873−0.14872−0.15039
MSE0.005580.024000.024000.023920.024470.024000.023920.02446
S ( x 0 ) Mean0.731480.643050.643050.643100.641980.643050.643090.64199
Bias0.00244−0.08560−0.08560−0.08594−0.08706−0.08560−0.08595−0.08705
MSE0.001340.008150.008150.008120.008340.008150.008120.00833
h ( x 0 ) Mean0.810880.983850.983850.983770.985790.983850.983760.98583
Bias−0.006980.165990.165990.165920.167940.165990.165910.16797
MSE0.005740.030150.030150.030030.030780.030150.030030.03079
Table 7. The ML and Bayesian estimates of the unknown parameters β , S ( x 0 ) , and h ( x 0 ) via the standard Bayes technique.
Table 7. The ML and Bayesian estimates of the unknown parameters β , S ( x 0 ) , and h ( x 0 ) via the standard Bayes technique.
nrParameters ML
Estimates
Bayes Estimates
SELINEXGE
τ = 0.001 τ = 2 τ = 5 q = 1 q = 3 q = 3
3024 β Mean4.001373.633793.633793.210652.750233.633793.377393.84007
Bias0.00137−0.36621−0.36621−0.78935−1.24977−0.36621−0.62261−0.15993
MSE0.524940.473240.473240.828901.667930.473240.684870.38669
S ( x 0 ) Mean0.730310.749930.749930.746330.743340.749930.743830.74885
Bias0.001260.020890.020890.017290.014290.020890.014790.01981
MSE0.001300.001300.001300.001210.001090.001300.001170.00121
h ( x 0 ) Mean0.816920.748930.748930.738740.716190.748930.708150.78359
Bias−0.00094−0.06892−0.06892−0.07911−0.10167−0.06892−0.10971−0.03426
MSE0.017290.016040.016040.016900.019020.016040.022250.01299
27 β Mean3.958863.734693.734693.242112.821723.734693.397403.83642
Bias−0.04114−0.26531−0.26531−0.75789−1.17828−0.26531−0.60260−0.16358
MSE0.433950.414950.414950.773791.507510.414950.641070.39905
S ( x 0 ) Mean0.732190.744650.744650.746810.743030.744650.744580.74802
Bias0.003150.016610.016610.017770.013980.016610.015540.01897
MSE0.001080.001110.001110.001160.001110.001110.001120.00121
h ( x 0 ) Mean0.809410.767510.767510.738810.722060.767510.711380.78332
Bias−0.00845−0.05034−0.05034−0.07904−0.09580−0.05034−0.10648−0.03453
MSE0.014320.013950.013950.016160.018380.013950.020890.01338
5040 β Mean3.993253.800543.800543.446233.071343.800543.585433.90027
Bias−0.00675−0.19947−0.19947−0.55368−0.92866−0.19947−0.41457−0.00973
MSE0.314440.345570.345570.492960.990450.345570.415680.32374
S ( x 0 ) Mean0.730170.740840.740840.740460.738400.740840.73876074216
Bias0.001130.011790.011790.011420.009350.011790.009710.01311
MSE0.000770.000900.000900.000820.000850.000900.000800.00090
h ( x 0 ) Mean0.815930.779970.779970.765050.749270.779970.744970.79609
Bias−0.00193−0.03789−0.03789−0.05281−0.06858−0.03789−0.07289−0.02177
MSE0.010340.011520.011520.011220.013020.011520.013550.01074
45 β Mean4.003623.775223.775223.472413.147403.775223.599343.91113
Bias0.00362−0.22479−0.22479−0.52759−0.85260−0.22479−0.40066−0.08887
MSE0.292720.322600.322600.548040.851180.322600.389370.29091
S ( x 0 ) Mean0.729600.741920.741920.740830.737390.741920.739320.74077
Bias0.000560.012880.012880.011790.000840.012880.010280.01173
MSE0.000720.000850.000850.000780.000760.000850.000760.00079
h ( x 0 ) Mean0.817860.775520.775520.765050.756230.775520.747190.79844
Bias2.80599 × 10 6 −0.04233−0.04233−0.05281−0.06163−0.04233−0.07067−0.01942
MSE0.009610.010800.010800.010670.011450.010800.012710.00963
10080 β Mean4.021043.894323.894323.686363.411833.894323.778203.92103
Bias0.02104−0.10568−0.10568−0.31364−0.58817−0.10568−0.22180−0.0790
MSE0.185770.204500.204500.250760.466170.204500.227670.21109
S ( x 0 ) Mean0.728470.735360.735360.735070.735450.735360.734120.73752
Bias−0.000580.006320.006320.006020.006400.006320.005080.00847
MSE0.000450.000510.000510.000510.000540.000510.000500.00056
h ( x 0 ) Mean0.821260.797730.797730.789710.775310.797730.778910.80144
Bias0.00340−0.02012−0.02012−0.02814−0.04255−0.02012−0.03895−0.01641
MSE0.006090.006760.006760.006810.007790.006760.007460.00699
90 β Mean4.025203.909553.909553.722223.478513.909553.805603.94136
Bias0.02520−0.09045−0.09045−0.27778−0.52149−0.09045−0.19440−0.05864
MSE0.170380.175020.175020.209350.381750.175020.190190.18087
S ( x 0 ) Mean0.728220.734470.734470.734240.734140.734470.733400.73599
Bias−0.000820.005420.005420.005200.005090.005420.004360.00695
MSE0.000420.0004390.0004390.000420.000460.0004390.000420.00047
h ( x 0 ) Mean0.822050.800620.800620.793330.781990.800620.783770.80538
Bias0.00419−0.01724−0.01724−0.02453−0.03586−0.01724−0.03409−0.01248
MSE0.005590.005780.005780.005710.006530.005780.006230.00597
Table 8. The ML and Bayesian estimates of the unknown parameters β , S ( x 0 ) , and h ( x 0 ) via the importance sampling technique.
Table 8. The ML and Bayesian estimates of the unknown parameters β , S ( x 0 ) , and h ( x 0 ) via the importance sampling technique.
nrParameters ML
Estimates
Bayes Estimates
SELINEXGE
τ = 0.001 τ = 2 τ = 5 q = 1 q = 3 q = 3
3024 β Mean4.001375.015285.015284.354813.706185.015284.680735.06025
Bias0.001371.015281.015280.35481−0.293821.015280.680731.06025
MSE0.524941.735841.735840.482660.292701.735841.008831.82667
S ( x 0 ) Mean0.730310.683940.683940.684970.684180.683940.682410.68929
Bias0.00126−0.04510−0.04510−0.04408−0.04487−0.04510−0.04663−0.03975
MSE0.001300.003410.003410.003180.003370.003410.003460.00285
h ( x 0 ) Mean0.816920.997700.997700.967020.933350.997700.943901.00351
Bias−0.000940.179840.179840.149160.115500.179840.126050.18565
MSE0.017290.054450.054450.039990.029750.054450.033390.05635
27 β Mean3.958864.880514.880514.217143.682354.880514.413884.57499
Bias−0.041140.880510.880510.21714−0.317650.880510.413880.57499
MSE0.433951.422621.422620.463460.304341.422620.758680.92222
S ( x 0 ) Mean0.732190.689240.689240.704140.7041 80.689240.703200.70628
Bias0.00315−0.03980−0.03980−0.02491−0.02486−0.03980−0.05840−0.02276
MSE0.001080.002880.002880.001960.001870.002880.002030.00173
h ( x 0 ) Mean0.809410.974330.974330.906140.8888180.974330.893880.91906
Bias−0.008450.156480.156480.088290.070330.156480.076020.10120
MSE0.014320.044850.044850.027100.021060.044850.024720.02903
5040 β Mean3.993254.930704.930704.294033.904634.930704.443504.63432
Bias−0.006750.930700.930700.29403−0.095370.930700.443500.63432
MSE0.314441.437521.437520.426350.234291.437520.632560.95327
S ( x 0 ) Mean0.730170.686720.686720.703820.700970.686720.703070.70263
Bias0.00113−0.04233−0.04233−0.02522−0.02807−0.04233−0.02597−0.02642
MSE0.000770.002930.002930.001640.001930.002930.001690.00182
h ( x 0 ) Mean0.815930.983460.983460.908160.905750.983460.898850.93003
Bias−0.001930.165600.165600.090310.087890.165600.080990.11218
MSE0.010340.045400.045400.022480.023190.045400.020640.03003
45 β Mean4.003624.304244.304244.229213.853494.304244.274844.11707
Bias0.003620.304240.304240.22921−0.146510.304240.274840.11707
MSE0.292720.531690.531690.444640.294070.531690.495910.41340
S ( x 0 ) Mean0.729600.715340.715340.715160.686140.724600.714940.72503
Bias0.00056−0.01370−0.01370−0.01388−0.00444−0.01370−0.01409−0.00402
MSE0.000720.001180.001180.001160.000960.001180.001170.00096
h ( x 0 ) Mean0.817860.871840.871840.870120.831950.871840.867230.83792
Bias2.80599 × 10 6 0.053980.053980.052260.014090.053980.049370.02006
MSE0.009610.017080.017080.016360.012970.017080.016050.01342
10080 β Mean4.021044.522754.522754.311803.977094.522754.330794.07013
Bias0.021040.522750.522750.31180−0.022910.522750.330790.07013
MSE0.185770.612020.612020.401220.222900.612020.420780.25983
S ( x 0 ) Mean0.728470.704620.704620.713010.726350.704620.712910.72651
Bias−0.00058−0.02442−0.02442−0.01603−0.00270−0.02442−0.01613−0.00253
MSE0.000450.001330.001330.000950.000610.001330.000950.00061
h ( x 0 ) Mean0.821260.911370.911370.878450.827620.911370.877220.82990
Bias0.003400.093520.093520.060590.009760.093520.059360.01204
MSE0.006090.019570.019570.013680.008370.019570.013540.00845
90 β Mean4.025203.712603.712603.725493.853243.712603.728113.87665
Bias0.02520−0.28740−0.28740−0.27451−0.14676−0.28740−0.27189−0.12335
MSE0.170380.268200.268200.291510.234230.268200.290720.23222
S ( x 0 ) Mean0.728220.744010.744010.743150.735800.744010.743130.73586
Bias−0.000820.014970.014970.014100.006760.014970.014080.00681
MSE0.000420.000700.000700.000750.000580.000700.000750.00058
h ( x 0 ) Mean0.822050.765090.765090.768200.794210.765090.767920.79492
Bias0.00419−0.05276−0.05276−0.04966−0.02365−0.05276−0.04994−0.02294
MSE0.005590.008960.008960.009660.007710.008960.009690.00768
Table 9. The ML and Bayesian estimates of the unknown parameters λ , S ( x 0 ) , and h ( x 0 ) via the standard Bayes technique.
Table 9. The ML and Bayesian estimates of the unknown parameters λ , S ( x 0 ) , and h ( x 0 ) via the standard Bayes technique.
nrParameters ML
Estimates
Bayes Estimates
SELINEXGE
τ = 0.001 τ = 2 τ = 5 q = 1 q = 3 q = 3
3024 λ Mean3.023502.889902.889902.800942.688372.889902.825932.92667
Bias0.02350−0.11010−0.11010−0.19906−0.31163−0.11010−0.17407−0.07333
MSE0.102610.121870.121870.136650.182710.121870.134980.11895
S ( x 0 ) Mean0.727330.693500.693500.689270.684030.693500.680350.70038
Bias−0.00171−0.03554−0.03554−0.03977−0.04501−0.03554−0.04870−0.02867
MSE0.003870.006050.006050.006390.007000.006050.007590.00534
h ( x 0 ) Mean0.891500.898450.898450.876060.841290.898450.844590.92167
Bias0.001640.080590.080590.058200.023490.080590.026730.10381
MSE0.022370.033860.033860.029300.025090.033860.028450.03753
27 λ Mean3.032972.776392.776392.685122.580702.776392.705132.79645
Bias0.03297−0.22361−0.22361−0.31488−0.41930−0.22361−0.29487−0.20355
MSE0.105380.123450.123450.167580.235140.123450.160400.11978
h ( x 0 ) Mean0.815310.952720.952720.936460.902210.952720.907940.98287
Bias−0.002550.134860.134860.118610.084350.134860.090090.16501
MSE0.022850.038880.038880.035250.027130.038880.030850.04856
S ( x 0 ) Mean0.729060.670940.670940.663740.657720.670940.654090.67485
Bias0.000016−0.05811−0.05811−0.06530−0.07132−0.05811−0.07495−0.05420
MSE0.003950.007060.007060.008220.009110.007060.009890.00664
5040 λ Mean3.005953.045943.045942.977492.902493.045942.997633.07076
Bias0.005950.045940.04594−0.02251−0.097510.04594−0.002370.07076
MSE0.052810.090920.090920.086490.084300.090920.090530.09503
S ( x 0 ) Mean0.726880.728750.728750.724800.723550.728750.720220.73277
Bias−0.00217−0.00029−0.00029−0.00424−0.00550−0.00029−0.008820.00373
MSE0.002110.003250.003250.003470.003340.003250.003700.00309
h ( x 0 ) Mean0.821740.815010.815010.805230.780350.815010.784250.82975
Bias0.00389−0.00284−0.00284−0.01262−0.03751−0.00284−0.033610.01189
MSE0.012090.018980.018980.019400.018720.018980.021170.01862
45 λ Mean3.005902.907532.907532.848202.788882.907532.864442.94098
Bias0.00590−0.09247−0.09247−0.15180−0.21112−0.09247−0.13556−0.05902
MSE0.051500.068030.068030.078860.097910.068030.076950.06776
S ( x 0 ) Mean0.726960.702640.702640.699160.699090.702640.694150.70881
Bias−0.00208−0.02641−0.02641−0.02989−0.02995−0.02641−0.03490−0.02023
MSE0.002020.003290.003290.003530.003700.003290.003980.00304
h ( x 0 ) Mean0.821580.878260.878260.866760.838160.878260.847770.88764
Bias0.003720.060400.060400.048910.020310.060400.029920.06979
MSE0.011610.018390.018390.016860.015000.018390.016020.02017
10080 λ Mean2.995863.145473.145473.115693.066873.145473.127953.16113
Bias−0.004140.145470.145470.115690.066870.145470.127950.16113
MSE0.010980.090930.090930.083840.065090.090930.088810.09265
S ( x 0 ) Mean0.727520.750510.750510.749980.748660.750510.747970.75305
Bias−0.001520.021470.021470.020930.019620.021470.018920.02400
MSE0.000440.002610.002610.002730.002540.002610.002710.00264
h ( x 0 ) Mean0.821230.763480.763480.755260.744680.763480.744050.76983
Bias0.00337−0.05438−0.05438−0.06560−0.07318−0.05438−0.07380−0.04803
MSE0.002500.015840.015840.017270.017380.015840.019050.01483
90 λ Mean2.998282.981792.981792.959652.920422.981792.969363.00405
Bias−0.00172−0.01821−0.01821−0.04035−0.07958−0.01821−0.030640.00405
MSE0.011960.037350.037350.036850.039980.037350.037080.03716
S ( x 0 ) Mean0.727940.721140.721140.721530.720510.721140.719280.72524
Bias−0.00111−0.00791−0.00791−0.00752−0.00853−0.00791−0.00977−0.00380
MSE0.000480.001520.001520.001500.001540.001520.001580.00143
h ( x 0 ) Mean0.820200.835150.835150.824410.812320.835150.814330.83768
Bias0.002350.017300.017300.00656−0.005540.01730−0.003530.01982
MSE0.002730.008700.008700.008190.008050.008700.008350.00864
Table 10. The ML and Bayesian estimates of the unknown parameters λ , S ( x 0 ) , and h ( x 0 ) via the importance sampling technique.
Table 10. The ML and Bayesian estimates of the unknown parameters λ , S ( x 0 ) , and h ( x 0 ) via the importance sampling technique.
nrParameters ML
Estimates
Bayes Estimates
SELINEXGE
τ = 0.001 τ = 2 τ = 5 q = 1 q = 3 q = 3
3024 λ Mean3.023504.107364.107364.033283.813244.107364.101854.16576
Bias0.023501.107361.107361.033280.813241.107361.101851.16576
MSE0.102611.490581.490581.347560.868031.490581.530821.65682
S ( x 0 ) Mean0.727330.870670.870670.872770.869820.870670.871010.87427
Bias−0.001710.141630.141630.143730.140770.141630.141960.14523
MSE0.003870.022350.022350.023280.022320.022350.022880.02339
h ( x 0 ) Mean0.891500.449800.449800.432130.425170.449800.410850.46252
Bias0.00164−0.36805−0.36805−0.38573−0.39268−0.36805−0.40701−0.35534
MSE0.022370.152920.152920.168070.171190.152920.185400.14449
27 λ Mean3.032973.891313.891313.780473.604963.891313.833323.90260
Bias0.032970.891310.891310.780470.604960.891310.833320.90260
MSE0.105381.023431.023430.818880.505811.023430.927421.02966
S ( x 0 ) Mean0.729060.849460.849460.846930.845660.849460.844810.85056
Bias0.0000160.120420.120420.117890.116620.120420.115770.12151
MSE0.003950.017090.017090.016810.016160.017090.016410.01718
h ( x 0 ) Mean0.815310.508060.508060.502680.491260.508060.483510.52755
Bias−0.00255−0.30980−0.30980−0.31518−0.32660−0.30980−0.33434−0.29030
MSE0.022850.114750.114750.119520.123640.114750.132730.10155
5040 λ Mean3.005953.359953.359953.358873.344643.359953.364143.38831
Bias0.005950.359950.359950.358870.344640.359950.364140.38831
MSE0.052810.252500.252500.254400.235300.252500.259590.27492
S ( x 0 ) Mean0.726880.785640.785640.786930.788880.785640.786390.79005
Bias−0.002170.056600.056600.057890.059830.056600.057340.06100
MSE0.002110.006180.006180.006320.006370.006180.006280.00649
h ( x 0 ) Mean0.821740.067630.067630.670580.662430.067630.667940.66911
Bias0.00389−0.14155−0.14155−0.14728−0.15543−0.14155−0.14991−0.14874
MSE0.012090.038640.038640.040220.041600.038640.041060.03969
45 λ Mean3.005903.260073.260073.285373.256343.260073.290453.29795
Bias0.005900.260070.260070.285370.256340.260070.290450.29795
MSE0.051500.161480.161480.182130.161700.161480.186210.19122
S ( x 0 ) Mean0.726960.770770.770770.776260.775290.770770.775670.77656
Bias−0.002080.041730.041730.047220.046240.041730.046620.04751
MSE0.002020.004270.004270.004870.004840.004270.004840.00493
h ( x 0 ) Mean0.821580.713910.713910.697510.696320.713910.694810.70328
Bias0.00372−0.10395−0.10395−0.12035−0.12153−0.10395−0.12305−0.11458
MSE0.011610.026300.026300.030710.031160.026300.031430.02963
10080 λ Mean2.995862.885662.885662.885832.889622.885662.886042.89182
Bias−0.00414−0.11434−0.11434−0.11417−0.11038−0.11434−0.11396−0.10818
MSE0.010980.059510.059510.057560.056760.059510.057530.05624
S ( x 0 ) Mean0.727520.702380.702380.702640.703690.702380.702580.70382
Bias−0.00152−0.02666−0.02666−0.02640−0.02535−0.02666−0.02646−0.02522
MSE0.000440.002770.002770.002720.002600.002770.002730.00259
h ( x 0 ) Mean0.821230.880260.880260.879430.876570.880260.879240.87720
Bias0.003370.062400.062400.061580.058710.062400.061380.05934
MSE0.002500.015510.015510.015190.014490.015510.015160.01460
90 λ Mean2.998282.807762.807762.814872.802152.807762.815062.80426
Bias−0.00172−0.19224−0.19224−0.18513−0.19785−0.19224−0.18494−0.19574
MSE0.011960.074590.074590.071040.078270.074590.070980.07743
S ( x 0 ) Mean0.727940.685840.685840.687610.684810.685840.687540.68495
Bias−0.00111−0.04320−0.04320−0.04143−0.04424−0.04320−0.04151−0.04409
MSE0.000480.003730.003730.003520.003920.003730.003530.00390
h ( x 0 ) Mean0.820200.919540.919540.915120.921280.919540.914930.92194
Bias0.002350.101680.101680.097260.103420.101680.097070.10408
MSE0.002730.020700.020700.019470.021560.020700.019430.02172
Table 11. The ML and Bayesian estimates of the unknown parameters β , λ , S ( x 0 ) , and h ( x 0 ) via the importance sampling technique.
Table 11. The ML and Bayesian estimates of the unknown parameters β , λ , S ( x 0 ) , and h ( x 0 ) via the importance sampling technique.
nrParameters ML
Estimates
Bayes Estimates
SELINEXGE
τ = 0.001 τ = 2 τ = 5 q = 1 q = 3 q = 3
3024 β Mean4.166075.282855.282855.069674.525765.282855.200785.04745
Bias0.166071.282851.282851.069670.525761.282851.200781.04745
MSE1.674641.723151.723151.194080.296291.723151.507111.15652
λ Mean3.034994.261974.261974.170483.950534.261974.216494.20813
Bias0.034991.261971.261971.170480.950531.261971.216491.20813
MSE0.252311.750371.750371.520031.017331.750371.641891.61824
S ( x 0 ) Mean0.730950.862710.862710.860640.857850.862710.859120.86186
Bias0.001910.133670.133670.131600.128810.133670.130080.13281
MSE0.003810.019490.019490.019090.018310.019490.018760.01924
h ( x 0 ) Mean0.805480.508120.508120.501850.486380.508120.478480.52651
Bias−0.01238−0.30973−0.30973−0.31601−0.33147−0.30973−0.33937−0.29135
MSE0.025860.110720.110720.114960.122780.110720.130280.09937
27 β Mean3.930865.206255.206255.028414.506135.206255.139474.96516
Bias−0.069141.206251.206251.028410.506131.206251.139470.96516
MSE0.996791.544811.544811.122310.284071.544811.379231.00406
λ Mean2.964854.214934.214934.128653.966554.214934.169674.20868
Bias−0.035161.214931.214931.128650.966561.214931.169671.20868
MSE0.225401.650261.650261.412921.044951.650261.517041.61560
S ( x 0 ) Mean0.725080.858820.858820.857560.860060.858820.856080.86378
Bias−0.003960.129780.129780.128520.131010.129780.127040.13473
MSE0.004150.018670.018670.018200.018770.018670.017870.01966
h ( x 0 ) Mean0.807820.517790.517790.510010.480250.517790.487210.51821
Bias−0.01003−0.30007−0.30007−0.30785−0.33761−0.30007−0.29964−0.28406
MSE0.023410.106510.106510.108940.126340.106510.123540.10339
5040 β Mean4.051594.661344.661344.580934.33284.661344.633204.69727
Bias0.051590.661340.661340.580930.433290.661340.633200.69727
MSE0.939070.541130.541130.433730.266930.541130.502200.59661
λ Mean2.99303.448723.448723.438083.429553.448723.443683.47011
Bias−0.006990.448720.448720.438080.429550.448720.443680.47011
MSE0.154510.283420.283420.274960.262210.283420.280920.30397
S ( x 0 ) Mean0.726590.777150.777150.776420.777880.777150.775350.78020
Bias−0.003280.048110.048110.047380.048830.048110.046300.05115
MSE0.002590.004610.004610.004590.004660.004610.004530.00484
h ( x 0 ) Mean0.815810.730210.730210.725260.712160.730210.716070.73270
Bias−0.00205−0.08764−0.08764−0.09259−0.10570−0.08764−0.10178−0.08516
MSE0.016640.024750.024750.025550.027240.024750.027380.02423
45 β Mean4.091124.459874.459874.390464.296154.459874.422714.47300
Bias0.091120.459870.459870.390460.296150.459870.422710.47300
MSE0.853130.302710.302710.242520.166490.302710.270750.31575
λ Mean3.013253.432713.432713.406643.376513.432713.413453.42173
Bias0.015320.432710.432710.406640.376510.432710.413450.42173
MSE0.151420.273120.273120.250980.224520.273120.257660.26411
S ( x 0 ) Mean0.728280.781360.781360.779340.777360.781360.778160.77985
Bias−0.000760.052310.052310.050300.048320.052310.049110.05081
MSE0.002630.005030.005030.004890.004630.005030.004810.00481
h ( x 0 ) Mean0.814540.710170.710170.707420.703280.710170.697740.72367
Bias−0.00332−0.10768−0.10768−0.11044−0.11458−0.10768−0.12011−0.09419
MSE0.015490.028200.028200.028790.028320.028200.031080.02514
10080 β Mean4.050714.598254.598254.563383.788554.598254.582273.86137
Bias0.050700.598250.598250.56338−0.211450.598250.58227−0.13863
MSE0.726170.474150.474150.424890.120620.474150.446980.09656
λ Mean2.991593.156423.156423.137283.104753.156423.138013.11545
Bias−0.008410.156420.156420.137280.104750.156420.138010.11545
MSE0.098570.065580.065580.059250.054180.065580.059470.05666
S ( x 0 ) Mean0.727030.729960.729960.726300.754050.729960.725840.75520
Bias−0.002010.000920.00092−0.002740.025010.00092−0.003200.026153
MSE0.001340.001690.001690.001810.002250.001690.001820.00228
h ( x 0 ) Mean0.817730.850160.850160.856490.740840.850160.853340.75007
Bias−0.000120.032310.032310.03864−0.077020.032310.03549−0.06778
MSE0.008530.013690.013690.014880.016340.013690.014590.01537
90 β Mean4.133164.126494.126494.110494.100474.126494.118534.14749
Bias0.133160.126490.126490.110490.100470.126490.118530.14749
MSE0.631630.090530.090530.086380.089430.090530.088110.10049
λ Mean3.036713.147053.147053.137273.124283.147053.138633.13364
Bias0.036710.147050.147050.137270.124280.147050.138630.13364
MSE0.086240.069070.069070.064790.059530.069070.065320.06197
S ( x 0 ) Mean0.731300.748220.748220.746850.744360.748220.746400.74528
Bias0.002250.019180.019180.017810.015320.019180.017350.01624
MSE0.001230.002070.002070.002040.001880.002070.002030.00189
h ( x 0 ) Mean0.814490.778040.778040.778740.782280.778040.775560.78958
Bias−0.00398−0.03982−0.03982−0.03911−0.03558−0.03982−0.04230−0.02828
MSE0.007690.013140.013140.013150.012280.013140.013420.01209
Table 12. Descriptive statistics of the fatigue life of 6061-T6 aluminum coupons’ cut data.
Table 12. Descriptive statistics of the fatigue life of 6061-T6 aluminum coupons’ cut data.
MeasureValueMeasureValue
n102Minimum223
Maximum560Mean397.88
Q1352Q3439
Median400Skewness−0.003
Kurtosis2.85Variance3884.30
Standard deviation62.32
Table 13. ML and Bayesian estimates of the unknown parameters β and λ via the standard Bayes technique for the fatigue life of the 6061-T6 aluminum coupons’ cut data.
Table 13. ML and Bayesian estimates of the unknown parameters β and λ via the standard Bayes technique for the fatigue life of the 6061-T6 aluminum coupons’ cut data.
ModelrParametersML
Estimates
Bayes Estimates
SELINEXGE
τ = 0.001 τ = 2 τ = 5 q = 1 q = 3 q = 3
BIED82 λ 411.27803.84233.84233.12722.49383.84233.34744.0774
β 1.52560.17210.17210.17170.17100.17210.16720.1745
92 λ 432.79104.32854.32853.51842.80214.32853.83184.5658
β 1.831220.19680.19680.19630.19550.19680.19170.1993
Table 14. ML estimates of the model parameters and the statistics of the AIC, BIC, CAIC, HQIC, and for the fatigue life of the 6061-T6 aluminum coupons’ cut data.
Table 14. ML estimates of the model parameters and the statistics of the AIC, BIC, CAIC, HQIC, and for the fatigue life of the 6061-T6 aluminum coupons’ cut data.
ModelsrML EstimatesAICBICCAICHQIC
α β λ
BIED8244.783920.9497150.0270−151.342308.684316.559308.929311.873
9243.279720.0114148.1150−172.933351.867359.742352.112355.056
IED82424.3190−248.584499.168501.793499.208500.231
92404.397−288.006578.013580.638578.053579.076
WIED8252.35130.67962540.0300−159.233324.467332.342324.712327.655
9253.50410.64382660.4100−180.559367.117374.992367.362370.306
IWD82138.18801.0897−293.464590.927596.177591.049593.053
92148.29701.2775−321.958647.916653.166648.037650.042
WED820.00040.00473.9449−151.652309.303317.178309.548312.492
920.00640.01061.1240−179.223364.447372.321364.691367.635
OFIED823.5426252.0410−157.345318.689323.939318.810320.815
923.7253250.2300−180.356364.713369.963364.834366.839
Table 15. Descriptive statistics of the survival times (in weeks) of the patients suffering from acute myelogenous leukemia data.
Table 15. Descriptive statistics of the survival times (in weeks) of the patients suffering from acute myelogenous leukemia data.
MeasureValueMeasureValue
n33Minimum1
Maximum156Mean40.88
Q14Q365
Median22Skewness1.16
Kurtosis3.12Variance2181.17
Standard deviation46.70
Table 16. ML and Bayesian estimates of the unknown parameters β and λ via the standard Bayes technique for the survival times (in weeks) of patients suffering from acute myelogenous leukemia data.
Table 16. ML and Bayesian estimates of the unknown parameters β and λ via the standard Bayes technique for the survival times (in weeks) of patients suffering from acute myelogenous leukemia data.
ModelrParametersML
Estimates
Bayes Estimates
SELINEXGE
τ = 0.001 τ = 2 τ = 5 q = 1 q = 3 q = 3
BIED26 λ 4.43131.75851.75851.51461.26421.75851.40251.9172
β 0.45240.36220.36220.35630.34810.36220.32980.3786
30 λ 4.82181.90921.90921.64831.37881.90921.56262.0654
β 0.51420.40580.40580.39930.39020.40580.37370.4220
Table 17. ML estimates of the model parameters and the statistics of the AIC, BIC, CAIC, HQIC, and for the patients suffering from acute myelogenous leukemia data.
Table 17. ML estimates of the model parameters and the statistics of the AIC, BIC, CAIC, HQIC, and for the patients suffering from acute myelogenous leukemia data.
ModelsrML EstimatesAICBICCAICHQIC
α β λ
BIED262.56990.46721.2091−40.40484.80887.80185.20885.815
305.59480.52820.5915−56.103116.205119.198116.605117.212
IED266.0189−46.88895.77797.27395.90696.280
306.0164−61.848125.695127.191125.824126.198
WIED260.37380.56735.1250−41.70189.40293.89290.23090.913
300.35350.57884.9131−57.983121.965126.455122.793123.476
IWD268.63920.6227−40.67287.34491.83488.17288.855
308.24900.6557−56.731119.462123.952120.29120.973
WED260.42690.44310.0478−44.42294.84499.33395.67196.354
300.49850.34790.0390−61.213128.426132.916129.254129.937
OFIED260.42023.6035−42.00788.01491.00788.41489.021
300.46003.3522−58.573121.146124.139121.546122.153
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Aldahlan, M.A.; Bakoban, R.A.; Alzahrani, L.S. On Estimating the Parameters of the Beta Inverted Exponential Distribution under Type-II Censored Samples. Mathematics 2022, 10, 506. https://0-doi-org.brum.beds.ac.uk/10.3390/math10030506

AMA Style

Aldahlan MA, Bakoban RA, Alzahrani LS. On Estimating the Parameters of the Beta Inverted Exponential Distribution under Type-II Censored Samples. Mathematics. 2022; 10(3):506. https://0-doi-org.brum.beds.ac.uk/10.3390/math10030506

Chicago/Turabian Style

Aldahlan, Maha A., Rana A. Bakoban, and Leena S. Alzahrani. 2022. "On Estimating the Parameters of the Beta Inverted Exponential Distribution under Type-II Censored Samples" Mathematics 10, no. 3: 506. https://0-doi-org.brum.beds.ac.uk/10.3390/math10030506

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