Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Reliability Estimation in Multicomponent Stress-strength Model based on Generalized Pareto Distribution

Parameshwar V Pandit*, Shubhashree Joshi

Department of Statistics, Bangalore University, Bangalore -560056

Abstract: The paper deals with the estimation of multicomponent system reliability where the system has k components with their strengths X₁, X₂, … Xk being independently and identically distributed random variables and each component is experiencing a random stress Y. The s-out-of-k system is said to function if atleast s out of k (1 ≤ s ≤ k) strength variables exceed the random stress. The reliability of such a system is derived when both strength and stress variables follow generalized Pareto distribution. The system reliability is estimated using maximum likelihood and Bayesian approaches. The maximum likelihood estimators are derived under both simple random sampling and ranked set sampling schemes. Lindley's approximation technique is used to get approximate Bayes estimators. The reliability estimators obtained from both the methods are compared by using mean squares error criteria and real data analysis is carried out to illustrate the procedure.

Keywords: generalized Pareto distribution, stress-strength reliability, ranked set sampling, simple random sampling, maximum likelihood estimator, Bayes estimator

Received June 14, 2018; Revised August 14, 2018; Accepted October 26, 2018

Citation: Parameshwar V Pandit, Shubhashree Joshi, Reliability Estimation in Multicomponent Stress-strength Model based on Generalized Pareto Distribution. American Journal of Applied Mathematics and Statistics, 6 (5). 210-217. January 2018. doi:10.12691/ajams-6-5-5

1. Introduction

The study of stress-strength model has received attention by the researchers because of its practical applications in the field of science and technology. The “stress-strength" reliability can be described as an assessment of reliability of a component in terms of random variable Y representing “stress” experienced by the component and X representing “strength” of a component available to overcome the possible stress, if the stress exceeds the strength, then the system will fail. The idea of stress-strength reliability R = P(X > Y) was introduced by Birnbaum [1] and developed by Birnbaum and McCarty [2]. Estimation of reliability in stress-strength model are considered in the literature when the stress and strength variables follow distributions such as exponential, Weibull, normal, gamma etc. Raqab et. al., [3] estimated the reliability of the model for a 3-parameter generalized exponential distribution. Wong [4] obtained confidence intervals for P(X > Y) when the underlying distribution is generalized Pareto. Angali et al. [5] considered Bayesian estimators of reliability for four parameter of bivariate exponential distribution under different loss functions.

The estimation of multicomponent stress-strength reliability for s-out-of-k and the other related systems have been extensively investigated by many authors in the literature. An s-out-of-k system functions when the system having k statistically independent and identically distributed components functions if s (1 ≤ s ≤ k) or more components with stand a common stress, which was first studied by Bhattacharyya and Johnson [6]. This type of systems can be seen in both industrial and military applications [7]. Estimation of multicomponent stress-strength reliability is considered for the log-logistic, generalized exponential, generalized inverted exponential, Rayleigh, Burr Type XII and generalized Rayleigh distributions respectively by Rao and Kantam [8], Rao [9,10,11,12] and Rao et al. [13]. Pandit and Kantu [14] considered estimation of multicomponent stress-strength reliability for parallel and series systems when strength and stress variables follow exponential distribution. Recently, Kizilaslan and Nadar [15] considered estimation of s-out-of-k stress-strength reliability using both classical and Bayesian approach when underlying distribution is Weibull.

In this paper, s-out-of-k system in stress-strength environment which has k independent and identical strength components and a common stress is studied. These kinds of situations may occur in real life. For example, in an electrical power station containing eight generating units which produce the electricity only if at least six units are working; the demand for the electricity of a district is fulfilled only if s-out-of-k wind rose are operating at all times and in a communication system for a navy can be successful only if six transmitters out of ten are operational to cover a district [15]. This paper considers the above problem when the underlying distribution are to follow generalized Pareto for the strength and stress variables.

