Next Article in Journal
An Amended Whale Optimization Algorithm for Optimal Bidding in Day Ahead Electricity Market
Next Article in Special Issue
Inferences of a Mixture Bivariate Alpha Power Exponential Model with Engineering Application
Previous Article in Journal
Some Common Fixed-Circle Results on Metric Spaces
Previous Article in Special Issue
Inference of Reliability Analysis for Type II Half Logistic Weibull Distribution with Application of Bladder Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian and Non-Bayesian Analysis of Exponentiated Exponential Stress–Strength Model Based on Generalized Progressive Hybrid Censoring Process

by
Manal M. Yousef
1,†,
Amal S. Hassan
2,†,
Huda M. Alshanbari
3,*,†,
Abdal-Aziz H. El-Bagoury
4,† and
Ehab M. Almetwally
5,6,†
1
Department of Mathematics, Faculty of Science, New Valley University, El-Khargah 72511, Egypt
2
Faculty of Graduate Studies for Statistical Research, Cairo University, Giza 12613, Egypt
3
Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
4
Basic Science Department, Higher Institute of Engineering and Technology, El-Mahala El-Kobra 6734723, Egypt
5
Faculty of Business Administration, Delta University of Science and Technology, Gamasa 11152, Egypt
6
The Scientific Association for Studies and Applied Research, Al Manzalah 35646, Egypt
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 28 July 2022 / Revised: 21 August 2022 / Accepted: 26 August 2022 / Published: 5 September 2022
(This article belongs to the Special Issue Computational Statistics & Data Analysis)

Abstract

:
In many real-life scenarios, systems frequently perform badly in difficult operating situations. The multiple failures that take place when systems reach their lower, higher, or extreme functioning states typically receive little attention from researchers. This study uses generalized progressive hybrid censoring to discuss the inference of R = P ( X < Y < Z ) for a component when it is exposed to two stresses, Y , Z , and it has one strength X that is regarded. We assume that both the stresses and strength variables follow an exponentiated exponential distribution with a common scale parameter. We obtain R’s maximum likelihood estimator and approximate confidence intervals. In addition, the Bayesian estimators for symmetric, such as squared error, and asymmetric loss functions, such as linear exponential, are developed. Credible intervals with the highest posterior densities are established. Monte Carlo simulations are used to evaluate and compare the effectiveness of the many proposed estimators. The process is then precisely described using an analysis of real data.

1. Introduction

The stress–strength system is one of the most widely used data analysis methods in a variety of fields, including industrial engineering, military applications, health, and applied sciences. The stress–strength system reliability is the evaluation of a component’s reliability in terms of the random variable X that represents the stress of the component that is exposed to Y that represents the component’s strength available to overcome the possible stress. When the stress exceeds the strength of the system, it fails. Due to its practical applications, Bhattacharyya and Johnson [1] were the first researchers to be interested in investigating and deriving the reliability of the stress–strength model. In recent years, a lot of works have been conducted on the problem of estimating the stress–strength model; see, for example, Ahmad et al. [2], Saraçoglu et al. [3], Hassan et al. [4], Kotb and Raqab [5], and Jana and Bera [6].
This paper is interested in the model R = P ( X < Y < Z ) that discusses the case where the strength Y should not only be greater than stress X but also be smaller than stress Z. The model R = P ( X < Y < Z ) is necessary when devices cease to function under extreme lower as well as extreme upper stress operating environments. The main idea of this model was first introduced by Chandra and Owen [7]. For example, a person’s blood pressure has two limits, systolic and diastolic, and their blood pressure should lie within these limits. As another example, some electrical components malfunction when positioned below and above a specific power generator. In a similar vein, many gadgets are non-functional in both high and low temperatures. The R = P ( X < Y < Z ) stress–strength models have a wide range of applications in engineering, psychology, genetics, clinical trials, and other fields; see Kotz et al. [8].
In the literature, Chandra and Owen [7] developed maximum likelihood estimators (MLEs) and uniform minimum unbiased estimators (UMVUEs) for R = P ( X < Y < Z ) . The lowest variance unbiased and the MLEs of R = P ( X < Y < Z ) were investigated by Singh [9], where X , Y , and Z are mutually independent random variables that follow the normal distribution. Dutta and Sriwastav [10] addressed the estimation of R when X , Y , and Z are exponentially distributed. Ahmad et al. [11] discussed the comparative inference on reliability estimation for a multi-component stress–strength model under power Lomax distribution. Almetwally et al. [12] developed an optimal censoring plan for multi-stress—strength reliability based on progressive first failure utilising Bayesian and non-Bayesian methods. Ivshin [13] presented the MLE and UMVUE of R when X , Y , and Z are either uniform or exponential random variables. The estimation of R = P ( X < Y < Z ) for the Weibull distribution in the presence of k outliers have been provided by Hassan et al. [14]. The estimation of the stress–strength reliability for a new exponential inverted Topp–Leone distribution has been provided by Metwally et al. [15]. Estimation by using various estimation methods of the stress–strength reliability of the exponentiated inverted Weibull distribution has been obtained by Abu El Azm et al. [16]. The inference of a fuzzy stress–strength reliability model for the inverse Rayleigh distribution has been discussed by Sabry et al. [17].
Under the assumption that the three samples were independent, Wang et al. [18] used nonparametric normal approximations and the Jackknife empirical likelihood to make a statistical inference for R. Patowary et al. [19] used the Monte Carlo simulation to study the technique of the reliability estimation for R = P ( X < Y < Z ) of n standby systems ( n = 1 , 2 ) . Yousif et al. [20] proposed the estimation reliability for R = P ( X < Y < Z ) using the exponentiated inverse Rayleigh distribution. Hameed et al. [21] discussed the estimation of R = P ( Y 1 < X < Y 2 ) using the inverse Kumaraswamy distribution. Attia and Karam [22] studied the Bayesian estimation of R = P ( X < Y < Z ) for the Dagum distribution. Some estimation methods of R = P ( X < Y < Z ) for the inverse Rayleigh distribution were considered by Raheem et al. [23] and Abd Elfattah and Taha [24]. Yousef and Almetwally [25] provided different estimates for a multi-stress–strength model when data are observed from the Kumaraswamy distribution.
Life testing studies, which can be characterized as mathematical and statistical models of survival analysis, are frequently employed in engineering, biology, mechanical, and other disciplines of research. In fact, we are unable to observe the failure time of all the units due to a variety of constraints, such as time and cost. Before all the observations fail, it is typical to stop in the middle of the process. Type-I and Type-II censoring schemes are the two most common censoring schemes among all censoring cases (see, for example, Meeker and Escobar [26]).
Before the final termination, an inevitable pause or loss of the experiment units is likely to occur. The constraint in those two censoring techniques, however, is that the units cannot be removed during the trial. Cohen [27] initially proposed a progressive censoring strategy to address this inflexibility. According to the progressive Type-II censored schemes defect, if the experimental units are highly trustworthy, this experiment will last a long period. As a result, Kundu and Joarder [28] proposed the progressively hybrid censoring scheme. Moreover, El-Sherpieny et al. [29] introduced progressive Type-I and Type-II hybrid censored schemes based on the maximum product spacing method as an alternative estimation method for this scheme. The implementation of n independent identical distributed units is employed for the censoring scheme. At min ( T , X m ) , the experimenter will terminate the operation. The time T as well as 1 m n are predetermined in this case. In the context of the progressive type-II censored approach, the experiment span will not be more than T. However, the observations we acquired would be insufficient if the predetermined termination time T is small. Cho et al. [30] presented a new censoring scheme called the generalized progressive hybrid censoring (GPHC) scheme, which allows us to obtain a specified series of failures.
The estimation of the parameters for the Gompertz distribution was obtained using the ML and the Bayesian methods under different loss functions based on Mohie El-Din et al. [31]. Tu and Gui [32] considered the estimation for the Kumaraswamy distribution under the GPHC. Nagy et al. [33] used a GPHC sample from the Burr XII distribution to estimate the unknown parameters, reliability, and hazard functions. Maswadah [34] improved the ML estimation method using the Runge–Kutta technique.
The goal of this paper is to establish the reliability inferences for the stress–strength variables that follow the exponentiated exponential distributions (EEDs) with a similar scale parameter. As far as we know, no earlier work has attempted to estimate R using an EED under the GPHC scheme. As a result, the estimations of R’s statistical inference based on data from a GPHC scheme will be discussed in this study. The lifetime distributions of the one strength Y and two stresses, X and Z, are considered to have independent EEDs with the same scale parameter. MLEs as well as Bayesian estimators under different loss functions of R are derived. The asymptotic confidence intervals as well as Bayesian credible intervals are created. To compare the different methodologies, extensive simulations were run, and dataset was evaluated for demonstration reasons.
This paper proceeds as follows: First, data description and reliability model are mentioned in Section 2. The MLEs as well as the approximate confidence intervals of R will be derived in Section 3. We calculate the Bayesian estimators with different loss functions in Section 4. Furthermore, we apply the Markov chain Monte Carlo (MCMC) algorithm to derive the Bayesian estimators and establish the highest posterior density (HPD) intervals under the sample generated by the MCMC algorithm in Section 5. Then, in Section 6, the data analysis is demonstrated. Finally, conclusions are arranged in the last Section 7.

2. Data Description and Reliability Model

