Next Article in Journal
Inerter-Based Eigenvector Orientation Approach for Passive Control of Supersonic Panel Flutter
Previous Article in Journal
Identification of Vital Genes for NSCLC Integrating Mutual Information and Synergy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Dependent Competing Risks Model with Baseline Proportional Hazards Models under Minimum Ranked Set Sampling

1
School of Mathematics, Yunnan Normal University, Kunming 650500, China
2
Department of Statistics, Tamkang University, Tamsui District, New Taipei City 251037, Taiwan
3
Department of Mathematics, Indian Institute of Technology Patna, Bihta 801103, India
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(6), 1461; https://doi.org/10.3390/math11061461
Submission received: 14 January 2023 / Revised: 12 March 2023 / Accepted: 13 March 2023 / Published: 17 March 2023
(This article belongs to the Section Computational and Applied Mathematics)

Abstract

:
The ranked set sampling (RSS) is an efficient and flexible sampling method. Based on a modified RSS named minimum ranked set sampling samples (MinRSSU), inference of a dependent competing risks model is proposed in this paper. Then, Marshall–Olkin bivariate distribution model is used to describe the dependence of competing risks. When the competing risks data follow the proportional hazard rate distribution, a dependent competing risks model based on MinRSSU sampling is constructed. In addition, the model parameters and reliability indices were estimated by the classical and Bayesian method. Maximum likelihood estimators and corresponding asymptotic confidence intervals are constructed by using asymptotic theory. In addition, the Bayesian estimator and highest posterior density credible intervals are established under the general prior. Furthermore, according to E-Bayesian theory, the point and interval estimators of model parameters and reliability indices are obtained by a sampling algorithm. Finally, extensive simulation studies and a real-life example are presented for illustrations.

1. Introduction

In order to enhance the efficiency of simple random sampling (SRS) scheme, the ranked set sampling (RSS) method was proposed by McIntyre [1] as an alternative scenario for making inferences about the target population using less sample resource. Since the RSS scheme could improve the inference accuracy by using informal ranking information, the RSS could be viewed as a more suitable sampling scheme to improve testing efficient and inferential accuracy when experiment is time and cost consuming. The traditional RSS scenario proposed by McIntyre [1] can be implemented as follows. Suppose that n 2 identical units are randomly selected from a population and are divided into n groups with n units in each testing group. The n units are conducted in each testing group simultaneously, the i-th failure times in group j is denoted as X ( i ) j , i , j = 1 , 2 , , n . Then, the i-th smallest ranked data X ( i ) i is recorded in the i-th group. Under this scheme, the sketch of the RSS can be described as follows:   
group testing data RSS data 1 : X ( 1 ) 1 X ( 2 ) 1 X ( n ) 1 X ( 1 ) 1 2 : X ( 1 ) 2 X ( 2 ) 2 X ( n ) 2 X ( 2 ) 2 n : X ( 1 ) n X ( 2 ) n X ( n ) n X ( n ) n
Therefore, observation { X ( 1 ) 1 , X ( 2 ) 2 , , X ( n ) n } is collected as the RSS sample. Inference with using RSS scheme has been investigated by many authors under various perspectives. Takahasi and Wakimoto [2] and Takahasi [3] have developed the mathematical background of the method (and its efficiency) based on RSS samples. Dell and Clutter [4] showed that mean estimation based on RSS samples gives better results than using SRS samples in terms of the mean square error and cost efficiency. For some recent contributions, one could also refer to the works of Basikhasteh et al. [5], Esemen and Gürler [6], Samawi et al. [7], and Zamanzade et al. [8] among others.
In practice, the ranking steps of traditional RSS are tedious, and ranking errors are easy to occur in the process. Moreover, n 2 units are needed to obtain the associated ranking failure times in this complete sampling process. From the perspective of time and cost, this sampling scheme is still less efficient for traditional RSS. Therefore, Muttlak [9] proposed the median RSS, and Samawi et al. [10] proposed the extreme RSS to reduce the influence of ranking errors on RSS. Further, many other modified RSS schemes have been proposed, such as Double-RSS (Al-Saleh and Al-Kadiri [11]), Multi-stage RSS (Al-Saleh and Al-Omari [12]), and quartile RSS (Muttlak [13]). In particular, Al-Odat and Al-Saleh [14] proposed an effective modified RSS scheme called minimum ranked set sampling with unequal samples (MinRSSU), which has attracted much attention in recent years. The MinRSSU sample is conducted with only n ( n + 1 ) / 2 units in its testing. Under MinRSSU scenario, n groups are conducted where i-th ( i = 1 , , n ) group contains i units in experiment. The smallest observation in each group is selected in the ranking process, and the MinRSSU data can be described as follows:
group testing data MinRSSU data 1 : X ( 1 ) 1 X ( 1 ) 1 2 : X ( 1 ) 2 X ( 2 ) 2 X ( 1 ) 2 n : X ( 1 ) n X ( 2 ) n X ( n ) n X ( 1 ) n
From the sketch of the MinRSSU sampling, it is noted that the size of MinRSSU sampling is smaller than that of the traditional RSS, and such a scenario also reduces the ranking error in consequence. Similarly, if the maximum lifetimes are selected from each testing group, one could obtain the maximum ranked set sampling with unequal samples (MaxRSSU) (Biradar and Santosha [15]). In the literature, MinRSSU and MaxRSSU samples sometimes are called together as the moving extreme ranking set sampling (MERSS). As another efficient sampling scenario, both MinRSSU and MaxRSSU have been extensively discussed by many authors. For instance, Yao et al. [16] verified that MERSS can provide more information than SRS by comparing Fisher information matrix. Qiu and Eftekharian [17] considered the information measurements of MinRSSU and MaxRSSU from the perspective of entropy. Rahmani and Razmkhah [18] also showed that the parameter estimation has higher accuracy in MinRSSU and MaxRSSU scenarios, even though there are ordering errors in the sampling process.
In reliability and survival analysis, the failure of an individual unit is often yielded by multiple failure causes that compete with each other throughout the life cycle of the unit. In the literature, this problem is referred to as the competitive risk model. In general, competing risk data includes the failure times and corresponding indicators showing the cause of failure. In theoretical and practical studies, inference of competing risks model has been analyzed under the assumption that there are independent failure causes. In this regard, various inferential works have been done under different perspectives. See, for example, the works of Ren and Gui [19], Wang et al. [20], Sarhan et al. [21], Mahto et al. [22], among others. However, due to the complexity of internal structure and external environment of units, such independence assumption sometimes inappropriate. For example, there are two main failure causes namely cancer recurrence and death in colon cancer, and these two causes are dependent and competing with each other in clinical analysis. In reliability application for a two-engine aeroplane, the breakdown of one aeroplane engine could cause additional pressure on the remaining functional engine and increase the risk of engine failure. In this regard, inference for competing risks model with dependent failure causes seems more reasonable. In particular, competing risks models involving two dependent causes of failure are frequently studied in theoretical and practical analysis, and different approaches are proposed for modelling dependent characteristic including multivariate lifetime distribution, copula model, and frailty model among others (e.g., Wang et al. [23], Gorfine and Hsu [24], Lo and Wilke [25], Fan et al. [26]). Among different approaches, Marshall and Olkin [27] proposed a popular bivariate distribution based on exponential and Weibull baseline model in reliability studies, which provided successful applications for competing risk analysis with dependent causes. Since then, the Marshall–Olkin lifetime model has been used extensively to analyze two dependent failure causes because of its simple structure and attractive statistical properties. For example, Guan et al. [28] considered the objective Bayesian analysis of dependent competing risks model using MOBE distribution. Feizjavadian and Hashemi [29] used Marshall–Olkin bivariate Weibull (MOBW) distribution to study diabetic retinopathy data. Paul et al. [30] considers a Bayesian analysis of a dependent competing risk model using Marshall–Olkin bivariate Pareto (MOBP) distribution. See also Kundu and Gupta [31], Mulinacci [32], Shen and Xu [33] for more discussions.
Let U be a random variable from a family of continuous distribution with a resilience parameter θ > 0 , the associated probability density function (PDF) and the cumulative distribution function (CDF) of U are, respectively, given by
f ( u ; θ ) = θ g ( u ) [ G ¯ ( u ) ] θ 1 , F ( u ; θ ) = 1 [ G ¯ ( u ) ] θ , u > 0 ,
where G ¯ ( u ) = 1 G ( u ) is an arbitrary continuous function with G ( 0 ) = 0 and G ( ) = 1 , and g ( u ) = d G ( u ) / d u . In addition, the survival function (SF) can be expressed as
S ( u ; θ ) = [ G ¯ ( u ) ] θ , u > 0 .
The model (1) is a subclass of a regular one parameter exponential family of distributions and refers as proportional hazard rate (PHR) model (e.g., Lawless [34]). By different choices of the baseline function G ( u ) , many well-known distributions are the specific cases of PHR model like exponential, Rayleigh, Pareto, among others (e.g., Contreras-Reyes et al. [35]). Inference for the PHR model has been studied by several authors. For example, some recent contributions of Chaturvedi et al. [36], Bai et al. [37], Kizilaslan [38], Popović et al. [39], among others. In this paper, it is assumed that the marginal distribution of the failure causes follows the PHR distribution, and then an Marshall–Olkin type bivariate PHR distribution (MOBPHR) is proposed for modelling the dependent competing risks.
Motivated by the former, this paper proposes inference for dependent competing risks model when the failure time is MOBPHR distributed under the MinRSSU scenario. The model parameters and reliability indices were estimated by the classical and Bayesian method. In addition, the expected-Bayesian (E-Bayesian) approach is also proposed for comparison. To the best of our knowledge, estimation for the dependent MinRSSU competing risks data has been not discussed in the literature.
The rest of the article is organized as follows. In Section 2, the MOBPHR model and the MinRSSU competing risks data are described, respectively. Section 3 presents the maximum likelihood estimators (MLEs) of the model parameters, and the associated approximate confidence intervals (ACIs) are also established in consequence. Bayesian and E-Bayesian estimations are presented in Section 4. Extensive simulation studies and a real-life example are conducted in Section 5 for illustration. Finally, some brief concluding remarks are given in Section 6.

2. Model and Data Description

2.1. Marshall–Olkin Bivariate Proportional Hazard Rate Distribution

Let U1, U3 and U3 be independent random variables following the PHR distributions (1) with parameters θ 1 ,   θ 2 and θ 3 , respectively. Define T1 = min(U1,U3) and T2 = min(U2,U3), then random vector ( T 1 , T 2 ) follow the MOBPHR distribution with parameters ( θ 1 , θ 2 , θ 3 ) , denoted by ( T 1 , T 2 ) MOBPHR ( θ 1 , θ 2 , θ 3 ) . For the sake of completeness, some basic properties of the MOBPHR ( θ 1 , θ 2 , θ 3 ) are provided in following theorems.
Theorem 1.
Suppose random vector ( T 1 , T 2 ) follows MOBPHR ( θ 1 , θ 2 , θ 3 ) . The joint SF of ( T 1 , T 2 ) is given by
S ( T 1 , T 2 ) ( t 1 , t 2 ) = S ( t 1 ; θ 1 ) S ( t 2 ; θ 2 + θ 3 ) , i f t 1 < t 2 ; S ( t 1 ; θ 1 + θ 3 ) S ( t 2 ; θ 2 ) , i f t 1 > t 2 ; S ( t ; θ 1 + θ 2 + θ 3 ) , i f t 1 = t 2 = t .
Proof. 
See Appendix A.    □
A direct calculation based on the relationship between SF and PDF yields the following results, with details omitted for simplicity.
Corollary 1.
Suppose random vector ( T 1 , T 2 ) follows MOBPHR ( θ 1 , θ 2 , θ 3 ) . The joint PDF of ( T 1 , T 2 ) can be expressed as
f ( T 1 , T 2 ) ( t 1 , t 2 ) = f ( t 1 ; θ 1 ) f ( t 2 ; θ 2 + θ 3 ) , t 1 < t 2 ; f ( t 1 ; θ 1 + θ 3 ) f ( t 2 ; θ 2 ) , t 1 > t 2 ; θ 3 θ 1 + θ 2 + θ 3 f ( t ; θ 1 + θ 2 + θ 3 ) , t 1 = t 2 = t ;
For the sake of simplicity, denote θ i j = θ i + θ j , i , j = 1 , 2 , 3 , i j , and θ 123 = θ 1 + θ 2 + θ 3 in rest of the paper. According to the above definition, one further has the following results.
Corollary 2.
Let ( T 1 , T 2 ) MOBPHR ( θ 1 , θ 2 , θ 3 ) , then
  • T = min { T 1 , T 2 } follows the PHR distribution with parameter θ 123 ;
  • The SF and the hazard rate function (HRF) of T at mission time t 0 are given by
    S ( t 0 ) = [ G ¯ ( t 0 ) ] θ 123 a n d H ( t 0 ) = θ 123 g ( t 0 ) [ G ¯ ( t 0 ) ] 1 .
