Abstract

In this work, we show how to estimate stress strength (SS) reliability when the stress (Y) and strength (X) distributions are generalized exponentials with a common scale parameter. The SS reliability estimator is considered in view of neoteric ranked set sampling (NRSS) and median ranked set sampling (MRRS). We acquire an estimate of the reliability (R) when such samples of the stress and strength random variables are gathered using the same NRSS technique. Furthermore, the reliability estimator is derived when the stress distribution data are in the pattern of MRSS with just an odd/even set size and the strength distribution data are derived from NRSS and vice versa. The simulation results are used to evaluate and understand the adequacy of a variety of estimators for the suggested schemes. Based on our simulated results, we found that NRSS-based stress strength reliability estimates are more efficient than MRSS-based stress strength reliability estimates. The analysis of real-world data is used to implement the recommended estimators.

1. Introduction

In the literature, the issue of inferring the SS model, represented by R = P (Y < X), has gotten a lot of attention. Engineering, statistics, and biostatistics are just a few of the disciplines where it is employed. When the stress surpasses the system’s strength, it is self-evident that the system breaks down; otherwise, it will operate normally. The SS reliability (SSR) is calculated in useful applications such as rocket-motor constructions and deterioration, fatigue failure of aircraft structures, and ageing concrete pressure vessels. The SS model was firstly proposed by Birnbaum [1], in which both X and Y follow specific distributions based on full samples, and several researchers have explored the estimate of the SSR model. The reader can refer to Sathe and Shah [2], Constantine et al. [3], Weerahandi and Johnson [4], Grimshaw [5], Saraçoyğlu and Kaya [6], and Krishnamoorthy and Lin [7].

To estimate the population mean, McIntyre [8] introduced the ranked set sampling (RSS) strategy for selecting a more informative sample than a simple random sample (SRS). A small number of units can be quickly ranked, either by eye in relation to the research variable or on the basis of an arbitrary function. Once precise sample measurement for a particular unit is complicated, costly, or time-consuming, RSS can indeed be employed, such as environmental management, ecology, sociology, and agriculture. The following is the fundamental principle of selecting a sample in RSS:(i)Divide the target population into n groups of size by allocating randomly chosen units.(ii)Despite knowing the values for the variable of interest, rank the units inside every set in order of importance. This could be influenced by personal professional evaluation or a variable that is connected to the variable of interest.(iii)Take the least ranked unit from the first set, the second smallest ranked unit from the second set, and so on until the biggest ranked unit from the last set is picked for real quantification.(iv)Repeat this technique p times to get a sample of size for the actual measurement.

In the literature, several researchers use RSS to estimate the parameters and SSR of particular distributions. Using RSS, Bhoj [9] estimated the parameters of the extreme value distribution. Fei et al. [10] examined the performance of RSS in the estimate process for more information about RSS. Under unbalanced RSS, Kvam and Samaniego [11] presented an estimate of the population mean and population distribution function. Lam et al. [12] utilized RSS to estimate the exponential distribution’s two parameters. Under various loss functions, Ali [13] provided the Bayesian estimators for the mean residual life function, Lorenz curve, and SSR. Hassan [14] used the RSS technique to explore maximum-likelihood and Bayesian estimators for the generalized exponential distribution (GED). The SSR estimates for Burr XII under some RSS designs were considered by Hassan et al. [15, 16], and Bantan et al. [17] handled the Zubair Lomax parameter estimators under RSS. Safariyan et al. [18] handled the point and interval SSR estimators using RSS. Al-Omari et al. [19] and Akgul et al. [20] derived, respectively, SSR estimators of exponentiated Pareto and generalized inverse Lindley distributions. Almarashi et al. [21] investigated the estimation of R = P[Y < X], in the case where Y and X are two independent Topp–Leone random variables. Bantan et al. [22] introduced and discussed parameter estimators of half logistic inverted Topp–Leone distribution from RSS scheme. Esemen et al. [23] proposed the SSR for the GED using maximum-likelihood (ML) and Bayesian methods via SRS and RSS designs.

