1. Introduction
Reliability estimation is a popular research topic that has recently received much attention. Many studies focus on the stress–strength system with two components, which can be regarded as two random variables
X and
Y, representing strength and stress. Its reliability is
. For this kind of system, extensive research has been carried out. Many scholars have discussed the estimation of this reliability under certain distributions, such as Weerahandi and Johnson [
1], Badr et al. [
2], Xu et al. [
3], Hassan et al. [
4], etc.
In the theory of reliability, the multicomponent stress–strength system is an extension of the classical stress–strength system, which contains one stress component and k independent strength components. The system will remain stable if at least s of the k strength values exceed the stress value. For instance, a bridge with k vertical cables that represent the strength of the structure is just a multicomponent stress–strength system. It will remain stable if the stress brought on by wind, high traffic volume, etc., does not exceed the values of at least s of its k strength components. Another example is an eight-cylinder vehicle engine, which operates as long as at least four of the eight cylinders are firing.
For
k strengths that are exposed to one stress, strengths
are independent random variables with the same cumulative distribution function (CDF)
, and stress
Y (independent of
) is a random variable with the CDF
. Then, the multicomponent stress–strength reliability (MSSR) can be described as
The estimation of the MSSR has been widely studied under complete sample data, where different distributions have been used. For example, Rao [
5] studied the maximum likelihood estimation and asymptotic confidence interval of the MSSR based on generalized exponential distribution. Jia et al. [
6] discussed the classical estimation of the MSSR based on generalized inverted exponential distribution. Nadar and Kizilaslan [
7] derived the estimation of the MSSR based on the Marshall–Olkin Bivariate Weibull distribution using frequentist and Bayesian methods. The Bayesian estimators were developed using Lindley’s approximation and Markov Chain Monte Carlo techniques. Kizilaslan and Nadar [
8] discussed the estimation of the MSSR based on the Bivariate Kumaraswamy distribution in a similar way.
Censoring is an effective tool that is often used in lifetime tests. In practical experiments, censoring plays an important role when there are not enough test units or when data for all test units cannot be collected because of the lack of time and resources. Type-I and Type-II censoring attract a lot of attention due to their mathematical simplicity. In Type-I censoring, if a pre-determined time is reached, the test will be stopped, whereas in Type-II censoring, if a pre-determined number of units fail, the test will be stopped. However, both censoring schemes may be inappropriate if the experimenter needs to intermittently remove units. Thus, progressive Type-II censoring is considered a better option and has been widely used in recent years. In this censoring, the intermittent removal of units is allowed. In addition, it saves time and expenses to a certain extent.
Progressive Type-II censoring in a lifetime test is as follows. Assume that
M units are tested in the experiment. When the first failure
happens,
units are randomly removed. When the second failure
happens,
units are randomly removed, and so on. Finally, when the
failure
happens, the remaining
active units are all removed. Note that
. In this way, the ordered lifetime data for
m elements are obtained using the censoring scheme
. This process is shown in
Figure 1. In particular, progressive Type-II censoring can be seen as an extension of Type-II censoring, which is derived by assuming
. Otherwise, the complete sample is derived by assuming
.
For samples that have been censored or are incomplete, there are few studies on the inference of the MSSR in the literature. Saini et al. [
9] studied the estimation of the MSSR for Burr XII distribution using progressive first-failure censoring data. Tsai et al. [
10] discussed the estimation of the MSSR for generalized exponential distribution using Type-I censoring data. Saini et al. [
11] obtained the estimation of the MSSR for Topp–Leone distribution using progressive censoring data.
Many probability distributions have been proposed to fit the lifetime data. Here, we select the Chen distribution proposed by Chen [
12]. Its probability density function (PDF), CDF, and hazard rate function (HRF) are, respectively, given by
where
and the shape parameters
. Hereafter, the Chen distribution can be represented by Chen
. Images of the PDFs and HRFs of the Chen distribution for certain cases are presented in
Figure 2. It can be seen that the Chen distribution has a flexible PDF and HRF when the shape parameters are taken at different values. When
, the PDFs have a significant peak; when
, the PDFs have a roughly decreasing trend; and when
, the PDFs are monotonically decreasing. Additionally, when
, the HRFs are monotonically increasing, and when
, the HRFs present as a bathtub-shaped curve. In the literature, the Chen distribution has received considerable attention and has been studied by several scholars, including Chen and Gui [
13], Rastogi et al. [
14], Sarhan et al. [
15], Mendez-Gonzalez et al. [
16], and the references therein.
Let
follow
independently, and let
Y (independent of
) follow
. They have a common second shape parameter and different first shape parameters. Now, based on Equations (
1)–(
3),
can be derived as
where
.
As far as we know, no research has been conducted on the inference of the MSSR for the Chen distribution using progressively censored data. Wang et al. [
17] discussed its classical estimation using Type-II censored data. So, our goal is to expand their research by carrying out its estimation from the Chen distribution using progressively censored data. In addition, the Bayesian estimation is also considered.
This article is organized as follows. In
Section 2, the maximum likelihood estimate (MLE) is studied. Then, the interval estimation is discussed in
Section 3. The asymptotic confidence interval (ACI) is derived based on the MLE. Additionally, the Bootstrap confidence interval (BootCI) is constructed. In
Section 4, the Bayes estimates based on gamma priors are derived using the Markov Chain Monte Carlo (MCMC) method. The Bayesian credible interval (BCI) and the highest posterior density credible interval (HPDCI) are also discussed. In
Section 5, Monte Carlo simulation studies are performed and a real data set is analyzed. Finally, the conclusions are reported in
Section 6.
2. Maximum Likelihood Estimation
In this section, we derive the MLE of
by adapting the approach of Wang et al. [
17] and extending its application to progressively Type-II censored data. Assume that
N systems with
K strength components are subjected to a lifetime test. We observe the lifetime data for
n systems with
k components using progressive Type-II censoring. The observed
and
are described as follows:
Observed strength values Observed stress values
where each row of values in
is the censored sample from
with the progressive Type-II censoring scheme
, and
is the censored sample from
with the progressive Type-II censoring scheme
. For simplicity, we denote the censoring scheme as
and
. Then, the likelihood function is
where
Now, using Equations (
2), (
3), and (
5), the likelihood function is
From Equation (
6), the log-likelihood function is
By taking the derivative of
with respect to
and
and setting them to 0, we obtain the following equations:
Using Equations (
8) and (
9), we derive
Theorem 1. For a given , the MLEs of and exist and are, respectively, given by and in Equation (10). Proof. Let
,
. Using inequality
, we obtain
By substituting Equations (
11) and (
12) into (
7), one obtains
From Equation (
10), one obtains
The equality holds if , where and . □
From Theorem 1, by substituting
and
into Equation (
7), the log-likelihood function of
is
By taking the derivative of
and setting it to 0, we obtain the following equation and its solution is the MLE
.
From Equation (
14), we can derive a nonlinear equation
, where
The above equation has a fixed-point solution for
. The MLE
can be derived using the fixed-point iterative approach as
, where
is the
iterative value of
. When
is very close to 0, the iteration process can be stopped. Then, according to Equation (
11), the MLEs
and
can be obtained using the MLE
.
Based on the invariance property of the MLEs, once we obtain the MLEs
,
, and
, the MLE
is
6. Conclusions
We study the estimation of the MSSR based on the Chen distribution using progressively censored data. The reliability inference of multicomponent stress–strength systems has attracted significant interest, with numerous scholars contributing to the field. Progressive Type-II censoring has been widely used in lifetime tests for nearly two decades. The Chen distribution is commonly used to model real-life data in the fields of lifetime analysis and reliability theory.
We begin by obtaining the MLE. Then, the ACI based on the MLE is derived, where the delta method is used. The percentile BootCI is also constructed. Additionally, the Bayes estimates, BCIs, and HPDCIs are obtained using the MCMC method. The stochastic simulations are performed using the R software. According to the results, the Bayes estimates under the gamma priors using the precautionary loss function exhibit the smallest MADs and MSEs among all the point estimates. Among all the interval estimates, the BCIs have the shortest lengths and the HPDCIs have the highest coverage probabilities. Finally, the applicability of the method is illustrated using real data.
Our study focuses solely on the estimation of the MSSR based on the Chen distribution using a common second shape parameter. In the future, we will explore the case of unequal shape parameters. Additionally, the multicomponent stress–strength system we study contains only one stress component and we do not consider scenarios involving multiple stress components. In the future, we will explore a more complex system and further extend the model to cases where there is more than one stress component. Furthermore, we assume that the strength random variables are independently and identically distributed, which may not accurately reflect certain real-world situations. Therefore, the study of multicomponent stress–strength systems with non-identical strength variables will be considered in the future.