From (3) with θ 3 = 0 , it is noted that variables T 1 and T 2 are statistical independent implying that parameter θ 3 can be regarded as the dependence structure between T 1 and T 2 . Furthermore, one can observe from (4) that the probability for events { T 1 < T 2 } , { T 1 > T 2 } and { T 1 = T 2 } are θ 1 θ 123 , θ 2 θ 123 and θ 3 θ 123 , respectively. In addition, it is also seen from Corollary 2 that the failure time of units can be described by variable T = min { T 1 , T 2 } following MOBPHR model.

2.2. Data Description and Notation

Without loss of generality, suppose n identical units are tested under the MinRSSU scenario, and their lifetimes are described by independent and identically distributed (i.i.d.) random variables T 1 , T 2 , , T n . For sake of simplicity, the following notations are used in this paper when the dependent MinRSSU competing risk data is available:
T i j :
failure time of ith unit under cause j = 1 , 2
T i :
observed lifetime of the i-th unit, i.e., T i = min ( T i 1 , T i 2 ) .
δ i :
indicator variable for the failure cause of ith unit with
δ i = 1 , unit fails due to cause 1 . 2 , unit fails due to cause 2 . 3 , unit fails due to both causes .
I i :
set of indices of case δ i = 1 , 2 and 3
| I i | :
cardinality of I i . It is assumed that | I i | = n i , i = 1 , 2 , 3
In this paper, the dependent lifetime variables of the competing risks are denoted by { ( T i 1 , T i 2 ) , i = 1 , 2 , , n } from population MOBPHR ( θ 1 , θ 2 , θ 3 ) . Therefore, the dependent competing risks sample is T i = min ( T i 1 , T i 2 ) , i = 1 , 2 , , n under MinRSSU scenario with failure cause δ i and the associated observations are given by
( t 1 , δ 1 ) , ( t 2 , δ 2 ) , , ( t n , δ n ) .
Theorem 2.
Suppose the observation (5) comes from MOBPHR model with parameters θ 1 , θ 2 and θ 3 . The likelihood function of parameters θ 1 , θ 2 and θ 3 can be expressed as
L ( θ 1 , θ 2 , θ 3 ) = i I 1 i θ 1 g ( t i ) [ G ¯ ( t i ) ] i θ 123 1 × i I 2 i θ 2 g ( t i ) [ G ¯ ( t i ) ] i θ 123 1 × i I 3 i θ 3 g ( t i ) [ G ¯ ( t i ) ] i θ 123 1 = n ! × θ 1 n 1 × θ 2 n 2 × θ 3 n 3 × i = 1 n g ( t i ) [ G ¯ ( t i ) ] i θ 123 1 ,
where n = n 123 = n 1 + n 2 + n 3 .
Proof. 
See Appendix B.    □

3. Classical Inference

In this section, the MLEs of model parameters and reliability indices are presented, and the correspondingly ACIs are also constructed based on asymptotic theory.

3.1. Maximum Likelihood Estimation

From (6), the log-likelihood function is given by
( θ 1 , θ 2 , θ 3 ) = ln n 123 ! + n 1 ln θ 1 + n 2 ln θ 2 + n 3 ln θ 3 + θ 123 ω + i = 1 n 123 ln g ( t i ) G ¯ ( t i )
where ω = i = 1 n 123 i ln G ¯ ( t i ) .
Theorem 3.
Suppose the observation (5) follow the MOBPHR distribution with parameters θ 1 , θ 2 , θ 3 . For n i > 0 , i = 1 , 2 , 3 , the MLE of θ i can be written as
θ ^ M i = n i ω , i = 1 , 2 , 3 .
Proof. 
See Appendix C.    □
Further, using invariance principle of maximum likelihood estimation, the MLEs of the SF and HRF can be constructed as
S ^ M ( t ) = [ G ¯ ( t 0 ) ] θ ^ M 123 and H ^ M ( t ) = θ ^ M 123 g ( t ) G ¯ ( t 0 ) ,
with θ ^ M 123 = θ ^ M 1 + θ ^ M 2 + θ ^ M 3 = n 123 ω .
Remark 1.
It is noted that the condition n i > 0 , i = 1 , 2 , 3 guarantees the existence of MLE θ ^ M i . If (without loss of generality) n i = 0 , sample (5) provides no information about θ i and the MLE of θ i does not exist in this manner.

3.2. Approximate Confidence Intervals

The expected Fisher information matrix and asymptotic theory in this subsection to obtain the ACIs of θ 1 , θ 2 , θ 3 . Let Θ = ( θ 1 , θ 2 , θ 3 ) , the expected Fisher information matrix I ( Θ ) can be constructed by
I ( Θ ) = [ I i j ] i , j = 1 , 2 , 3 = E 2 ( θ 1 , θ 2 , θ 3 ) θ i θ j i , j = 1 , 2 , 3
where I i i = n i θ i 2 for i = 1 , 2 , 3 and I i j = 0 if i j . Under mild regularity conditions, the asymptotic distribution of the MLE Θ ^ M = ( θ ^ M 1 , θ ^ M 2 , θ ^ M 3 ) is
Θ ^ M Θ N 0 , I 1 ( Θ ^ M ) ,
where I 1 ( Θ ^ M ) is the inverse of the expected Fisher information matrix given by
I 1 ( Θ ^ ) = V a r ( θ ^ M 1 ) C o v ( θ ^ M 1 , θ ^ M 2 ) C o v ( θ ^ M 1 , θ ^ M 3 ) C o v ( θ ^ M 1 , θ ^ M 2 ) V a r ( θ ^ M 2 ) C o v ( θ ^ M 2 , θ ^ M 3 ) C o v ( θ ^ M 1 , θ ^ M 3 ) C o v ( θ ^ M 2 , θ ^ M 3 ) V a r ( θ ^ M 3 ) .
For arbitrary 0 < γ < 1 , a 100 ( 1 γ ) % ACI of θ i is given by
θ ^ M i z γ / 2 V a r ( θ ^ M i ) , θ ^ M i + z γ / 2 V a r ( θ ^ M i ) , i = 1 , 2 , 3 ,
where z γ is the upper γ -th quantile of the standard normal distribution.
Sometimes, the previous obtained ACIs may have a negative lower bound which further yields worse performance for interval estimation of positive parameter. In order to avoid such drawback, one can use logarithmic transformation method and delta technique to construct the ACI in consequence. The asymptotic normality distribution of ln θ ^ M i can be obtained as
ln θ ^ M i ln θ i V a r ( ln θ ^ M i ) N ( 0 , 1 ) , i = 1 , 2 , 3 ,
where notation V a r ( ln θ ^ M i ) = V a r ( θ ^ M i ) / θ ^ M i 2 , i = 1 , 2 , 3 . Therefore, a 100 ( 1 γ ) % alternative ACI of θ i can be expressed as
θ ^ M i exp z γ / 2 V a r ( ln θ ^ M i ) , θ ^ M i exp z γ / 2 V a r ( ln θ ^ M i ) , i = 1 , 2 , 3 .
Let g ( Θ ) be arbitrary function of parameter Θ , the asymptotic distribution of g ( Θ ) can be constructed by using delta method as
g ^ M ( Θ ) g ( Θ ) N ( 0 , V ( Θ ^ M ) ) ,
where notations g ^ M ( Θ ) = g ( Θ ^ M ) is the MLE of g ( Θ ) and
V ( Θ ^ M ) = [ g T ( Θ ^ M ) ] I 1 ( Θ ^ M ) [ g ( Θ ^ M ) ] , g ( Θ ^ M ) = g ( Θ ) λ 1 , g ( Θ ) θ 2 , g ( Θ ) θ 3 | Θ = Θ ^ M .
Therefore, the 100 ( 1 γ ) % confidence interval of g ( Θ ) can be further constructed following similar approach as previous. Under this manner, the ACIs of the SF S ( t 0 ) and the HRF H ( t 0 ) can be obtained directly and detailed expressions are omitted for saving space.

4. Bayesian Inference

In inferential approach, Bayesian estimation has received much attention as an alternative to classical likelihood approach where the capability of incorporating prior information makes it very valuable in analysis. In this section, estimation problems based on traditional Bayesian and E-Bayesian approaches are discussed under different loss functions, respectively, and corresponding algorithms are given to solve Bayesian point estimation and credible interval.

4.1. Prior Information and Posterior Analysis

Assume that parameters θ 1 , θ 2 and θ 3 are statistically independent, and the gamma conjugate prior for θ i , i = 1 , 2 , 3 is implemented with hyper-parameters a i and b i as
π 0 ( θ i | a i , b i ) = b i a i Γ ( a i ) θ i a i 1 e b i θ i , a i > 0 , b i > 0 , θ i > 0 .
Therefore, the joint prior of ( θ 1 , θ 2 , θ 3 ) can be written as
π ( θ 1 , θ 2 , θ 3 ) = π 0 ( θ 1 | a 1 , b 1 ) π 0 ( θ 2 | a 2 , b 2 ) π 0 ( θ 3 | a 3 , b 3 ) = b 1 a 1 b 2 a 2 b 3 a 3 Γ ( a 1 ) Γ ( a 2 ) Γ ( a 3 ) θ 1 a 1 1 e b 1 θ 1 θ 2 a 2 1 e b 2 θ 2 θ 3 a 3 1 e b 3 θ 3 .
Based on (6) and (9), the joint posterior density of θ 1 , θ 2 , θ 3 can be expressed as
π ( θ 1 , θ 2 , θ 3 | d a t a ) i = 1 3 ( b i ω ) a i + n i Γ ( a i + n i ) θ i a i + n i 1 exp ( b i ω ) θ i ,
and the marginal posterior density of θ i is also gamma distribution given by
π ( θ i | d a t a ) = ( b i ω ) a i + n i Γ ( a i + n i ) θ i a i + n i 1 exp ( b i ω ) θ i .

4.2. Bayesian Estimation under Symmetry and Asymmetry Losses

Loss function plays an important role in Bayesian analysis and decision theory. The symmetric loss gives equal weighted to overestimation and underestimation, whereas asymmetric loss may be appropriate for estimation when overestimation and underestimation have different consequences. Suppose δ is an estimator of parameter θ , the following symmetric squared error loss and asymmetric linear-exponential (Linex) loss are used in this paper:
  • L 1  (squared error loss):  L ( δ , θ ) = ( δ θ ) 2 ,
  • L 2  (Linex loss):  L ( δ , θ ) = e d ( δ θ ) d ( δ θ ) 1 , d 0 ,
In addition, Bayesian estimator with squared error and Linex losses are the posterior expectation and 1 d ln E ( e d θ | d a t a ) , respectively. Through straightforward computation, the Bayes estimators of model parameters and reliability indices can be obtained as follows and the detailed proofs are omitted for simplicity and saving space.
Theorem 4.
Under squared error loss L 1 , one has
  • The Bayes estimator of θ i is given by
    θ ^ B i L 1 ( a i , b i ) = a i + n i b i ω , i = 1 , 2 , 3 .
  • The Bayes estimators of the SF and HRF can be expressed, respectively, as
    S ^ B L 1 ( t 0 ) = i = 1 3 b i ω b i ω ln G ¯ ( t 0 ) a i + n i and H ^ B L 1 ( t 0 ) = η ( t 0 ) i = 1 3 a i + n i b i ω ,
    with η ( t 0 ) = g ( t 0 ) G ¯ ( t 0 ) .
