Abstract

We focus on estimating the stress-strength reliability model when the strength variable is subjected to the step-stress partially accelerated life test. Based on the assumption that both stress and strength random variables follow Weibull distribution with a common first shape parameter, the inferences for this reliability system are constructed. The maximum likelihood, two parametric bootstraps, and Bayes estimates are obtained. Moreover, approximate confidence intervals, asymptotic variance-covariance matrix, and highest posterior density credible intervals are derived. A simulation study and application to real-life data are conducted to compare the proposed estimation methods developed here and also check the accuracy of the results.

1. Introduction

Stress-strength models have attracted many statisticians for many years due to their applicability in different and diverse areas such as engineering, economics, and quality control, and, in the last years, there have been numerous applications to medical and engineering problems.

In the last ten years, many authors have been interested in studying the application of the simple stress-strength reliability model, which is more handled theoretically and at the same time is more simple and applicable to implement in practice. This model of reliability contains a strength variable and a stress variable , which is exposed to it. Such a system will properly function when exceeds ; namely, denotes the reliability system. Many estimation studies of reliability system are considered by several statistician researchers under both complete and censored samples from different models, for example, exponential distribution under progressive type-II censoring by Saraçoğlu et al. [1], Weibull distribution under complete samples by Kundu and Gupta [2], Kumaraswamy distribution under upper record values by Nadar and Kızılaslan [3], Kumaraswamy distribution under progressive type-II censoring by Nadar et al. [4], Lomax distribution under record values by Mahmoud et al. [5], Burr X distribution under complete samples by Surles and Padgett [6], inverse Lindley distribution under complete samples by Sharma et al. [7], exponential and Weibull by Kumar and Siju [8], Weibull-Gamma distribution under progressively type-II censored samples by Mahmoud et al. [9], general exponential form distribution under complete samples by Mokhlis et al. [10], Rayleigh distribution under complete samples by Afshin [11], modified Weibull model under progressively type-II censored samples by Soliman et al. [12], Lindley distribution using progressively first-failure censoring by Kumar et al. [13], generalized inverted exponential distribution under progressively first-failure censoring by Krishna et al. [14], Kumaraswamy exponential distribution under progressively type-II censored samples by El-Sagheer and Mansour [15], Burr XII distribution under progressively first-failure censored samples by Saini et al. [16], and generalized Maxwell failure distribution under progressive first-failure censoring by Saini et al. [17].

In previous reliability studies, it is evident that it is difficult to observe the lifetime of highly reliable components because few failures occur in a limited test time due to the very long lifetimes under normal test conditions. Therefore, to overcome this problem, we are looking for a catalyst for early failure of the components. Since testing under normal conditions takes a long time, then the development of accelerated life testing (ALT) or partially accelerated life test (PALT) is needed, where units are subjected to a more severe environment (increased or decreased stress levels) than the normal operating environment so that failures can be induced quickly. In this case, ALT or PALT allows experimenters to control higher stress levels to be used in the test.

In PALT, only part of the test components run under a higher stress level than the normal level, while all the test components run under a higher stress level in ALT. We use PALT when the acceleration factor is unknown, where items are examined at both normal and accelerated conditions. According to Nelson [18], there are three types of stress in PALT: constant stress, step stress, and progressive stress. In step-stress partially accelerated life test (S-SPALT), items are tested at a normal level; if it does not fail, then the stress will be changed at a certain time. This type allows the experimenter to select multiple stress factors, for instance, temperature, voltage stress, thermal and electrical cycling and shock, vibration, mechanical stress, and radiation in life testing.

Many authors have studied the inference based on the S-SPALT models for different probability distributions under various cases for censored or complete data, including Weibull by Zhang et al. [19] and Ismail [20], extended Weibull by Zhang and Shi [21], Burr type XII by Abd-Elfattah et al. [22], exponentiated exponential distribution by Abdel-Hamid and Al-Hussaini [23], Lomax by El-Sagheer and Ahsanullah [24], Gompertz by Ismail [25], modified Weibull by Mahmoud et al. [26], Kumaraswamy Inverse Weibull by El-Sagheer and Mohamed [27], and Weibull-Gamma by El-Sagheer et al. [28].