Rezaei et al. [16] proposed generalized Pareto distribution by assuming X and Y as independent generalized Pareto distribution with common scale parameter and different shape parameter.

The probability density function and cumulative function of generalized Pareto distribution are given by:

and

where α and λ are the shape and scale parameters, respectively.

Here, reliability in a multicomponent stress-strength based on X, Y are two independent random variables, which follow generalized Pareto distributions with shape parameters α and β and with common scale parameter λ.

Let the random variables Y, X₁, X₂, … Xk be independent, G(y) be the cummulative distribution function of Y and F(x) be the common cummulative distribution function of X₁, X₂, … Xk. The reliability in a multicomponent stress-strength model developed by Bhattacharyya and Johnson [6] is

[picture] | (1)

The estimators of multicomponent system reliability are derived using maximum likelihood method under simple random sampling (SRS) and ranked set sampling (RSS) schemes. Also the Bayesian estimates are obtained using Lindley's approximation. The RSS method was introduced by McIntyre [17] and several authors are interested to study statistical inference related to reliability under RSS, as it has applications in different fields such as, reliability [18,19], environment [20,21,22].Estimation of reliability based on RSS is considered by Sengupta and Mukhati [23]. Muttlak et al. [24] estimated reliability when X and Y follow exponential distribution. Hussian [25] discussed estimation of stress-strength model for generalized inverted exponential distribution based on RSS and SRS, using maximum likelihood method to estimate R. Hassan et al. [26] studied the estimation of R when X and Y are independently distributed Burr XII random variables under different sampling schemes.

In section 2, reliability of the system is derived. Section 3 and 4, deals with Maximum likelihood estimation of Rs,k based on Simple random sampling (SRS) and ranked set sampling(RSS) respectively. In section 5, Bayes estimators of Rs,k is obtained using Lindley's approximation. Section 6 is devoted to simulation study in which the comparison of estimators of reliability are studied. In section 7, real data application is given and the conclusions are given in section 8.

The multicomponent stress-strength reliability

In this section multicomponent system reliability is considered when X and Y follow generalized Pareto distribution with parameters (α₁, λ) and (α₂, λ) respectively,

Where 1+λy = t

[picture] | (2)

2. MLE of Rs,k under SRS

Let and be two ordered random samples of size n, m respectively. Here, Strength and stress variables follow generalized Pareto distribution with shape parameters α₁, α₂ and scale parameter λ, then the likelihood function is given by

Thus, the log-likelihood function is

The likelihood equations for estimating α₁, α₂ and λ are

[picture] | (3)

[picture] | (4)

[picture] | (5)

From (3), (4) and (5), the MLE of α₁, α₂ and λ is

[picture] | (6)

[picture] | (7)

where can be obtained as the solution of non-linear equation of the form,

[picture] | (8)

here, can be obtained by using any iterative scheme

where, is the iterate of The iteration procedure should be stopped when is sufficiently small. After obtaining, the MLEs of α₁ and α₂ are obtained from (6) and (7). Hence, the MLE of is obtained by using the invariance property of MLEs, that is,

[picture] | (9)

3. MLE of Rs,k under RSS

Let be a ranked set samples with sample size where is the set size and is the number of cycles drawn from generalized Pareto distribution with parameters and λ and be a ranked set samples with sample size where is the set size and is the number of cycles drawn from generalized Pareto distribution with parameters and λ respectively. For convenience, denote and asand respectively.

The pdf of the random variables and are given by

[picture] | (10)

[picture] | (11)

The likelihood function of α₁, α₂ and λ is given by

Thus, the log-likelihood function of α₁, α₂ and λ is

where is a constant. The likelihood equations are

[picture] | (12)

[picture] | (13)

[picture] | (14)

A closed form expression for equations (12)-(14) is difficult to obtain analytically. Hence, one can use any iterative technique to solve these equations. The MLE of is obtained by using the invariance property of MLEs, that is,

[picture] | (15)

4. Bayes Estimation of Rs,k