Several RSS modifications have been proposed in the literature; two important ranking systems are discussed in the following subsections.

1.1. NRSS Design

The NRSS system, which varies from the original RSS method, was presented by Zamanzade and Al-Omari [24]. If compared with SRS and RSS methods, this method has been demonstrated to give more accurate estimators for the population mean and variance. It is employed in circumstances where sorting sample observations are considerably easier than getting their precise values. The NRSS system is described in the following way:(i)Assign n2 randomly selected units from the target population and rank the sample units according to the predetermined ordering criteria.(ii)Select as the ranking unit for i = 1, …, n if n is odd. If n is even, choose the ranked unit, where if i is an even and if i is an odd for i = 1, …, n.(iii)Repeat this process p times to investigate a final sample of size N = np.

Taconeli and Cabral [25] developed several two-stage sampling strategies based on NRSS and evaluated their performance. Koyuncu and Karagoz [26] used the NRSS architecture in statistical quality control to develop robust control charts. The inverse Weibull distribution parameters were estimated using the NRSS method by Sabry and Shaaban [27].

1.2. MRSS Design

The MRSS was proposed by Muttlak [28] as a modification of the RSS to reduce RSS lost efficiency due to ranking errors and to improve the effectiveness of the population mean estimator. The approach is described in the following manner:(i)Distribute the randomly chosen units from the target population into sets, each of size .(ii)Despite knowing the values for the variable of interest, rank the units within each set in order of importance. It might have been defined as a subjective professional evaluation or a variable that is connected to the variable of interest.(iii)If the set size is odd, use , the lowest rank unit from each sample for real measurements.(iv)If the number of pieces in the set is even, choose the smallest ranked unit from the first samples, the for the measurement, and the smallest ranked unit from the second samples, the for the measurement.(v)To acquire a final sample of size , repeat this method p times.

The asymptotic distribution of the parameter estimators of the simple linear regression model using MRSS was studied by Alodat et al. [29]. Al-Nasser and Gogah [30] presented a single sampling approach for ranking data under the assumption of a GED. Using one and two auxiliary variables, Koyuncu [31] developed novel regression estimators in MRSS and NRSS schemes. Hassan et al. [32] discussed the SSR when X and Y are generalized inverted exponential distributions.

To date, no research has been published that deal with estimating SSR using the NRSS technique. Here, we address SSR-based estimation when the stress Y and the strength X are both independent generalized exponential random variables using NRSS and MRSS schemes. As a result, we do the following:(i)A comparison of stress strength reliability estimates under the NRSS and MRSS systems is taken into account.(ii)When the stress and strength random variables are chosen from the same sample design (NRSS/MRSS), this problem is examined.(iii)When the stress and strength random variables are chosen from various sample designs, this issue is taken into account.(iv)A simulation study is conducted to validate the theoretical study and two applications to real data are provided.

The following sections are arranged in the following order: Section 2 discusses the significance of GED and its reliability parameter. The SSR estimator utilizing NRSS is shown in Section 3. When the observed data for X are obtained from MRSS and the observed data for Y are received from NRSS and the opposite, Section 4 produces the SSR estimator. Sections 5 and 6 give the experimental results and data analysis, respectively. Finally, concluding remarks are handled in Section 7.

2. Generalized Exponential Model

Since 1995, exponentiated distributions have been extensively researched in statistics, with a number of authors developing distinct classes of these distributions. Gupta and Kundu [33] proposed the GED, which is a very appealing extension of the exponential distribution. In the literature, the GED is also known as the exponentiated exponential distribution. It is often used in evaluating lifetime or survival statistics and is of considerable interest. The cumulative distribution function (CDF) and the probability density function of the GED are given, respectively, as follows:where and are shape and scale parameters, respectively. Many authors have looked into the properties of the GED; for example, see Gupta and Kundu [34], Raqab [35], Zheng and Park [36], Shirke et al. [37], Srivastava et al. [38], and Gupta and Kundu [39].

