Estimation of Reliability in Multicomponent Stress-Strength Based on Generalized Rayleigh Distribution

Similar documents
System in Weibull Distribution

Reliability estimation in Pareto-I distribution based on progressively type II censored sample with binomial removals

PROBABILITY AND STATISTICS Vol. III - Analysis of Variance and Analysis of Covariance - V. Nollau ANALYSIS OF VARIANCE AND ANALYSIS OF COVARIANCE

PARAMETER ESTIMATION IN WEIBULL DISTRIBUTION ON PROGRESSIVELY TYPE- II CENSORED SAMPLE WITH BETA-BINOMIAL REMOVALS

Estimation in Step-stress Partially Accelerated Life Test for Exponentiated Pareto Distribution under Progressive Censoring with Random Removal

BAYESIAN AND NON BAYESIAN ESTIMATION OF ERLANG DISTRIBUTION UNDER PROGRESSIVE CENSORING

XII.3 The EM (Expectation-Maximization) Algorithm

BAYESIAN CURVE FITTING USING PIECEWISE POLYNOMIALS. Dariusz Biskup

Statistical analysis of Accelerated life testing under Weibull distribution based on fuzzy theory

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Case Study of Cascade Reliability with weibull Distribution

Bayesian estimation using MCMC approach based on progressive first-failure censoring from generalized Pareto distribution

Excess Error, Approximation Error, and Estimation Error

Statistical inference for generalized Pareto distribution based on progressive Type-II censored data with random removals

Research Article Maximum Likelihood Estimator of AUC for a Bi-Exponentiated Weibull Model

Inference for the Rayleigh Distribution Based on Progressive Type-II Fuzzy Censored Data

Accelerated Life Testing in Interference Models with Monte-Carlo Simulation

Change Point Estimation for Pareto Type-II Model

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Chapter 11: Simple Linear Regression and Correlation

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Least Squares Fitting of Data

COS 511: Theoretical Machine Learning

Interval Estimation of Stress-Strength Reliability for a General Exponential Form Distribution with Different Unknown Parameters

Multipoint Analysis for Sibling Pairs. Biostatistics 666 Lecture 18

Designing Fuzzy Time Series Model Using Generalized Wang s Method and Its application to Forecasting Interest Rate of Bank Indonesia Certificate

Test for Intraclass Correlation Coefficient under Unequal Family Sizes

Modeling and Simulation NETW 707

Double Acceptance Sampling Plan for Time Truncated Life Tests Based on Transmuted Generalized Inverse Weibull Distribution

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation

Economics 130. Lecture 4 Simple Linear Regression Continued

A new Approach for Solving Linear Ordinary Differential Equations

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Applied Mathematics Letters

Least Squares Fitting of Data

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)

Estimation of the Mean of Truncated Exponential Distribution

Research Article Green s Theorem for Sign Data

DERIVATION OF THE PROBABILITY PLOT CORRELATION COEFFICIENT TEST STATISTICS FOR THE GENERALIZED LOGISTIC DISTRIBUTION

1.3 Hence, calculate a formula for the force required to break the bond (i.e. the maximum value of F)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

On the Eigenspectrum of the Gram Matrix and the Generalisation Error of Kernel PCA (Shawe-Taylor, et al. 2005) Ameet Talwalkar 02/13/07

LINEAR REGRESSION ANALYSIS. MODULE VIII Lecture Indicator Variables

Chapter 13: Multiple Regression

2 Complement Representation PIC. John J. Sudano Lockheed Martin Moorestown, NJ, 08057, USA

Outline. Prior Information and Subjective Probability. Subjective Probability. The Histogram Approach. Subjective Determination of the Prior Density

Determination of the Confidence Level of PSD Estimation with Given D.O.F. Based on WELCH Algorithm

Parameters Estimation of the Modified Weibull Distribution Based on Type I Censored Samples

Several generation methods of multinomial distributed random number Tian Lei 1, a,linxihe 1,b,Zhigang Zhang 1,c

Comparative Analysis of Bradley-Terry and Thurstone-Mosteller Paired Comparison Models for Image Quality Assessment

Computing MLE Bias Empirically

Durban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications

2nd International Workshop on Materials Engineering and Computer Sciences (IWMECS 2015)

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

Assignment 5. Simulation for Logistics. Monti, N.E. Yunita, T.

BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS. M. Krishna Reddy, B. Naveen Kumar and Y. Ramu

