1. Introduction
The reliability function, denoted by
, is defined as the probability of a failure-free operation until time
t. Thus, if the random variable
X denotes the lifetime of an item or a system, then
. Another measure of reliability under a stress–strength setup is the probability
, which represents the reliability of an item or a system of random strength
X subject to random stress
Y, under a bivariate setting. These two reliability measures are frequently used in many applications. There is a substantial body of literature on the estimation and testing of the parameters
and
under progressive censoring. For instance, ref. [
1] proposed shrinkage estimators of
for the one-parameter exponential distribution. The authors of [
2,
3] estimated
under type-I and type-II censoring, while for estimating
, they used the complete sample case. This strategy’s exceptional benefits have attracted extensive academic research, as evidenced by studies such as [
4,
5,
6,
7,
8,
9,
10,
11].
There are many scenarios in life-testing and reliability experiments in which units are lost or removed from the experiment before failure. This removal may be unintentional, as in the case of an unexpected breakdown of an experimental unit, or it may be designed into the study. Unintentional loss may occur due to unforeseen circumstances such as lack of funding or lack of access to testing facilities. However, units are often deliberately removed from the test for reasons such as freeing up testing facilities for other experiments, saving time and money, or meeting ethical considerations in experiments involving human subjects. Among the various censoring schemes, the Type-II progressive censoring scheme has become very common in recent years because it allows the experimenter to remove active units during the experiment. The Maximum Likelihood Estimators (MLEs) are commonly used to estimate the reliability parameters under Type-II progressive censoring ([
12]). However, the efficiency of MLEs can be enhanced by integrating non-sample information, often available in the form of prior hypotheses regarding the parameter under consideration. Shrinkage estimation is one approach to amalgamate this non-sample information with existing estimators, and a comprehensive exploration of shrinkage and similar estimation methodologies is provided by [
13,
14].
Motivated by the potential improvements in estimation accuracy, our study focuses on devising shrinkage estimators for reliability and stress–strength parameters within a specific family of lifetime distributions. These parameters play crucial roles in assessing the performance and reliability of various systems, ranging from mechanical components to medical devices. Our research goes beyond mere theoretical formulation; we introduce novel shrinkage-type estimators that incorporate out-of-sample information, thereby offering a practical and robust solution to estimation challenges encountered in real-world scenarios. By addressing limitations inherent in traditional models, our proposed estimators exhibit superior performance, as demonstrated through rigorous simulation studies. Furthermore, our study underscores the practical relevance of these estimators by showcasing their efficacy in industrial contexts by applying them to real-world data. Through real-world applications, we illustrate how our shrinkage estimators can facilitate more accurate and reliable assessments of system reliability and stress–strength parameters.
The rest of this manuscript is organized as follows. In
Section 2, we give a brief introduction to existing results on the MLEs and then outline and develop the proposed improved estimation strategies for the reliability measures. The exact biases and efficiencies of the proposed estimators are also derived in this section.
Section 3 contains an extensive simulation study to evaluate the performance of the proposed estimators.
Section 4 is devoted to the application of the methods to real data sets. We give some concluding remarks in
Section 5.
3. Simulation Study
Here we conduct a Monte Carlo simulation study with a small sample size to assess the performance of the estimators proposed in this paper. The simulation setting and assumptions are as follows:
: The true value of reliability is taken to be at ;
m: number of observed failures is taken to be 10;
n: sample size was taken to be 100;
: obtained from (the prior guesses of );
: progressive Type-II censoring schemes which are taken to be
Scheme 1 (): ; Scheme 2 (): ; and Scheme 3 (): .
For each combination of
R and
, 500 Monte Carlo samples of size
were generated from the distribution given in (
1), taking
. The proposed estimators for
are calculated under progressive Type-II censoring and their CIs are computed according to Algorithms 1 and 2. Let
be a CI of
and
, observed values of lower and upper bounds of the proposed CI. The empirical expected length and coverage probability of the intervals are, respectively, computed as
Table 1 represents the coverage probability (CP) and expected length (EL) of the estimators for
under the three progressive Type-II censoring schemes (
) from an exponential model. From this table, we can see that the proposed estimators have, in most of the cases, higher coverage probabilities and shorter expected lengths than the usual estimators. This is more evident near the prior information
, as the S and PT estimators have CPs closer to the target 0.95 with shorter ELs. In general, the asymptotic intervals are conservative while the intervals based on bootstrap are liberal.
As for the estimation of , the simulation setup was as follows:
: the true value of is taken to be ;
: the prior guesses of are taken to be ;
: number of censored observations of X is taken to be 10;
: number of censored observations of Y is taken to be 8;
: progressive Type-II censoring schemes for X, assumed to be Scheme 1 (): ; Scheme 2 (): ; and Scheme 3 (): ;
: progressive Type-II censoring schemes of Y, assumed to be Scheme 1 (): ; Scheme 2(): ; and Scheme 3(): .
For each combination of
and
, 500 samples of size
were generated for
X from the distribution given in (
1), taking
and
and 500 samples of size
were generated for
Y from the same distribution with
and
. The proposed estimators for
are calculated under progressive Type-II censoring and their CIs are computed by using Algorithm 2.
The results of these simulations are presented in
Table 2.
The shrinkage-type estimators, PT and S, outperform the MLEs both in their CPs and Els and this is more pronounced than the case of the estimation of . The PT estimator has shorter expected lengths in the asymptotic CIs as compared to the bootstrap CIs. In general, for the S estimator, the asymptotic CIs are conservative while the bootstrap CIs are liberal. As we expected all confidence intervals are close to their theoretical coverage and are shortest near the true values of .