Inference of fuzzy reliability model for inverse Rayleigh distribution

: In this paper, the question of inference of the reliability parameter of fuzzy stress strength is attached to the difference between stress and strength values when and are independently distributed from inverse Rayleigh random variables. Including fuzziness in the stress-strength interference enables researchers to make more sensitive and precise analyses about the underlying systems. The maximum product of the spacing method for the reliability of fuzzy stress intensity inference has been introduced. As classical estimation methods and Bayesian estimation methods are used to estimate the reliability parameter , the maximum product of spacing and maximum likelihood estimation methods is used. The maximum product of spacing under fuzzy reliability of stress strength model is introducing in this paper. Markov Chain Monte Carlo approach is used to obtain Bayesian estimators of traditional and fuzzy reliability of stress strength for inverse Rayleigh model by using the Metropolis-Hastings algorithm. Using an extensive Monte Carlo simulation analysis, the outputs of the fuzzy reliability and traditional reliability estimators are contrasted. Finally, for example, and to verify the efficiency of the proposed estimators, a genuine data application is used


Introduction
The stress-strength reliability model is significant in reliability literature, medical, economic, and engineering applications where systems or specialized units may be exposed to randomly occurring environmental stresses such as pressure, temperature, and humidity.In this case, the survival of the system depends on its resistance.During World War II, it was found that some of the equipment such as radar and communication systems failed to work efficiently when they were used in an environment different from the environment for which they were designed.For that, experts started to consider the effects of environmental conditions while evaluating the reliability of equipment.
The computations of fuzzy reliability have been addressed by [1], where if and are independent but not identical random variables in distribution.The idea behind fuzzy reliability is that including fuzziness in the stress-strength interference enables researchers to make more sensitive and precise analyses about the underlying systems of life reliability and the system becomes more stable and reliable when the difference gets larger.The advantages of the fuzzy stress-strength reliability model over the traditional stress-strength reliability model are in considering the randomness in reliability engineering and the fuzziness of operating time.For more information see [1][2][3][4][5].
The fuzzy reliability is defined as (1) where is a fuzzy set and is an appropriate membership function on ; that is assumed increasing on the difference (readers are encouraged to read [5,6] who used the definition of the fuzzy stress-strength model to estimate , when and were independent inverse exponential random variables).
The probability that the system is strong enough to overcome the stress imposed on it is defined as system reliability.Traditional reliability may be equally explained as the region under the receiver operating characteristic (ROC) for diagnostic test or biomarkers, see reference [7].The ROC curve is exceedingly used in medical, biological, economic, and health service research, to evaluate the reliability and distinguish between two groups of subjects, generally non-satisfied and satisfied subjects.The research conducted on the traditional stress-strength reliability model focuses on computing, calculating, and estimating the reliability of different stress and strength distributions.For example, reference [8] estimated the traditional reliability of the stress-strength model for a generalized exponential distribution with three parameters.Confidence intervals estimation of traditional reliability of stress-strength model for generalized Pareto distribution has been discussed by [9].The stressstrength model of a generalized logistic distribution has been studied by [10].Reference [11] estimated the when and independent Lindley populations.In 2020, [12] and others discussed the estimation of when and are independent exponentiated Pareto random variables when samples are selected using some ranked set sampling designs.Reference [13] presented a comprehensive review of the traditional reliability of the stress-strength model.In 2021, [14] estimated the traditional stress strength reliability by the use of the MPS estimation method.
In our study, we used failure times in insulating fluid between two electrodes subjected to a voltage of 34 kV and 36 kV as an application and for illustrative purposes.These failure times were randomly observed and there was no reliable information available.Moreover, the failure times model may be difficult to measure due to the complexity of the action of electrodes.So, we used stress-strength reliability model in the presence of fuzziness.
In this article, estimation of fuzzy stress-strength reliability model , when and are independent but not symmetrically distributed inverse Rayleigh random variables, is discussed.The product of the spacing method was presented to infer the reliability of fuzzy stress strength by using different methods.The proposed estimators are obtained using the maximum likelihood estimation method (MLE) and the maximum product of the spacing estimation method (MPS) as well as Bayesian estimation when prior distributions are assumed exponential.Besides, a Monte Carlo simulation study is made to analyze and compare the performance of the different estimators.A real data application is conducted for illustration purposes and to test the estimated functions of the reliability parameter .Finally, the paper is concluded.