LECTURE :FACTOR ANALYSIS

Statistics for Economics & Business

Convexity preserving interpolation by splines of arbitrary degree

A Comparative Study for Estimation Parameters in Panel Data Model

ON WEIGHTED ESTIMATION IN LINEAR REGRESSION IN THE PRESENCE OF PARAMETER UNCERTAINTY

SEMI-EMPIRICAL LIKELIHOOD RATIO CONFIDENCE INTERVALS FOR THE DIFFERENCE OF TWO SAMPLE MEANS

A Hybrid Variational Iteration Method for Blasius Equation

Small-Sample Equating With Prior Information

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Hydrological statistics. Hydrological statistics and extremes

/ n ) are compared. The logic is: if the two

Chapter 8 Indicator Variables

Statistics for Business and Economics

Interval Estimation in the Classical Normal Linear Regression Model. 1. Introduction

BERNSTEIN-TYPE APPROXIMATION OF SMOOTH FUNCTIONS

Quantum Particle Motion in Physical Space

MD. LUTFOR RAHMAN 1 AND KALIPADA SEN 2 Abstract

ECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics

ASYMPTOTIC PROPERTIES OF ESTIMATES FOR THE PARAMETERS IN THE LOGISTIC REGRESSION MODEL

ISSN X Reliability of linear and circular consecutive-kout-of-n systems with shock model

Bayesian Estimation and Prediction of Generalized Pareto Distribution Based on Type II Censored Samples

ANALYSIS OF SIMULATION EXPERIMENTS BY BOOTSTRAP RESAMPLING. Russell C.H. Cheng

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore

International Journal of Mathematical Archive-9(3), 2018, Available online through ISSN

Journal of Global Research in Computer Science A MARKOV CHAIN MODEL FOR ROUND ROBIN SCHEDULING IN OPERATING SYSTEM

STAT 3008 Applied Regression Analysis

Joint Statistical Meetings - Biopharmaceutical Section

1 Review From Last Time

A Generalization of the Weibull Distribution with Applications

DISCRETE GENERALIZED RAYLEIGH DISTRIBUTION

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6

VERIFICATION OF FE MODELS FOR MODEL UPDATING

First Year Examination Department of Statistics, University of Florida

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

Lecture 16 Statistical Analysis in Biomaterials Research (Part II)

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors

Elastic Collisions. Definition: two point masses on which no external forces act collide without losing any energy.

RELIABILITY ASSESSMENT

Numerical Solution of Ordinary Differential Equations

ON A CLASS OF RENEWAL QUEUEING AND RISK PROCESSES

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010

x i1 =1 for all i (the constant ).

Manpower Planning Process with Two Groups Using Statistical Techniques

Transcription:

Journal of Modern Appled Statstcal Methods Volue 13 Issue 1 Artcle 4 5-1-014 Estaton of Relablty n Multcoponent Stress-Strength Based on Generalzed Raylegh Dstrbuton Gadde Srnvasa Rao Unversty of Dodoa, Dodoa, Tanzana, gaddesrao@yahoo.co Follow ths and addtonal wors at: http://dgtalcoons.wayne.edu/jas Part of the Appled Statstcs Coons, Socal and Behavoral Scences Coons, and the Statstcal Theory Coons Recoended Ctaton Rao, Gadde Srnvasa (014) "Estaton of Relablty n Multcoponent Stress-Strength Based on Generalzed Raylegh Dstrbuton," Journal of Modern Appled Statstcal Methods: Vol. 13 : Iss. 1, Artcle 4. DOI: 10.37/jas/1398918180 Avalable at: http://dgtalcoons.wayne.edu/jas/vol13/ss1/4 Ths Regular Artcle s brought to you for free and open access by the Open Access Journals at DgtalCoons@WayneState. It has been accepted for ncluson n Journal of Modern Appled Statstcal Methods by an authorzed edtor of DgtalCoons@WayneState.