Ahuja and Nash [35] presented the two-parameter EED which was subsequently investigated by Gupta and Kundu [36]. The shape and scale parameters of this distribution are similar to those of the gamma and Weibull distributions. In many circumstances, it has a better fit than the Weibull and Gamma distributions (Raja and Mir [37]). The applications of the EED were substantial, and we point out: models to determine bout standards for analysis of animal conduct (Yeates et al. [38]); software reliability increase fashions for vital best metrics (Subburaj et al. [39]); models for episode height and period for eco-hydro-climatic packages (Biondi et al. [40]); and therapy price modeling (Kannan et al. [41]). The probability density function (PDF) and cumulative distribution function (CDF) of EED are defined as
f ( x ) = α λ e λ x ( 1 e λ x ) α 1 ; x > 0 , α , λ > 0 ,
and
F ( x ) = ( 1 e λ x ) α ,
respectively, for x > 0 , α > 0 and λ > 0 . We write X EED( λ , α ), and here, α is the shape parameter and λ is the scale parameter.
Let X EED( λ , α 1 ), Y EED( λ , α 2 ), and Z EED( λ , α 3 ), and they are independent. The reliability formula of the stress–strength model that the probability of a component strength falling in between two stresses is given by:
R = 0 x x z f ( x ) f ( y ) f ( z ) d y d z d x , = α 2 α 3 ( α 1 + α 2 ) ( α 1 + α 2 + α 3 ) .
The reliability R of stress–strength model with different values of parameters, see Figure 1.
Figure 1 shows the reliability R of the stress–strength model where the higher the value of parameter α 3 , the higher the reliability value. We also note that with the change in the parameter values, the reliability rating changes, and in most matters, it gives high results and covers most of the values.
Assume that our research group is made of n independent units with the same lifetime distribution, where X 1 , X 2 , , X n represent the corresponding lifetime. The integers k and m, k < m , as well as R 1 , R 2 , , R m , which can satisfy the equation i = 1 m R i + m = n function as preplanned integers, have been under predetermination between zero and n. On the arrival of the first failure X 1 , we withdraw R 1 units. When the second failure, X 2 , occurs, we remove R 2 units at random from the n 2 R 1 survivors. With the rest of the survival units removed, the procedure is repeated and ended at T * = max ( min ( T , X m ) , X k ) . It vastly improved prior approaches by allowing us to choose whether or not to continue the experiment if the sample size is insufficient at the predetermined cut-off time T. Researchers would prefer to obtain m failures under the GPHC scheme, but they can alternatively choose k failures, which are considered the bare minimum. The GPHC scheme is referred to as R 1 , R 2 , , R m . Let D be the observed failure times before arriving at the predefined time T.
The GPHC scheme can be classified into the following categories:
Case 1: X 1 , , X d , , X k for T < X k < X m ,
Case 2: X 1 , , X k , , X d for X k < T < X m ,
Case 3: X 1 , , X k , , X m for X k < X m < T ,
According to the GPHC scheme, the joint density function for three different cases is as follows:
f X ( x ) = A * i = 1 D f ( x i ; m ; n ) [ F ¯ ( x i ; m ; n ) ] R i * [ F ¯ ( T ) ] R τ * ,
where
D = k If T < X k , m , n < X m , m , n d If X k , m , n < T < X m , m , n m If X k , m , n < X m , m , n < T , A * = i = 1 D j = i m ( R j * + 1 ) ,
R * = ( R 1 , , R d , 0 , , 0 , R k * = n k j = 1 d R j ) If T < X k , m , n < X m , m , n ( R 1 , , R d ) If X k , m , n < T < X m , m , n ( R 1 , , R m ) If X k , m , n < X m , m , n < T
with R τ * the number of surviving units that are removed at T, given by
R τ * = 0 If T < X k , m , n < X m , m , n n d j = 1 d R j If X k , m , n < T < X m , m , n 0 If X k , m , n < X m , m , n < T .

3. Maximum Likelihood Estimation

The ML procedure is a popular and effective strategy used by statisticians when dealing with reliability issues and survival analysis. This section provides the MLE of R and the approximate CI is established.

3.1. Maximum Likelihood Estimator of R

Here, the MLE of R is determined, first by obtaining the MLEs of α 1 , α 2 , α 3 , and λ . Plugging the PDF and CDF of the EED, i.e., (1) and (2), into the likelihood formula (4), the likelihood function of α 1 , α 2 , α 3 , and λ is expressed as:
l ( α 1 , α 2 , α 3 , λ ) = A 1 * A 2 * A 3 * α 1 D 1 α 2 D 2 α 3 D 3 λ D 1 + D 2 + D 3 × j = 1 D 1 e λ x j ( 1 e λ x j ) α 1 1 ( 1 ( 1 e λ x j ) α 1 ) R j 1 ( 1 ( 1 e λ T 1 ) α 1 ) R d 1 + 1 * × j = 1 D 2 e λ y j ( 1 e λ y j ) α 2 1 ( 1 ( 1 e λ y j ) α 2 ) R j 2 ( 1 ( 1 e λ T 2 ) α 2 ) R d 2 + 1 * × j = 1 D 3 e λ z j ( 1 e λ z j ) α 3 1 ( 1 ( 1 e λ z j ) α 3 ) R j 3 ( 1 ( 1 e λ T 3 ) α 3 ) R d 3 + 1 * .
Taking the logarithm for (5), we obtain the log-likelihood function
L = ln l ( α 1 , α 2 , α 3 , λ ) = ln A 1 * A 2 * A 3 * + D 1 ln α 1 + D 2 ln α 2 + D 3 ln α 3 + ( D 1 + D 2 + D 3 ) ln λ λ [ j = 1 D 1 x j + j = 1 D 2 y j + j = 1 D 3 z j ] + ( α 1 1 ) j = 1 D 1 ln ψ 1 ( x j , λ ) + ( α 2 1 ) j = 1 D 2 ln ψ 2 ( y j , λ ) + ( α 3 1 ) j = 1 D 3 ln ψ 3 ( z j , λ ) + j = 1 D 1 R j 1 ln ( 1 ( ψ 1 ( x j , λ ) ) α 1 ) + j = 1 D 2 R j 2 ln ( 1 ( ψ 2 ( y j , λ ) ) α 2 ) + j = 1 D 3 R j 3 ln ( 1 ( ψ 3 ( z j , λ ) ) α 3 ) + R d 1 + 1 * ln ( 1 ( ψ 1 ( T 1 , λ ) ) α 1 ) + R d 2 + 1 * ln ( 1 ( ψ 2 ( T 2 , λ ) ) α 2 ) + R d 3 + 1 * ln ( 1 ( ψ 3 ( T 3 , λ ) ) α 3 )
where
ψ 1 ( x j , λ ) = 1 e λ x j , ψ 2 ( y j , λ ) = 1 e λ y j , ψ 3 ( z j , λ ) = 1 e λ z j ,
and ψ 1 ( T 1 , λ ) is as given by (7) with x j = T 1 . Similarly for ψ 2 ( T 2 , λ ) and ψ 3 ( T 3 , λ ) . We take the partial derivatives of (6) for α 1 , α 2 , α 3 , and λ , respectively, and obtain a set of likelihood equations as follows:
L α 1 = D 1 α 1 + j = 1 D 1 [ ln ψ 1 ( x j , λ ) R j 1 E 1 ( ψ 1 ( x j , λ ) , α 1 ) ] R d 1 + 1 * E 1 ( ψ 1 ( T 1 , λ ) , α 1 ) , L α 2 = D 2 α 2 + j = 1 D 2 [ ln ψ 2 ( y j , λ ) R j 2 E 2 ( ψ 2 ( y j , λ ) , α 2 ) ] R d 2 + 1 * E 2 ( ψ 2 ( T 2 , λ ) , α 2 ) , L α 3 = D 3 α 3 + j = 1 D 3 [ ln ψ 3 ( z j , λ ) R j 3 E 3 ( ψ 3 ( z j , λ ) , α 3 ) ] R d 3 + 1 * E 3 ( ψ 3 ( T 3 , λ ) , α 3 ) , L λ = D 1 + D 2 + D 3 λ [ j = 1 D 1 x j + j = 1 D 2 y j + j = 1 D 3 z j ] + ( α 1 1 ) j = 1 D 1 Q 1 ( ψ 1 ( x j , λ ) , λ ) + ( α 2 1 ) j = 1 D 2 Q 2 ( ψ 2 ( y j , λ ) , λ ) + ( α 3 1 ) j = 1 D 3 Q 3 ( ψ 3 ( z j , λ ) , λ ) j = 1 D 1 R j 1 S 1 ( ψ 1 ( x j , λ ) , λ ) j = 1 D 2 R j 2 S 2 ( ψ 2 ( y j , λ ) , λ ) j = 1 D 3 R j 3 S 3 ( ψ 3 ( z j , λ ) , λ ) R d 1 + 1 * S 1 ( ψ 1 ( T 1 , λ ) , λ ) R d 2 + 1 * S 2 ( ψ 2 ( T 2 , λ ) , λ ) R d 3 + 1 * S 3 ( ψ 3 ( T 3 , λ ) , λ ) ,
where
E i ( ψ i ( w j , λ ) , α i ) = ( ψ i ( w j , λ ) ) α i 1 ( ψ i ( w j , λ ) ) α i ln ψ i ( w j , λ ) , E i ( ψ i ( T i , λ ) , α i ) = ( ψ i ( T i , λ ) ) α i 1 ( ψ i ( T i , λ ) ) α i ln ψ i ( T i , λ ) , Q i ( ψ i ( w j , λ ) , λ ) = w j e λ w j ψ i ( w j , λ ) , S i ( ψ i ( w j , λ ) , λ ) = α i ( ψ i ( w j , λ ) ) α i 1 ( ψ i ( w j , λ ) ) α i Q i ( ψ i ( w j , λ ) , λ ) ,
where w = x , y , z , i = 1 , 2 , 3 .
There are no analytical solutions for the unknown parameters α 1 , α 2 , α 3 , and λ from (8). As a result, (8) may be maximized using a fairly straightforward iterative approach such as the Newton–Raphson (NR) procedure to obtain the appropriate MLEs of α 1 , α 2 , α 3 , and λ . The corresponding precise distribution (or exact confidence intervals (CIs) ) of α 1 , α 2 , α 3 , and λ cannot be obtained because the MLEs of α 1 , α 2 , α 3 , and λ cannot be obtained in closed form. For this purpose, we recommend implementing the ‘maxLik’ package for any given dataset ( α 1 , α 2 , α 3 , λ ) and then performing ML calculations using the NR iterative method. We can substitute these estimates to obtain the MLE of R as:
R ^ = α ^ 2 α ^ 3 ( α ^ 1 + α ^ 2 ) ( α ^ 1 + α ^ 2 + α ^ 3 ) .

3.2. Asymptotic Confidence Interval