Here we assume that all the parameters α₁, α₂ and λ are unknown and independent random variables with gamma priors

The pdf of gamma random variables X with parameters is

where .

Thus the joint prior α₁, α₂ and λ is

substituting and , the corresponding joint posterior distribution is given by

where

Then the Bayes estimator of under squared error (SE) loss function is given by

[picture] | (16)

It can be seen that, equation (16) cannot be reduced to a closed form. Hence, one can use Lindley's approximation method.

The simplest method to approximate is Lindley's [27] approximation method which approaches the ratio of the integrals as a whole and produces a single numerical result. If n is large, according to Lindley approximation, any ratio of the integrals of the form

where

is a function of and , is log of likelihood, is log of joint prior of and , can be written as

where and are the MLEs of and respectively,

and subscripts 1,2,3 on right-hand sides refer to and respectively, and

and is the element of the inverse of the matrix having elements .

In our case, and

we have,

and can be obtained as follows for i, j = 1,2,3

and the values of for i, j, k = 1, 2, 3

Since, ,

and . Hence

The quantities and Lijk, i, j, k = 1, 2, 3 are evaluated by replacing by .

Then, the Bayes estimator of is

5. Simulation Study

Simulation study consists of estimating multicomponent stress-strength reliability when the sample is generated from generalized Pareto distribution under SRS and RSS using ML and Bayesian approaches. The comparison of the estimates are done though mean squared error criteria based on 100000 random samples of size n, m, mx, my, rx and ry for the strength and stress populations, with different parameter values. The values of used for comparison are (1,1,0.6), (1,1.5,0.5) and (0.5,0.2,0.5). The corresponding true values of stress-strength reliability for s-out-of-k system with (s, k) =(1,3) are 0.75, 0.8476, 0.4747 and that for (s, k) = (2,4) are 0.6, 0.7229, 0.3315. The Bayesian estimates under squared error loss function using gamma prior are a₁ = 9, a₂ = 4, a₃ = 1, b₁ = 3, b₂ = 2, b₃ = 1 (prior1) and a₁ = 1, a₂ = 1, a₃ = 1, b₁ = 1, b₂ = 1, b₃ = 1 (prior2).

From the simulation study, it is observed the MSEs for the estimates decreases as the sample size increases in all the cases. The Bayes estimates of the under the squared error loss function have the smaller MSEs. It is also be seen that when comparing the maximum likelihood estimates under SRS and RSS, RSS performs better as it has smaller MSE when compared to SRS.

Table 1. MLEs under SRS and RSS, Bayes estimators and MSEs for the estimators of Rs,k

Table 2. MLEs under SRS and RSS, Bayes estimators and MSEs for the estimators of Rs,k

Table 3. MLEs under SRS and RSS, Bayes estimators and MSEs for the estimators of Rs,k

6. Data Analysis

In this section, we present a real data which was originally reported by Badar and Priest [28]. The data represents the strength measured in GPA for single carbon fibers and impregnated 1000-carbon fiber tows. Single fibers were tested under tension at gauge lengths of 20 mm (Data sets I) and 10 mm (Data set II), with sample sizes n = 69 and m = 63 respectively.

Data set I:

1.312, 1.314, 1.479, 1.552, 1.700, 1.803, 1.861, 1.865, 1.944, 1.958, 1.966, 1.997, 2.006, 2.021, 2.027, 2.055, 2.063, 2.098, 2.140, 2.179, 2.224, 2.240, 2.253, 2.270, 2.272, 2.274, 2.301, 2.301, 2.359, 2.382,2.382, 2.426, 2.434, 2.435, 2.478, 2.490, 2.511, 2.514, 2.535, 2.554, 2.566, 2.570, 2.586, 2.629, 2.633, 2.642, 2.648, 2.684, 2.697, 2.726, 2.770, 2.773, 2.800, 2.809, 2.818, 2.821, 2.848, 2.880, 2.954, 3.012, 3.067, 3.084, 3.090, 3.096, 3.128, 3.233, 3.433, 3.585, 3.585