Journal of Modern Appled Statstcal Methods May 014, Vol. 13, No. 1, 367-379. Copyrght 014 JMASM, Inc. ISSN 1538 947 Estaton of Relablty n Multcoponent Stress-Strength Based on Generalzed Raylegh Dstrbuton Gadde Srnvasa Rao Unversty of Dodoa Dodoa, Tanzana A ultcoponent syste of coponents havng strengths followng - ndependently and dentcally dstrbuted rando varables x 1, x,, x and each coponent experencng a rando stress Y s consdered. The syste s regarded as alve only f at least s out of (s < ) strengths exceed the stress. The relablty of such a syste s obtaned when strength and stress varates are gven by a generalzed Raylegh dstrbuton wth dfferent shape paraeters. Relablty s estated usng the axu lelhood (ML) ethod of estaton n saples drawn fro strength and stress dstrbutons; the relablty estators are copared asyptotcally. Monte-Carlo sulaton s used to copare relablty estates for the sall saples and real data sets llustrate the procedure. Keywords: Generalzed Raylegh dstrbuton, relablty estaton, stress-strength, ML estaton, confdence ntervals Introducton Surles and Padgett (1998, 001) ntroduced the two-paraeter Burr Type X dstrbuton and naed t the generalzed Raylegh dstrbuton. Note that the twoparaeter generalzed Raylegh dstrbuton s a partcular eber of the generalzed Webull dstrbuton, orgnally proposed by Mudholar and Srvastava (1993). The two-paraeter Burr Type X dstrbuton s referred to as the generalzed Raylegh dstrbuton (GRD). For > 0 and > 0, the twoparaeter GRD has the densty functon; 1 ( x) ( x) f ( x;, ) xe 1 e for x 0 (1) Dr. G. Srnvasa Rao s a Professor of Statstcs. Eal at: gaddesrao@yahoo.co. 367

ESTIMATION OF RELIABILITY IN STRESS-STRENGTH and the dstrbuton functon s gven by ( x) F( x;, ) 1 e for x 0. () Here and are the shape and scale paraeters respectvely. The GRD has been studed extensvely by Kundu and Raqab (005) and Raqab and Kundu (005). The two-paraeter GRD s denoted by GR( )., Surles and Padgett (001) showed that the two-paraeter GR dstrbuton can be used effectvely n odelng strength as well as general lfete data. Ths artcle studes relablty n a ultcoponent stress-strength based on X, Y, two ndependent rando varables, where X and Y fallow generalzed Raylegh dstrbutons wth shape paraeters and respectvely and wth coon scale paraeter. Let the rando saples y, x1, x,... x be ndependent, G(y) be the contnuous dstrbuton functon of Y and F(x) be the coon contnuous dstrbuton functon of x1, x,... x. The relablty n a ultcoponent stressstrength odel developed by Bhattacharyya and Johnson (1974) s R =P at least s of the( x, x,... x )exceed Y s, 1 s [1 G( y)] [ G( y)] df( y), (3) where x1, x,... x are ndependently and dentcally dstrbuted (d) wth coon dstrbuton functon F(x), ths syste s subjected to coon rando stress Y. The probablty n (3) s called relablty n a ultcoponent stress-strength odel (Bhattacharyya & Johnson, 1974). The survval probablty of sngle coponent stress-strength versons have been consdered by several authors assung varous lfete dstrbutons for the stress-strength rando varates (Ens & Gesser, 1971; Downtown, 1973; Awad & Gharraf, 1986; McCool, 1991; Nand & Ach, 1994; Surles & Padgett, 1998; Raqab & Kundu, 005; Kundu & Gupta, 005, 006; Raqab, et al., 008; Kundu & Raqab, 009). Relablty n a ultcoponent stress-strength was developed by Bhattacharyya and Johnson (1974) and Pandey and Borhan Uddn (1985) and the references theren cover the study of estatng P( Y X ) n any standard dstrbutons assgned to one or both of stress, strength varates. Recently Srnvasa Rao and Kanta (010) 368

GADDE SRINIVASA RAO studed estaton of relablty n ultcoponent stress-strength for log-logstc dstrbuton. Suppose a syste, wth dentcal coponents, functons f s(1 s ) or ore of the coponents sultaneously operate. In ts operatng envronent, the syste s subjected to a stress Y whch s a rando varable wth dstrbuton functon G(.). The strengths of the coponents, that s the nu stresses to cause falure, are ndependent and dentcally dstrbuted rando varables wth dstrbuton functon F(.). Then the syste relablty, whch s the probablty that the syste does not fal, s the functon R gven n (3). The estaton of survval probablty n a ultcoponent stress-strength syste when the stress, strength varates are followng Raylegh dstrbuton s not pad uch attenton. Therefore, ths artcle studes the estaton of relablty n ultcoponent stress-strength odel wth reference to Raylegh dstrbuton. Maxu Lelhood Estator of R Let X ~ GR(, ) and Y ~ GR(, ) wth unnown shape paraeters, and coon scale paraeter, where X and Y are ndependently dstrbuted. The relablty n ultcoponent stress- strength for generalzed Raylegh dstrbuton usng (3) results n: 1 ( y) ( y) ( y) ( y) Rs, 1 1 e 1 e ye 1 e dy s 0 s 1 1 [1 ] [ ] t t t dt where t 1 e y s 0 1 1 z [1 z] dz f z t, 0 s (, 1). After splfcaton ths reduces to ( )! Rs, ( j) ( )! s j0 1 (4) 369