The GED has a wide range of applications. The analysis of Los Angeles rainfall data has been studied by Madi and Raqab [40]. Subburaj et al. [41] suggested software reliability growth models for critical quality measures. Sarhan [42] published GED parameter estimates for a competing risk model based on censored and incomplete data. In addition, the censoring times of 36 small electrical equipments subjected to an automated life test were considered. Biondi et al. [43] investigated models for episode peak and length for eco-hydro-climatic applications. Cota-Felix et al. [44] used the GED to estimate the mean life of power system equipment with limited end-of-life failure data. Kannan et al. [45] investigated the performance of the cure rate model based on GED. They used data from a study on drug treatment programs to demonstrate the model’s performance. Rao [46] suggested a multicomponent SSR estimator of GED with an application to data from an airplane’s air conditioning system. Sadeghpour et al. [47] proposed an SSR estimator based on a lower record ranked set sampling scheme from the GED, and analysis to rocket-motor data was provided.

Let X be an GED with parameters and Y be another GED with parameters , where X and Y are independent. The SSR is calculated assuming that the scale parameter in both models is the same but with different shape parameters and ; that is, X ∼ GED and Y ∼ GED . Therefore, the SSR is given as follows:

Expression (3) is a function of parameters and .

3. Reliability Estimator Based on NRSS

In this section, we obtain the SSR when both X and Y are independent GEDs based on NRSS. Suppose that , is selected NRSS from GED with sample size N = n. Let , be the chosen NRSS from GED with sample size M = m. The observed NRSS design (e.g., n = 3) is illustrated in Table 1.

The likelihood function (LF) of the observed samples from strength X via NRSS is given as follows:where (0) = 0, (n + 1) =  + 1,  = n2, and Similarly, the LF of the observed samples from stress Y via NRSS is as follows:where (0) = 0, (m + 1) =  + 1,  = m2, and . The log LF based on NRSS is offered as follows:

The partial derivatives with respect to , and are computed as follows:where for simplified we write and

To produce ML estimators of the population parameters, equations (7) and (8), and (9) are quantitatively determined using an iterative approach. Consequently, the SSR estimator is obtained using (3).

4. Reliability Estimator Using NRSS versus MRSS

Herein, we estimate the SSR where both X and Y are independent GEDs viz. MRSS and NRSS. Four different estimators are derived: the first and second SSR estimators are derived when the collected samples of stress random variable are generated from NRSS and strength random variable is generated from MRSS for odd set size (MRSSO) and vice versa. The third and fourth SSR estimators are calculated where the collected data of stress are taken from NRSS and strength data are drawn from MRSS for even set size (MRSSE) and conversely.

4.1. ML Estimators of

Suppose that , is a MRSSO observed from GED distribution, with sample size . The observed MRSSO is illustrated in Table 2.

Let be a NRSS observed from GED distribution, with sample size M = m. The LF of observed data would be as described as follows:where The log LF of the GED and GED distributions based on NRSS and MRSSO schemes is as follows:

The partial derivative of is derived in (9). The partial derivatives of and are as follows:

The ML estimators of , and are derived by equating (9), (12), and (13) by zero. Despite the fact that the suggested estimators cannot be stated in closed forms, they may be easily produced using a numerical method. As a result, we can obtain the SSR estimator from (3) using the invariance property.

4.2. ML Estimator of

Suppose that , is a NRSS observed from GED , with sample size N = npx. Let be a MRSSO observed from GED , with size . The log LF of observed data, say L3, would be as described as follows:where The log LF of the X and Y distributions that depend on NRSS and MRSSO schemes is as follows:

The partial derivative with respect to is derived in (8). The partial derivatives with respect to and are as follows:

Equating (8), (16), and (17) by zero, the ML estimators of , and are obtained. Although the suggested estimators are not in closed forms, they may be easily produced using a numerical approach. Hence, we obtain the SSR estimator by inserting the ML estimators of and in (3).