The stress strength model
An increase in the values of can be thought of equivalently as the increase in the difference of .With such consideration, the membership function can be redefined as , where Let and be two independent inverse Rayleigh random variables with scale parameters and , respectively.The inverted Rayleigh [ distribution has the following cumulative distribution function (CDF) and probability density function (PDF) for : and respectively, where is a scale parameter.The traditional reliability of the stress-strength model for inverse Rayleigh distribution was studied and calculated to be (see Kotz et al. [9]).Therefore, the fuzzy reliability of stress-strength is given by (4) The traditional reliability is always greater than the fuzzy reliability , and as .Figure 1 shows different values for when and changes simultaneously, Figure 2 shows fuzzy reliably values for different values of the constant and when also and changes simultaneously.

Inference of stress-strength model
In this section, the two methods (MLE and MPS) of estimation are used to estimate the fuzzy reliability parameter .Let and be two independent random samples from IR ( ) and IR ( ), respectively.

Likelihood estimation
The joint likelihood function of the IR distribution for the stress-strength model can be written as (5) and the log-likelihood function is given as (6) The normal equations for unknown parameters , are obtained by differentiating ( 6) partially concerning the parameters , and equate them to zero.The estimators can be obtained as a solution to the following equation: where , .The estimate by using MLE is .By using the invariance property of MLE, traditional reliability and the fuzzy reliability of IR distribution for the stressstrength model are obtained by using MLE's as following

Maximum product of spacing estimation
The maximum product spacing for stress-strength model is denoted as following.
)7) such that , where The natural logarithm of the product spacing function of the exponential distribution for stressstrength model is denoted as following . )9) To obtain the normal equations for the unknown parameters, we partially differentiate Eq (8) concerning the parameters , and equate them with zero.The estimators can be obtained as a solution of the following equations.
The above nonlinear equations can't be solved analytically to find of .So, by using optimization algorism as conjugate-gradient or Newton-Raphson optimization methods, the estimators of are obtained.Using the invariance property of MPS estimators of the MPS's which have been discussed by [13][14][15][16][17][18] and have concluded that it is the same as that of MLE, traditional reliability and fuzzy reliability for IR distribution for the stress-strength model can be computed.

Bayesian estimation
Assume that the parameters and are random variables with exponential prior distributions [11] with density function given by and , where and are hyperparameters of the prior distributions of and .The posterior distributions of and are obtained by combining the likelihood function with the prior distributions of both and and are found to be and Also, and Now, let and .After some manipulations and simplifications, we have that and respectively.The joint posterior distribution of and is then given by Let and since , then using standard transformation techniques, the joint posterior pdf of and will be Under squared error loss function, the Bayes estimate is the expected value of and is given by, )10) where .
The above integrals are hard to obtain, so numerical technique is used to evaluate the Bayes estimate.Markov Chain Monte Carlo (MCMC) approach is used to obtain Bayesian estimators.An important subclass of the MCMC techniques is Gibbs's sampling and more general Metropolis within Gibbs samplers see [24,25].The Metropolis-Hastings algorithm, together with the Gibbs sampling, are the two most popular examples of an MCMC method.It's similar to acceptance-rejection sampling, and the Metropolis-Hastings algorithm considers that to each iteration of the algorithm, a candidate value can be generated from the IR distributions.We use the Metropolis-Hastings within Gibbs sampling steps to generate random samples from conditional posterior densities of , and .For more information, see [19][20][21][22][23].

