Abstract

We consider the estimation of stress-strength reliability based on lower record values when and are independently but not identically inverse Rayleigh distributed random variables. The maximum likelihood, Bayes, and empirical Bayes estimators of are obtained and their properties are studied. Confidence intervals, exact and approximate, as well as the Bayesian credible sets for are obtained. A real example is presented in order to illustrate the inferences discussed in the previous sections. A simulation study is conducted to investigate and compare the performance of the intervals presented in this paper and some bootstrap intervals.

1. Introduction

The inverse Rayleigh distribution is an important lifetime distribution in survival analysis that has many applications in the area of reliability studies. Voda [1] presented some properties of the maximum likelihood estimator for inverse Rayleigh distribution and mentioned that the distribution of lifetimes of several types of experimental units can be approximated by the inverse Rayleigh distribution. The cumulative distribution function (cdf) of the inverse Rayleigh distribution with scale parameter is given by The probability density function (pdf) corresponding to (1) is given by We will denote inverse Rayleigh distribution with scale parameter by .

In many real life applications such as meteorology, hydrology, sports, life-tests and so on, we are dealing with record values. In industry and reliability studies, many products may fail under stress. For example, a wooden beam breaks when sufficient perpendicular force is applied to it, an electronic component ceases to function in an environment of too high temperature, and a battery dies under the stress of time. But the precise breaking stress or failure point varies even among identical items. Hence, in such experiments, measurements may be made sequentially and only values larger (or smaller) than all previous ones are recorded. Data of this type are called record data. Thus, the number of measurements made is considerably smaller than the complete sample size. This measurement saving can be important when the measurements of these experiments are costly if the entire sample was destroyed.

Let be a sequence of independent and identically distributed (iid) random variables with an absolutely continuous cumulative distribution function (cdf) and probability density function (pdf) . An observation is called an upper record if its value exceeds all previous observations; that is, is an upper record if   for every . An analogous definition can be given for lower records. Records were first introduced and studied by Chandler [2]. Interested readers may refer to Arnold et al. [3], Ahmadi [4], and Gulati and Padgett [5] for more details and applications of record values.

The problem of estimating arises in the context of mechanical reliability of a system with strength and stress and is chosen as a measure of system reliability. The system fails if and only if, at any time, the applied stress is greater than its strength. This type of reliability model is known as the stress-strength model. This problem also arises in situations where and represent lifetimes of two devices and one wants to estimate the probability that one fails before the other. For example, in biometrical studies, the random variable may represent the remaining lifetime of a patient treated with a certain drug while represent the remaining lifetime when treated by another drug. The estimation of stress-strength reliability is very common in the statistical literature. The reader is referred to Kotz et al. [6] for other applications and motivations for the study of the stress-strength reliability.

The problem of estimating the stress-strength reliability in the inverse Rayleigh distribution was considered by Rao et al. [7] for ordinary samples. Soliman et al. [8] discussed different methods of estimation for the inverse Rayleigh distribution based on lower record values. Sindhu et al. [9] and Feroze and Aslam [10] considered the Bayesian estimation for the parameter of the inverse Rayleigh distribution under left censored data and under singly and doubly type II censored data, respectively. In this paper, we consider the problem of estimating the stress-strength reliability in the inverse Rayleigh distribution based on lower record values.

The rest of the paper is organized as follows. In Section 2, we discussed likelihood inference for the stress-strength reliability, while in Section 3 we considered Bayesian inference. In Section 4, we presented a real example. A simulation study is described in Section 5. Finally conclusion of the paper is provided in Section 6.

2. Likelihood Inference