Recently, in parallel with progress in engineering, technology, and manufacturing, the experimenters may want to investigate the stress-strength reliability in case the strength component of the reliability system is exposed to an ALT. In this paper, we study a simple stress-strength model when the component strength exposes to S-SPALT. This system can be described as follows: such a system starts with the normal use condition of the strength variable and stress variable . If the system does not fail before the prespecified time , then the strength variable runs at an acceleration factor . This model will help us to evaluate when induced early failures to . Moreover, force to failure on strength may help us to see the effect of change on due to not only stress variable but also exposed stress by accelerating on strength variable . For this reason, we consider the S-SPALT model introduced by DeGroot and Goel [29] for strength variable .

The outline of the paper is as follows. In Section 2, assumptions of S-SPALT for the stress-strength reliability model are provided. Section 3 deals with the maximum likelihood estimate and asymptotic confidence intervals. Two parametric bootstrap methods are proposed in Section 4. In a Bayes paradigm, estimation techniques have been assayed in Section 5. In Section 6, a simulation study is conducted to compare the proposed procedures. In Section 7, a real-life data example is presented to illustrate the application of the proposed inference procedures. Finally, a conclusion is furnished in Section 8.

2. Assumptions of S-SPALT for Reliability System

Suppose that denote the lifetime of a test item as strength under S-SPALT can be determined, according to DeGroot and Goel [29], by the relationwith probability density function (PDF)where is the lifetime under normal use condition, is the time when stress is changed, and is the acceleration factor as . Suppose and are independent random variables following Weibull distribution (WD) with parameters and , respectively, that is, and , where the parameter is common and known, considering strength under S-SPALT with the PDF and CDF and primary stress with PDF and CDF . According to Çetinkaya [30], a partially accelerated life test implemented stress-strength reliability estimation can be written as

The PDF of S-SPALT implemented strength variable as suggested by DeGroot and Goel [29], which is given as follows:and CDF is given as follows:

Also, the PDF and CDF of primary stress are given by

Then, by using equations (4) and (6) in equation (3), the reliability of such a system can be obtained as

Then, if we put and , , the reliability of such a system, , devolves to one-parameter exponential distribution. If , equation (7) becomes the reliability for a simple stress-strength system without any acceleration. From Figures 1 and 2, we notice the following: (i) The reliability of the system increases with increasing stress change time , when the acceleration factor is fixed. (ii) Increasing on acceleration factor reduces the reliability quickly.

In equation (7), is common and known, is the predetermined stress change time, and , , and are unknown and need to be estimated; then

3. Maximum Likelihood Inference

Let be a random sample of strength from and be a random sample of stress from . Then, by considering equations (4) and (6), the likelihood function of observed samples in this reliability system is given by (see Çetinkaya [30])and equally

Hence, the logarithm of the likelihood function may then be written as

Taking the first partial derivatives of the log-likelihood in (11) with respect to , , and , we getwhere and . To get the MLEs of the unknown parameters, denoted by , , and , we should equate , , and to zero; thus,

Then, we use the Newton–Raphson iteration method to solve (15). Therefore, the MLE of , denoted by , can be obtained by considering the invariance property of the MLEs by replacing the parameters with their estimates as follows:where the parameter is common and known.

3.1. Asymptotic Confidence Interval

In this subsection, we construct an asymptotic confidence interval (ACI) for based on the asymptotic normal property of MLEs. Let be the MLEs of ; according to Cohen [31], the observed Fisher information matrix, denoted by , is defined bywhere

Also, the variance of is obtained by using the delta method as follows:where the first partial derivatives included in (20) can be easily obtained and is the element of the inverse of the information matrix as given bywhere , , , , and . Therefore, the ACI of is constructed aswhere is upper percentile of standard normal variate .

4. Parametric Bootstrap

In this section, we propose a resampling technique, the bootstrap procedure, to obtain a more widely used confidence interval. DiCiccio and Efron [32] introduced the bootstrap method and showed that the bootstrap method can improve the accuracy of the confidence intervals, especially when the sample is small such that the normal approximation is inappropriate. Besseris [33] showed that the bootstrap method can provide tighter confidence intervals. Reiser et al. [34] compared difference bootstrap confidence intervals by applying Monte Carlo simulation. Here, we study two bootstrap methods: bootstrap-p and bootstrap-t. These bootstrap confidence intervals work as follows.