ESTIMATION OF RELIABILITY IN STRESS-STRENGTH because and are ntegers. The probablty n (4) s called relablty n a ultcoponent stress-strength odel. If and are not nown, t s necessary to estate and to estate R. In ths artcle and are estated usng the ML ethod. The estates are substtuted n to obtan an estate of R usng equaton (4). It s nown that the ethod of Maxu Lelhood Estaton (MLE) has nvarance property. In ths drecton, ths artcle proposes the ML estator for the relablty of a ultcoponent stress-strength odel by consderng the estators of the paraeters of stress, strength dstrbutons by ML ethod of estaton n a generalzed Raylegh dstrbuton. x x... x ; y y... y be two ordered rando saples of Let 1 n 1 sze n, respectvely on strength, stress varates each followng GRD wth shape paraeters and, coon scale paraeter. The log-lelhood functon of the observed saple s L(,, ) ( n )ln nln ln n n ( n )ln x y j ln x ln y j 1 j1 1 j1 n ( x ) ( y ) e e j ( 1) ln 1 ( 1) ln 1 1 j1 (5) The MLEs of, and, for exaple, ˆ, ˆ and ˆ, respectvely can be obtaned as the teratve soluton of 1 n L n 0 ln 1 ( x ) e 0 (6) j1 L ( y ) 0 ln 1 j e 0 (7) 370

GADDE SRINIVASA RAO L ( n) 0 ( 1) j1 ( y j ) j y j ye ( 1) 0 1 e n n ( x ) xe x yj ( x ) 1 j1 1 1 e (8) and then fro (6), (7) and (8) ˆ ˆ n 1 j1 n ln 1 ln 1 ˆ ( x ) e ˆ ( y j ) e (9) (10) where ˆ can be obtaned as the soluton of non-lnear equaton g( ) 0 ( ) ( ) j y n x xe ye j n ( x ) ( y j ) n 1 1 e j1 1 e n ( x) ( y) e e ln 1 ln 1 1 1 n x y j 0 ( x) ( yj) 11 e j11 e (11) Therefore, ˆ s sple teratve soluton of non-lnear equaton g( ) 0. Once ˆ s nown, ˆ and ˆ can be obtaned fro (9) and (10) respectvely. Therefore, the MLE of R becoes! ˆ R ˆ ˆ ˆ ( j) where. ( )! ˆ s j0 1 (1) 371

ESTIMATION OF RELIABILITY IN STRESS-STRENGTH The asyptotc confdence nterval for R, s calculated as: Frst, the asyptotc varance of the MLE s gven by 1 1 L V( ˆ) and V( ˆ L E ) E n (13) The asyptotc varance (AV) of an estate of R whch a functon of two ndependent statstcs, for exaple,, s gven by Rao (1973). ˆ Rs, ˆ Rs, AV(R ˆ )=V( ) V( ) (14) Fro the asyptotc optu propertes of MLEs (Kendall & Stuart, 1979) and of lnear unbased estators (Davd, 1981), t s nown that MLEs are asyptotcally equally effcent havng the Craer-Rao lower bound as ther asyptotc varance as gven n (13). Thus, fro Equaton (14), the asyptotc varance of Rˆs, can be obtaned. To avod the dffculty of dervaton of R, the dervatves of R are obtaned for (s,)=(1,3) and (,4) separately, they are gven by R1,3 3 R1,3 3ˆ and. (3 ˆ ) (3 ˆ ) R,4 1(7 ˆ ) R,4 1 ˆ (7 ˆ ) and. (3 ˆ )(4 ˆ ) (3 ˆ )(4 ˆ ) Thus, ˆ 9 ˆ 1 1 AV( R1,3 )= 4 (3 ˆ ) n and AV( Rˆ )= ˆ ˆ. (3 ˆ )(4 ˆ ) n 144 (7 ) 1 1,4 4 As n,, Rˆ s, Rs, d N(0,1), AV( Rˆ ) 37