Theorem 5.
Under Linex loss L 2 , one has
  • The Bayes estimator of θ i is given by
    θ ^ B i L 2 ( a i , b i ) = 1 d ln b i ω b i ω + d a i + n i , i = 1 , 2 , 3 .
  • The Bayes estimators of the SF and HRF can be expressed, respectively, as
    S ^ B L 2 ( t 0 ) = 1 d ln m = 0 ( d ) m m ! i = 1 3 b i ω b i ω m ln G ¯ ( t 0 ) a i + n i ,
    and
    H ^ B L 2 ( t 0 ) = 1 d · i = 1 3 ln b i ω b i ω + η ( t 0 ) d a i + n i .
For arbitrary 0 < γ < 1 , a 100 ( 1 γ ) % Bayesian highest posterior density interval (HPDI) of θ i can be also constructed. A simple constrained optimization approach is provided as follows. Let C ( γ ) = [ θ i l ( γ ) , θ i r ( γ ) ] denote the 100 ( 1 γ ) % credibility interval of θ i for i = 1 , 2 , 3 under the posterior (11) satisfying
θ i l ( γ ) θ i r ( γ ) π ( θ i | d a t a ) d θ i = 1 γ .
Thus, the 100 ( 1 γ ) % HPDI of θ i denoted by C * ( γ ) = [ θ i l * ( γ ) , θ i r * ( γ ) ] can be obtained by solving following constrained optimization problem as
Minimize θ i r ( γ ) θ i l ( γ ) Subject to θ i l ( γ ) θ i r ( γ ) π ( θ i | d a t a ) d θ i = 1 γ , 0 < θ i l ( γ ) < θ i r ( γ ) .
which is equivalent to minimize the Lagrangian function defined by
h i = θ i r ( γ ) θ i l ( γ ) + τ i θ i l ( γ ) θ i r ( γ ) π ( θ i | d a t a ) d θ i 1 + γ ,
where τ i is the Lagrangian multiplier and solutions θ i l * ( γ ) and θ i r * ( γ ) can be further obtained numerically from following nonlinear equations as
π θ i l ( γ ) | d a t a = π θ i r ( γ ) | d a t a , θ i l ( γ ) θ i r ( γ ) π θ i | d a t a d θ i = 1 γ .

4.3. E-Bayesian Estimation

The E-Bayesian estimation was firstly proposed by Han [40], which is a special hierarchical Bayesian estimation scenario, and then was further discussed by Han [41,42]. Recently, the E-Bayesian approach has attracted wide attention and has been discussed by many authors. See, for example, some recent works of Basheer et al. [43], Wu et al. [44], Okasha and Wang [45] for comprehensive discussions.
Following the idea of Han [40], the hyper-parameters a i and b i should be selected to guarantee the prior density function π 0 ( θ i | a i , b i ) decreases in θ i under E-Bayesian scenario, which implies that 0 < a i < 1 and b i > 0 when π 0 ( θ i | a i , b i ) θ i < 0 . Under such requirements, it is assumed that the hyper-parameters a i and b i are statistically independent and have the following three priors as
π i 1 ( a i , b i ) = 1 c i , 0 < a i < 1 , 0 < b i < c i , i = 1 , 2 , 3 , π i 2 ( a i , b i ) = 2 b i c i 2 , 0 < a i < 1 , 0 < b i < c i , i = 1 , 2 , 3 , π i 3 ( a i , b i ) = 2 ( c i b i ) c i 2 , 0 < a i < 1 , 0 < b i < c i , i = 1 , 2 , 3 .
Theorem 6.
Under squared error loss L 1 , it is seen that
  • For prior π i 1 ( a i , b i ) , the E-Bayesian estimator of θ i are given by
    θ ^ E i L 1 π i 1 = 2 n i + 1 2 c i ln 1 c i ω , i = 1 , 2 , 3 .
  • For prior π i 2 ( a i , b i ) , the E-Bayesian estimator of θ i can be written as
    θ ^ E i L 1 π i 2 = 2 n i + 1 c i 2 ω ln 1 c i ω + c i i = 1 , 2 , 3 .
  • For prior π i 3 ( a i , b i ) , the E-Bayesian estimator of θ i can be expressed as
    θ ^ E i L 1 π i 3 = 2 n i + 1 c i 2 ( c i ω ) ln 1 c i ω c i i = 1 , 2 , 3 .
Proof. 
See Appendix D.    □
Theorem 7.
Under squared error loss L 1 , one has that
  • For prior π i 1 ( a i , b i ) , the E-Bayesian estimators of S ( t 0 ) and H ( t 0 ) are given by
    S ^ E L 1 π 1 ( t 0 ) = i = 1 3 ln G ¯ ( t 0 ) c i 0 c i ( b i ω ) n i b i ω ln G ¯ ( t 0 ) n i + 1 ln b i ω b i ω ln G ¯ ( t 0 ) 1 d b i
    and
    H ^ E L 1 π 1 ( t 0 ) = η ( t 0 ) i = 1 3 2 n i + 1 2 c i ln 1 c i ω .
  • For prior π i 2 ( a i , b i ) , the E-Bayesian estimators of S ( t 0 ) and H ( t 0 ) can be written as
    S ^ E L 1 π 2 ( t 0 ) = i = 1 3 2 ln G ¯ ( t 0 ) c i 2 0 c i b i ( b i ω ) n i b i ω ln G ¯ ( t 0 ) n i + 1 ln b i ω b i ω ln G ¯ ( t 0 ) 1 d b i
    and
    H ^ E L 1 π 2 ( t 0 ) = η ( t 0 ) i = 1 3 2 n i + 1 c i 2 c i + ω ln 1 c i ω .
  • For prior π i 3 ( a i , b i ) , the E-Bayesian estimators of S ( t 0 ) and H ( t 0 ) can be expressed as
    S ^ E L 1 π 3 ( t 0 ) = i = 1 3 2 ln G ¯ ( t 0 ) c i 2 0 c i ( c i b i ) ( b i ω ) n i b i ω ln G ¯ ( t 0 ) n i + 1 ln b i ω b i ω ln G ¯ ( t 0 ) 1 d b i
    and
    H ^ E L 1 π 3 ( t 0 ) = η ( t 0 ) i = 1 3 2 n i + 1 c i 2 ( c i ω ) ln 1 c i ω c i .
Proof. 
See Appendix E.    □
Analogously to Theorems 6 and 7, the E-Bayesian estimators of parameters and reliability indices can be obtained under Linex loss, and the proofs are omitted for concision.
Theorem 8.
Under Linex loss L 2 , it is seen that
  • For prior π i 1 ( a i , b i ) , the E-Bayesian estimator of θ i are given by
    θ ^ E i L 2 π i 1 = 2 n i + 1 2 c i d ( c i ω ) ln 1 + d c i ω + ω ln 1 d ω + d ln 1 + c i d ω , i = 1 , 2 , 3 .
  • For prior π i 2 ( a i , b i ) , the E-Bayesian estimator of θ i can be written as
    θ ^ E i L 2 π i 2 = 2 n i + 1 2 c i 2 d c i 2 ω 2 ln 1 + d c i ω + 2 ω d d 2 ln 1 + c i d ω + ω 2 ln 1 d ω + c i d , i = 1 , 2 , 3 .
  • For prior π i 3 ( a i , b i ) , the E-Bayesian estimator of θ i can be expressed as
    θ ^ E i L 2 π i 3 = 2 n i + 1 2 c i 2 d ( c i ω ) 2 ln 1 + d c i ω + d 2 + 2 c i d 2 ω d ln 1 + c i d ω + 2 c i ω ω 2 ln 1 d ω c i d , i = 1 , 2 , 3 .
Theorem 9.
Let D be the range of ν = ( a 1 , b 1 , a 2 , b 2 , a 3 , b 3 ) , it is noted under Linex loss that
  • For prior π i 1 ( a i , b i ) , the E-Bayesian estimators of S ( t 0 ) and H ( t 0 ) are given by
    S ^ E L 2 π 1 ( t 0 ) = 1 c 1 c 2 c 3 d D ln m = 0 ( d ) m m ! i = 1 3 b i ω b i ω m ln G ¯ ( t 0 ) a i + n i d ν
    and
    H ^ E L 2 π 1 ( t 0 ) = i = 1 3 2 n i + 1 2 c i d ( c i ω ) ln 1 + η ( t 0 ) d c i ω + ω ln 1 η ( t 0 ) d ω + η ( t 0 ) d ln 1 + c i η ( t 0 ) d ω .
  • For prior π i 2 ( a i , b i ) , the E-Bayesian estimators of S ( t 0 ) and H ( t 0 ) can be written as
    S ^ E L 2 π 2 ( t 0 ) = 1 d D ln m = 0 ( d ) m m ! k = 1 3 b i ω b i ω m ln G ¯ ( t 0 ) a i + n i i = 1 3 2 b i c i 2 d ν
    and
    H ^ E L 2 π 2 ( t 0 ) = i = 1 3 2 n i + 1 2 c i 2 d c i 2 ω 2 ln 1 + η ( t 0 ) d c i ω + ω 2 ln 1 η ( t 0 ) d ω + 2 ω η ( t 0 ) d η ( t 0 ) d 2 ln 1 + c i η ( t 0 ) d ω + c i η ( t 0 ) d .
  • For prior π i 3 ( a i , b i ) , the E-Bayesian estimators of S ( t 0 ) and H ( t 0 ) can be expressed as
    S ^ E L 2 π 3 ( t 0 ) = 1 d D ln m = 0 ( d ) m m ! i = 1 3 b i ω b i ω m ln G ¯ ( t 0 ) a i + n i i = 1 3 2 ( c i b i ) c i 2 d ν
    and
    H ^ E L 2 π 3 ( t 0 ) = i = 1 3 2 n i + 1 2 c i 2 d ( c i ω ) 2 ln 1 + η ( t 0 ) d c i ω + 2 c i ω ω 2 ln 1 η ( t 0 ) d ω + η ( t 0 ) d 2 + 2 c i η ( t 0 ) d 2 ω η ( t 0 ) d ln 1 + c i η ( t 0 ) d ω c i η ( t 0 ) d .
In order to construct HPDI estimates under E-Bayesian scenario, a sampling based approach is proposed as follows. Based on (6), (9) and (12), the posterior density function of θ i , a i , b i , i = 1 , 2 , 3 with respect to hyper-prior π i 1 , π i 2 , π i 3 can be expressed, respectively, as
π E B π i 1 ( θ i , a i , b i | d a t a ) 1 c i b i a i Γ ( a i ) θ i a i + n i 1 e ( b i ω ) θ i , π E B π i 2 ( θ i , a i , b i | d a t a ) 2 c i 2 b i a i + 1 Γ ( a i ) θ i a i + n i 1 e ( b i ω ) θ i , π E B π i 3 ( θ i , a i , b i | d a t a ) 2 ( c i b i ) c i 2 b i a i Γ ( a i ) θ i a i + n i 1 e ( b i ω ) θ i .
Further, the conditional marginal posterior distribution of θ i , a i , b i , i = 1 , 2 , 3 can be obtained from the posterior distribution (13), respectively, as follows:
θ i π i 1 , θ i π i 2 , θ i π i 3 G a m m a ( a i + n i , b i ω ) i = 1 , 2 , 3 ,
a i π i 1 , a i π i 2 , a i π i 3 ( b i θ i ) a i Γ ( a i ) , i = 1 , 2 , 3 ,
and
b i π i 1 b i a i e θ i b i , b i π i 2 b i a i + 1 e θ i b i , b i π i 3 ( c i b i ) b i a i e θ i b i .
It can be seen from (15) and (16) that the full conditional posterior distribution of a i π i 1 , a i π i 2 , a i π i 3 and b i π i 1 , b i π i 2 , b i π i 3 cannot be analyzed by using some well-known models. Therefore, the Markov-Chain Monte Carlo (MCMC) method is used here for complex computation. To be specific, the Metropolis-Hastings (M-H) algorithm via Gibbs sampling is implemented for posterior sampling under the E-Bayesian scenario. Then, a MCMC approach called Algorithm 1 is provided to obtain the E-Batesian estimates of θ 1 , θ 2 , θ 3 , SF S ( t ) and HRF  H ( t ) .