Data set II:

1.901, 2.132, 2.203, 2.228, 2.257, 2.350, 2.361, 2.396, 2.397, 2.445, 2.454, 2.474, 2.518, 2.522, 2.525, 2.532, 2.575, 2.614, 2.616, 2.618, 2.624, 2.659, 2.675, 2.738, 2.740, 2.856, 2.917, 2.928, 2.937, 2.937, 2.977, 2.996, 3.030, 3.125, 3.139, 3.145, 3.220, 3.223, 3.235, 3.243, 3.264, 3.272, 3.294, 3.332, 3.346, 3.377, 3.408, 3.435, 3.493, 3.501, 3.537, 3.554, 3.562, 3.628, 3.852, 3.852, 3.871, 3.886, 3.971, 4.024, 4.027, 4.225, 4.395, 5.020

For the above data sets, we fit the generalized Pareto model and also checked the validity of the model using Kolmogorov-Smirnov (K-S) test for each data set. It was found that for data set I and II, the K-S distanced are 0.0479 and 0.0986 with the corresponding p value are 0.5638 and 0.6425 respectively. From the result, it shows that generalized Pareto distribution fits better for the data sets. The maximum likelihood estimate and Bayes estimate, based on the parameters and are obtained as under prior1 and under prior 2. For s =2 and k = 4 MLE and Bayes estimators are under prior1 and under prior2.

7. Conclusion

The main aim of this article, is to study the multicomponent system reliability which has k independent and identical strength components and each component exposed to a common random stress by assuming both strength and stress variables follow generalized Pareto distribution. The reliability of the system is estimated using maximum Likelihood under SRS and RSS scheme and Bayes approaches. The performance of these estimates are compared using MSEs, the results show that the MLE has greater MSE when compared to Bayes estimates. The ML estimates under RSS have lesser MSE than SRS. However, as sample size increases, MSEs of both the approaches, i.e., SRS and RSS are close to each other.

References