CIs, which contain the population parameter with a given high probability, show the level of uncertainty in an estimate derived from sample data. Utilizing the large-sample normal distribution of the MLE is the most popular technique for establishing confidence bounds for the parameters. For a large sample size, MLEs are consistent and asymptotically normally distributed.
Here, the asymptotic CI of R is determined using the asymptotic distribution of R ^ , which was obtained from the asymptotic distribution of α 1 , α 2 , α 3 , and λ . The observed Fisher information matrix is denoted by
I ( θ ) = [ I i j ] = 2 L θ i θ j , i , j = 1 , 2 , 3 .
The elements in lines of the I ( θ ) matrix can be obtained by differentiating from (8) two times with respect to α 1 , α 2 , α 3 , and λ :
I 11 = D 1 α 1 2 + j = 1 D 1 R j 1 E 1 ( ψ 1 ( x j , λ ) , α 1 ) α 1 + R d 1 + 1 * E 1 ( ψ 1 ( T 1 , λ ) , α 1 ) α 1 , I 22 = D 2 α 2 2 + j = 1 D 2 R j 2 E 2 ( ψ 2 ( y j , λ ) , α 2 ) α 2 + R d 2 + 1 * E 2 ( ψ 2 ( T 2 , λ ) , α 2 ) α 2 , I 33 = D 3 α 3 2 + j = 1 D 3 R j 3 E 3 ( ψ 3 ( z j , λ ) , α 3 ) α 3 + R d 3 + 1 * E 3 ( ψ 3 ( T 3 , λ ) , α 3 ) α 3 , I 14 = j = 1 D 1 Q 1 ( ψ 1 ( x j , λ ) R j 1 E 1 ( ψ 1 ( x j , λ ) , α 1 ) λ R d 1 + 1 * E 1 ( ψ 1 ( T 1 , λ ) , α 1 ) λ , I 24 = j = 1 D 2 Q 2 ( ψ 2 ( y j , λ ) R j 2 E 2 ( ψ 2 ( y j , λ ) , α 2 ) λ R d 2 + 1 * E 2 ( ψ 2 ( T 2 , λ ) , α 2 ) λ , I 34 = j = 1 D 3 Q 3 ( ψ 3 ( z j , λ ) R j 3 E 3 ( ψ 3 ( z j , λ ) , α 3 ) λ R d 3 + 1 * E 3 ( ψ 3 ( T 3 , λ ) , α 3 ) λ ,
I 44 = D 1 + D 2 + D 3 λ 2 + ( α 1 1 ) j = 1 D 1 Q 1 ( ψ 1 ( x j , λ ) , λ ) λ + ( α 2 1 ) j = 1 D 2 Q 2 ( ψ 2 ( y j , λ ) , λ ) λ , + ( α 3 1 ) j = 1 D 3 Q 3 ( ψ 3 ( z j , λ ) , λ ) λ j = 1 D 1 R j 1 S 1 ( ψ 1 ( x j , λ ) , λ ) λ j = 1 D 2 R j 2 S 2 ( ψ 2 ( y j , λ ) , λ ) λ , j = 1 D 3 R j 3 S 3 ( ψ 3 ( z j , λ ) , λ ) λ R d 1 + 1 * S 1 ( ψ 1 ( T 1 , λ ) , λ ) λ R d 2 + 1 * S 2 ( ψ 2 ( T 2 , λ ) , λ ) λ R d 3 + 1 * S 3 ( ψ 3 ( T 3 , λ ) , λ ) λ , I 12 = I 21 = I 13 = I 31 = I 32 = I 23 = 0 ,
where
E i ( ψ i ( w j , λ ) , α i ) α i = ( E i ( ψ i ( w j , λ ) , α i ) ) 2 ( ψ i ( w j , λ ) ) α i , E i ( ψ i ( T i , λ ) , α i ) α i = ( E i ( ψ i ( T i , λ ) , α i ) ) 2 ( ψ i ( T i , λ ) ) α i , E i ( ψ i ( w j , λ ) , α i ) λ = S i ( ψ i ( w j , λ ) , λ ) [ 1 α i + ln ψ i ( w j , λ ) + E i ( ψ i ( w j , λ ) , α i ) ] , Q i ( ψ i ( w j , λ ) , λ ) λ = Q i ( ψ i ( w j , λ ) , λ ) ( ( w j + Q i ( ψ i ( w j , λ ) , λ ) ) , S i ( ψ i ( w j , λ ) , λ ) λ = S i ( ψ i ( w j , λ ) , λ ) [ ( α i 1 ) Q i ( ψ i ( w j , λ ) , λ ) + S i ( ψ i ( w j , λ ) , λ ) w j ] ,
where w = x , y , z , i = 1 , 2 , 3 .
As D 1 , D 2 , D 3 , then
[ n 1 ( α ^ 1 α 1 ) , n 2 ( α ^ 2 α 2 ) , n 3 ( α ^ 3 α 3 ) ] N ( 0 , I 1 ( α 1 , α 2 , α 3 ) ) ,
where
I 1 ( α 1 , α 2 , α 3 ) = 1 I 11 0 0 0 1 I 22 0 0 0 1 I 33 .
The asymptotic distribution of R ^ is found as
[ ( R ^ R ) ] N ( 0 , σ R 2 ) ,
where σ R 2 = M T I 1 M is the asymptotic variance of R, and
M = R α 1 , R α 2 , R α 3 ,
where
R α 1 = α 2 α 3 ( 2 ( α 1 + α 2 ) + α 3 ) ( α 1 + α 2 ) 2 ( α 1 + α 2 + α 3 ) 2 , R α 2 = α 3 ( α 1 ( α 1 + α 3 ) α 2 2 ) ( α 1 + α 2 ) 2 ( α 1 + α 2 + α 3 ) 2 , R α 3 = α 2 ( α 1 + α 2 + α 3 ) 2 .
Therefore, a 100 ( 1 γ ) % asymptotic CI of R can be constructed as,
( R ^ z 1 γ 2 σ ^ R 2 , R ^ + z 1 γ 2 σ ^ R 2 ) ,
where z γ 2 is the 100 γ th standard normal percentile and σ ^ R is the standard deviation of the MLE of R.

4. Bayesian Estimation

In this section, the Bayesian estimation of R is obtained when data are observed using the GPHC based on the squared error loss function (SELF) and linear exponential (Linex) loss function, which are defined, respectively, by
L 1 = ( ϑ , ϑ ˇ ) = ( ϑ ˇ ϑ ) 2 , L 2 = ( ϑ , ϑ ˇ ) = e g ( ϑ ˇ ϑ ) g ( ϑ ˇ ϑ ) 1 ,
where ϑ ˇ is an estimator of ϑ . Denote the prior and posterior distributions of ϑ by π ( ϑ ) and π * ( ϑ x ̲ ) , respectively. Under the SELF and Linex loss function, the Bayesian estimation of any function B ( ϑ ) of ϑ is given by
B ( ϑ ) S E = E [ B ( ϑ ) x ̲ ] = 0 B ( ϑ ) π * ( ϑ x ̲ ) d ϑ , B ( ϑ ) L N = 1 c ln E ( e c B ( ϑ ) ) = 1 c ln 0 e c B ( ϑ ) π * ( ϑ x ̲ ) d ϑ .
Prior distribution is important for the development of Bayes estimators.
Under the assumption of gamma prior distributions, we investigate at this estimate problem. The gamma family of distributions is particularly adaptable and can be thought of as an appropriate prior for any unknown parameter(s) because it has a wide variety of shapes depending on its parameter values; for further information, see references Kundu [42] and Dey et al. [43]. Furthermore, the independent gamma priors are simple and transparent, which may help to avoid many challenging inferential problems.
Therefore, it is assumed here that α 1 , α 2 , α 3 , and λ follow independent gamma distributions with α 1 G ( η 1 , ζ 1 ) , α 2 G ( η 2 , ζ 2 ) , α 3 G ( η 3 , ζ 3 ) , and λ G ( η 4 , ζ 4 ) , with probability densities given by, respectively,
π ( α i ) = α i η i 1 Γ ( η i ) ζ i η i e α i ζ i , π ( λ ) = λ η 4 1 Γ ( η 4 ) ζ 4 η 4 e λ ζ 4 , α i > 0 , η i , ζ i > 0 , i = 1 , 2 , 3 .
Using the informative prior (10) and the likelihood function (5), the joint posterior density can be derived as follows:
π * ( α 1 , α 2 , α 3 , λ ) = i = 1 3 A i * α i D i + η i 1 Γ ( η i ) ζ i η i λ η 4 1 + i = 1 3 D i Γ ( η 4 ) ζ 4 η 4 e λ ζ 4 × e i = 1 D 1 λ x i + α 1 ζ 1 ( α 1 1 ) ln ψ 1 ( x i , λ ) R i 1 ln ( 1 ( ψ 1 ( x i , λ ) ) α 1 ) + R d 1 + 1 * ln ( 1 ( ψ 1 ( T 1 , λ ) ) α 1 ) × e i = 1 D 2 λ y i + α 2 ζ 2 ( α 2 1 ) ln ψ 2 ( y i , λ ) R i 2 ln ( 1 ( ψ 2 ( y i , λ ) ) α 2 ) + R d 2 + 1 * ln ( 1 ( ψ 2 ( T 2 , λ ) ) α 2 ) × e i = 1 D 3 λ z i + α 3 ζ 3 ( α 3 1 ) ln ψ 3 ( z i , λ ) R i 3 ln ( 1 ( ψ 3 ( z i , λ ) ) α 3 ) + R d 3 + 1 * ln ( 1 ( ψ 3 ( T 3 , λ ) ) α 3 ) .
Then, the Bayes estimates of R under the SELF and Linex loss function, say R ˜ S E and R ˜ L N are given by
R ˜ S E = 0 0 0 R π * ( α 1 , α 2 , α 3 , λ ) d α 1 d α 2 d α 3 , R ˜ L N = 1 c 0 0 0 e c R π * ( α 1 , α 2 , α 3 , λ ) d α 1 d α 2 d α 3 .
Obviously, it is not possible to compute (11) analytically. The MCMC approaches can be applied to approximate (11).
The marginal posterior densities of the parameters α 1 , α 2 , α 3 , and λ can be derived as
π * ( α 1 ) α 1 D 1 + η 1 1 e i = 1 D 1 α 1 ( 1 ζ 1 ln ψ 1 ( x i , λ ) ) R i 1 ln ( 1 ( ψ 1 ( x i , λ ) ) α 1 ) + R d 1 + 1 * ln ( 1 ( ψ 1 ( T 1 , λ ) ) α 1 ) π * ( α 2 ) α 2 D 2 + η 2 1 e i = 1 D 2 α 2 ( 1 ζ 2 ln ψ 2 ( y i , λ ) ) R i 2 ln ( 1 ( ψ 2 ( y i , λ ) ) α 2 ) + R d 2 + 1 * ln ( 1 ( ψ 2 ( T 2 , λ ) ) α 2 ) π * ( α 3 ) α 3 D 3 + η 3 1 e i = 1 D 3 α 3 ( 1 ζ 3 ln ψ 3 ( z i , λ ) ) R i 3 ln ( 1 ( ψ 3 ( z i , λ ) ) α 3 ) + R d 3 + 1 * ln ( 1 ( ψ 3 ( T 3 , λ ) ) α 3 ) π * ( λ ) λ η 4 1 + i = 1 3 D i e λ ζ 4 e i = 1 D 1 λ x i ( α 1 1 ) ln ψ 1 ( x i , λ ) R i 1 ln ( 1 ( ψ 1 ( x i , λ ) ) α 1 ) + R d 1 + 1 * ln ( 1 ( ψ 1 ( T 1 , λ ) ) α 1 ) × e i = 1 D 2 λ y i ( α 2 1 ) ln ψ 2 ( y i , λ ) R i 2 ln ( 1 ( ψ 2 ( y i , λ ) ) α 2 ) + R d 2 + 1 * ln ( 1 ( ψ 2 ( T 2 , λ ) ) α 2 ) × e i = 1 D 3 λ z i ( α 3 1 ) ln ψ 3 ( z i , λ ) R i 3 ln ( 1 ( ψ 3 ( z i , λ ) ) α 3 ) + R d 3 + 1 * ln ( 1 ( ψ 3 ( T 3 , λ ) ) α 3 ) ,
The marginal posterior densities in (12) are not well-known distributions, so we will use the Metropolis–Hastings (MH) sampler to generate the values of α 1 , α 2 , α 3 , and λ with normal proposal distribution to generate samples from it in (12)
Furthermore, the approach of Chen and Shao [44] was extensively used to construct HPD intervals with unknown benefit distribution parameters for Bayesian estimate. For example, a 95 % HPD interval can be created using two endpoints from the MCMC sample outputs: 2.5 % and 97.5 % percentiles, respectively. The Θ parameters’ Bayes, trustworthy intervals are calculated as follows:
  • Sorted parameters as α ˜ l [ 1 ] < α ˜ l [ 2 ] < < α ˜ l [ N ] ; l = 1 , 2 , 3 , λ ˜ [ 1 ] < λ ˜ [ 2 ] < < λ ˜ [ N ] , and R [ 1 ] < R [ 2 ] < < R [ N ] , and N is the length of MCMC generated.
  • The 95 % symmetric credible intervals of α ˜ 1 , α ˜ 2 , α ˜ 3 , λ ˜ , and R ˜ become α ˜ l L 25 1000 , α ˜ l L 975 1000 , λ ˜ L 25 1000 , λ ˜ L 975 1000 , and R ˜ L 25 1000 , R ˜ L 975 1000 .

5. Simulation Study

We give some simulation results in this section to show how the various strategies presented in this paper perform in practice. Different situations have been used as:
  • Case I: α 1 = 0.5 , α 2 = 4 , α 3 = 20 , λ = 0.5 and T 1 = 14 , T 2 = 6 , T 3 = 18 .
  • Case II: α 1 = 0.8 , α 2 = 5 , α 3 = 12 , λ = 1.5 and T 1 = 1.6 , T 2 = 2 , T 3 = 2.5 .
  • Case III: α 1 = 0.3 , α 2 = 2 , α 3 = 15 , λ = 0.5 and T 1 = 3.5 , T 2 = 3 , T 3 = 6 .
  • Case IV: α 1 = 0.7 , α 2 = 3 , α 3 = 8 , λ = 2 and T 1 = 1 , T 2 = 1.5 , T 3 = 1.9 .
We have considered different sample sizes ( n ) for each sample as n 1 = 20 , n 2 = 25 , n 3 = 15 , and n 1 = 30 , n 2 = 40 , n 3 = 30 ; different effective sample sizes ( m ) for each sample as m 1 = 15 , m 2 = 18 , m 3 = 12 , and m 1 = 17 , m 2 = 22 , m 3 = 14 ; different k values for each sample as k 1 = 12 , k 2 = 16 , k 3 = 10 , and k 1 = 15 , k 2 = 20 , k 3 = 12 ; and two different progressive censoring schemes, namely Scheme-I R 1 = ( n 1 m 1 , r e p ( 0 , m 1 1 ) ) , R 2 = ( n 2 m 2 , r e p ( 0 , m 2 1 ) ) , R 3 = ( n 3 m 3 , r e p ( 0 , m 3 1 ) ) and Scheme-II R 1 = ( 2 , r e p ( 0 , ( m 1 3 ) / 2 ) , n 1 m 1 3 , r e p ( 0 , ( m 1 3 ) / 2 ) , 1 ) , R 2 = ( 1 , r e p ( 0 , ( m 2 2 ) / 2 ) , n 2 m 2 2 , r e p ( 0 , ( m 2 4 ) / 2 ) , 1 ) , R 3 = ( 1 , r e p ( 0 , ( m 3 3 ) / 2 ) , n 3 m 3 2 , r e p ( 0 , ( m 3 1 ) / 2 ) , 1 ) , where r e p points to replicate censored scheme. In each case, we compute and construct the 95% exact and the credible CIs for both loss functions. The methods are repeated 5000 times, and the average estimators (AvE), mean squared errors (MSE), average length of CIs (L.CI) with related coverage percentages (CP), and average length of HPD credible CIs (L.CCI) are reported.
For prior distribution, the hyperparameters are chosen using elective hyperparameters based on mean and variance of gamma prior distribution.
To find out how to elicit hyperparameters of the independent joint prior, we can utilise the likelihood method’s estimate and variance–covariance matrix. The mean and variance of gamma priors can be used to represent the derived hyperparameters.
η j = 1 N i = 1 N ϑ j i ^ 2 1 N 1 i = 1 N ϑ j i ^ 1 N i = 1 N ϑ j i ^ 2 ; j = 1 , , p ,
ζ j = 1 N i = 1 N ϑ j i ^ 1 N 1 i = 1 N ϑ j i ^ 1 N i = 1 N ϑ j i ^ 2 ; j = 1 , , p ,
where N is the number of simulation iterations.
For MCMC techniques, we replicate the process 10,000 times of MH algorithms. Figure 2 and Figure 3 are explained by Table 1 that discusses different label of these figures as MSE.i.1 is MLE, MSE.i.2 is SELF, MSE.i.3 is Linex c = −0.5, and MSE.i.4 is Linex c = 1.5, where i = 1 , , 8 . The X label indicates different parameters as param1 is α 1 , param2 is α 2 , param3 is α 3 , and param4 is λ . On the left side of Table 2, Table 3, Table 4 and Table 5, the average of the frequentist estimates and the average of the Bayesian MCMC estimates of α 1 , α 2 , α 3 , and λ are determined. Additionally, on the right side of Table 2, Table 3, Table 4 and Table 5, we report the L.CI and L.CCI of α 1 , α 2 , α 3 , and λ that were obtained. Two highly regarded R software packages are used to conduct extensive evaluations: the ‘CODA’ package for computing Bayes estimates using MCMC techniques and the ‘maxLik’ tool for computing ML estimates using the NR algorithm.
The heatmap plots of the MSE of simulation results for parameters and R of Cases I, II, III, and IV are provided in Figure 2 and Figure 3, respectively. Using one variable on each axis, heatmaps are used to display relationships between two variables (MSE and methods). You can determine if there are any trends in the values for one or both variables by monitoring how cell colors vary across each axis. MSE was represented by cell coloring, with the darker color representing a higher MSE value than the lighter color. The results are presented in Table 2, Table 3, Table 4 and Table 5, and the following observations can be made:
  • ML and Bayesian estimates of population parameters are quite good based on AvE where they tend to actual values.
  • As the sample size increases, the MSE decreases as expected for ML and Bayesian estimations.
  • For a given sample size, the MSE also declines with m.
  • If n and m are held constant, the MSEs decrease as the acceptable bare minimum of failures, k, rises.
  • Based on AvE, MSE, and length of CI, we note scheme II is better than scheme I in same times.
  • Bayesian estimates outperform ML estimates in terms of AvE, MSE, and length of CI because they incorporate prior knowledge based on a gamma informative prior.
  • When the symmetric and asymmetric loss functions are compared, Bayes estimates under the asymmetric loss function are found to be more accurate than others under the symmetric loss function.
  • The average length of HPD credible intervals for Bayesian estimation is better than length of asymptotic confidence interval.

6. Application of Transformer Insulation Data

This section examines an actual dataset in order to demonstrate the methodologies proposed in the previous sections. This dataset is also used to display the EED based on the GPHC method. The findings of a stress–strength life test of transformer insulation were published in chapter three of Nelson’s book [45]. The test included three levels of voltage, which were 35:4 kv, 42:4 kv, and 46:7 kv, respectively, with a normal voltage of 14:4 kv. When 42:4 kv, the data are 0.6, 13.4, 15.2, 19.9, 25.0, 30.2, 32.8, 44.4, and 56.2. When 46:7 kv, the data are 3.1, 8.3, 8.9, 9.0, 13.6, 14.9, 16.1, 16.9, 21.3, and 48.1. When 35:4 kv, the data are 40.1, 59.4, 71.2, 166.5, 204.7, 229.7, 308.3, and 537.9. We want to calculate the reliability of the stress–strength model P ( X < Y < Z ) . Firstly, Table 6 provides the values of the estimate with different measures as the Kolmogorov–Smirnov (KS) statistic along with its p-value, Akaike information criterion (AIC), Bayesian information criterion (BIC), corrected AIC (CAIC), and Hannan–Quinn information criterion (HQIC) for the EED based on each sample of the dataset. The EED fits each dataset according to the KS test and the accompanying p-value. We have also included two plots based on the estimated model. On the left, the estimated and the empirical CDF of the EED are displayed, and on the right, the P-P plot estimated of the EED is displayed, see Figure 4, Figure 5 and Figure 6 for each dataset, respectively. Figure 4, Figure 5 and Figure 6 prove the fitting of each sample to show how each sample fits into the EED. The ML and Bayesian estimators via the complete sample are listed in Table 7. We notice that the standard error (SE) of the Bayesian estimates is lower than others for the ML estimates. Thus, the Bayesian estimation is the optimum estimation of the parameters for the EED based on the GPHC. The reliability of P ( X < Y < Z ) of the Bayesian estimation method is higher than that of the ML process, proving the stated conclusion.
It is first assumed that X E E D ( λ , α 1 ) , Y E E D ( λ , α 2 ) , and Z E E D ( λ , α 3 ) . Following that, these are the ML estimates for the unknown parameters: α 1 = 0.6089, α 2 = 0.4964, α 3 = 3.5237, and λ = 0.0131; and L 0 = 129.9995 for the corresponding log-likelihood value. Second, assume that X E E D ( λ 1 , α 1 ) , Y E E D ( λ 2 , α 2 ) , and Z E E D ( λ 3 , α 3 ) . The following are the ML estimates for the unknown parameters: α 1 = 1.3049, α 2 = 4.7632, α 3 = 1.8156, λ 1 = 0.0443, λ 2 = 0.1768, and λ 3 = 0.0071; and L 1 = 116.6313 for the corresponding log-likelihood value. The following tests of the hypotheses are conducted by us:
H 0 : λ 1 = λ 2 = λ 3 ; H 1 : λ 1 λ 2 λ 3
Then, the likelihood ratio tests value is 2 ( L 0 L 1 ) = −26.73644 in this instance, with degrees of freedom as 2 (difference in the number of parameters for the two models). Moreover, the p-value of the Chi-squared test is one; therefore, the null hypothesis cannot be rejected. As a result, H 0 is a valid assumption in this situation.
Figure 7 depicts the trace and normal curve of the posterior distribution for the MCMC estimate of the stress–strength model for the EED based on the GPHC. The MCMC samples are shown as a pairs plot in Figure 8, which displays the pairwise relationship between the parameters in the top plot, the correlation coefficients in the bottom plot, and the marginal frequency distribution for each parameter on the diagonal. Moreover, as shown in Figure 9, for the reliability stress–strength estimate for the EED based on the complete sample, convergence begins at 2000 iterations or fewer. The MCMC samples are displayed as a pairs plot in Figure 8, which depicts the pairwise relationship between the parameters as independent, with the scatter plots matrix in the top plot, the correlation coefficients in the bottom plot, and the marginal frequency distribution for each parameter on the diagonal. The parameters P3 and P4 are shown to be medially connected in this diagram, where p1 is an α 1 , p2 is an α 2 , p3 is an α 3 , and p4 is a λ .
For each component of this model, we suggested using the following GPHC sample: X = (0.6, 13.4, 15.2, 19.9, 25.0, 30.2, 32.8, 44.4), Y = (3.1, 8.3, 8.9, 9.0, 13.6, 14.9, 16.1), Z = (40.1, 59.4, 71.2, 166.5, 204.7, 229.7, 308.3), R 1 = (0, 0, 0, 0, 0, 0, 0, 1, 0), R 2 = (0, 0, 0, 0, 0, 0, 0, 0, 1), and R 3 = (0, 0, 0, 0, 0, 0, 0, 2).
The ML and Bayesian estimations of the EED parameters based on the GPHC sample of the stress–strength model are discussed in Table 8. The trace and normal curve of the posterior distribution for the MCMC estimate of the stress–strength for the EED based on the GPHC are shown in Figure 10. The MCMC samples are given as a pairs plot in Figure 11, which illustrates the pairwise relationship between the parameters in the top plot, the correlation coefficients in the bottom plot, and the marginal frequency distribution for each parameter on the diagonal. Moreover, as shown in Figure 12, for the reliability stress–strength estimate of the stress–strength model for the EED based on the censored sample, convergence begins at 2000 iterations or less.

7. Summary and Conclusions

In this study, the statistical inference of R = P ( X < Y < Z ) for a component which has a strength that is independent on opposite lower and upper bound stresses, when the stresses and strength follow EED is discussed. We assume that both stresses and strength random variables are independent, having EED with a common scale parameter. Using the generalized progressive hybrid censoring design, various point and interval estimators for the reliability model R are obtained in view of the ML and Bayesian approaches. Bayesian estimators are achieved by the MCMC method, as well as the MH algorithm, based on SELF and Linex loss functions, which are all conducted in light of informative priors. The CIs are derived using asymptotic distribution theory as well as Bayes credible intervals being constructed. The Monte Carlo simulation is conducted for the comparison of the effectiveness of the proposed estimates through some measures, such as average values, MSE, and CIs lengths. The outcomes of the study showed that the Bayes estimates produce a lower MSE for four parameters and stress–strength reliability based on a generalized progressive hybrid censoring scheme. Finally, a progressively censored real dataset is presented for illustration purposes. More in-depth studies, such as the application of the Bayesian estimation for the exponentiated exponential distribution based on various schemes, should be discussed in the future.

Author Contributions

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

Funding

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R299), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Datasets are available in the application section.

Acknowledgments

The authors thank the support from the Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R299), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. We thank the referees for their valuable suggestions which improved the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bhattacharyya, G.K.; Johnson, R.A. Estimation of reliability in a multi-component stress-strength model. J. Am. Stat. Assoc. 1974, 69, 966–970. [Google Scholar] [CrossRef]
  2. Ahmad, K.E.; Jaheen, Z.F.; Yousef, M.M. Inference on Pareto distribution as stress- strength model based on generalized order statistics. J. Appl. Stat. Sci. 2008, 17, 247–257. [Google Scholar]
  3. Saraçoglu, B.; Kaya, M.F.; Abd-Elfattah, A.M. Comparison of estimators for stress-strength reliability in the Gompertz case. Hacet. J. Mathmatics Stat. 2009, 38, 339–349. [Google Scholar]
  4. Hassan, A.S.; Nagy, H.F.; Muhammed, H.Z.; Saad, M.S. Estimation of multicomponent stress-strength reliability following Weibull distribution based on upper record values. J. Taibah Univ. Sci. 2020, 14, 244–253. [Google Scholar] [CrossRef]
  5. Kotb, M.S.; Raqab, M.Z. Statistical inference for modified Weibull distribution based on progressively type-II censored data. Math. Comput. Simul. 2019, 162, 233–248. [Google Scholar] [CrossRef]
  6. Jana, N.; Bera, S. Estimation of parameters of inverse Weibull distribution and application to multi-component stress-strength model. J. Appl. Stat. 2022, 49, 169–194. [Google Scholar] [CrossRef]
  7. Chandra, S.; Owen, D.B. On estimating the reliability of a component subjected to several different stresses (strengths). Nav. Res. Logist. Q. 1975, 22, 31–39. [Google Scholar]
  8. Kotz, S.; Lumelskii, Y.; Pensky, M. The Stress–Strength Model and Its Generalizations: Theory and Applications; World Scientific: New York, NY, USA, 2003. [Google Scholar]
  9. Singh, N. On the estimation of P(Y1<X<Y2). Commun. Stat. Theory Methods 1980, 9, 1551–1561. [Google Scholar]
  10. Dutta, K.; Sriwastav, G. An n-standby system with P(X<Y<Z). IAPQR Trans. 1986, 12, 95–97. [Google Scholar]
  11. Ahmad, H.H.; Almetwally, E.M.; Ramadan, D.A. A comparative inference on reliability estimation for a multi-component stress-strength model under power Lomax distribution with applications. AIMS Math. 2022, 7, 18050–18079. [Google Scholar] [CrossRef]
  12. Almetwally, E.M.; Alotaibi, R.; Mutairi, A.A.; Park, C.; Rezk, H. Optimal Plan of Multi-Stress–Strength Reliability Bayesian and Non-Bayesian Methods for the Alpha Power Exponential Model Using Progressive First Failure. Symmetry 2022, 14, 1306. [Google Scholar] [CrossRef]
  13. Ivshin, V.V. On the estimation of the probabilities of a double linear inequality in the case of uniform and two-parameter exponential distributions. J. Math. Sci. 1998, 88, 819–827. [Google Scholar] [CrossRef]
  14. Hassan, A.S.; Elsherpieny, E.A.; Shalaby, R.M. On the estimation of P(Y<X<Z) for Weibull distribution in the presence of k outliers. Int. J. Eng. Res. Appl. 2013, 3, 1728–1734. [Google Scholar]
  15. Metwally, A.S.M.; Hassan, A.S.; Almetwally, E.M.; Kibria, B.M.; Almongy, H.M. Reliability analysis of the New exponential inverted Topp–Leone distribution with applications. Entropy 2021, 23, 1662. [Google Scholar] [CrossRef]
  16. Abu El Azm, W.S.; Almetwally, E.M.; Alghamdi, A.S.; Aljohani, H.M.; Muse, A.H.; Abo-Kasem, O.E. Stress-strength reliability for exponentiated inverted Weibull distribution with application on breaking of Jute fiber and Carbon fibers. Comput. Intell. Neurosci. 2021, 2021, 4227346. [Google Scholar] [CrossRef]
  17. Sabry, M.A.; Almetwally, E.M.; Alamri, O.A.; Yusuf, M.; Almongy, H.M.; Eldeeb, A.S. Inference of fuzzy reliability model for inverse Rayleigh distribution. AIMS Math. 2021, 6, 9770–9785. [Google Scholar] [CrossRef]
  18. Wang, Z.; Xiping, W.; Guangming, P. Nonparametric statistical inference for P(X<Y<Z). Sankya A Indian J. Stat. 2013, 75, 118–138. [Google Scholar]
  19. Patowary, A.N.; Sriwastav, G.L.; Hazarika, J. Inference of R=P(X<Y<Z) for n-standby system: A Monte Carlo simulation approach. IOSR-J. Math. 2016, 12, 18–22. [Google Scholar]
  20. Yousif, S.M.; Karam, N.S.; Karam, G.S.; Abood, Z.M. Stress strength reliability estimation for R=P(T<X<Z) using exponentiated inverse Rayleigh distribution. In AIP Conference Proceedings; AIP Publishing LLC: Melville, NY, USA, 2020; Volume 2307, p. 020013-1–9. [Google Scholar]
  21. Hameed, B.A.; Salman, A.N.; Kalaf, B.A. On estimation of (Y1<X<Y2) in cased inverse Kumaraswamy distribution. Iraqi J. Sci. 2020, 61, 845–853. [Google Scholar]
  22. Attia, A.M.; Karam, N.S. Bayesian estimation of stress-strength R=P(X<Y<Z) Dagum distribution. J. Phys. Conf. Ser. 2021, 1963, 012041. [Google Scholar]
  23. Raheem, S.H.; Kalaf, B.A.; Salman, A.N. Comparison of some of estimation methods of stress-strength model:R=P(X<Y<Z). Baghdad Sci. J. 2021, 18, 1103–1111. [Google Scholar]
  24. Abd Elfattah, A.M.; Taha, M.A. On the estimation ofP(y<x<z) for inverse Rayleigh distribution in the presence of outliers. J. Stat. Appl. Probab. Lett. 2021, 8, 181–189. [Google Scholar]
  25. Yousef, M.M.; Almetwally, E.M. Multi stress-strength reliability based on progressive first failure for Kumaraswamy model: Bayesian and non-Bayesian estimation. Symmetry 2021, 13, 2120. [Google Scholar] [CrossRef]
  26. Meeker, W.Q.; Escobar, L.A. Statistical Methods for Reliability Data; Wiley: New York, NY, USA, 1998. [Google Scholar]
  27. Cohen, A.C. Progressively censored samples in life testing. Technometrics 1963, 5, 327–339. [Google Scholar] [CrossRef]
  28. Kundu, D.; Joarder, A. Analysis of type-II progressively hybrid censored data. Comput. Stat. Data Anal. 2006, 50, 2509–2528. [Google Scholar] [CrossRef]
  29. El-Sherpieny, E.S.A.; Almetwally, E.M.; Muhammed, H.Z. Progressive Type-II hybrid censored schemes based on maximum product spacing with application to Power Lomax distribution. Phys. A Stat. Mech. Its Appl. 2020, 553, 124251. [Google Scholar] [CrossRef]
  30. Cho, Y.; Sun, H.; Lee, K. An estimation of the entropy for a Rayleigh distribution based on doubly-generalized type-II hybrid censored samples. Entropy 2014, 16, 3655–3669. [Google Scholar] [CrossRef]
  31. Mohie El-Din, M.M.; Nagy, M.; Abu-Moussa, M.H. Estimation and prediction for Gompertz distribution under the generalized progressive hybrid censored data. Ann. Data Sci. 2019, 6, 673–705. [Google Scholar] [CrossRef]
  32. Tu, J.; Gui, W. Bayesian inference for the Kumaraswamy distribution under generalized progressive hybrid censoring. Entropy 2020, 22, 1032. [Google Scholar] [CrossRef]
  33. Nagy, M.; Sultan, K.S.; Abu-Moussa, M.H. Analysis of the generalized progressive hybrid censoring from Burr Type-XII lifetime model. AIMS Math. 2021, 6, 9675–9704. [Google Scholar] [CrossRef]
  34. Maswadah, M. Improved maximum likelihood estimation of the shape-scale family based on the generalized progressive hybrid censoring scheme. J. Appl. Stat. 2022, 49, 2825–2844. [Google Scholar] [CrossRef]
  35. Ahuja, J.C.; Nash, S.W. The generalized gompertz verhulst family of distributions. Sankhya Indian J. Stat. 1967, 29, 144–156. [Google Scholar]
  36. Gupta, R.D.; Kundu, D. Generalized exponential distributions. Aust. N. Z. J. Stat. 1999, 41, 173–188. [Google Scholar] [CrossRef]
  37. Raja, T.A.; Mir, A.H. On extension of some exponentiated distributions with application. J. Contemp. Math. Sci. 2011, 6, 393–400. [Google Scholar]
  38. Yeates, M.P.; Tolkamp, B.J.; Allcroft, D.J.; Kyriazakis, I. The use of mixed distribution models to determine bout criteria for analysis of animal behaviour. J. Theor. Biol. 2001, 213, 413–425. [Google Scholar] [CrossRef]
  39. Subburaj, R.; Gopal, G.; Kapur, P.K. A software reliability growth model for vital quality metrics. S. Afr. J. Ind. Eng. 2007, 18, 93–108. [Google Scholar] [CrossRef]
  40. Biondi, F.; Kozubowski, T.J.; Panorska, A.K.; Saito, L. A new stochastic model of episode peak and duration for eco-hydro-climatic applications. Ecol. Model. 2008, 211, 383–395. [Google Scholar] [CrossRef]
  41. Kannan, N.; Kundu, D.; Nair, P.; Tripathi, R.C. The generalized exponential cure rate model with covariates. J. Appl. Stat. 2010, 37, 1625–1636. [Google Scholar] [CrossRef]
  42. Kundu, D. Bayesian inference and life testing plan for the Weibull distribution in presence of progressive censoring. Technometrics 2008, 50, 144–154. [Google Scholar] [CrossRef]
  43. Dey, S.; Singh, S.; Tripathi, Y.M.; Asgharzadeh, A. Estimation and prediction for a progressively censored generalized inverted exponential distribution. Stat. Methodol. 2016, 32, 185–202. [Google Scholar] [CrossRef]
  44. Chen, M.H.; Shao, Q.M. Monte Carlo estimation of Bayesian credible and HPD intervals. J. Comput. Graph. Stat. 1999, 8, 69–92. [Google Scholar]
  45. Nelson, W.B. Accelerated Testing: Statistical Models, Test Plans, and Data Analysis; John Wiley & Sons: Hoboken, NJ, USA, 2009; Volume 344, pp. 1625–1636. [Google Scholar]
Figure 1. The reliability R of stress–strength model with different values of parameters.
Figure 1. The reliability R of stress–strength model with different values of parameters.
Axioms 11 00455 g001
Figure 2. Heatmaps of MSE for parameters and R based on different estimation methods: Case I and II.
Figure 2. Heatmaps of MSE for parameters and R based on different estimation methods: Case I and II.
Axioms 11 00455 g002
Figure 3. Heatmaps of MSE for parameters and R based on different estimation methods: Case III and IV.
Figure 3. Heatmaps of MSE for parameters and R based on different estimation methods: Case III and IV.
Axioms 11 00455 g003
Figure 4. Empirical CDF and P-P plots for the EED for dataset 1.
Figure 4. Empirical CDF and P-P plots for the EED for dataset 1.
Axioms 11 00455 g004
Figure 5. Empirical CDF and P-P plots for the EED for dataset 2.
Figure 5. Empirical CDF and P-P plots for the EED for dataset 2.
Axioms 11 00455 g005
Figure 6. Empirical CDF and P-P plots for the EED for dataset 3.
Figure 6. Empirical CDF and P-P plots for the EED for dataset 3.
Axioms 11 00455 g006
Figure 7. Trace and normal curve of posterior distribution for MCMC estimation of stress–strength for EED based on complete sample.
Figure 7. Trace and normal curve of posterior distribution for MCMC estimation of stress–strength for EED based on complete sample.
Axioms 11 00455 g007
Figure 8. Pairs plots of the MCMC samples for parameter estimates of stress–strength for EED based on complete sample.
Figure 8. Pairs plots of the MCMC samples for parameter estimates of stress–strength for EED based on complete sample.
Axioms 11 00455 g008
Figure 9. Trace and normal curve of posterior distribution for MCMC estimation of reliability stress–strength for EED based on complete sample.
Figure 9. Trace and normal curve of posterior distribution for MCMC estimation of reliability stress–strength for EED based on complete sample.
Axioms 11 00455 g009
Figure 10. Convergence plots of MCMC for parameter estimates of the EED for censored sample.
Figure 10. Convergence plots of MCMC for parameter estimates of the EED for censored sample.
Axioms 11 00455 g010
Figure 11. Pairs plot of the MCMC samples for parameter estimates of stress–strength for EED based on complete sample.
Figure 11. Pairs plot of the MCMC samples for parameter estimates of stress–strength for EED based on complete sample.
Axioms 11 00455 g011
Figure 12. Trace and normal curve of posterior distribution for MCMC estimation of reliability stress–strength for EED based on censored sample.
Figure 12. Trace and normal curve of posterior distribution for MCMC estimation of reliability stress–strength for EED based on censored sample.
Axioms 11 00455 g012
Table 1. Labels of MSE results for heatmaps plots.
Table 1. Labels of MSE results for heatmaps plots.
scheme1
n 1 , n 2 , n 3 20, 25, 1530, 40, 30
m 1 , m 2 , m 3 15, 18, 1217, 22, 1425, 30, 2227, 35, 26
k 1 , k 2 , k 3 12, 16, 1015, 20, 1220, 24, 2023, 30, 23
labelMSE1MSE2MSE3MSE4
scheme2
n 1 , n 2 , n 3 20, 25, 1530, 40, 30
m 1 , m 2 , m 3 15, 18, 1217, 22, 1425, 30, 2227, 35, 26
k 1 , k 2 , k 3 12, 16, 1015, 20, 1220, 24, 2023, 30, 23
labelMSE5MSE6MSE7MSE8
Table 2. AvE and MSE for estimation methods: Case I.
Table 2. AvE and MSE for estimation methods: Case I.
α 1 = 0.5 , α 2 = 4 , α 3 = 20 , λ = 0.5
T 1 = 14 , T 2 = 6 , T 3 = 18 ML SELF Linex c = −0.5 Linex c = 1.5 ML SELF c = −0.5 c = 1.5
Scheme n 1 , n 2 , n 3 m 1 , m 2 , m 3 k 1 , k 2 , k 3 AvE MSE AvE MSE AvE MSE AvE MSE L.CI CP L.CCI L.CCI L.CCI
120, 25, 1515, 18, 1212, 16, 10 α 1 0.53780.02040.50730.00980.50540.00970.50160.00950.540195.80%0.37810.36990.3701
α 2 4.13540.80593.99700.02343.99420.02343.98870.02353.482294.40%0.57160.57150.5701
α 3 20.05130.546019.99540.021619.99260.021819.98720.02212.892595.80%0.56760.56390.5553
λ 0.50580.00210.50750.00190.50690.00180.50570.00180.177694.40%0.15840.15830.1581
R0.71530.00090.72450.00040.72490.00030.72560.00030.108995.80%0.07280.07280.0723
17, 22, 1415, 20, 12 α 1 0.53110.01760.50160.00720.50040.00710.49810.00710.505995.20%0.31510.31540.3159
α 2 4.06150.75334.00260.01254.00110.01253.99800.01253.397295.80%0.42280.42210.4160
α 3 20.06480.224519.99610.012819.99450.012819.99120.01291.841995.20%0.41680.41730.4176
λ 0.50100.00170.50390.00140.50350.00140.50260.00140.162795.80%0.14650.14650.1460
R0.71740.00080.72550.00030.72580.00030.72620.00030.103895.20%0.06130.06120.0612
30, 40, 3025, 30, 2220, 24, 20 α 1 0.51510.01120.50530.00700.50370.00700.50050.00680.411395.40%0.32280.32070.3156
α 2 4.08740.43184.00560.01924.00300.01913.99770.01902.555795.60%0.54810.54830.5488
α 3 20.00910.033719.99210.022019.98940.022119.98410.02230.719395.40%0.57630.57510.5755
λ 0.50190.00100.50330.00090.50300.00090.50230.00090.124495.60%0.11340.11350.1138
R0.72080.00040.72470.00030.72510.00030.72570.00020.080795.40%0.06180.06090.0595
27, 35, 2623, 30, 23 α 1 0.52480.01030.51280.00670.51170.00670.50950.00650.385295.20%0.30920.30570.3035
α 2 4.05180.38753.99270.01303.99110.01303.98790.01312.434395.40%0.43930.44300.4463
α 3 20.02080.030519.99520.015019.99350.015019.99010.01520.688495.20%0.46500.46430.4668
λ 0.49990.00080.50180.00080.50150.00080.50090.00080.112795.40%0.10710.10700.1068
R0.71940.00040.72350.00020.72370.00020.72410.00020.076195.20%0.06040.06020.0593
220, 25, 1515, 18, 1212, 16, 10 α 1 0.52350.01550.50640.01050.50470.01040.50120.01020.479294.40%0.38280.38190.3814
α 2 4.11420.74234.00010.02153.99750.02153.99230.02153.351095.20%0.56420.56290.5569
α 3 20.04420.596920.00320.023820.00040.023819.99480.02383.026594.40%0.57830.57870.5797
λ 0.50180.00180.50290.00160.50240.00160.50130.00160.164895.20%0.15440.15450.1547
R0.71870.00070.72480.00040.72510.00040.72580.00040.099394.40%0.07360.07360.0734
17, 22, 1415, 20, 12 α 1 0.53020.01480.50000.00840.49880.00830.49640.00820.451994.80%0.34700.34590.3470
α 2 4.13040.59434.00280.01424.00110.01423.99760.01422.981594.40%0.47770.47730.4754
α 3 20.05810.444320.00360.012620.00190.012619.99840.01262.605494.80%0.44900.44770.4501
λ 0.50260.00150.50540.00150.50490.00150.50400.00150.152594.40%0.14140.14110.1411
R0.71740.00060.72590.00030.72620.00030.72660.00030.091094.80%0.06610.06590.0664
30, 40, 3025, 30, 2220, 24, 20 α 1 0.52350.01550.50640.01050.50470.01040.50120.01020.479294.40%0.38280.38190.3814
α 2 4.11420.74234.00010.02153.99750.02153.99230.02153.351095.20%0.56420.56290.5569
α 3 20.04420.596920.00320.023820.00040.023819.99480.02383.026594.40%0.57830.57870.5797
λ 0.50180.00180.50290.00160.50240.00160.50130.00160.164895.20%0.15440.15450.1547
R0.71870.00070.72480.00040.72510.00040.72580.00040.099394.40%0.07360.07360.0734
27, 35, 2623, 30, 23 α 1 0.53020.01480.50000.00840.49880.00830.49640.00820.451994.80%0.34700.34590.3470
α 2 4.13040.59434.00280.01424.00110.01423.99760.01422.981594.40%0.47770.47730.4754
α 3 20.05810.444320.00360.012620.00190.012619.99840.01262.605494.80%0.44900.44770.4501
λ 0.50260.00150.50540.00150.50490.00150.50400.00150.152594.40%0.14140.14110.1411
R0.71740.00060.72590.00030.72620.00030.72660.00030.091094.80%0.06610.06590.0664
Table 3. AvE and MSE for estimation methods: Case II.
Table 3. AvE and MSE for estimation methods: Case II.
α 1 = 0.8 , α 2 = 5 , α 3 = 12 , λ = 1.5
T 1 = 1 . 6 , T 2 = 2 , T 3 = 2 . 5 MLSELFLinex c = −0.5Linex c = 1.5MLSELF c = −0.5c = 1.5
Scheme n 1 , n 2 , n 3 m 1 , m 2 , m 3 k 1 , k 2 , k 3 AvEMSEAvEMSEAvEMSEAvEMSEL.CICPL.CCIL.CCIL.CCI
120, 25, 1515, 18, 1212, 16, 10 α 1 0.85590.05690.80050.01200.79840.01190.79430.01190.909895.60%0.40970.41190.4129
α 2 5.49732.80464.99910.01844.99640.01844.99110.01846.275095.20%0.52090.52080.5209
α 3 12.71088.275512.00990.025712.00690.025512.00090.025312.265295.60%0.60730.60750.6097
λ 1.53200.03221.50220.00881.50040.00871.49690.00870.693195.20%0.35070.34790.3437
R0.56970.00320.58150.00020.58180.00020.58230.00020.217495.60%0.05550.05550.0561
17, 22, 1415, 20, 12 α 1 0.85110.04820.79760.00980.79630.00980.79360.00980.837995.20%0.37940.37990.3837
α 2 5.57101.89044.99930.01414.99760.01414.99410.01426.283596.00%0.43700.43890.4459
α 3 12.94785.206112.00250.014412.00080.014411.99740.014410.310095.20%0.45410.45430.4611
λ 1.54170.03161.49970.00761.49850.00761.49610.00760.672996.00%0.32860.32950.3314
R0.57380.00240.58180.00020.58200.00020.58230.00020.191495.20%0.05200.05210.0523
30, 40, 3025, 30, 2220, 24, 20 α 1 0.82290.03080.81070.01260.80870.01250.80470.01220.682495.00%0.42150.42320.4261
α 2 5.23411.24924.99580.02004.99320.02004.98800.02024.288495.00%0.54200.54000.5509
α 3 12.22603.706311.99740.021211.99460.021311.98910.02147.502095.00%0.56330.56570.5616
λ 1.50960.01751.50540.00691.50400.00691.50100.00680.518195.00%0.31880.31820.3167
R0.57460.00160.58000.00020.58030.00020.58080.00020.152795.00%0.05600.05540.0552
27, 35, 2623, 30, 23 α 1 0.83500.02980.80130.00890.80000.00880.79750.00880.663495.40%0.35350.35160.3475
α 2 5.32220.94195.00560.01225.00390.01225.00060.01224.150095.40%0.42980.43210.4286
α 3 12.24971.781011.99890.014311.99730.014311.99400.01446.566895.40%0.47540.47280.4689
λ 1.51760.01681.50400.00591.50290.00591.50070.00590.502595.40%0.29070.29180.2929
R0.57230.00130.58110.00020.58120.00020.58160.00020.139495.40%0.04820.04780.0480
220, 25, 1515, 18, 1212, 16, 10 α 1 0.87420.05140.80410.01330.80210.01320.79820.01310.840295.60%0.44720.44780.4468
α 2 5.44322.39054.99160.01994.98910.01994.98400.02015.812395.40%0.55180.55630.5516
α 3 12.93129.111412.00840.024412.00550.024311.99960.024311.266795.60%0.60500.60270.6014
λ 1.54690.03231.49740.00961.49570.00961.49220.00960.680495.40%0.39290.39240.3888
R0.57300.00230.58120.00030.58140.00030.58200.00020.186295.60%0.06390.06370.0636
17, 22, 1415, 20, 12 α 1 0.85600.05010.79860.01020.79720.01010.79450.01010.835094.80%0.38810.38870.3886
α 2 5.38972.17904.99940.01294.99780.01294.99470.01295.586894.20%0.43660.43560.4377
α 3 12.66665.164212.00190.013612.00030.013511.99700.01359.587194.80%0.43290.43480.4344
λ 1.52700.02701.50030.00641.49910.00641.49680.00640.635794.20%0.30100.30240.3051
R0.57180.00220.58160.00020.58180.00020.58220.00020.179094.80%0.05330.05280.0525
30, 40, 3025, 30, 2220, 24, 20 α 1 0.84690.03090.80560.01110.80360.01100.79970.01090.665095.40%0.39440.39350.3911
α 2 5.26951.12454.98840.02004.98560.02014.98020.02044.024494.80%0.55690.55730.5509
α 3 12.28753.879511.99320.021311.99030.021411.98450.02157.645995.40%0.57720.57850.5787
λ 1.51970.01501.50780.00691.50630.00681.50340.00680.474894.80%0.30960.30920.3105
R0.57220.00140.58070.00020.58100.00020.58150.00020.141395.40%0.05350.05360.0540
27, 35, 2623, 30, 23 α 1 0.83820.03020.80750.00830.80620.00820.80370.00810.658194.60%0.35040.34630.3435
α 2 5.25281.10205.00780.01255.00610.01245.00270.01243.927295.60%0.44390.44240.4393
α 3 12.16361.122311.99980.014011.99810.014011.99460.01406.756594.60%0.45030.44440.4456
λ 1.51690.01461.49890.00521.49790.00521.49580.00520.469395.60%0.28510.28610.2842
R0.57610.00130.58020.00020.58040.00010.58070.00010.137794.60%0.04770.04740.0469
Table 4. AvE and MSE for estimation methods: Case III.
Table 4. AvE and MSE for estimation methods: Case III.
α 1 = 0.3 , α 2 = 2 , α 3 = 15 , λ = 0.5
T 1 = 3 . 5 , T 2 = 3 , T 3 = 6 ML SELF Linex c = −0.5 Linex c = 1.5 ML SELF c = −0.5 c = 1.5
Scheme n 1 , n 2 , n 3 m 1 , m 2 , m 3 k 1 , k 2 , k 3 AvE MSE AvE MSE AvE MSE AvE MSE L.CI CP L.CCI L.CCI L.CCI
120, 25, 1515, 18, 1212, 16, 10 α 1 0.31660.00710.31480.00680.31340.00670.31060.00640.323295.40%0.29740.29500.2912
α 2 2.11130.27751.99280.01931.99020.01941.98510.01972.020496.20%0.53800.53950.5349
α 3 15.18040.929515.00860.023715.00570.023514.99990.02333.716295.40%0.58480.58590.5892
λ 0.50890.00350.50810.00280.50720.00270.50550.00270.229896.20%0.19080.19040.1895
R0.74610.00100.74910.00090.74950.00090.75040.00090.120295.40%0.10970.10960.1087
17, 22, 1415, 20, 12 α 1 0.31190.00560.30610.00440.30510.00430.30320.00430.288796.20%0.24710.24600.2461
α 2 2.14960.28022.00040.01251.99890.01251.99580.01251.992595.20%0.45000.44940.4458
α 3 15.06550.693214.99440.013314.99270.013314.98930.01353.257096.20%0.43730.44000.4415
λ 0.50550.00290.50210.00250.50150.00250.50020.00250.208795.20%0.19770.19760.1959
R0.74670.00080.75220.00060.75250.00060.75310.00060.106796.20%0.09070.09010.0897
30, 40, 3025, 30, 2220, 24, 20 α 1 0.31180.00400.31390.00430.31280.00420.31080.00410.243595.20%0.24010.23880.2377
α 2 2.08370.18502.00010.02051.99750.02051.99230.02051.655695.40%0.55070.55220.5557
α 3 15.11721.152015.00280.021714.99990.021714.99410.02184.186595.20%0.54530.54590.5549
λ 0.50200.00180.50220.00140.50160.00140.50060.00140.168695.40%0.14890.14840.1467
R0.74810.00060.74910.00060.74940.00060.75010.00050.092495.20%0.08860.08840.0886
27, 35, 2623, 30, 23 α 1 0.30810.00390.30390.00360.30320.00350.30160.00350.243294.40%0.22440.22370.2223
α 2 2.05000.15202.00080.01221.99920.01221.99590.01231.517195.20%0.42660.42880.4296
α 3 15.06180.904715.00490.012815.00330.012715.00010.01273.724394.40%0.41930.41920.4187
λ 0.49870.00150.50140.00140.50100.00140.50020.00140.151895.20%0.15030.15050.1510
R0.74930.00060.75290.00050.75310.00050.75370.00050.093194.40%0.08240.08250.0821
220, 25, 1515, 18, 1212, 16, 10 α 1 0.31890.00710.31210.00470.31080.00460.30840.00450.298395.80%0.24140.23990.2371
α 2 2.08280.24771.98380.02221.98090.02241.97520.02271.925895.40%0.56880.56590.5681
α 3 15.10290.554114.99620.020314.99350.020314.98800.02042.251095.80%0.56430.56090.5610
λ 0.50420.00290.50560.00240.50480.00230.50330.00230.210195.40%0.18470.18420.1828
R0.74410.00090.74960.00060.75000.00060.75080.00060.110595.80%0.09240.09250.0925
17, 22, 1415, 20, 12 α 1 0.31800.00610.30180.00460.30090.00450.29900.00440.232294.60%0.24710.24690.2440
α 2 2.10920.22762.00180.01362.00010.01361.99680.01361.822294.60%0.47060.46930.4651
α 3 15.02120.329514.99730.013714.99560.013714.99230.01372.089394.60%0.44690.44770.4434
λ 0.50300.00250.50360.00230.50290.00240.50160.00240.197494.60%0.19270.19210.1904
R0.74560.00090.75380.00060.75410.00060.75470.00060.101794.60%0.09190.09210.0917
30, 40, 3025, 30, 2220, 24, 20 α 1 0.31070.00370.30960.00350.30870.00340.30670.00340.236695.80%0.22270.22190.2201
α 2 2.03200.13281.98920.02211.98660.02221.98150.02261.424595.20%0.56940.56540.5668
α 3 15.10210.702114.98790.020414.98530.020514.98010.02073.263495.80%0.52990.53130.5318
λ 0.50050.00140.50270.00120.50220.00120.50130.00120.147495.20%0.13360.13330.1332
R0.74890.00060.75020.00050.75050.00050.75120.00050.090195.80%0.08200.08170.0806
27, 35, 2623, 30, 23 α 1 0.31170.00380.30860.00350.30780.00340.30620.00330.236895.60%0.21510.21360.2131
α 2 2.08030.12502.00170.01202.00020.01201.99720.01201.408696.00%0.40760.40620.4030
α 3 15.04940.404214.99780.013114.99610.013114.99280.01322.487295.60%0.43940.44120.4450
λ 0.50650.00140.50550.00130.50510.00120.50430.00130.143796.00%0.14060.14060.1401
R0.74820.00050.75100.00050.75130.00050.75190.00050.088295.60%0.08050.07920.0776
Table 5. AvE and MSE for estimation methods: Case IV.
Table 5. AvE and MSE for estimation methods: Case IV.
α 1 = 0.7 , α 2 = 3 , α 3 = 8 , λ = 2
T 1 = 1 , T 2 = 1 . 5 , T 3 = 1 . 9 ML SELF Linex c = −0.5 Linex c = 1.5 ML SELF c = −0.5 c = 1.5
Scheme n 1 , n 2 , n 3 m 1 , m 2 , m 3 k 1 , k 2 , k 3 AvE MSELF AvE MSE AvE MSE AvE MSE L.CI CP L.CCI L.CCI L.CCI
120, 25, 1515, 18, 1212, 16, 10 α 1 0.76350.04170.70330.01230.70140.01230.69750.01220.761195.00%0.42890.42890.4266
α 2 3.38701.34373.00050.02102.99780.02102.99250.02114.287695.80%0.55080.54800.5441
α 3 9.36169.30147.99960.02277.99670.02277.99100.022810.276295.00%0.58670.59110.5934
λ 2.08070.10862.00150.01371.99940.01361.99530.01351.253995.80%0.44730.44720.4465
R0.55230.00330.55420.00050.55460.00050.55530.00050.225195.00%0.08550.08510.0846
17, 22, 1415, 20, 12 α 1 0.76470.04040.70220.00920.70080.00920.69800.00910.718596.00%0.37310.37250.3693
α 2 3.43651.03652.99490.01292.99330.01292.99010.01304.252094.80%0.43270.43210.4262
α 3 9.30347.59587.99830.01457.99660.01457.99330.01459.344596.00%0.46640.46580.4637
λ 2.10140.10451.99550.00911.99420.00911.99170.00911.204394.80%0.38270.38300.3834
R0.55070.00310.55440.00040.55460.00040.55510.00040.213196.00%0.07350.07330.0726
30, 40, 3025, 30, 2220, 24, 20 α 1 0.73360.02330.70890.01240.70690.01230.70280.01200.584695.40%0.40940.41290.4123
α 2 3.24290.67142.99390.02132.99130.02142.98610.02153.070794.60%0.55880.56270.5669
α 3 8.88946.92827.99360.02397.99080.02397.98520.02409.720995.40%0.61590.61380.6050
λ 2.07190.07052.00180.01071.99990.01071.99610.01061.002794.60%0.41260.41230.4110
R0.55500.00210.55310.00050.55340.00050.55420.00050.179495.40%0.08210.08210.0823
27, 35, 2623, 30, 23 α 1 0.72850.02280.70300.00850.70170.00850.69910.00840.582495.80%0.34470.34380.3420
α 2 3.16710.48013.00040.01472.99870.01472.99530.01472.638595.40%0.48200.47780.4715
α 3 8.82466.84638.00270.01388.00090.01377.99750.01379.743995.80%0.46010.46160.4614
λ 2.05960.05422.00900.00772.00770.00772.00510.00760.883095.40%0.33850.34010.3398
R0.55610.00200.55420.00030.55450.00030.55490.00030.176295.80%0.06930.06940.0691
220, 25, 1515, 18, 1212, 16, 10 α 1 0.74240.04400.69940.01400.69720.01390.69280.01380.806295.00%0.44810.44680.4455
α 2 3.31671.13182.99270.02312.98980.02322.98400.02343.985393.80%0.56160.56550.5752
α 3 9.04368.03468.00550.02408.00250.02407.99650.023912.374795.00%0.59220.59760.6043
λ 2.07100.10912.00110.01151.99910.01141.99520.01141.265693.80%0.41000.41170.4110
R0.54970.00390.55540.00060.55580.00060.55660.00060.245795.00%0.09020.09050.0910
17, 22, 1415, 20, 12 α 1 0.75700.04050.68840.00930.68700.00930.68420.00940.785497.00%0.37880.37660.3752
α 2 3.39831.11823.00620.01463.00440.01463.00070.01453.843995.00%0.45590.45550.4568
α 3 9.06105.52698.01050.01518.00880.01508.00540.014912.975097.00%0.47610.47330.4721
λ 2.10420.09832.00070.00941.99930.00941.99660.00941.160595.00%0.38200.38080.3814
R0.55880.00390.55720.00040.55750.00040.55800.00040.245297.00%0.07590.07570.0755
30, 40, 3025, 30, 2220, 24, 20 α 1 0.74160.02430.70760.01100.70570.01090.70190.01070.590095.40%0.40900.40820.4065
α 2 3.22920.71062.99570.02022.99310.02022.98780.02033.183294.80%0.54590.54660.5530
α 3 8.79516.38678.00280.02517.99970.02517.99350.025110.198195.40%0.59560.60370.6049
λ 2.04790.06931.99910.01091.99720.01091.99360.01091.016094.80%0.40720.40570.4051
R0.55100.00240.55340.00040.55370.00040.55440.00040.191995.40%0.07870.07920.0793
27, 35, 2623, 30, 23 α 1 0.74430.02280.69700.00870.69560.00870.69280.00870.563494.80%0.36350.36100.3605
α 2 3.21760.64042.99130.01402.98950.01412.98600.01433.021996.20%0.45730.46040.4599
α 3 7.93324.15157.99810.01447.99640.01447.99310.014410.587894.80%0.45550.45520.4541
λ 2.06550.06742.00350.00792.00230.00791.99980.00790.986296.20%0.35120.35040.3508
R0.55330.00230.55540.00030.55570.00030.55620.00040.189594.80%0.07060.07050.0707
Table 6. ML and Bayesian estimates for reliability in stress–strength model based on complete sample.
Table 6. ML and Bayesian estimates for reliability in stress–strength model based on complete sample.
EstimatesSEKSp-ValueAICCAICBICHQIC
X α 1 1.30440.57290.23960.599480.580482.580480.974979.7292
λ 0.04430.0173
Y α 2 4.76662.71560.19270.832260.949162.949161.343560.0979
λ 0.17690.0537
Z α 3 1.81490.93480.19010.8864103.7330106.1330103.8919102.6614
λ 0.00700.0028
Table 7. ML and Bayesian estimates for stress–strength model based on complete sample.
Table 7. ML and Bayesian estimates for stress–strength model based on complete sample.
MLBayesian
EstimatesSELFLowerUpperEstimatesSELFLowerUpper
α 1 0.61700.22160.18271.05130.65060.16470.28961.0329
α 2 0.50310.17810.15400.85230.54160.17160.21890.8249
α 3 3.68451.74660.26127.10774.01671.65911.09437.0266
λ 0.01360.00350.00670.02040.01450.00330.00820.0202
R0.34450.3503
Table 8. ML and Bayesian estimates for reliability of stress–strength model.
Table 8. ML and Bayesian estimates for reliability of stress–strength model.
MLBayesian
T 1 , T 2 , T 3 EstimatesSELFLowerUpperEstimatesSELFLowerUpper
25,60,550 α 1 0.55770.20120.16330.95200.58170.18330.25930.9484
α 2 0.50390.17970.15170.85620.56180.17390.23270.8495
α 3 2.72591.31520.14815.30373.42251.55950.93345.6519
λ 0.01050.00310.00440.01660.01180.00300.00620.0168
R0.34160.3683
20,50,500 α 1 0.57260.20670.16740.97780.60690.20070.24820.9376
α 2 0.51160.18260.15370.86950.57400.17200.21010.8999
α 3 2.90211.40570.14685.65733.60671.36331.04575.1271
λ 0.01110.00320.00470.01740.01240.00230.00670.0169
R0.34350.3662
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yousef, M.M.; Hassan, A.S.; Alshanbari, H.M.; El-Bagoury, A.-A.H.; Almetwally, E.M. Bayesian and Non-Bayesian Analysis of Exponentiated Exponential Stress–Strength Model Based on Generalized Progressive Hybrid Censoring Process. Axioms 2022, 11, 455. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms11090455

AMA Style

Yousef MM, Hassan AS, Alshanbari HM, El-Bagoury A-AH, Almetwally EM. Bayesian and Non-Bayesian Analysis of Exponentiated Exponential Stress–Strength Model Based on Generalized Progressive Hybrid Censoring Process. Axioms. 2022; 11(9):455. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms11090455

Chicago/Turabian Style

Yousef, Manal M., Amal S. Hassan, Huda M. Alshanbari, Abdal-Aziz H. El-Bagoury, and Ehab M. Almetwally. 2022. "Bayesian and Non-Bayesian Analysis of Exponentiated Exponential Stress–Strength Model Based on Generalized Progressive Hybrid Censoring Process" Axioms 11, no. 9: 455. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms11090455

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