Let and be independent random variables. Let be the stress-strength reliability. Then, We will consider estimating based on lower record values on both variables. Let be a set of lower records from and let be an independent set of lower records from . The likelihood functions are given by (Ahsanullah [11]) where and are the pdf and cdf of , respectively, and and are the pdf and cdf of , respectively. Substituting , , , and in the likelihood functions and using (4), we obtain It can be easily shown that the log-likelihood functions are given by Differentiate above equations with respect to the parameters and and equating with zero, the maximum likelihood estimators (MLEs) of and based on the lower record values are given by Therefore using the invariance properties of the maximum likelihood estimation, the MLE of is given by To study the distribution of we need the distributions of and . Consider first ; the pdf of the th lower record value is given by (Ahsanullah [11]) Consequently, using standard procedure of transformation of random variables, the pdf of is given by This is recognized as the inverted gamma distribution; that is, . Similarly, the pdf of is given by Thus . Therefore we can find the pdf of Consider . Note that, by the properties of the inverted gamma distribution and its relation with the gamma distribution, we have and . Hence and . Note that, by the independence of two random quantities, we have Hence, has a scaled distribution. It follows that the distribution of is that of which can be obtained using simple transformation techniques. This fact can be used to construct the following confidence interval for : Records are rare in practice (Arnold et al. [3]) and sample sizes are often very small; however, intervals based on the asymptotic normality of MLEs can be of interest in cases when the number of records is sufficiently large. This is because of their optimal asymptotic properties under very general conditions (Lehmann [12]). Note that as , where denotes convergence in distribution and is the asymptotic variance given by the reciprocal of the Fisher information: Similarly, as , where . Let and such that where ; it follows that . Since we have where A approximate confidence interval for based on this asymptotic result is given by where is obtained by substituting for and the MLEs and in the asymptotic standard deviation . In these calculations we assumed that is the smaller sample in size and is the larger. However, if this is not the case then the formula for the asymptotic variance in the asymptotic interval should be modified accordingly.

3. Bayesian Inference

Consider the likelihood functions of and based on the two sets of lower record values from the inverse Rayleigh distribution mentioned in previous section. These suggest that the conjugate family of prior distributions for and is the Gamma family of probability distributions: where , , , and are the parameters of the prior distributions of and , respectively. Combining the prior distributions in (21), with the likelihood functions in (5), we obtain the posterior distributions of and as It follows that It follows that the posterior distribution of is equal to that of , where and The Bayes estimator under squared error loss is the mean of this posterior distribution which may be approximated. A Bayesian confidence interval for is given by The case of a noninformative prior can be treated similarly. One of the well-known noninformative priors is the Jeffreys prior that say, . This suggests that prior densities for and are proportional to and , respectively. Using direct arguments one can show that and . It follows that the posterior distribution of is equal to that of , where and . Therefore a Bayesian confidence interval for is given by

Now consider the case when the parameters of prior distributions are themselves unknown. We consider the conjugate prior distributions for and above when the parameters and are unknown. In the empirical Bayes model, we must estimate them. So we calculate the marginal distribution of lower records, with densities Using (5) and (21), we obtain It can be shown that the maximum likelihood estimators (MLEs) of and based on the marginal distributions (28) are With substitution and for and in prior distributions and using similar arguments above, one can show that It follows that It follows that ; the empirical posterior distribution of is equal to that of , where and . A Bayesian confidence interval for is given by The construction of highest posterior density (HPD) regions requires finding the set , where is the largest constant such that . This often requires numerical optimization techniques. Chen and Shao [13] presented a simple Monte Carlo technique to approximate the HPD region.

4. A Real Example

In order to illustrate the inferences discussed in the previous sections, in this section, we present a data analysis for two data sets reported by Stone [14]. He reports an experiment in which specimens of solid epoxy electrical-insulation were studied in an accelerated voltage life test. For each of two voltage levels 52.5 and 57.5 kV, 20 specimens were tested. Failure times, in minutes, for the insulation specimens are given as follows.

Data Set 1. It belongs to voltage level 52.5: 4690, 740, 1010, 1190, 2450, 1390, 350, 6095, 3000, 1458, 6200, 550, 1690, 745, 1225, 1480, 245, 600, 246, 1805.

Data Set 2. It belongs to voltage level 57.5: 510, 1000, 252, 408, 528, 690, 900, 714, 348, 546, 174, 696, 294, 234, 288, 444, 390, 168, 558, 288.