4.1. Bootstrap-p

(1)Generate random samples from and from in (5) and (6), respectively. Calculate the MLEs of , , and .(2)Use , , and to generate independent bootstrap samples from and from . Calculate the MLEs of unknown parameters based on the bootstrap samples, denoted by , , and .(3)Calculate the bootstrap estimate of in (16), and denote by .(4)Repeat Steps 2 and 3 times; then we have .(5)Let be the CDF of . Define for given . Then, two-side percentile confidence intervals of are given by

4.2. Bootstrap-t

(1)The same as the bootstrap-p.(2)The same as the bootstrap-p.(3)The same as the bootstrap-p.(4)Obtain the -statistics , where given in (20).(5)Repeat Steps 2, 3, and 4 times; then we have .(6)Let be the CDF of . Define for given . Then, two-side bootstrap-t confidence intervals of are given by

5. Bayes Estimation

Bayes estimation is quite different from MLE and bootstrap methods because it takes into consideration both the information from observed sample data and the prior information. It can characterize the problems more rationally and reasonably. Assume that both parameters and have independent Gamma priors, while the parameter has usual noninformative prior; see Carlin and Louis [35]:

Here, , , , and are the hyperparameters that reflect the prior knowledge about the unknown parameters. The joint prior of the unknown parameters , , and is then given by

Via Bayes’ theorem, based on the considered joint prior (27) and the likelihood (10), the posterior distribution of , , and given data takes the form

It is clear that the conditional posterior densities of , , and can be written as

Thus, under the squared error loss function, the Bayes estimate of , denoted by , can be obtained as the mean of the posterior function as given in the following:

From (29) and (30), the full conditional posterior densities of and are and , respectively. Thus, the samples of and can be generated by using any Gamma routine. On the other hand, the expression of cannot be written as any well-known distribution. One can use the method proposed by Devroye [36] to generate sample data from this distribution.

However, the Metropolis–Hastings (MH) with the Gibbs sampling scheme by using normal proposal can be effectively used to simulate random samples from (25)–(27). The MH algorithm and Gibbs sampler work as follows:(1)Use the MLEs as the initial value, denoted by (2)Set (3)Generate from (4)Generate from (5)Using MH algorithm,(i)Generate from the proposal normal distribution .(ii)Evaluate the acceptance probabilities(iii)Generate from Uniform distribution.(iv)If , accept the proposal and set ; else set .(6)Compute at , , and (7)Set (8)Repeat Steps times and obtain ,

Then, the Bayesian estimators of under the squared error loss function are given bywhere is burn-in to guarantee the convergence and to remove the affection of the selection of initial values. Therefore, the highest posterior density (HPD) Bayes credible interval is given by

6. Simulation Study

In this section, we apply a Monte Carlo simulation to assess the performance of MLEs, bootstrap, and Bayes estimator methods for the stress-strength reliability model with component strength under PALT, along with their ACIs, bootstrap CIs, and HPD credible intervals. Furthermore, we study the variations on reliability with the different cases for both acceleration factor and stress change time . The performance of estimators is evaluated in terms of mean square error (MSE) for the point estimates, also coverage probability (CP), and average lengths (ALs) for interval estimates (asymptotic, bootstrap, and HPD). We consider five sample sizes such as , , , , and for eight cases of the true values of the parameters, stress change times, acceleration factor, and corresponding actual values of , when the common parameter . These cases are as follows:

Case 1. , , , , and .

Case 2. , , , , and .

Case 3. , , , , and .

Case 4. , , , , and .

Case 5. , , , , and .

Case 6. , , , , and .

Case 7. , , , , and .

Case 8. , , , , and .
The first four cases provide us with the upper (around 0.81477) actual values of and their outcomes presented in Tables 14, while the second four cases provide us with the lower (around 0.3395) actual values of and their results are given in Tables 58. The study is performed for replicates. For each replication, 1000 bootstrap samples are used. The Bayes estimates and the credible intervals are computed based on 12000 MCMC samples and discard the first values 2000 as “burn-in.” For Bayesian analysis, we consider informative prior with hyperparameters values , , , and . The results of this study are reported in Tables 18. In the first rows of these tables, the estimates and their ACIs are given, respectively, for all sample sizes . Besides, the MSEs, ALs for ACIs, and their corresponding CPs are reported in the second rows in all sample sizes.