[1]   | Birnbaum, Z. W. (1956). ’On a use of Mann-Whitney statistics’. Proceeding Third Berkley Symposium on Mathematical Statistics and Probability, 1, 13-17.
[2]   | Birnbaum, Z. W. and McCarty, B.C. (1958). ‘A distribution-free upper confidence bounds for Pr(Y < X) based on independent samples of X and Y’. The Annals of Mathematical Statistics, 29(2), 558-562.
[3]   | Raqab, M. Z., Madi, M. T. and Kundu, D (2008). Estimation of P(Y < X) for the 3-parameter generalized exponential distribution. Communications in Statistics-Theory and Methods, 37(18), 2854-2864.
[4]   | Wong, A. (2012). Interval estimation of P(Y < X) for generalized Pareto distribution. Journal of Statistical Planning and Inference, 142, 601-607.
[5]   | Angali, K. A, Latifi S. M, Hanagal, D. D. (2014). Bayesian estimation of bivariate exponential distributions based on LINEX and Quadratic loss functions: a survival approach with censored samples. Communications in Statistics-Theory and Methods, 43(1), 31-44.
[6]   | Bhattacharyya, G. K. and Johnson, R. A. (1974). Estimation of reliability in multicomponent stress-strength model. Journal of American Statistical Association, 69, 966-970.
[7]   | W. Kuo and M. J. Zuo. (2003). Optimal Reliability Modeling, Principles and Applications. New York, NY, USA: Wiley.
[8]   | G. S. Rao and R.R.L. Kantam. (2010). Estimation of reliability in multicomponent stress-strength model: Log-logistic distribution. Electronic Journal of Applied Statistical Analysis, 3(2), 75-84.
[9]   | G. S. Rao. (2012). Estimation of reliability in multicomponent stress-strength model based on generalized exponential distribution. Colombian Journal of statistics, 35(1), 67-76.
[10]   | G. S. Rao. (2012). Estimation of reliability in multicomponent stress-strength model based on Rayleigh distribution. Probability Statistics forum, 5, 150-161.
[11]   | G. S. Rao (2014). Estimation of reliability in multicomponent stress-strength model based on generalized Rayleigh distribution. Journal of Modern Applied Statistical Methods, 13(1), 367-379.
[12]   | G. S. Rao, R. R. L. Kantam, K. Rosaiah, and J. P. Reddy (2013). Estimation of reliability in multicomponent stress-strength model based on inverse Rayleigh distribution. Journal of Statistics and Applied Probability, 3, 261-267.
[13]   | G. S. Rao, Aslam, M and Kundu, D (2015). Burr Type XII distribution parametric estimation and estimation of reliability in multicomponent stress-strength. Communications in Statistics-Theory and Methods, 44(23), 4953-4961.
[14]   | Pandit, P. V and Kantu, Kala, J (2013). System reliability estimation in multicomponent exponential stress-strength models. International Journal of Reliability and Applications, 14(2), 97-105.
[15]   | Nadar, M. and Kizilaslan, F. (2015). ‘Classical and Bayesian estimation of Reliability in Multicomponent Stress-Strength Model based on Weibull distribution’, Revista Clombiana de Estadistica, 38(2), 467-484.
[16]   | S. Rezaei, R. Tahmasbi and M. Mahmoodi (2010). Estimation of Pr(X > Y) for generalized Pareto distribution. Journal of Statistical Planning and Inference, 140, 480-494.
[17]   | G. A. McIntyre (1952). A method for unbiased selective sampling using ranked sets. Australian Journal of Agricultural Research, 3, 385-390.
[18]   | Kvam, P. H, Samaniego, F, J. (1993). On maximum likelihood estimation based on ranked set sampling with application to reliability. In Basu: A., ed, Advanced in Reliability. Amsterdam: North- Holland, 215-229.
[19]   | Kvam, P. H, Samaniego, F, J. (1994). Nonparametric maximum likelihood estimation based on ranked set sampling. American Journal of statistical Association, 89, 526-537.
[20]   | Abu-Dayyeh, W. A., Muttlak, H. A. (1996). Using ranked set sampling for testing hypothesis on the scale parameter for exponential and uniform distribution. Pakistan Journal of statistics, 12, 131-138.
[21]   | H. A. Muttlak (1997). Median ranked set sampling. Journal of Applied Statistical Science, 6, 245-255.
[22]   | Yu, P. L. H., Lam. K (1997). Regression estimator in ranked set sampling. Biometrics, 53, 1070-1081.
[23]   | S. Sengupta and S. Mukhuti (2008). Unbiased estimation of Pr(X > Y) using ranked set sample data. Statistics, 42, 223-230.
[24]   | H. A. Muttlak, W. A. Abu-Dayyah, M. F. Saleh, and E. Al-Sawi (2010). Estimating P(Y < X) using ranked set sampling in case of the exponential distribution. Communications in Statistics: Theory and Methods, 39, 1855-1868.
[25]   | M. A. Hussian (2014). Estimation of stress-strength model for generalized inverted exponential distribution using ranked set sampling. International Journal of Advance in Engineering and Technology, 6, 2354-2362.
[26]   | A. S. Hassan, S. M. Assar, and M. Yahya (2015). Estimation of P(Y < X) for Burr distribution under several Modifications for ranked set sampling. Australian Journal of Basic and Applied Sciences, 9(1), 124-140.
[27]   | D. V. Lindley (1980). Approximate Bayes method. Trabajos de Estadistica, 3, 281-288.
[28]   | Badar, M. G. & Priest, A. M. (1982). Statistical aspects of fibre and bundle strength in hybrid composites. In T. Hayashi, K. Kawata, and S. Umekawa (eds.), Progress in Science and Engineering Composites, (pp. 1129-1136). Tokyo: ICCM-IV.