4.4. Some Results of Bayesian and E-Bayesian Estimation

In this section, the posterior risk (PR) of Bayesian estimators and the E-posterior risk (EPR) of E-Bayesian estimators are further presented for complementary which generally can be used to measure the associated estimated risk.
Theorem 10.
For Bayes estimators
  • under squared error loss L 1 , the PR of θ ^ B i L 1 can be expressed as
    PR θ ^ B i L 1 ( a i , b i ) = a i + n i ( b i ω ) 2 , i = 1 , 2 , 3 .
  • under Linex loss L 2 , the PR of θ ^ B i L 2 can be expressed as
    PR θ ^ B i L 2 ( a k , b i ) = d ( a i + n i ) b i ω ( a i + n i ) ln 1 + d b i ω , i = 1 , 2 , 3 .
Proof. 
See Appendix F.    □
Following Han [46], the EPR of various E-Bayesian estimators are also obtained.
Theorem 11.
For E-Bayesian estimators,
  • under SE loss, the EPR of θ ^ E i L 1 π i 1 , θ ^ E i L 1 π i 2 , θ ^ E i L 1 π i 3 with respect to priors π k 1 , π k 2 , π i 3 can be expressed, respectively, as
    E PR θ ^ E i L 1 π i 1 ( a i , b i ) = 2 n i + 1 2 c i 1 ω c i 1 ω , E PR θ ^ E i L 1 π i 2 ( a i , b i ) = 2 n i + 1 c i 2 ln 1 c i ω ω 1 c i ω + 1 ω , E PR θ ^ E i L 1 π i 3 ( a i , b i ) = 2 n i + 1 c i 2 ln 1 c i ω c i ω .
  • under Linex loss, the EPR of θ ^ E i L 2 π i 1 , θ ^ E i L 2 π i 2 , θ ^ E i L 2 π i 3 with respect to priors π i 1 , π i 2 , π i 3 can be expressed, respectively, as
    E PR θ ^ E i L 2 π i 1 ( a i , b i ) = 2 n i + 1 2 c i d d ln 1 c i ω ( c i ω ) ln 1 + d c i ω ω ln 1 d ω d ln 1 + c i d ω , E PR θ ^ E i L 2 π i 2 ( a i , b i ) = 2 n i + 1 2 c i 2 d 2 ω d ln 1 c i ω + c i d c i 2 ω 2 ln 1 + d c i ω 2 ω d d 2 ln 1 + c i d ω ω 2 ln 1 d ω , E PR θ ^ E i L 2 π i 3 ( a i , b i ) = 2 n i + 1 2 c i 2 d ( 2 c i d 2 ω d ) ln 1 c i ω c i d ( c i ω ) 2 ln 1 + d c i ω d 2 + 2 c i d 2 ω d ln 1 + c i d ω 2 c i ω ω 2 ln 1 d ω .
Proof. 
See Appendix G.    □
In addition, some order results among E-Bayesian estimators are also presented as follows.
Theorem 12.
Let 0 < c i ω < 1 , the E-Bayesian estimators under squared error loss feature following results
  • θ ^ E i L 1 π i 2 < θ ^ E i L 1 π i 1 < θ ^ E i L 1 π i 3 ,
  • lim ω θ ^ E i L 1 π i 1 = lim ω θ ^ E i L 1 π i 2 = lim ω θ ^ E i L 1 π i 3 .
Proof. 
See Appendix H.    □
Theorem 13.
Let 0 < c i ω < 1 , the E-Bayesian estimators from Linex loss have following results
  • for d < ω c i , θ ^ E i L 2 π i 2 < θ ^ E i L 2 π i 1 < θ ^ E i L 2 π i 3 .
  • lim ω θ ^ E i L 2 π i 1 = lim ω θ ^ E i L 2 π i 2 = lim ω θ ^ E i L 2 π i 3 .
Proof. 
See Appendix I.    □
Algorithm 1: E-Bayesian interval estimation based M-H algorithm via Gibbs sampling
Step 1
For prior π i 1 in (12), given c i , set initial value a i ( 0 ) = 0.5 , b i ( 0 ) = c i / 2 , θ i ( 0 ) = θ ^ M i , i = 1 , 2 , 3 .
Step 2
Set s = 1 .
Step 3
For given ( b i ( s 1 ) , θ i ( s 1 ) ) and (15), generate a i ( s ) in following steps.
  Step 3.1
Propose a conditional transition probability function q ( a i ( * ) | a i ( o l d ) ) for a next a i ( * ) , given the current a i ( o l d ) . The transition probability function could be a normal distribution with the current status a i ( o l d ) as the center location and a reasonable variance or a uniform distribution over a reasonable range.
  Step 3.2
Let a i ( o l d ) = a i ( s 1 ) , and generate a i ( * ) from the proposed conditional transition distribution q ( a i ( * ) | a i ( o l d ) ) and u U ( 0 , 1 ) .
  Step 3.3
Calculate the acceptance probability
r = min 1 , π ( a i ( * ) | d a t a ) q ( a i ( o l d ) | a i ( * ) ) π ( a i ( o l d ) | d a t a ) q ( a i ( * ) | a i ( o l d ) ) .
  Step 3.4
If u r , set a i ( s ) = a i ( * ) . Otherwise, set a i ( s ) = a i ( o l d ) .
Step 4
For ( a i ( s ) , θ i ( s 1 ) ) , generate b i ( s ) from (16) using similar method of Step 3.
Step 5
For given ( a i ( s ) , b i ( s ) ) , generate θ i ( s ) in (14) from G a m m a ( a i + n i , b i ω ) .
Step 6
Set s = s + 1 .
Step 7
Repeat Steps 3–6 N times and obtain N samples of θ i , i = 1 , 2 , 3 as θ i ( 1 ) , θ i ( 2 ) , , θ i ( N ) .
Step 8
Let η = η ( θ 1 , θ 2 , θ 3 ) be any function of ( θ 1 , θ 2 , θ 3 ) , η ( 1 ) , η ( 2 ) , , η ( N ) can be obtained based N samples of ( θ 1 , θ 2 , θ 3 ) .
Step 9
Set K being burn-in period, and arrange last N K estimates of η in ascending
order as η * [ 1 ] , η * [ 2 ] , , η * [ M ] with M = N K . Then, for arbitrary 0 < γ < 1 ,
a 100 ( 1 γ ) % credible interval of η can be obtained as
η [ l ] , η [ l + M [ M γ + 1 ] ] , l = 1 , 2 , , [ M γ ] ,
where [ y ] denotes the greatest integer less than or equal to y. Therefore, the 100 ( 1 γ ) % HPDI can be constructed as the l * -th one satisfying
η [ l * + M [ M γ + 1 ] ] η [ l * ] = min l = 1 [ M γ ] η [ l + M [ M γ + 1 ] ] η [ l ] .
Step 10
For prior π i 2 and π i 3 , i = 1 , 2 , 3 in (12), the results could be obtained similarly as above procedure.

5. Numerical Illustration

5.1. Simulation Studies

Extensive simulation studies are conducted for investigating the performance of the proposed methods when dependent MinRSSU competing risks samples are available, the associated point estimates are evaluated in terms of absolute bias (AB) and mean squared error (MSE), and the interval estimates are investigated by average length (AL) and coverage probability (CP). In this illustration, the following approach entitled Algorithm 2 is provided to generate dependent MinRSSU competing risks samples.
Algorithm 2: Generation dependent MinRSSU competing risks samples
Step 1
For given c 1 , c 2 and c 3 , generate a i and b i , i = 1 , 2 , 3 from priors (12).
Step 2
Based on generated a i and b i , generate θ i from gamma prior (8).
Step 3
For ( θ 1 , θ 2 , θ 3 ) and n, generate n ( n + 1 ) 2 random data { t 1 , t 2 , , t n } from the PHR distribution with parameters θ 123 .
Step 4
Generate a group random samples of 1 , 2 , 3 each has the size of n ( n + 1 ) 2 from a multinomial distribution with parameters θ 1 θ 123 , θ 2 θ 123 and θ 3 θ 123 . When 1 appears that the sample { t i , i = 1 , 2 , , n } generated in Step 3 is based on cause one; when 2 appears that the sample { t i , i = 1 , 2 , , n } is treated as the failure times due to cause two; otherwise it is assumed that both failure causes occur at t i simultaneously for i = 1 , 2 , , n . Then, a complete dependent competing risks sample of size n ( n + 1 ) 2 is obtained.
Step 5
For the data presented in Step 3, selecting n samples by using MinRSSU methods. Therefore, a group of MinRSSU competing risks sample is generated.
Since the proposed MOBPHR model is a general family, the following three distributions are selected as its special cases in the study:
Exponential baseline model : F E ( t ; θ ) = 1 e θ t with G E ( t ) = 1 e t , Pareto baseline model : F P ( t ; θ ) = 1 1 t θ with G P ( t ) = 1 1 t , Rayleigh baseline model : F R ( t ; θ ) = 1 e t 2 2 θ with G R ( t ) = 1 e t 2 2 .
In the simulation study, two sets of random chosen values of the parameters ( θ 1 , θ 2 , θ 3 ) = ( 0.67 , 0.50 , 0.83 ) and ( 0.30 , 0.40 , 0.50 ) , correspondingly the associated hyper-parameters ( c 1 , c 2 , c 3 ) are taken to be ( 1.5 , 2 , 1.2 ) and ( 3 , 2.5 , 2 ) satisfying the prior means are the same as the original means. Therefore, the criteria quantities ABs, MSEs, CPs and ALs are carried out based on 10,000 repetitions. All simulation results are tabulated in Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9 and Table 10, where the significance level is 0.95 for interval estimates, and the results of model parameters are just presented for θ 1 and θ 3 for concision. In addition, it is worth mentioning that the Bayesian interval estimates of reliability indices S ( t 0 ) and H ( t 0 ) are obtained by using sampling approach like Algorithm 1 based on the posterior function (11).
From the simulated results in Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9 and Table 10, one can observe the following conclusions for point and interval estimates:
  • For point estimation, AB and MSEs of point estimates from E-Bayesian methods under π 1 , π 2 and π 3 , as well as Bayesian and maximum likelihood methods, decrease with increase of the sample size n.
  • For a given n, AB and MSEs of Bayes and E-Bayesian estimates under π 1 , π 2 and π 3 are smaller than those of MLEs. In addition, three E-Bayesian estimates are superior to the Bayesian results, where estimates under prior π 2 perform best among all E-Bayesian estimates.
  • For interval estimation, when sample size n increases, ALs of ACIs, HPDI and intervals from E-Bayesian scenario decrease in all cases.
  • Under fixed sample size n, ALs of Bayesian HPDIs are smaller than those of ACIs, but slightly larger than those of credible intervals obtained from E-Bayesian approaches in general.
  • The CPs of different interval estimates are all close to nominal significance level.
  • The PRs of E-Bayesian are all smaller than those of Bayesian estimates, where the PRs under the prior π 2 seems to perform best among all estimates.

5.2. Real Data Illustrations