4.3. ML Estimator of

Suppose that , is a MRSSE observed from GED , with sample size . The observed MRSSE is illustrated in Table 3.

Let , be a NRSS observed from GED distribution, with sample size M = m The LF of sample, denoted by L4, would be as described as follows:where The log LF of the GED and GED based on MRSSE and NRSS schemes is as follows:

The partial derivative of is derived in (9). The partial derivatives of and are as follows:

Equating (9), (20), and (21) by zero, the ML estimators of , and are obtained. Despite the fact that the suggested estimators cannot be stated in closed forms, they may be easily produced using a numerical method. Consequently, we can get the reliability estimator using (3).

4.4. ML Estimator of

Suppose that , is a NRSS observed from GED , with sample size N = n. Let , be a MRSSE observed from GED , with sample size . The observed LF, say L5, would be as described as follows:where The log LF of X and Y distributions based on NRSS and MRSSE is as follows:

The partial derivative of is given in (8). The partial derivatives of and are as follows:

After setting (8), (24), and (25) to zero and solving them numerically, we get the ML estimators of , and It is straightforward to create a reliability estimator of R using the invariance property.

5. Numerical Study

We devised a simulation scenario to investigate and compare the behavior of SSR estimates derived from various sampling techniques. Measures such as the absolute bias (), the mean squared error (MSE), and the relative efficiency (RE) criteria are utilized to check the validity of estimates. The set sizes of X and Y are taken as (n, m) = () = (2, 2), (4, 4), (6, 6), (3, 3), (5, 5), and (7, 7), and number of cycles are chosen as  =  = 5. Consequently, the sizes of the samples are  =  = (10, 10), (20, 20), (30, 30), (15, 15), (25, 25), and (35, 35). The values of parameters are selected as  = (0.5, 0.4), (4, 0.5), and (5, 0.4) and  = 1, and therefore, the exact value of R is equal to 0.5556, 0.8889, and 0.925. 1000 random samples from X ∼ GED and Y ∼ GED are generated. The relative efficiency of SSR estimates of the RSS with respect to NRSS is defined as follows:

The values of ABs, MSEs, and REs of the SSR estimate are listed in Tables 4 and 5 for specific values of (n, m), , and distribution parameters.

Some conclusions, based on Tables 4-5 and Figures 19, are summarized as follows:(i)The SSR estimates when X and Y random variables are selected from NRSS are more efficient than the others under different ranking schemes, that is, and () for chosen values of sample sizes and true value of R (see, e.g., Figures 12 and Tables 45).(ii)The MSE of reliability estimates decreases with increasing set size, in all schemes (see Figures 12 and Tables 45).(iii)The MSE of SSR estimates decreases with increased R values in all ranked schemes (see Figure 3 and Tables 45).(iv)The MSE of SSR estimates takes the smallest values based on NRSS for all set sizes (see Figures 43 and Tables 45).(v)The MSE of is smaller than the corresponding for at all set sizes (see Figures 56 and Tables 45).(vi) is more efficient than at all set sizes (see Figures 78 and Tables 45).(vii)The MSE of the in case of odd set size is smaller than the corresponding in case of even set size at all cases (see Figure 9, Tables 4 and 5).(viii)The values of AB decrease as the sample size increases in most sampling schemes (see Tables 4 and 5).(ix)For all even set sizes, the AB of decreases as the value of R grows (see Table 5).(x)At set size (5, 5), the AB of decreases as the value of R increases (see Table 4).(xi)The AB of takes the smallest values in comparison with the others, i.e., for odd set size (see Table 4).(xii)For all true values of R at (2, 2), (4, 4), and (6, 6), the AB of receives the lowest values compared with the other sampling designs (see Table 5).

6. Data Analysis

In this section, two applications were utilized to verify the proposed estimators explored in the simulation investigation.

6.1. First Application