We fit the inverse Rayleigh distribution to the two data sets separately. We used the Kolmogorov-Smirnov (K-S) tests for each data set to fit the inverse Rayleigh model. It is observed that for data sets 1 and 2, the K-S distances are 0.2119 and 0.2280 with the corresponding values 0.2878 and 0.2496, respectively. Therefore, it is clear that inverse Rayleigh model fits well to both the data sets. We plot the empirical distribution functions and the fitted distribution functions in Figures 1 and 2. These Figures show that the empirical and fitted models are very close for each data set.

Next, we consider the lower record values from the above observed data as follows:: ,: .Based on the above data, we obtain the MLEs of and as 240100 and 112896, respectively. Therefore, the MLE of becomes . The corresponding 95% confidence interval based on (15) is equal to (0.3242, 0.9041). Letting , , and , the Bayesian 95% confidence interval based on (25) is equal to (0.3371, 0.8453). So, the Bayesian 95% confidence interval based on (26) is equal to (0.3242, 0.9041). Finally, the Bayesian 95% confidence interval based on (32) is equal to (0.3935, 0.8745).

5. A Simulation Study

In this section, a simulation study is conducted to investigate and compare the performance of the confidence intervals presented in this paper and some bootstrap intervals. There are several bootstrap based intervals discussed in the literature (Efron and Tibshirani [15]). It is important here to note that all inference procedures in this paper depend only on the smallest records, and . Therefore we will use the parametric bootstrap based on the marginal distribution of as given in (9). In what follows we describe the bootstrapping procedure.(1)Calculate , , and , the maximum likelihood estimators of , , and based on and .(2)Generate from the distribution given in (9) with replaced by and generate similarly.(3)Calculate , , and using the and obtained in Step (2).(4)Repeat Steps (2) and (3) times to obtain .Then we can calculate the following bootstrap intervals.

Normal Interval. The simplest bootstrap interval is the normal interval where is the bootstrap estimate of the standard error based on .

Basic Pivotal Interval. The bootstrap basic pivotal confidence interval is where is the quantile of .

Percentile Interval. The bootstrap percentile interval is defined by that is, just use the and quantiles of the bootstrap sample.

Interested readers may refer to DiCiccio and Efron [16] and the references contained therein to observe more details.

In the simulation design we used all combinations of and . We used and . The value of is determined by the choice of and . The confidence level taken is and . For each combination of the simulation indices we generated 2000 samples of lower records from the distributions of and . For each generated pair of samples we calculated the following intervals:(1)ML: the interval based on the MLE given in (15),(2)Bayes: the interval based on the Bayes estimator given in (25),(3)J.B: the interval based on the Bayes estimator given in (26),(4)E.B: the interval based on the empirical Bayes estimator given in (32),(5)Norm: the normal interval,(6)Basic: the basic pivotal interval,(7)Perc: the percentile interval.

The empirical coverage probability and expected lengths of intervals are obtained by using the 2000 replications. In the interval based on the Bayes estimators we used , , , and wherever we need them. For bootstrap intervals we used 1000 bootstrap samples. The results of our simulations are given in Tables 1 and 2.

6. Conclusion and Discussion

Based on simulation results in Tables 1 and 2, we observe that the length of the intervals is maximized when and gets shorter and shorter as we move away to the extremes. Increasing the sample size on either variable also results in shorter intervals. The performance of both basic pivotal interval and percentile interval is similar in terms of expected length but in terms of coverage rate percentile interval has the better performance. The percentile interval appears to be the best among bootstrap intervals. The interval based on the MLE and the interval based on the Bayes estimator given in (26) appear to perform almost as well as the percentile interval. The interval based on the Bayes estimator given in (25) has best coverage rate and the short expected length between the other intervals when , but it has the low coverage rate and the long expected length for small values of since it is dependent on choice of and values. Furthermore, the interval based on the empirical Bayes estimator has the shortest expected length between the other intervals but it has the low coverage rate.

Based on the above discussion, we can conclude that, between the intervals obtained in this paper, the intervals based on the MLE, and the Bayes estimator given in (26) and between the bootstrap intervals, percentile interval simultaneously has the short expected length and very good coverage rate in comparison with the other intervals. Hence, we recommend using this confidence interval in all.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors express their sincere thanks to the editor and the referees for their constructive criticisms and excellent suggestions which led to a considerable improvement in the presentation of the paper.