A real data example is presented for showing the applications of the proposed methods. In this data set, the time (in years) from human immunodeficiency virus (HIV) infection to the development of acquired immune deficiency syndrome (AIDS) or switch to syncytium inducing (SI) phenotype in 329 homosexual men was discussed. The data set was also original reported by Putter et al. [47] and Geskus [48], where events AIDS and SI can be regarded as competing risks denoted as δ i = 1 and δ i = 2 , respectively. In order to obtain dependent competing risks data, we randomly select some data from the overall data and treated these data as failure times due to causes 1 and 2 simultaneous, which is denoted as δ i = 3 . Before further proceeding, the MOBPHR model is implemented with exponential, Pareto and Rayleigh baseline distributions for fitting the data set. Since there is no general goodness-of-fit approach for bivariate model, here we provide the Kolmogorov–Smirnov (KS) distance and associated p-value for three baseline models by ignoring the failure causes and the results are presented in Table 11. In addition, empirical cumulative distributions via theoretical distribution plot, Probability-Probability (P-P) and Quantile-Quantile (Q-Q) plots under three baseline distributions are presented in Figure 1. It is conducted that the Rayleigh distribution can be used as a proper model to fit the real data set, and then we use the MOBPHR model with Rayleigh baseline marginal distribution for further illustration.
Furthermore, a group of MinRSSU competing risk samples with a sample size of n = 25 were obtained by sampling AIDSSI data through MinRSSU sampling strategy. The sample data set is shown in Table 12. Since there is no historical information for model parameters, the hyper-parameters are set to be ( c 1 , c 2 , c 3 ) = ( 0.01 , 0.01 , 0.01 ) in this illustration indicating that the priors are proper but almost non-informative. The point and interval estimates are presented in Table 13 and Table 14 with significance level 0.95 for interval estimation and d = 1 in Linex loss, where the associated interval lengths are also provided in squared brackets.
Based on the listed results shown in Table 13 and Table 14, it is observed that the classical and Bayes estimates appear close to each other, whereas the E-Bayesian HPDIs and ACIs have shortest and longest interval length, respectively. For further illustration, CDF and PDF of Rayleigh distribution for the minimum of the failure causes are plotted in Figure 2 by using different parameter estimates, where the notations E-Bayes-1, E-Bayes-2 and E-Bayes-3 represent the E-Bayesian estimates obtained under priors π 1 , π 2 and π 3 , respectively. From Figure 2, it is shown that the plots of Rayleigh PDF and CDF via classical and Bayes estimates also appear similar performance. In addition, Figure 3 shows the trend diagram of the M-H sampling algorithm for the parameters θ 1 , θ 2 and θ 3 , it is also noted that the E-Bayesian sampling has good mixing performance.
Finally, Figure 4 also shows the plots of HRF H ( t ) via different point estimates in mission time t ( 0 , 15 ) , and the plots of interval lengths of HRF are presented as well. It can be seen that there are certain differences between classical and Bayesian point estimates of HRF, and the associated Bayesian and E-Bayesian point estimates are similar. In addition, it is also seen that the ACI of HRF at different mission time features longest interval length, whereas the E-Bayesian HPDI of HRF features shortest interval length.

6. Conclusions

This paper studies estimation of dependent competing risk model under the MinRSSU scenario. When the MOBPHR distribution is proposed for modeling failure causes, point and interval estimation of the unknown parameters as well as the reliability indices are estimated by classical, Bayesian and E-Bayesian methods, respectively. Simulation studies and a real-life example are carried out to investigate the performance of different methods, the results indicate that all classical and Bayesian estimates work satisfactory, and the proposed E-Bayesian method performs better than other results in general. Although the MOBPHR distribution is constructed in this paper to model the lifetime of dependent failure causes, this work can be extended to other bivariate models when the MinRSSU competing risks samples are available. In addition, there are many types of modified RSS sampling strategies. In this study, MinRSSU sampling method is selected according to the characteristics of life-time data. In future studies, RSS or modified RSS sampling method should be selected according to different characteristics of life-time data for corresponding statistical analysis. Secondly, the PHR model considered in this paper is a single-parameter distribution model. Due to the complexity of failure causes, the competing risks data may come from more complex distribution models. Therefore, the multi-parameter distribution models, such as generalized gamma, beta prime and generalized inverse Gaussian, will be considered in future studies.

Author Contributions

Funding acquisition, L.W.; Methodology, Y.Z. and L.W.; Writing-original draft, T.-R.T. and Y.M.T. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by Yunnan Fundamental Research Projects (No. 202101AT070103).

Data Availability Statement

The numerical data used to support the findings of this study are included within the article.

Acknowledgments

Authors would like to thank the Editor and reviewers for their valuable comments and suggestions which improve this paper significantly.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Proof of Theorem 1

Using the independent property of U 1 , U 2 and U 3 , one has
S ( T 1 , T 2 ) ( t 1 , t 2 ) = P ( T 1 > t 1 , T 2 > t 2 ) = P ( U 1 > t 1 , U 2 > t 2 , U 3 > max t 1 , t 2 ) = S ( t 1 ; θ 1 ) S ( t 2 ; θ 2 ) S ( max t 1 , t 2 ) = S ( t 1 ; θ 1 ) S ( t 2 ; θ 2 + θ 3 ) , if t 1 < t 2 ; S ( t 1 ; θ 1 + θ 3 ) S ( t 2 ; θ 2 ) , if t 1 > t 2 ; S ( t ; θ 1 + θ 2 + θ 3 ) , if t 1 = t 2 = t .
Then the assertion is proved.

Appendix B. The Proof of Theorem 2

For MinRSSU competing risks data ( t i , δ i ) , i = 1 , , n , a brief explanation of the probability contribution of ( t i , δ i ) is shown here. For ( t i , δ i ) = ( t i , 1 ) , it means that the first failure of unit in i-th group is yielded due to cause one, then one equivalently has t i = T i 1 and t i < T i 2 . Using distribution property of order statistic, the likelihood contribution of data ( t i , δ i ) = ( t i , 1 ) can be established as
P ( t i < T i 1 < t i + d t i , T i 2 > t i ) = S T 1 , T 2 ( t 1 , t 2 ) t 1 | ( t 1 , t 2 ) = ( t i , t i ) · i · S t i ; θ 123 i 1 d t i = i θ 1 g ( t i ) [ G ¯ ( t i ) ] i θ 123 1 .
Similarly, the likelihood contribution of ( t i , δ i ) = ( t i , 2 ) can be obtained as well.
For ( t i , δ i ) = ( t i , 3 ) , failure is yielded by causes one and two simultaneously implying that T i 1 = T i 2 = t i . Then the likelihood contribution is given by
P ( t i < T i 1 = T i 2 < t i + d t i ) = θ 3 θ 123 f t i ; θ 123 i · S t i ; θ 123 i 1 d t i = i θ 3 g ( t i ) [ G ¯ ( t i ) ] i θ 123 1 .
Therefore, the assertion is completed.

Appendix C. The Proof of Theorem 3

By taking derivatives of ( θ 1 , θ 2 , θ 3 ) with respect to θ 1 , θ 2 and θ 3 , one can obtain MLE θ ^ M i , i = 1 , 2 , 3 directly. Using the inequality ln x x 1 , x > 0 , let x = θ i θ ^ M i , then one has
n i ln θ i θ i ω n i + n i ln θ ^ M i , i = 1 , 2 , 3 .
Using the above inequality, it can be obtained that
( θ 1 , θ 2 , θ 3 ) ln n 123 ! + n 1 ln θ ^ M 1 + n 2 ln θ ^ M 2 + n 3 ln θ ^ M 3 n 1 n 2 n 3 + i = 1 n 123 ln g t i G ¯ t i
Since n i = θ ^ M i ω , and n 1 + n 2 + n 3 = ( θ ^ M 1 + θ ^ M 2 + θ ^ M 3 ) ω = θ ^ M 123 ω . Then it is observed that
( θ 1 , θ 2 , θ 3 ) ln n 123 ! + n 1 ln θ ^ M 1 + n 2 ln θ ^ M 2 + n 3 ln θ ^ M 3 + θ ^ M 123 ω + i = 1 n 123 ln g t i G ¯ t i = ( θ ^ M 1 , θ ^ M 2 , θ ^ M 3 ) ,
where equality holds if and only if θ i = θ ^ M i , i = 1 , 2 , 3 . Therefore, the assertion is completed.

Appendix D. The Proof of Theorem 6

Under squared error loss, the E-Bayesian estimator of θ i , i = 1 , 2 , 3 can be obtained by using priors π i 1 ( a i , b i ) , π i 2 ( a i , b i ) and π i 3 ( a i , b i ) as follows
θ ^ E i L 1 π i 1 ( a i , b i ) = 0 c i 0 1 θ ^ B i L 1 ( a i , b i ) π i 1 ( a i , b i ) d a i d b i , θ ^ E i L 1 π i 2 ( a i , b i ) = 0 c i 0 1 θ ^ B i L 1 ( a i , b i ) π i 2 ( a i , b i ) d a i d b i ,
and
θ ^ E i L 1 π i 3 ( a i , b i ) = 0 c i 0 1 θ ^ B i L 1 ( a i , b i ) π i 3 ( a i , b i ) d a i d b i .
Therefore, the assertion is completed.

Appendix E. The Proof of Theorem 7

Under squared error loss, the E-Bayesian estimator of reliability indices S ( t 0 ) and H ( t 0 ) with respect to prior π i 1 ( a i , b i ) , i = 1 , 2 , 3 can be expressed, respectively, as
S ^ E L 1 π 1 ( t 0 ) = D S ^ B L 1 ( t 0 ) π 1 1 ( a 1 , b 1 ) π 2 1 ( a 2 , b 2 ) π 3 1 ( a 3 , b 3 ) d ν
and
H ^ E L 1 π 1 ( t 0 ) = D H ^ B L 1 ( t 0 ) π 1 1 ( a 1 , b 1 ) π 2 1 ( a 2 , b 2 ) π 3 1 ( a 3 , b 3 ) d ν .
When the priors π i 2 ( a i , b i ) and π i 3 ( a i , b i ) are available, the associated E-Bayesian estimator of S ( t 0 ) and H ( t 0 ) could be also obtained similarly. The details are omitted for simplicity. Therefore, the assertion is completed.

Appendix F. The Proof of Theorem 10

Under squared error and Linex loss functions, the PR of the Bayes estimator can be obtained with respect to posterior density (11) as
PR θ ^ B i L 1 a i , b i = 0 θ ^ B i L 1 θ i 2 π θ i | d a t a d θ i ,
and
PR θ ^ B i L 2 a i , b i = 0 e d ( θ ^ B i L 2 θ i ) d θ ^ B i L 2 θ i 1 π θ i | d a t a d θ i ,
for i = 1 , 2 , 3 . Therefore, the assertion is completed.

Appendix G. The Proof of Theorem 11

Following Han [46], the EPR of θ ^ E i L 1 π i 1 and θ ^ E i L 2 π i 1 under squared error and Linex loss functions can be expressed, respectively, as
E PR θ ^ E i L 1 π i 1 = 0 c i 0 1 0 θ ^ B i θ i 2 π θ i | d a t a π i 1 a i , b i d θ i d a i d b i ,
and
E PR θ ^ E i L 2 π i 1 = 0 c i 0 1 0 e d ( θ ^ B i L 2 θ i ) d θ ^ B i L 2 θ i 1 π θ i | d a t a π i 1 a i , b i d θ i d a i d b i .
In addition, the EPR of θ ^ E i L 1 π i 2 , θ ^ E i L 1 π i 3 and θ ^ E i L 2 π i 2 , θ ^ E i L 2 π i 3 can be obtained under both loss functions similarly, and the details are omitted for saving space. Therefore, the assertion is completed.

Appendix H. The Proof of Theorem 12

Using expressions of θ ^ E i L 1 π i 1 , θ ^ E i L 1 π i 2 and θ ^ E i L 1 π i 3 , it is seen that
θ ^ E i L 1 π i 3 θ ^ E i L 1 π i 1 = θ ^ E i L 1 π i 1 θ ^ E i L 1 π i 2 = 2 n i + 1 c i 1 2 ω c i ln 1 c i ω 1 .
For 0 < x = c i ω < 1 and using Taylor expansion ln ( 1 x ) = m = 1 x m m , one has
1 2 1 x ln 1 x 1 = m = 1 m 1 2 m ( m + 1 ) x m > 0 .
Thus, the order relation among θ ^ E i L 1 π i 1 , θ ^ E i L 1 π i 2 and θ ^ E i L 1 π i 3 is proved. Whereas the second result can be shown by taking direct limitation. Therefore, the assertion is completed.

Appendix I. The Proof of Theorem 13