For illustration purposes, we take two data sets and illustrate all of the information. The two data sets were first published via Xia et al. [48], which represent the breaking strength of jute fibers at 5 mm, 10 mm, 15 mm, and 20 mm gauge lengths. We used the data sets about the (Tables 6, 7) breaking strengths of jute fiber at gauge lengths of 5 and 15 mm with sizes of the two samples equal to 30. The data sets are scheduled as follows.

We verify the effectiveness of the fitted model using the Kolmogorov–Smirnov (KS) goodness-of-fit test and its P value (PV). The KS distance between the empirical distribution and the fitted distribution is 0.1375 with PV = 0.622 for the first data set and 0.1828 with PV = 0.269 for the second data set, indicating satisfactory fits. The estimated densities, estimated CDF, PP plots, and empirical survival for the considered data appear in Figures 10 and 11.(1)Let and the estimates of , and are as follows:  = 29.464,  = 0.987,  = 2.102, and  = 0.246, and the log-likelihood value is  = 402.381.(2)Let and the estimates of , and are as follows: 11.91, 3.651, and 0.990, and the associated log-likelihood value is  = 402.6. The following hypotheses are examined.and in this case,  = 0.238, and hence, the null hypothesis is accepted; i.e., the assumption of is justified.

The NRSS and MRSS designs are applied to chosen data sets based on previous theoretical outcomes. Table 8 shows SSR estimates from the GED using suggested sampling techniques for some selected values of , and m using 5 repeated cycles. The MATHEMATICA program is used to create NRSS and MRSS samples from the data sets that have been evaluated.

6.2. Second Application

Here, we examine two real data sets and show how the recommended approaches may be applied in practice. Lawless [49] provided (Table 9) the first data set, while the second was provided (Table 10) by Linhart and Zucchini [50]. For different lifetime distributions, both are analyzed and fitted. Separately, we fitted the GED to each data set. The first data set reported here is from deep groove ball bearing endurance testing. For each of the 23 ball bearings in the life test, the numbers shown are the number of million revolutions before failure.

We use the KS goodness-of-fit test and its PV to validate the effectiveness of the fitted model. For the first data, the KS distance between the empirical and fitted distributions is 0.10583, with PV = 0.959. The KS distance in the second data set is 0.17197, with PV = 0.338. Figures 12 and 13 show the estimated densities, estimated CDF, PP plots, and empirical survival for both the data under consideration.

Let and the estimates of , and are as follows:  = 22.879,  = 0.9364,  = 0.1675, and , and the log-likelihood value is .

Let and the estimates of , and are as follows:  = 5.327,  = 0.3681, and  = 0.9584, and the associated log-likelihood value is  = −265.235. The following hypotheses are examined:and in this case,  = 0.116, and hence, the null hypothesis is accepted; i.e., the assumption of is justified.

Based on earlier theoretical results, the NRSS and MRSS designs are applied to selected data sets. Table 11 gives the SSR estimates from the GED via proposed sampling methods for some selected values of , and m using 5 replicated cycles. The MATHEMATICA program is used to generate NRSS and MRSS samples from considered data sets.

7. Concluding Remarks

In this study, we look at R = P[Y < X], for SSR estimation, assuming that strength X and stress Y are independent generalized exponential random variables. Under the MRSS and NRSS systems, the SSR estimators are considered in five cases. When both the stress and strength random variables have the same sample design NRSS, one of them is considered, whereas the other is used when the X and Y distributions have dissimilar sampling designs. We use the simulation research to evaluate and compare various behavior estimates using a set of criteria. In all issues, the MSE of SSR estimates using NRSS yielded the least values compared with the others based on MRSS data, according to the study’s outcomes. In general, SSR estimates that use both stress and strength data from the NRSS are more efficient than those that use strength and stress data from the MRSS. Furthermore, the SSR estimates when X is MRSS and Y is NRSS data are superior to the other in the reversed case (when X is NRSS and Y is MRSS data). Real data applications illustrate these results.

Data Availability

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

Conflicts of Interest

The authors declare that there are no conflicts of interest.