GADDE SRINIVASA RAO and the asyptotc 100(1 )% confdence nterval for R s gven by Rˆ Z AV( Rˆ ) s, (1 ) s,. The asyptotc 100(1 )% confdence nterval for R 1,3 s gven by Rˆ Z 1,3 (1 /) 3 ˆ 1 1 (3 ˆ ) n, where ˆ ˆ / ˆ. The asyptotc 100(1 )% confdence nterval for R,4 s gven by Rˆ Z 1 ˆ (7 ˆ ) 1 1, where ˆ ˆ / ˆ. (3 ˆ )(4 ˆ ) n,4 (1 /) Z where (1 /) s the (1 / ) th percentle of the standard noral dstrbuton. Sulaton Study and Data Analyss Sulaton Study Results based on Monte Carlo sulatons to copare the perforance of the R usng dfferent saple szes are presented. 3,000 rando saple of sze 10(5)35 each fro stress populaton, strength populaton were generated for (, ) = (3.0,1.0), (.5,1.0), (.0,1.0), (1.5,1.0), (1.0,1.0), (1.5,.0),(1.5,.5) and (1.5,3.0) on lnes of Bhattacharyya and Johnson (1974). The ML estators of and were then substtuted n to obtan the relablty n a ultcoponent stressstrength for (s, ) = (1, 3), (, 4). The average bas and average ean square error (MSE) of the relablty estates over the 3,000 replcatons are gven n Tables 1 and. Average confdence length and coverage probablty of the sulated 95% confdence ntervals of R are gven n Tables 3 and 4. The true value of relablty n ultcoponent stress- strength wth the gven cobnatons of (, ) for (s, ) = (1, 3) are 0.563, 0.600, 0.643, 0.69, 0.750, 0.800, 0.833, 0.857, 0.875 and for (s, ) = (, 4) are 0.355, 0.400, 0.454, 0.519, 0.600, 0.674, 0.75, 373

ESTIMATION OF RELIABILITY IN STRESS-STRENGTH 0.76, 0.790. Thus the true value of relablty n ultcoponent stress- strength ncreases as ncreases for a fxed whereas relablty n ultcoponent stress- strength decreases as ncreases for a fxed n both the cases of (s, ). Therefore, the true value of relablty s ncreases as decreases and vce-versa. The average bas and average MSE are decreases as saple sze ncreases for both (s, ). It verfes the consstency property of the MLE of R. Also the bas s negatve n both stuatons of (s, ). Whereas, aong the paraeters the absolute bas and MSE are ncreases as ncreases for a fxed n both the cases of (s, ) and the absolute bas and MSE are decreases as ncreases for a fxed n both the cases of (s, ). The average length of the confdence nterval also decreases as the saple sze ncreases. The coverage probablty s close to the nonal value n all cases but slghtly less than 0.95 n ost of the cobnatons. Overall, the perforance of the confdence nterval s good for all cobnatons of paraeters. Whereas, aong the paraeters observed, the sae phenoenon for average length and average coverage probablty were observed n the case of average bas and MSE. Table 1. Average bas of the sulated estates of R (, ) (s,) (n,) (3.5,1.5) (3.0,1.5) (.5,1.5) (.0,1.5) (1.5,1.5) (1.5,.0) (1.5,.5) (1.5,3.0) (1.5,3.5) (10,10) -0.0036-0.0187-0.01650-0.01357-0.00984-0.00649-0.0049-0.00360-0.0049 (15,15) -0.01517-0.0141-0.0168-0.01075-0.0085-0.0059-0.0048-0.00311-0.007 (1,3) (0,0) -0.00860-0.00773-0.00669-0.00548-0.00367-0.0077-0.00179-0.00114-0.00101 (5,5) -0.00851-0.00766-0.00657-0.0051-0.00357-0.0015-0.001-0.00060-0.00017 (30,30) -0.00679-0.00613-0.0058-0.0041-0.0090-0.00175-0.00098-0.00046-0.00010 (35,35) -0.00655-0.00610-0.00517-0.00413-0.0055-0.00147-0.00073-0.0001-0.00008 (10,10) -0.01113-0.01143-0.0114-0.0107-0.00819-0.00657-0.00539-0.0047-0.00348 (15,15) -0.00908-0.00945-0.00950-0.00903-0.00776-0.00601-0.00447-0.0034-0.009 (,4) (0,0) -0.00601-0.00599-0.00571-0.00508-0.00400-0.00176-0.0010-0.00089-0.00138 (5,5) -0.0051-0.00505-0.00473-0.00405-0.0093-0.00168-0.00073-0.00054-0.00045 (30,30) -0.0040-0.00416-0.00391-0.00338-0.0047-0.00144-0.00063-0.00043-0.00040 (35,35) -0.00386-0.00390-0.00377-0.0037-0.00160-0.0015-0.00058-0.00019-0.00035 374