Using expressions of θ ^ E i L 2 π i 1 , θ ^ E i L 2 π i 2 and θ ^ E i L 2 π i 3 , it is seen that
θ ^ E i L 1 π i 3 θ ^ E i L 1 π i 1 = θ ^ E i L 1 π i 1 θ ^ E i L 1 π i 2 = 2 n i + 1 2 d ω d c i 2 ω d c i ln 1 c i ω d ω c i 2 ω c i ln 1 c i ω d c i = J ( d ) .
Following similar line as Theorem 12 and denoting x 1 = c i ω ( 0 , 1 ) and x 2 = c i ω d ( 0 , 1 ) for d < ω c i , then using Taylor expansion ln ( 1 x ) = m = 1 x m m , one has
J ( d ) = 2 n i + 1 2 d m = 1 1 m ( m + 1 ) c i ω d m 1 c i ω m 1 ,
implying J ( d ) > 0 , the first result is proved. The second result is established by taking direct limitation. Therefore, the assertion is completed.

References

  1. McIntyre, G.A. A method for unbiased selective sampling using ranked sets. Aust. J. Agric. Res. 1952, 3, 385–390. [Google Scholar] [CrossRef]
  2. Takahasi, K.; Wakimoto, K. On unbiased estimates of the population mean based on the sample stratified by means of ordering. Ann. Inst. Stat. Math. 1968, 20, 1–31. [Google Scholar] [CrossRef]
  3. Takahasi, K. Practical note on estimation of population means based on samples stratified by means of ordering. Ann. Inst. Stat. Math. 1970, 22, 421–428. [Google Scholar] [CrossRef]
  4. Dell, T.R.; Clutter, J.L. Ranked set sampling theory with order statistics background. Biometrics 1972, 28, 545–555. [Google Scholar] [CrossRef]
  5. Basikhasteh, M.; Lak, F.; Afshari, M. Bayesian estimation of stress-strength reliability for two-parameter bathtub-shaped lifetime distribution based on maximum ranked set sampling with unequal samples. J. Stat. Comput. Simul. 2020, 90, 1–16. [Google Scholar] [CrossRef]
  6. Esemen, M.; Gürler, S. Parameter estimation of generalized Rayleigh distribution based on ranked set sample. J. Stat. Comput. Simul. 2018, 88, 615–628. [Google Scholar] [CrossRef]
  7. Samawi, H.M.; Rochani, H.; Linder, D.; Chatterjee, A. More efficient logistic analysis using moving extreme ranked set sampling. J. Appl. Stat. 2015, 44, 753–766. [Google Scholar] [CrossRef]
  8. Zamanzade, E.; Mahdizadeh, M.; Samawi, H.M. Efficient estimation of cumulative distribution function using moving extreme ranked set sampling with application to reliability. AStA-Adv. Stat. Anal. 2020, 104, 485–502. [Google Scholar] [CrossRef]
  9. Muttlak, H.A. Median ranked set sampling. J. Appl. Stat. Sci. 1997, 6, 577–586. [Google Scholar]
  10. Samawi, H.M.; Ahmed, M.S.; Abu-Dayyeh, W. Estimating the population mean using extreme ranked set sampling. Biom. J. 1996, 38, 577–586. [Google Scholar] [CrossRef]
  11. Al-Saleh, M.F.; Al-Kadiri, M.A. Double-ranked set sampling. Stat. Probab. Lett. 2000, 48, 205–212. [Google Scholar] [CrossRef]
  12. Al-Saleh, M.F.; Al-Omari, A.I. Multistage ranked set sampling. J. Stat. Plan. Inference 2002, 102, 273–286. [Google Scholar] [CrossRef]
  13. Muttlak, H.A. Investigating the use of quartile ranked set samples for estimating the population mean. Appl. Math. Comput. 2003, 146, 437–443. [Google Scholar] [CrossRef]
  14. Al-Odat, M.; Al-Saleh, M.F. A variation of ranked set sampling. J. Appl. Stat. Sci. 2001, 10, 137–146. [Google Scholar]
  15. Biradar, B.S.; Santosha, C.D. Estimation of the mean of the exponential distribution using maximum ranked set sampling with unequal samples. Open J. Stat. 2014, 4, 641–649. [Google Scholar] [CrossRef] [Green Version]
  16. Yao, D.; Chen, W.; Yang, R.; Long, C. Fisher information in moving extreme ranked set sampling with application to parameter estimation. J. Syst. Sci. Complex. 2022, 35, 361–372. [Google Scholar] [CrossRef]
  17. Qiu, G.; Eftekharian, A. Extropy information of maximum and minimum ranked set sampling with unequal samples. Commun. Stat.-Theory Methods 2021, 50, 2979–2995. [Google Scholar] [CrossRef]
  18. Rahmani, H.; Razmkhah, M. Perfect ranking test in moving extreme ranked set sampling. Stat. Pap. 2017, 58, 855–875. [Google Scholar] [CrossRef]
  19. Ren, J.; Gui, W. Statistical analysis of adaptive type-II progressively censored competing risks for Weibull models. Appl. Math. Model. 2021, 98, 323–342. [Google Scholar] [CrossRef]
  20. Wang, L.; Tripathi, Y.M.; Lodhi, C. Inference for Weibull competing risks model with partially observed failure causes under generalized progressive hybrid censoring. J. Comput. Appl. Math. 2020, 368, 112537. [Google Scholar] [CrossRef]
  21. Sarhan, A.M.; Hamilton, D.C.; Smith, B. Statistical analysis of competing risks models. Reliab. Eng. Syst. Saf. 2010, 95, 953–962. [Google Scholar] [CrossRef]
  22. Mahto, A.K.; Lodhi, C.; Tripathi, Y.M.; Wang, L. Inference for partially observed competing risks model for Kumaraswamy distribution under generalized progressive hybrid censoring. J. Appl. Stat. 2021, 49, 2064–2092. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, Y.C.; Emura, T.; Fan, T.H.; Lo, S.M.S.; Wilke, R.A. Likelihood-based inference for a frailty-copula model based on competing risks failure time data. Qual. Reliab. Eng. Int. 2020, 36, 1622–1638. [Google Scholar] [CrossRef]
  24. Gorfine, M.; Hsu, L. Frailty-based competing risks model for multivariate survival data. Biometrics 2011, 67, 415–426. [Google Scholar] [CrossRef]
  25. Lo, S.M.S.; Wilke, R.A. A copula model for dependent competing risks. J. R. Stat. Soc. C-Appl. 2010, 59, 359–376. [Google Scholar] [CrossRef] [Green Version]
  26. Fan, T.H.; Wang, Y.F.; Ju, S.K. A competing risks model with multiply censored reliability data under multivariate Weibull distributions. IEEE. Trans. Reliab. 2019, 68, 462–475. [Google Scholar] [CrossRef]
  27. Marshall, A.W.; Olkin, I. A multivariate exponential distribution. J. Am. Stat. Assoc. 1967, 62, 30–44. [Google Scholar] [CrossRef]
  28. Guan, Q.; Tang, Y.C.; Xu, A.C. Objective Bayesian analysis for bivariate Marshall-Olkin exponential distribution. Comput. Stat. Data Anal. 2013, 64, 299–313. [Google Scholar] [CrossRef]
  29. Feizjavadian, S.H.; Hashemi, R. Analysis of dependent competing risks in the presence of progressive hybrid censoring using Marshall-Olkin bivariate Weibull distribution. Comput. Stat. Data Anal. 2015, 82, 19–34. [Google Scholar] [CrossRef]
  30. Paul, B.; Dey, A.K.; Kundu, D. Bayesian analysis of three parameter absolute continuous Marshall-Olkin bivariate Pareto distribution. Commun. Stat. Case Stud. Data Anal. Appl. 2018, 4, 57–68. [Google Scholar] [CrossRef]
  31. Kundu, D.; Gupta, A.K. Bayes estimation for the Marshall-Olkin bivariate Weibull distribution. Comput. Stat. Data Anal. 2013, 57, 271–281. [Google Scholar] [CrossRef]
  32. Mulinacci, S. Archimedean-based Marshall-Olkin distributions and related dependence structures. Methodol. Comput. Appl. Probab. 2018, 20, 205–236. [Google Scholar] [CrossRef] [Green Version]
  33. Shen, Y.; Xu, A. On the dependent competing risks using Marshall-Olkin bivariate Weibull model: Parameter estimation with different methods. Commun. Stat.-Theory Methods 2018, 47, 5558–5572. [Google Scholar] [CrossRef]
  34. Lawless, J.F. Statistical Models and Methods for Lifetime Data; Wiley: New York, NY, USA, 2003. [Google Scholar]
  35. Contreras-Reyes, J.E.; Gallardo, D.I.; Kharazmi, O. Time-dependent residual Fisher information and distance for some special continuous distributions. Commun. Stat.-Simul. Comput. 2022, 1–21. [Google Scholar] [CrossRef]
  36. Chaturvedi, A.; Pathak, A.; Kumar, N. Statistical inferences for the reliability functions in the proportional hazard rate models based on progressive type-II right censoring. J. Stat. Comput. Simul. 2019, 89, 2187–2217. [Google Scholar] [CrossRef]
  37. Bai, X.; Shi, Y.; Liu, Y.; Liu, B. Reliability inference of stress-strength model for the truncated proportional hazard rate distribution under progressively type-II censored samples. Appl. Math. Model. 2018, 65, 377–389. [Google Scholar] [CrossRef]
  38. Kizilaslan, F. E-Bayesian estimation for the proportional hazard rate model based on record values. Commun. Stat.-Simul. Comput. 2019, 48, 350–371. [Google Scholar] [CrossRef]
  39. Popovicć, B.V.; Ristić, M.M.; Genc, A.I. Dependence properties of multivariate distributions with proportional hazard rate marginals. Appl. Math. Model. 2020, 77, 182–198. [Google Scholar] [CrossRef]
  40. Han, M. The structure of hierarchical prior distribution and its applications. Oper. Res. Manag. Sci. 1997, 63, 31–40. [Google Scholar]
  41. Han, M. E-Bayesian estimation of the reliability derived from Binomial distribution. Appl. Math. Model. 2011, 35, 2419–2424. [Google Scholar] [CrossRef]
  42. Han, M. The E-Bayesian and hierarchical Bayesian estimations of Pareto distribution parameter under different loss functions. J. Stat. Comput. Simul. 2017, 87, 577–593. [Google Scholar] [CrossRef]
  43. Basheer, A.M.; Okasha, H.M.; EI-Baz, A.H.; Tarabia, A.M.K. E-Bayesian and hierarchical Bayesian estimations for the inverse Weibull distribution. Ann. Data Sci. 2021, 1–13. [Google Scholar] [CrossRef]
  44. Wu, M.; Shi, Y.; Wang, Y. E-Bayesian estimation for competing risk model under progressively hybrid censoring. J. Syst. Eng. Elect. 2016, 27, 936–944. [Google Scholar]
  45. Okasha, H.M.; Wang, J. E-Bayesian estimation for the geometric model based on record statistics. Appl. Math. Model. 2016, 40, 658–670. [Google Scholar] [CrossRef]
  46. Han, M. E-Bayesian estimation and its E-posterior risk of the exponential distribution parameter based on complete and type-I censored samples. Commun. Stat.-Theory Methods 2020, 49, 1858–1872. [Google Scholar] [CrossRef]
  47. Putter, H.; Fiocco, M.; Geskus, R.B. Tutorial in biostatistics: Competing risks and multistate models. Stat. Med. 2007, 26, 2389–2430. [Google Scholar] [CrossRef] [PubMed]
  48. Geskus, R.B. Data Analysis with Competing Risks and Intermediate States; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