Simulation study
In this section, we provide a numerical comparison using the Monte Carlo simulation algorithm.We explain our algorithm through an application of fuzzy and traditional stress-strength models by different estimation methods.In this current simulation, we will compare MLE, MPS, and Bayesian estimation methods based on traditional and fuzzy stress-strength measures for estimating the parameter of IR distribution.The comparison is made through bias and mean squared errors (MSE) of the different measures.The simulations are made using the program for several combinations of the parameters and , and .Simulation Algorithm: We build our model by creating all simulation controls.At this point, we must follow the following steps in order: Step 1: Suppose different values of the parameters vector of IR distribution.
Step 2: Choose the different sample size of strength and different sample size of stress respectively.
Step 3: Generate the sample random values of IR distribution by using quantile function in equation and .
Step 4: Solve differential equations for each estimation method.To obtain the estimators of the parameters for IR distribution, we calculate , when k = 1; and when k = 5.
Step 5: Repeat this experiment (L-1) iterations.In each experiment, the parameter values are the same.
The generated random values are certainly varying from experiment to experiment even though the sample size is not changed.In the end, we have L-values of mean and MSE, and we restricted the number of repeats in this experiment to 10,000.Take the averages of these values and call them Monte Carlo estimates.
After completing the treatment stage, simulated outcomes are listed in Tables 1-4, Figure 3, and the following observations were observed:  The Bias and MSE decrease as sample sizes increase for all estimates.


For fixed values of , the biases and MSE's of estimates of parameters are increasing with increasing.


For fixed values of , the biases and MSE's of estimates of parameters are increasing with increasing, but the estimate of R for fuzzy stress strength and tradition are decreasing in approximately most situations.


For fixed values of , the bias and MSE of estimates of , R for fuzzy stress strength and tradition are decreasing with increasing, but the estimates of are increasing, in approximately most situations.


The MPS method is found to be superior to the MLE and Bayesian methods in most cases.


From the observed results of reliability, we note the efficiency of the fuzzy stress strength is over traditional stress strength in most situations according to Bias and MSE.


In fuzzy stress strength, the efficiency is better with decreased values of k according to Bias and MSE.

Application of real data
The numerical results of tradition and fuzzy stress-strength reliability estimation of the IR distribution for real data are presented in this section.
Two real stress and strength data sets contained times to breakdown down an insulating fluid between electrodes recorded at different voltages; these data have been discussed by [26].Data I and data II as presented in Table 4, are the failure times (in minutes) are presented, which are for an insulating fluid between two electrodes subject to a voltage of 34 kV (data set I) and 36 kV (data set II). Table 6 provides information about the estimated parameters of the IR model and the corresponding traditional and fuzzy reliability measures The graph of MCMC estimates for and using the MH algorithm are the plotting of estimates, histogram of estimates, and convergence of estimates are shown in Figure 3.In Figure 4, we note the convergence of MCMC estimates for and in the first quartile iteration.

Conclusions
The new approach of estimating fuzzy stress-strength reliability is getting much attention because of the properties of , which makes the analysis more sensitive and more reliable.Also, when the study results are not known completely, the use of traditional methods may be misleading, and the need for new approaches that can handle such situations is very important.In this paper, the stress and strength variables were distributed as inverted Rayleigh distribution.It can be noted that different membership functions will provide different measures of .It is also noted that the MPS method is superior to the MLE and Bayesian methods in most cases.

Figure 1 .
Figure 1.Values for traditional reliability parameter when and changes simultaneously.

Figure 2 .
Figure 2. Values for fuzzy reliability parameter for different values of when and changes simultaneously.

Figure 3 .
Figure 3. Relative Efficiency for different measures.

Figure 4 .
Figure 4. Convergence of MCMC estimation of and .

Table 1 .
Bias and MSE for MLE, MPS and Bayesian of IR parameters under the stress-

Table 2 .
Bias and MSE for MLE, MPS, and Bayesian of IR parameters under the stress-

Table 3 .
Bias and MSE for MLE, MPS and Bayesian of IR parameters under the stress-

Table 4 .
Bias and MSE for MLE, MPS, and Bayesian of IR parameters under the stress-

Table 5 .
Data sets of times to breakdown down an insulating fluid between electrodes recorded at different voltages.

Table 6 .
Tradition and the fuzzy stress-strength estimation of the IR distribution.