GADDE SRINIVASA RAO Table. Average MSE of the sulated estates of R (, ) (s,) (n,) (3.5,1.5) (3.0,1.5) (.5,1.5) (.0,1.5) (1.5,1.5) (1.5,.0) (1.5,.5) (1.5,3.0) (1.5,3.5) (10,10) 0.0159 0.01437 0.01300 0.01109 0.00854 0.0067 0.00477 0.00375 0.00303 (15,15) 0.0105 0.00988 0.00894 0.00764 0.00590 0.00436 0.00334 0.0064 0.0014 (1,3) (0,0) 0.00713 0.00666 0.00599 0.00509 0.0039 0.0089 0.001 0.00175 0.00143 (5,5) 0.0059 0.00551 0.00494 0.00418 0.003 0.0036 0.00181 0.00144 0.00117 (30,30) 0.00460 0.0048 0.00383 0.0033 0.0048 0.0018 0.00139 0.00110 0.00090 (35,35) 0.0040 0.00374 0.00337 0.0086 0.000 0.0016 0.0015 0.00099 0.00081 (10,10) 0.01801 0.01877 0.01900 0.01831 0.01611 0.0130 0.01077 0.00889 0.00744 (15,15) 0.0185 0.01337 0.01351 0.0198 0.01140 0.0093 0.0076 0.00630 0.0059 (,4) (0,0) 0.00906 0.00936 0.00938 0.00895 0.00781 0.00635 0.00518 0.0049 0.00361 (5,5) 0.00754 0.00778 0.00778 0.00740 0.00643 0.005 0.0046 0.0035 0.0096 (30,30) 0.00594 0.0061 0.00610 0.00578 0.00500 0.00404 0.0038 0.0071 0.008 (35,35) 0.005 0.00538 0.00538 0.0051 0.00445 0.00361 0.0095 0.0044 0.0005 Table 3. Average confdence length of the sulated 95% confdence ntervals of R (, ) (s,) (n,) (3.5,1.5) (3.0,1.5) (.5,1.5) (.0,1.5) (1.5,1.5) (1.5,.0) (1.5,.5) (1.5,3.0) (1.5,3.5) (10,10) 0.4091 0.401 0.3880 0.363 0.34 0.763 0.400 0.113 0.1884 (15,15) 0.3399 0.3334 0.310 0.999 0.658 0.79 0.1981 0.1747 0.1559 (1,3) (0,0) 0.975 0.911 0.794 0.601 0.96 0.1961 0.1701 0.1497 0.1335 (5,5) 0.675 0.617 0.51 0.338 0.063 0.176 0.159 0.1346 0.101 (30,30) 0.453 0.398 0.300 0.140 0.1887 0.161 0.1398 0.13 0.1099 (35,35) 0.77 0.6 0.135 0.1986 0.1753 0.1498 0.1301 0.1146 0.103 (10,10) 0.4569 0.4719 0.480 0.4760 0.4498 0.4055 0.3637 0.374 0.966 (15,15) 0.3834 0.3953 0.401 0.3968 0.3739 0.3366 0.3019 0.719 0.466 (,4) (0,0) 0.3396 0.349 0.3533 0.3479 0.360 0.90 0.610 0.346 0.13 (5,5) 0.3056 0.3143 0.3180 0.3131 0.934 0.68 0.350 0.11 0.191 (30,30) 0.813 0.891 0.93 0.875 0.691 0.409 0.153 0.1935 0.175 (35,35) 0.611 0.684 0.714 0.670 0.500 0.40 0.003 0.1801 0.1631 375