7. Application to Real-Life Data

In this section, for illustrative purposes, real-life data are presented to inspect the inference procedures. We used two sets of data introduced by Badar and Priest [37] and used by Kundu and Gupta [2]. Set 1, denoted by given in Table 9, is strength measured in GPA for single carbon fibers tested under tension at gauge lengths of 20 mm. Set 2, denoted by given in Table 10, is strength measured in GPA for single carbon fibers tested under tension at gauge lengths of 10 mm. Both Set 1 and Set 2 are of size 63. For the purpose of the goodness of fit test, the Kolmogorov–Smirnov (KS) distance between the empirical and the fitted distribution functions has been computed. It is 0.087588 and the associated value is 0.7192 for Set 1, while it is 0.050055 and the associated value is 0.9975 for Set 2. Hence, the value for KS has the highest value for Sets 1 and 2. This leads us to conclude that WD is the best fit for the two real data sets. Empirical, , and plots are shown in Figures 3 and 4, which make it clear that the WD fits the data very well.

Under the S-SPALT implementation on the strength variable , the corresponding estimates of parameters change depending on the stress change time . For example, at and , MLEs of , , and . Thus, using (16), the MLE of and its corresponding ACI is calculated as and with length 0.18936, respectively. Also, by implementing the iterative algorithms mentioned in Section 4, the boot-p and boot-t CIs are computed as and with lengths 0.18536 and 0.18026.

Now, to compute the Bayesian estimate of , the prior distributions of the parameters and are needed to specify. Because we have no prior information, we used the noninformative Gamma prior for and that is, the hyper-parameters are close to or equal to zero and , . Under the MCMC technique, the posterior analysis was done across combining MH algorithm within Gibbs sampler. To conduct the MCMC algorithm, which was described in Section 5, the initial values for the parameters , , and were taken to be their MLEs. In addition, 12000 MCMC samples were generated. To avoid the effect of the initial values, we expunge the first samples as “burn-in.” Thus, the Bayesian estimate of and its corresponding CRI is calculated as and with length 0.17441, respectively. Figure 5 displays chain values for . The histogram and the kernel density estimate of are shown in Figure 6. Also, we obtained more results with different values of as shown in Table 11.

8. Conclusions

In this paper, we consider the estimation of the stress-strength reliability model when the strength variable is subjected to the S-SPALT. By exposing the strength variable to the acceleration factor, a normal stress-strength model is forced to early failure; this design of the reliability has an effect on the level of reliability. We can illustrate it as follows: (i) In a specific stress change time, the increasing on acceleration loadings decreases the stress-strength reliability. (ii) Delaying stress change time helps maintain a level of reliability.

The main aim of this paper is to study the effect of external stress loading on the strength variable. This external stress includes, for instance, temperature, voltage stress, thermal and electrical cycling and shock, vibration, mechanical stress, and radiation. Thus, the estimations of the stress-strength reliability model with the corresponding ACIs using the maximum likelihood, two parametric bootstrap, and Bayesian estimation methods are obtained. For illustrative purposes, we have applied a real-life example. A simulation study is computerized to inspect and compare the rendition of the proposed methods for different sample sizes , different acceleration factor , and different stress change time . From the results, we observe the following:(1)It is clear that, from Tables 18, as sample sizes increase, the MSEs and average interval lengths decrease.(2)MLE, bootstrap, and Bayesian methods have very close estimates, and their ACIs have quite high CPs (around 0.95).(3)The ALs are decreasing with parallel to increasing on the actual value of .(4)Bayes estimates perform better than the MLEs in terms of MSEs.(5)The ACIs of the MLEs have the smallest ALs, and the ACIs of the bootstrap-t have the largest ALs. At the same time, the Bayesian HPD intervals have the highest CPs.(6)Finally, we can conclude that the proposed inference methods give consistent results.(7)Sometimes, it is worth noting that the available data may be affected by uncertainties and/or inaccuracies. Then, strictly speaking, it would be necessary to carry out a fuzzy preprocessing of the data; see [38] as a future work.

Data Availability

The article contains data generated from the statistical models used in the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the Taif University Researchers Supporting Project (no. TURSP-2020/318), Taif University, Taif, Saudi Arabia.