Figure 1. Empirical cumulative distribution plot, Q-Q and P-P plots for the exponential, Pareto and Rayleigh distribution under AIDSSI data.
Figure 1. Empirical cumulative distribution plot, Q-Q and P-P plots for the exponential, Pareto and Rayleigh distribution under AIDSSI data.
Mathematics 11 01461 g001
Figure 2. The CDF and PDF plots of Rayleigh distribution for real data.
Figure 2. The CDF and PDF plots of Rayleigh distribution for real data.
Mathematics 11 01461 g002
Figure 3. Plots of θ 1 , θ 2 and θ 3 under different prior distributions.
Figure 3. Plots of θ 1 , θ 2 and θ 3 under different prior distributions.
Mathematics 11 01461 g003
Figure 4. Point estimate and interval length of HRF at different mission time under real data.
Figure 4. Point estimate and interval length of HRF at different mission time under real data.
Mathematics 11 01461 g004
Table 1. The AB and MSEs (within bracket) of point estimates with ( c 1 , c 2 , c 3 ) = ( 1.5 , 2 , 1.2 ) .
Table 1. The AB and MSEs (within bracket) of point estimates with ( c 1 , c 2 , c 3 ) = ( 1.5 , 2 , 1.2 ) .
nMLEBayesianE-Bayesian
π 1 π 2 π 3
θ 1 θ 3 θ 1 θ 3 θ 1 θ 3 θ 1 θ 3 θ 1 θ 3
ModelExponential Case
200.2200.2420.2190.2420.2030.2270.1970.2210.2100.234
[0.082][0.099][0.081][0.099][0.069][0.087][0.063][0.081][0.076][0.094]
350.1580.1740.1570.1740.1510.1680.1480.1660.1540.171
[0.040][0.051][0.041][0.050][0.037][0.047][0.035][0.045][0.039][0.049]
500.1350.1460.1350.1460.1310.1430.1290.1410.1330.145
[0.030][0.035][0.029][0.034][0.028][0.033][0.027][0.032][0.029][0.034]
ModelPareto Case
200.2200.2440.2190.2440.2030.2290.1970.2230.2100.236
[0.083][0.105][0.083][0.105][0.070][0.091][0.064][0.085][0.077][0.098]
350.1640.1820.1640.1810.1570.1750.1540.1720.1600.178
[0.044][0.054][0.044][0.053][0.040][0.050][0.038][0.048][0.042][0.052]
500.1280.1500.1270.1490.1240.1460.1230.1440.1260.148
[0.027][0.037][0.027][0.036][0.025][0.035][0.024][0.034][0.026][0.036]
ModelRayleigh Case
200.2160.2340.2150.2340.2000.2200.1930.2140.2060.227
[0.079][0.093][0.079][0.093][0.067][0.082][0.062][0.076][0.073][0.088]
350.1610.1750.1610.1740.1540.1680.1510.1660.1570.171
[0.042][0.050][0.042][0.050][0.038][0.047][0.037][0.045][0.041][0.049]
500.1330.1470.1330.1460.1290.1430.1270.1420.1310.145
[0.029][0.035][0.029][0.035][0.027][0.033][0.026][0.032][0.028][0.034]
Table 2. The ALs and CPs (within bracket) of interval estimates with ( c 1 , c 2 , c 3 ) = ( 1.5 , 2 , 1.2 ) .
Table 2. The ALs and CPs (within bracket) of interval estimates with ( c 1 , c 2 , c 3 ) = ( 1.5 , 2 , 1.2 ) .
nMLEBayesianE-Bayesian
π 1 π 2 π 3
θ 1 θ 3 θ 1 θ 3 θ 1 θ 3 θ 1 θ 3 θ 1 θ 3
ModelExponential Case
201.1672.8061.0711.3300.9711.1010.9471.0810.9951.122
[0.952][0.928][0.941][0.942][0.937][0.940][0.938][0.935][0.940][0.939]
351.0341.3990.8251.3170.7410.8390.7310.8300.7520.847
[0.958][0.951][0.952][0.966][0.945][0.949][0.942][0.948][0.944][0.950]
500.6751.0681.0211.2390.6240.7040.6190.6990.6310.709
[0.947][0.959][0.960][0.973][0.939][0.948][0.937][0.948][0.941][0.950]
ModelPareto Case
201.1682.7981.0771.3210.9711.0990.9481.0790.9961.119
[0.955][0.917][0.937][0.939][0.933][0.933][0.933][0.933][0.938][0.938]
350.8261.4661.0351.3170.7420.8390.7320.8310.7530.848
[0.946][0.944][0.950][0.959][0.938][0.943][0.938][0.943][0.937][0.946]
500.6741.0821.0161.2410.6230.7050.6170.7000.6300.710
[0.957][0.959][0.962][0.972][0.945][0.943][0.946][0.946][0.944][0.946]
ModelRayleigh Case
201.1602.9201.0651.3200.9651.0950.9411.0760.9901.116
[0.958][0.922][0.928][0.944][0.932][0.944][0.932][0.945][0.938][0.947]
350.8241.3891.0381.3080.7400.8330.7300.8250.7510.842
[0.952][0.951][0.942][0.963][0.933][0.946][0.933][0.946][0.937][0.946]
500.6731.0651.0171.2360.6230.7020.6170.6970.6300.707
[0.952][0.955][0.958][0.971][0.939][0.945][0.940][0.943][0.938][0.945]
Table 3. The AB and MSEs (within bracket) of point estimates with ( c 1 , c 2 , c 3 ) = ( 1.5 , 2 , 1.2 ) .
Table 3. The AB and MSEs (within bracket) of point estimates with ( c 1 , c 2 , c 3 ) = ( 1.5 , 2 , 1.2 ) .
nMLEBayesianE-Bayesian
π 1 π 2 π 3
S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 )
ModelExponential Case    t 0 = 0.5
200.0660.3860.0640.3860.0600.3570.0590.3370.0610.380
[0.007][0.263][0.006][0.262][0.006][0.220][0.005][0.190][0.006][0.256]
350.0490.2760.0480.2750.0460.2640.0460.2560.0470.273
[0.004][0.128][0.004][0.128][0.003][0.116][0.003][0.108][0.003][0.127]
500.0420.2320.0410.2320.0400.2250.0400.2200.0400.231
[0.003][0.090][0.003][0.090][0.002][0.084][0.002][0.080][0.003][0.090]
ModelPareto Case    t 0 = 2
200.0620.1960.0610.1960.0570.1810.0570.1720.0570.193
[0.006][0.070][0.006][0.070][0.005][0.058][0.005][0.050][0.005][0.068]
350.0460.1400.0460.1400.0440.1340.0440.1300.0440.139
[0.003][0.032][0.003][0.032][0.003][0.029][0.003][0.027][0.003][0.032]
500.0380.1140.0380.1140.0370.1110.0370.1080.0370.114
[0.002][0.022][0.002][0.022][0.002][0.021][0.002][0.019][0.002][0.022]
ModelRayleigh Case    t 0 = 1
200.0660.3810.0640.3810.0590.3530.0590.3330.0610.376
[0.007][0.257][0.006][0.257][0.005][0.215][0.005][0.186][0.006][0.250]
350.0490.2740.0480.2740.0460.2620.0460.2560.0470.271
[0.004][0.124][0.004][0.123][0.003][0.113][0.003][0.105][0.003][0.123]
500.0410.2310.0410.2300.0400.2240.0390.2190.0400.229
[0.003][0.088][0.003][0.088][0.002][0.082][0.002][0.078][0.003][0.087]
Table 4. The ALs and CPs (within bracket) of interval estimates with ( c 1 , c 2 , c 3 ) = ( 1.5 , 2 , 1.2 ) .
Table 4. The ALs and CPs (within bracket) of interval estimates with ( c 1 , c 2 , c 3 ) = ( 1.5 , 2 , 1.2 ) .
nMLEBayesianE-Bayesian
π 1 π 2 π 3
S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 )
ModelExponential Case    t 0 = 0.5
200.3261.9130.3021.8000.2951.7580.2951.7140.2941.794
[0.954][0.941][0.931][0.947][0.938][0.956][0.949][0.957][0.927][0.954]
350.2451.3870.2331.3360.2311.3160.2301.2980.2311.326
[0.952][0.946][0.938][0.947][0.947][0.954][0.950][0.955][0.941][0.951]
500.2051.1440.1971.1040.1961.1010.1961.0900.1961.112
[0.943][0.947][0.935][0.940][0.943][0.947][0.944][0.945][0.939][0.945]
ModelPareto Case     t 0 = 2
200.3110.9580.2800.8980.2730.8790.2760.8570.2700.901
[0.948][0.941][0.926][0.940][0.934][0.949][0.943][0.949][0.926][0.947]
350.2330.6960.2180.6660.2150.6600.2160.6510.2140.670
[0.957][0.952][0.944][0.948][0.947][0.957][0.951][0.958][0.938][0.956]
500.1940.5730.1840.5570.1830.5510.1840.5460.1820.553
[0.952][0.947][0.935][0.942][0.941][0.950][0.947][0.951][0.936][0.950]
ModelRayleigh Case    t 0 = 1
200.3261.9040.3021.7910.2961.7490.2951.7060.2951.784
[0.952][0.947][0.934][0.943][0.943][0.952][0.947][0.953][0.937][0.950]
350.2451.3810.2331.3290.2311.3110.2311.2940.2311.321
[0.953][0.954][0.943][0.948][0.949][0.954][0.952][0.951][0.947][0.956]
500.2051.1420.1971.1090.1961.0980.1961.0880.1961.101
[0.956][0.955][0.946][0.957][0.954][0.956][0.953][0.954][0.950][0.949]
Table 5. The PR and EPR of Bayes estimates for ( θ 1 , θ 3 ) with ( c 1 , c 2 , c 3 ) = ( 1.5 , 2 , 1.2 ) .
Table 5. The PR and EPR of Bayes estimates for ( θ 1 , θ 3 ) with ( c 1 , c 2 , c 3 ) = ( 1.5 , 2 , 1.2 ) .
nBayesianE-Bayesian
π 1 π 2 π 3
PR( θ 1 )PR( θ 3 )EPR( θ 1 )EPR( θ 3 )EPR( θ 1 )EPR( θ 3 )EPR( θ 1 )EPR( θ 3 )
ModelExponential Case
200.0790.0970.0720.0910.0680.0870.0760.095
350.0410.0520.0390.0500.0380.0490.0410.051
500.0280.0350.0270.0340.0270.0340.0280.035
ModelPareto Case
200.0790.0970.0720.0910.0680.0870.0760.095
350.0420.0520.0400.0500.0390.0490.0410.051
500.0280.0360.0270.0350.0270.0340.0280.035
ModelRayleigh Case
200.0780.0960.0710.0900.0670.0860.0750.093
350.0410.0510.0390.0490.0380.0480.0400.050
500.0280.0350.0270.0340.0270.0340.0280.035
Table 6. The AB and MSEs (within bracket) of point estimates with ( c 1 , c 2 , c 3 ) = ( 3 , 2.5 , 2 ) .
Table 6. The AB and MSEs (within bracket) of point estimates with ( c 1 , c 2 , c 3 ) = ( 3 , 2.5 , 2 ) .
nMLEBayesianE-Bayesian
π 1 π 2 π 3
θ 1 θ 3 θ 1 θ 3 θ 1 θ 3 θ 1 θ 3 θ 1 θ 3
ModelExponential Case
200.1180.1460.1180.1450.1070.1360.1030.1330.1120.141
[0.023][0.037][0.023][0.037][0.019][0.032][0.017][0.030][0.021][0.035]
350.0900.1080.0890.1080.0850.1040.0840.1030.0870.106
[0.013][0.020][0.013][0.019][0.012][0.018][0.011][0.017][0.012][0.019]
500.0740.0890.0740.0880.0720.0860.0710.0850.0730.087
[0.009][0.013][0.009][0.013][0.008][0.012][0.008][0.012][0.009][0.013]
ModelPareto Case
200.1170.1500.1170.1500.1070.1410.1030.1370.1110.145
[0.023][0.038][0.023][0.038][0.019][0.033][0.017][0.031][0.021][0.036]
350.0890.1080.0890.1080.0840.1050.0830.1030.0860.107
[0.013][0.020][0.013][0.019][0.012][0.018][0.011][0.017][0.012][0.019]
500.0720.0920.0720.0920.0700.0900.0690.0890.0710.091
[0.008][0.014][0.009][0.014][0.008][0.013][0.007][0.013][0.008][0.014]
ModelRayleigh Case
200.1190.1440.1180.1440.1080.1350.1040.1310.1130.139
[0.024][0.035][0.024][0.035][0.020][0.031][0.018][0.028][0.022][0.033]
350.0900.1110.0890.1110.0850.1070.0830.1060.0870.109
[0.013][0.021][0.013][0.021][0.012][0.019][0.011][0.019][0.013][0.020]
500.0730.0900.0730.0900.0700.0880.0690.0870.0710.089
[0.009][0.013][0.008][0.013][0.008][0.013][0.008][0.012][0.008][0.013]
Table 7. The ALs and CPs (within bracket) of interval estimates with ( c 1 , c 2 , c 3 ) = ( 3 , 2.5 , 2 ) .
Table 7. The ALs and CPs (within bracket) of interval estimates with ( c 1 , c 2 , c 3 ) = ( 3 , 2.5 , 2 ) .
nMLEBayesianE-Bayesian
π 1 π 2 π 3
θ 1 θ 3 θ 1 θ 3 θ 1 θ 3 θ 1 θ 3 θ 1 θ 3
ModelExponential Case
200.6592.4260.5520.7510.5270.6650.5090.6500.5450.680
[0.960][0.945][0.937][0.946][0.947][0.941][0.940][0.943][0.951][0.943]
350.4622.6690.4550.7720.4050.5100.3970.5040.4130.517
[0.950][0.965][0.936][0.959][0.930][0.946][0.932][0.946][0.936][0.947]
500.3780.8370.4790.7400.3440.4270.3390.4240.3490.431
[0.951][0.967][0.950][0.968][0.936][0.945][0.936][0.945][0.938][0.948]
ModelPareto Case
200.6612.3660.5530.7600.5290.6680.5110.6540.5470.683
[0.964][0.942][0.932][0.944][0.940][0.941][0.938][0.937][0.947][0.945]
350.4661.2520.4760.7670.4090.5100.4010.5040.4170.517
[0.955][0.964][0.952][0.958][0.945][0.947][0.945][0.945][0.947][0.947]
500.3770.8430.4750.7430.3430.4280.3380.4250.3480.432
[0.957][0.971][0.954][0.968][0.941][0.941][0.940][0.941][0.944][0.940]
ModelRayleigh Case
200.6592.4070.5520.7520.5280.6640.5100.6500.5450.679
[0.955][0.948][0.936][0.950][0.942][0.948][0.940][0.946][0.948][0.947]
350.4631.2710.4680.7760.4060.5120.3980.5060.4140.519
[0.952][0.963][0.937][0.955][0.931][0.939][0.927][0.937][0.934][0.937]
500.3780.8250.4780.7400.3440.4270.3400.4240.3490.431
[0.951][0.972][0.948][0.965][0.934][0.940][0.938][0.940][0.937][0.945]
Table 8. The AB and MSEs (within bracket) of point estimates with ( c 1 , c 2 , c 3 ) = ( 3 , 2.5 , 2 ) .
Table 8. The AB and MSEs (within bracket) of point estimates with ( c 1 , c 2 , c 3 ) = ( 3 , 2.5 , 2 ) .
nMLEBayesianE-Bayesian
π 1 π 2 π 3
S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 )
ModelExponential Case    t 0 = 1
200.0640.2370.0630.2370.0580.2190.0580.2080.0590.233
[0.006][0.102][0.006][0.102][0.005][0.085][0.005][0.074][0.005][0.099]
350.0480.1710.0470.1710.0450.1630.0450.1590.0460.169
[0.004][0.049][0.003][0.049][0.003][0.044][0.003][0.041][0.003][0.048]
500.0400.1410.0390.1410.0380.1370.0380.1340.0380.140
[0.003][0.033][0.002][0.033][0.002][0.031][0.002][0.030][0.002][0.033]
ModelPareto Case     t 0 = 5
200.0480.0480.0480.0480.0450.0440.0460.0420.0450.047
[0.004][0.004][0.003][0.004][0.003][0.003][0.004][0.003][0.003][0.004]
350.0350.0340.0350.0340.0340.0330.0340.0310.0340.034
[0.002][0.002][0.002][0.002][0.002][0.002][0.002][0.001][0.002][0.002]
500.0300.0280.0300.0290.0300.0280.0300.0270.0300.028
[0.001][0.001][0.001][0.001][0.001][0.001][0.001][0.001][0.001][0.001]
ModelRayleigh Case     t 0 = 1.5
200.0600.3490.0600.3490.0560.3230.0560.3070.0560.342
[0.006][0.212][0.005][0.213][0.005][0.178][0.005][0.156][0.005][0.207]
350.0470.2630.0460.2630.0450.2520.0450.2450.0450.261
[0.004][0.117][0.003][0.117][0.003][0.107][0.003][0.098][0.003][0.116]
500.0380.2100.0380.2090.0370.2030.0370.1990.0370.208
[0.002][0.074][0.002][0.073][0.002][0.069][0.002][0.065][0.002][0.073]
Table 9. The ALs and CPs (within bracket) of interval estimates with ( c 1 , c 2 , c 3 ) = ( 3 , 2.5 , 2 ) .
Table 9. The ALs and CPs (within bracket) of interval estimates with ( c 1 , c 2 , c 3 ) = ( 3 , 2.5 , 2 ) .
nMLEBayesianE-Bayesian
π 1 π 2 π 3
S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 ) S ( t 0 ) H ( t 0 )
ModelExponential Case    t 0 = 1
200.3211.1720.2931.1040.2861.0760.2871.0460.2831.101
[0.951][0.946][0.931][0.944][0.941][0.954][0.949][0.953][0.931][0.954]
350.2400.8550.2270.8240.2240.8110.2250.7990.2230.819
[0.953][0.946][0.940][0.940][0.940][0.953][0.948][0.955][0.936][0.950]
500.2000.7070.1910.6880.1900.6800.1910.6720.1900.682
[0.956][0.948][0.941][0.946][0.944][0.953][0.949][0.952][0.941][0.951]
ModelPareto Case     t 0 = 5
200.2590.2350.2220.2210.2150.2160.2190.2100.2090.221
[0.945][0.941][0.925][0.938][0.934][0.944][0.940][0.948][0.923][0.945]
350.1890.1720.1710.1660.1690.1630.1710.1600.1660.164
[0.956][0.949][0.940][0.954][0.941][0.958][0.948][0.956][0.934][0.954]
500.1560.1410.1450.1370.1440.1360.1450.1340.1420.136
[0.951][0.946][0.932][0.941][0.938][0.951][0.945][0.950][0.933][0.950]
ModelRayleigh Case     t 0 = 1.5
200.3131.7530.2821.6430.2751.6090.2771.5660.2711.652
[0.953][0.944][0.932][0.945][0.938][0.951][0.948][0.957][0.929][0.950]
350.2331.2860.2181.2390.2151.2190.2161.2020.2141.230
[0.949][0.943][0.933][0.946][0.935][0.944][0.943][0.951][0.931][0.941]
500.1941.0600.1841.0220.1831.0190.1841.0080.1821.031
[0.954][0.951][0.942][0.948][0.948][0.954][0.951][0.952][0.942][0.951]
Table 10. The PR and EPR of Bayes estimators for parameters ( θ 1 , θ 3 ) with ( c 1 , c 2 , c 3 ) = ( 3 , 2.5 , 2 ) .
Table 10. The PR and EPR of Bayes estimators for parameters ( θ 1 , θ 3 ) with ( c 1 , c 2 , c 3 ) = ( 3 , 2.5 , 2 ) .
nBayesianE-Bayesian
π 1 π 2 π 3
PR( θ 1 )PR( θ 3 )EPR( θ 1 )EPR( θ 3 )EPR( θ 1 )EPR( θ 3 )EPR( θ 1 )EPR( θ 3 )
ModelExponential Case
200.0240.0360.0220.0330.0200.0320.0230.035
350.0130.0190.0120.0190.0120.0180.0120.019
500.0090.0130.0080.0130.0080.0120.0090.013
ModelPareto Case
200.0240.0360.0220.0340.0200.0320.0230.035
350.0130.0190.0120.0190.0120.0180.0130.019
500.0090.0130.0080.0130.0080.0120.0090.013
ModelRayleigh Case
200.0240.0360.0220.0330.0200.0320.0230.035
350.0130.0190.0120.0190.0120.0180.0130.019
500.0090.0130.0080.0130.0080.0120.0090.013
Table 11. Goodness-of-fit results for real data.
Table 11. Goodness-of-fit results for real data.
Baseline ModelMLELog-LikelihoodKS Distancep-Value
Exponential θ = 0.144644−965.1150.19985< 0.05
Pareto θ = 0.585787−1066.5870.32379< 0.05
Rayleigh θ = 0.032023−899.5500.071860.1551
Table 12. MinRSSU samples of HIV infection competing risks data.
Table 12. MinRSSU samples of HIV infection competing risks data.
i12345678910111213
t i 5.1207.3021.3034.9665.7002.5652.1772.5134.2143.9532.3160.8241.462
δ i 2132221231322
i141516171819202122232425
t i 2.5330.9691.4401.4560.1122.3220.1482.8141.4370.8620.1370.142
δ i 121322323323
Table 13. Point estimates for unknown parameters under HIV infection competing risk data.
Table 13. Point estimates for unknown parameters under HIV infection competing risk data.
ParameterMLEBayesian E-Bayesian
Squared Error LossLinex Loss Squared Error LossLinex Loss
θ 1 0.006948170.008322310.00831653 π 1 0.007642450.00763715
π 2 0.007642280.00763700
π 3 0.007642630.00763729
θ 2 0.016675600.017227740.01721578 π 1 0.017369210.01735715
π 2 0.017368810.01735683
π 3 0.017369610.01735748
θ 3 0.011117070.011276600.01126878 π 1 0.011811060.01180287
π 2 0.011810790.01180264
π 3 0.011811340.01180309
Table 14. Interval estimates for unknown parameters under HIV infection competing risk data.
Table 14. Interval estimates for unknown parameters under HIV infection competing risk data.
ParameterACIBayesian-HPD E-Bayesian HPD
θ 1 (0.00289202,0.01669318)(0.00243423,0.01507565) π 1 (0.00201062,0.01368888)
[0.01380116][0.01264143] [0.01167826]
π 2 (0.00209560,0.01388392)
[0.01178832]
π 3 (0.00189844,0.01360728)
[0.01170883]
θ 2 (0.00947024,0.02936312)(0.00830230,0.02694340) π 1 (0.00824659,0.02670847)
[0.01989289][0.01864110] [0.01846188]
π 2 (0.008535855,0.02714980)
[0.01861395]
π 3 (0.00826488,0.02625957)
[0.01799469]
θ 3 (0.00366841,0.03369012)(0.00466253,0.01989804) π 1 (0.00418741,0.01908602)
[0.03002171][0.015235515] [0.01489861]
π 2 (0.00426602,0.01903927)
[0.01477325]
π 3 (0.00441722,0.01915059)
[0.01473338]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, Y.; Wang, L.; Tsai, T.-R.; Tripathi, Y.M. Estimation of Dependent Competing Risks Model with Baseline Proportional Hazards Models under Minimum Ranked Set Sampling. Mathematics 2023, 11, 1461. https://doi.org/10.3390/math11061461

AMA Style

Zhou Y, Wang L, Tsai T-R, Tripathi YM. Estimation of Dependent Competing Risks Model with Baseline Proportional Hazards Models under Minimum Ranked Set Sampling. Mathematics. 2023; 11(6):1461. https://doi.org/10.3390/math11061461

Chicago/Turabian Style

Zhou, Ying, Liang Wang, Tzong-Ru Tsai, and Yogesh Mani Tripathi. 2023. "Estimation of Dependent Competing Risks Model with Baseline Proportional Hazards Models under Minimum Ranked Set Sampling" Mathematics 11, no. 6: 1461. https://doi.org/10.3390/math11061461

APA Style

Zhou, Y., Wang, L., Tsai, T. -R., & Tripathi, Y. M. (2023). Estimation of Dependent Competing Risks Model with Baseline Proportional Hazards Models under Minimum Ranked Set Sampling. Mathematics, 11(6), 1461. https://doi.org/10.3390/math11061461

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