ESTIMATION OF RELIABILITY IN STRESS-STRENGTH Table 4. Average coverage probablty of the sulated 95% confdence ntervals of R (, ) (s,) (n,) (3.5,1.5) (3.0,1.5) (.5,1.5) (.0,1.5) (1.5,1.5) (1.5,.0) (1.5,.5) (1.5,3.0) (1.5,3.5) (1,3) (,4) (10,10) 0.9090 0.9140 0.9193 0.980 0.9317 0.9303 0.967 0.960 0.967 (15,15) 0.910 0.9150 0.9187 0.913 0.950 0.973 0.990 0.987 0.953 (0,0) 0.9303 0.9347 0.9370 0.9390 0.9390 0.9377 0.9380 0.9347 0.9303 (5,5) 0.97 0.967 0.9317 0.9360 0.9383 0.9373 0.9357 0.9333 0.977 (30,30) 0.9353 0.9403 0.943 0.9463 0.9490 0.9463 0.9433 0.9400 0.9390 (35,35) 0.9317 0.9330 0.9353 0.9387 0.9387 0.9400 0.9363 0.9337 0.997 (10,10) 0.9113 0.9153 0.903 0.943 0.973 0.977 0.957 0.950 0.943 (15,15) 0.9103 0.9160 0.9190 0.9193 0.933 0.963 0.997 0.960 0.93 (0,0) 0.9310 0.9340 0.9370 0.9380 0.9367 0.9367 0.9360 0.933 0.987 (5,5) 0.937 0.987 0.9313 0.9353 0.9370 0.9350 0.933 0.9307 0.957 (30,30) 0.9357 0.9397 0.9437 0.9437 0.9477 0.9457 0.940 0.9403 0.9397 (35,35) 0.997 0.9317 0.9367 0.9360 0.9407 0.9393 0.9367 0.937 0.9300 Data Analyss Strength data, whch was orgnally reported by Badar and Prest (198), represents the strength easured n GPA for sngle carbon fbers and pregnated 1,000-carbon fber tows. Sngle fbers were tested under tenson at gauge lengths of 0 (Data Set I) and 10 (Data Set II), wth saple szes n = 69 and = 63 respectvely (see Data sets I and II). Data Set I (gauge lengths of 0 ). 1.31, 1.314, 1.479, 1.55, 1.700, 1.803, 1.861, 1.865, 1.944, 1.958, 1.966, 1.997,.006,.01,.07,.055,.063,.098,.140,.179,.4,.40,.53,.70,.7,.74,.301,.301,.359,.38,.38,.46,.434,.435,.478,.490,.511,.514,.535,.554,.566,.570,.586,.69,.633,.64,.648,.684,.697,.76,.770,.773,.800,.809,.818,.81,.848,.880,.809,.818,.81,.848,.880,.954, 3.01, 3.067, 3.084, 3.090, 3.096, 3.18, 3.33, 3.433, 3.585, 3.585. Data Set II (gauge lengths of 10 ). 1.901,.13,.03,.8,.57,.350,.361,.396,.397,.445,.454,.474,.518,.5,.55,.53,.575,.614,.616,.618,.64,.659,.675,.738,.740,.856,.917,.98,.937,.937,.977,.996, 3.030, 3.15, 3.139, 3.145, 3.0, 3.3, 3.35, 3.43, 3.64, 3.7, 3.94, 3.33, 3.346, 3.377, 3.408, 3.435, 3.493, 3.501, 3.537, 3.554, 3.56, 3.68, 3.85, 3.871, 3.886, 3.971, 4.04, 4.07, 4.5, 4.395, 5.00. 376

GADDE SRINIVASA RAO Surles and Padgett (1998, 001) observed that generalzed Raylegh wors well for strength data. Raqab and Kundu (005) analyzed the data by subtractng 1.0 and 1.8 fro the frst and second data set respectvely. The transfored data sets correspond to 0 and 10 gauge lengths are assued to follow GR(, ) and GR(, ) respectvely. The obtaned fnal estates for these two data sets are ˆ =.441, ˆ = 1.416, and ˆ = 0.8598. Also they checed the valdty of the odels usng the Kologorov-Srnov (K-S) tests for each data set. It was observed that for Data Sets I and II, the K-S dstances are 0.09 and 0.1 wth the correspondng p values of 0.6069 and 0.845 respectvely. It ndcates that the GR odel provdes reasonable ft to the transfored data sets. Based on estates of and the MLE of R becoe ˆR 1,3 = 0.63588 and ˆR,4 = 0.44484. The 95% confdence ntervals for R 1,3 becoe (0.55680, 0.71496) and for R,4 becoe (0.34387, 0.54581). Conclusons Ths artcle used real data sets to nvestgate ultcoponent stress-strength relablty for a generalzed Raylegh dstrbuton when both stress, strength varates follow the sae populaton. Asyptotc confdence ntervals for ultcoponent stress-strength relablty were estated usng the ML ethod. Sulaton results ndcate that the average bas and average MSE decreases as saple sze ncreases n both cases of (s, ). Aong the paraeters the absolute bas and MSE are ncreases (decreases) as ncreases ( ncreases) n both the cases of (s, ). The length of the confdence nterval also decreases as the saple sze ncreases and coverage probablty s close to the nonal value n all sets of paraeters consdered. References Awad, M. & Gharraf, K. (1986). Estaton of p(y<x) n Burr case: A coparatve study. Councatons n Statstcs - Sulatons & Cop., 15: 389-403. Badar, M. G. & Prest, A. M. (198). Statstcal aspects of fber and bundle strength n hybrd copostes. In T. Hayash, K. Kawata, and S. Ueawa (eds.), Progress n Scence and Engneerng Copostes, (pp. 119-1136). Toyo: ICCM-IV. 377

ESTIMATION OF RELIABILITY IN STRESS-STRENGTH Bhattacharyya, G. K. & Johnson, R. A. (1974). Estaton of relablty n ultcoponent stress strength odel. JASA, 69: 966-970. Davd, H. A. (1981). Order Statstcs. New Yor: John Wley and Sons. Downtown, F. (1973). The estaton of p(x>y) n the noral case. Technoetrcs, 15: 551-558. Ens, P. & Gesser, S. (1971). Estaton of the probablty that Y<X. JASA, 66: 16-168. Kendall, M. G. & Stuart, A. (1979). The Advanced Theory of Statstcs, (Vol. ). London: Charles Grffn and Copany Lted. Kundu, D. & Gupta, R. D. (005). Estaton of p(y<x) for the generalzed exponental dstrbuton. Metra, 61(3): 91-308. Kundu, D. & Gupta, R. D. (006). Estaton of p(y<x) for Webull dstrbuton, IEEE Transactons on Relablty, 55(): 70-80. Kundu, D. & Raqab, M. Z. (005). Generalzed Raylegh dstrbuton: dfferent ethods of estaton. Coputatonal Statstcs and Data Analyss, 49: 187-00. Kundu, D. and Raqab, M. Z. (009). Estaton of R= p(y<x) for threeparaeter Webull dstrbuton, Statstcs and Probablty Letters, 79: 1839-1846. McCool, J. I. (1991). Inference on p(y<x) n the Webull case, Councatons n Statstcs - Sulatons & Cop., 0: 19-148. Mudholar, G. S. & Srvastava, D. K. (1993). Exponentated Webull faly for analyzng bathtub falure data, IEEE Transactons on Relablty, 4: 99-30. Nand, S. B. & Ach, A. B. (1994). A note on estaton of p(x>y) for soe dstrbutons useful n lfe- testng, IAPQR Transactons, 19(1): 35-44. Pandey, M. & Uddn, B. (1985). Estaton of relablty n ultcoponent stress strength odel followng Burr dstrbuton. Proceedngs of the Frst Asan congress on Qualty and Relablty, (pp. 307 31). New Delh, Inda. Rao, C. R. (1973). Lnear Statstcal Inference and ts Applcatons. Inda: Wley Eastern Lted. Raqab, M. Z. & Kundu, D. (005). Coparson of dfferent estators of p(y<x) for a scaled Burr type X dstrbuton, Councatons n Statstcs Sulaton and Coputaton, 34(): 465-483. Raqab, M. Z., Mad, M. T. & Kundu, D. (008). Estaton of p(y<x) for the 3-paraeter generalzed exponental dstrbuton, Councatons n Statstcs - Theory and Methods, 37(18): 854-864. 378

GADDE SRINIVASA RAO Srnvasa Rao, G. & Kanta, R. R. L. (010). Estaton of relablty n ultcoponent stress strength odel: log-logstc dstrbuton, Electronc Journal of Appled Statstcal Analyss, 3(): 75-84. Surles, J. G. & Padgett, W. J. (1998). Inference for p(y<x) n the Burr Type X odel, Journal of Appled Statstcal Scences, 7: 5-38. Surles, J. G. & Padgett, W. J. (001). Inference for relablty and stressstrength for a scaled Burr Type X dstrbuton, Lfete Data Analyss, 7: 187-00. 379