Estimation of Stress- Strength Reliability model using finite mixture of exponential distributions

Similar documents
Functions of Random Variables

Lecture 2 - What are component and system reliability and how it can be improved?

Bayes Estimator for Exponential Distribution with Extension of Jeffery Prior Information

Analysis of Variance with Weibull Data

CHAPTER VI Statistical Analysis of Experimental Data

A New Family of Transformations for Lifetime Data

Comparison of Parameters of Lognormal Distribution Based On the Classical and Posterior Estimates

LECTURE - 4 SIMPLE RANDOM SAMPLING DR. SHALABH DEPARTMENT OF MATHEMATICS AND STATISTICS INDIAN INSTITUTE OF TECHNOLOGY KANPUR

Point Estimation: definition of estimators

BAYESIAN INFERENCES FOR TWO PARAMETER WEIBULL DISTRIBUTION

BAYESIAN ESTIMATOR OF A CHANGE POINT IN THE HAZARD FUNCTION

Summary of the lecture in Biostatistics

Comparing Different Estimators of three Parameters for Transmuted Weibull Distribution

VOL. 3, NO. 11, November 2013 ISSN ARPN Journal of Science and Technology All rights reserved.

Bayesian Inferences for Two Parameter Weibull Distribution Kipkoech W. Cheruiyot 1, Abel Ouko 2, Emily Kirimi 3

STK4011 and STK9011 Autumn 2016

Chapter 14 Logistic Regression Models

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

The Generalized Inverted Generalized Exponential Distribution with an Application to a Censored Data

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Multivariate Transformation of Variables and Maximum Likelihood Estimation

A Combination of Adaptive and Line Intercept Sampling Applicable in Agricultural and Environmental Studies

Introduction to local (nonparametric) density estimation. methods

Chapter 5 Properties of a Random Sample

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

Chapter 4 Multiple Random Variables

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Chapter 3 Sampling For Proportions and Percentages

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Analysis of a Repairable (n-1)-out-of-n: G System with Failure and Repair Times Arbitrarily Distributed

CHAPTER 3 POSTERIOR DISTRIBUTIONS

Confidence Intervals for Double Exponential Distribution: A Simulation Approach

Estimation of the Loss and Risk Functions of Parameter of Maxwell Distribution

A new Family of Distributions Using the pdf of the. rth Order Statistic from Independent Non- Identically Distributed Random Variables

Solving Constrained Flow-Shop Scheduling. Problems with Three Machines

ENGI 4421 Joint Probability Distributions Page Joint Probability Distributions [Navidi sections 2.5 and 2.6; Devore sections

A NEW MODIFIED GENERALIZED ODD LOG-LOGISTIC DISTRIBUTION WITH THREE PARAMETERS

Unimodality Tests for Global Optimization of Single Variable Functions Using Statistical Methods

Module 7. Lecture 7: Statistical parameter estimation

Some Notes on the Probability Space of Statistical Surveys

Comparison of Dual to Ratio-Cum-Product Estimators of Population Mean

Analysis of Lagrange Interpolation Formula

Dynamic Analysis of Axially Beam on Visco - Elastic Foundation with Elastic Supports under Moving Load

(Monte Carlo) Resampling Technique in Validity Testing and Reliability Testing

Modified Moment Estimation for a Two Parameter Gamma Distribution

Accelerated Life Test Sampling Plans under Progressive Type II Interval Censoring with Random Removals

Continuous Distributions

1 Solution to Problem 6.40

X ε ) = 0, or equivalently, lim

( ) = ( ) ( ) Chapter 13 Asymptotic Theory and Stochastic Regressors. Stochastic regressors model

Lecture 3. Sampling, sampling distributions, and parameter estimation

Maximum Likelihood Estimation

ECONOMETRIC THEORY. MODULE VIII Lecture - 26 Heteroskedasticity

THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA

Special Instructions / Useful Data

Multi Objective Fuzzy Inventory Model with. Demand Dependent Unit Cost and Lead Time. Constraints A Karush Kuhn Tucker Conditions.

LINEAR REGRESSION ANALYSIS

Bayes Interval Estimation for binomial proportion and difference of two binomial proportions with Simulation Study

1. The weight of six Golden Retrievers is 66, 61, 70, 67, 92 and 66 pounds. The weight of six Labrador Retrievers is 54, 60, 72, 78, 84 and 67.

Chapter 4 (Part 1): Non-Parametric Classification (Sections ) Pattern Classification 4.3) Announcements

BAYESIAN ESTIMATION OF GUMBEL TYPE-II DISTRIBUTION

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

Bootstrap Method for Testing of Equality of Several Coefficients of Variation

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

Some Statistical Inferences on the Records Weibull Distribution Using Shannon Entropy and Renyi Entropy

On generalized fuzzy mean code word lengths. Department of Mathematics, Jaypee University of Engineering and Technology, Guna, Madhya Pradesh, India

Lecture 8: Linear Regression

Optimal Strategy Analysis of an N-policy M/E k /1 Queueing System with Server Breakdowns and Multiple Vacations

ESTIMATION OF PARAMETERS OF THE MARSHALL-OLKIN EXTENDED LOG-LOGISTIC DISTRIBUTION FROM PROGRESSIVELY CENSORED SAMPLES

A NEW GENERALIZATION OF ERLANG DISTRIBUTION WITH BAYES ESTIMATION

Lecture 3 Probability review (cont d)

To use adaptive cluster sampling we must first make some definitions of the sampling universe:

Simulation Output Analysis

Analyzing Fuzzy System Reliability Using Vague Set Theory

IJOART. Copyright 2014 SciResPub.

Likelihood and Bayesian Estimation in Stress Strength Model from Generalized Exponential Distribution Containing Outliers

Econometric Methods. Review of Estimation

Lecture Notes Types of economic variables

Non-uniform Turán-type problems

MAX-MIN AND MIN-MAX VALUES OF VARIOUS MEASURES OF FUZZY DIVERGENCE

ENGI 3423 Simple Linear Regression Page 12-01

Chapter 13, Part A Analysis of Variance and Experimental Design. Introduction to Analysis of Variance. Introduction to Analysis of Variance

Estimation and Testing in Type-II Generalized Half Logistic Distribution

STATISTICAL PROPERTIES OF LEAST SQUARES ESTIMATORS. x, where. = y - ˆ " 1

Lecture 1. (Part II) The number of ways of partitioning n distinct objects into k distinct groups containing n 1,

Mean is only appropriate for interval or ratio scales, not ordinal or nominal.

Module 7: Probability and Statistics

Simple Linear Regression

The Mathematical Appendix

Lecture 9: Tolerant Testing

Bounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy

Miin-Jye Wen National Cheng Kung University, City Tainan, Taiwan, R.O.C. Key words: general moment, multivariate survival function, set partition

Chapter 8. Inferences about More Than Two Population Central Values

On the Bayesian Estimation for two Component Mixture of Maxwell Distribution, Assuming Type I Censored Data

Parameter, Statistic and Random Samples

Dependence In Network Reliability

A Study of the Reproducibility of Measurements with HUR Leg Extension/Curl Research Line

Chapter Business Statistics: A First Course Fifth Edition. Learning Objectives. Correlation vs. Regression. In this chapter, you learn:

ρ < 1 be five real numbers. The

TESTS BASED ON MAXIMUM LIKELIHOOD

Transcription:

Iteratoal Joural of Computatoal Egeerg Research Vol, 0 Issue, Estmato of Stress- Stregth Relablty model usg fte mxture of expoetal dstrbutos K.Sadhya, T.S.Umamaheswar Departmet of Mathematcs, Lal Bhadur college, P.G. Cetre, Waragal Departmet of Mathematcs, Kaatya Uversty, Waragal. ABSTRACT: I ths paper cosdered a stuato where stress ad stregth follow fte mxture of expoetal dstrbutos to fd the relablty of a system. It has bee studed whe stress follow expoetal dstrbuto ad stregth follow fte mxture of expoetal dstrbutos ad both stress-stregth follow the fte mxture of expoetal dstrbutos. The geeral expresso for the relablty of a system s obtaed. The relablty s computed umercally for dfferet values of the stress ad stregth parameters. We estmate the parameters of the relablty stress-stregth models by the method of maxmum lelhood estmato. The role of fte mxture of expoetal dstrbutos s llustrated usg a real lfe data o tme to death of two groups of leuaema patets. KEYWORDS: Expoetal dstrbuto, Fte mxture of expoetal dstrbutos, Maxmum Lelhood Estmato, Relablty, Stress-Stregth model.. INTRODUCTION Relablty of a system s the probablty that a system wll adequately perform ts teded purpose for a gve perod of tme uder stated evrometal codtos []. I some cases system falures occur due to certa type of stresses actg o them. Thus system composed of radom stregths wll have ts stregth as radom varable ad the stress appled o t wll also be a radom varable. A system fals wheever a appled stress exceeds stregth of the system. I a fte mxture model, the dstrbuto of radom quatty of terest s modelled as a mxture of a fte umber of compoet dstrbutos varyg proportos []. The flexblty ad hgh degree of accuracy of fte mxture models have bee the ma reaso for ther successful applcatos a wde rage of felds the bologcal physcal ad socal sceces. The estmato of relablty based o fte mxture of pareto ad beta dstrbutos was studed by Maya, T. Nar (007)[]. I relablty theory, the mxture dstrbutos are used for the aalyss of the falure tmes of a sample of tems of coheret laser used telecommucato etwor. I a expermet, oe hudred ad three laser devces were operated at a temperature of 70 degree Celsus utl all had faled. The expermet was ru loger tha oe year before all the devces had faled, because most of the devces were extremely relable. The sample thus cossts of two dstct populatos, oe wth a very short mea lfe ad oe wth a much loger mea lfe. Ths ca be cosdered as a example of a mxture of two expoetal dstrbutos wth probablty desty fucto of the form f ( x) p exp( x) ( p) exp( x), 0 p, 0,, The above model wll be useful to predct how log all maufactured lasers should be lfe tested to assure that the fal product cotaed o devce from the fat mortalty populato. I the preset paper we dscuss the statstcal aalyss of fte mxture of expoetal dstrbutos the cotext of relablty theory. We gve the defto ad propertes of the fte mxture of expoetal dstrbutos. We derve the relablty, whe the stregth X follows fte mxture of expoetal ad the stress Y taes expoetal ad fte mxture of expoetal. We dscuss estmato procedure for fte mxture of expoetal dstrbutos by the method of maxmum lelhood estmato ad also estmato of stress-stregth relablty. We llustrate the method for a real data o survval tmes of leuaema patets ad fally gve the cocluso. Iss 50-005 November 0 Page 9

Estmato Of Stress- Stregth Relablty Model... II. STATISTICAL MODEL: The assumptos tae ths model are () The radom varables X ad Y are depedet. () The values of stress ad stregth are o-egatve. If X deotes the stregth of the compoet ad Y s the stress mposed o t, the the relablty of the compoet s gve by [], x R P( X Y) g( y) dy f ( x) dx 00 () where f(x) ad g(y) are probablty desty fuctos of stregth ad stress respectvely. A fte mxture of expoetal dstrbuto wth -compoets ca be represeted the form f ( x) p f ( x) p f ( x)... p f ( x) () The r th where p 0,,,..., p momet of the mxture of two expoetal dstrbutos ( r r E x ) x p exp( x) p exp( x) dx Whe r =, Whe r =, 0 ( r) ( r) p ( p ) r r Ex ( ) p p p ( p ) Ex ( ) 0 p, 0,, Thus the varace s gve by p ( p ) ( p )( p ) p ( p ) V( x) I ths paper we are cosderg two cases. They are () Stress follows expoetal dstrbuto ad stregth follows fte mxture of expoetal dstrbutos. () Stress ad stregth follows fte mxture of expoetal dstrbutos. III. RELIABILITY COMPUTATIONS: Let X be the stregth of the -compoets wth probablty desty fuctos f (x); =,,..,. Stregth X follows fte mxture of expoetal dstrbuto wth pdf f ( x) p exp( x), 0, x 0, p 0,,,..., ; p..case() The stress Y follows expoetal dstrbuto: Whe the stress Y follows expoetal wth pdf g( y) exp( y), y 0, 0 As a specal case of () wth =, we have f ( x) p exp( x) ( p ) exp( x), 0, x 0,(,) Ad f X ad Y are depedet, the the relablty R from () x exp( ) exp( ) ( ) exp( ) R y p x p x dydx 00 () Iss 50-005 November 0 Page 40

Estmato Of Stress- Stregth Relablty Model... R p ( p ) (4) As a specal case of () wth =, we have f ( x) p exp( x) p exp( x) p exp( x), 0, x 0,(,,); p x exp( ) exp( ) exp( ) exp( ) R y p x p x p x dydx 00 R p p p (5) f ( x) p f ( x) p f ( x)... p f ( x) I geeral from (), we get R where p 0,,,..., p p (6) From table ad table ad fgs. ad, t s observed that f stress parameter creases the the value of relablty creases, f stregth parameter creases the the value of relablty decreases...case () The stress Y follows fte mxture of expoetal dstrbutos: For =, we have f ( x) p exp( x) ( p ) exp( x),, 0 g( y) p exp( x) ( p ) exp( x),, 0 4 4 4 Ad f X ad Y are depedet, the the relablty R from () x exp( ) ( ) exp( ) exp( ) ( ) exp( ) R p y p y p x p x dydx 4 4 00 p p p ( p ) ( p ) p ( p )( p ) For =, we have 4 4 f ( x) p exp( x) p exp( x) p exp( x), 0, x 0,(,,); p The x R p exp( y) p exp( y) p exp( y) 4 4 4 5 5 5 6 6 6 00 p exp( x) p exp( x) p exp( x) dydx R p p p p p p 4 4 4 4 4 4 (7) Iss 50-005 November 0 Page 4

R Estmato Of Stress- Stregth Relablty Model... p p p p p p p p p p p p 5 5 5 6 6 6 5 5 5 6 6 5 6 pj p j j I geeral from (), The R pj p j j f ( x) p f ( x) p f ( x)... p f ( x) where p 0,,,..., p From table ad table 4 ad fgs. ad 4, t s observed that f stress parameter creases the the value of relablty decreases, f stregth parameter creases the the value of relablty creases... Hazard Rate: Let t deotes lfe tme of a compoet wth survval fucto S(t). The the survval fucto of the model s obtaed as S( t) p exp( t) ( p ) exp( t), 0, t 0, 0 p For the model the hazard rate h(t) s gve by f() t ht () st () I geeral, p exp( t ) ( p ) exp( t) p exp( t) ( p )exp( t) ht () p exp( t) p exp( t) () Fte mxture of expoetal possesses decreasg hazard rate ad costat hazard rate depedg upo the values of the parameters. Fg 5, show the behavour of hazard rate at varous tme pots. IV. ESTIMATION OF PARAMETERS: We estmate the parameters of the models by the method of maxmum lelhood estmato.cosder the stuato whe there are oly two sub populatos wth mxg proportos p &(-p ) ad f (x) ad f (x) are expoetal destes wth parameters & respectvely. The lelhood fucto s gve by L(,, p / y) p exp( y j) ( p) exp( y j) j () Where y j deoted the falure tme of the j th ut belogg to the th sub populato j=,,... ; =, ad y y, y,... y ; y, y,... y Maxmzato of log lelhood fucto of () w.r.t the parameters yelds the followg equato L p ( p ) exp( y ) exp( y ) j j j j (8) (9) (0) Iss 50-005 November 0 Page 4

Estmato Of Stress- Stregth Relablty Model... L c p ( p ) exp y j y j j j log( L) log c log( p ) log( p ) y y j j j j 0, 0, 0 yj yj j j p p j y j j y j p where The above results ca be geeralzed for ay, gvg the followg estmators p & j y j () (4) (5) (6) (7) 4.. Estmato of Stress Stregth relablty: () Whe the stregth X follows fte mxture of expoetal dstrbutos wth parameters ad p ad the stress Y follows expoetal dstrbuto wth parameter, the the M.L.E of R s gve as R p (8) () Whe the stregth X follows fte mxture of expoetal dstrbutos wth parameters ad p ad the stress Y follows fte mxture of expoetal dstrbuto wth parameter j ad p j the the M.L.E of R s gve as R pj p j (9) V. DATA ANALYSIS: We cosder a data o tme to death of two groups of leuaema patets whch s gve Table 5 (see Fegl ad Zele, 965) to llustrate the procedure of estmato. We the estmate the parameters usg M.L.E techque. Table 6 provdes the values of the estmates by M.L.E method. Table 7 provdes the maxmum lelhood estmate of survval fucto at varous tme pots. Table 8 provdes the hazard rate fucto at varous tme pots. Iss 50-005 November 0 Page 4

Estmato Of Stress- Stregth Relablty Model... Case () Stress has expoetal dstrbuto ad Stregth has mxture two of expoetal dstrbutos: Table Table P = R P R 0. 0. 0. 0.7 0. 0.875 0. 0. 66667 0. 0.7 0. 0.777778 0. 0. 0.75 0. 0.7 0. 0.7 0. 0.4 0.8 0. 0.7 0.4 664 0. 0.8 0. 0.7 8 0. 0.8574 0. 0.7 846 0. 0.7 0.875 0. 0.7 0.7 0. 0.8 0.888889 0. 0.7 0.8 0.466667 0. 0.9 0.9 0. 0.7 0.9 0.475 0. 0.90909 0. 0.7 0.4765 P =P Case() Stress-Stregth has mxture two of expoetal dstrbutos: Table Table 4 4 R P =P = 4 R 0. 0. 0.7 0.875 0. 0. 0. 0. 0.7 0.777778 0. 0. 66667 0. 0. 0.7 0.7 0. 0. 0.75 0. 0.4 0.7 664 0. 0.4 0.8 0. 0.7 8 0. 0.8 0. 0.7 846 0. 0.8574 0. 0.7 0.7 0. 0.7 0.875 0. 0.8 0.7 0.466667 0. 0.8 0.888889 0. 0.9 0.7 0.475 0. 0.9 0.9 0. 0.7 0.4765 0. 0.90909 Table 5 Survval tmes of leuaema patets AG +ve AG ve 4 56 6 4 6 65 56 00 9 5 65 4 08 8 7 6 4 56 65 7 4 0 4 4 Table 6 Estmates of parameters of survval tmes of leuaema patets. 0.06 0.055 p 5 Iss 50-005 November 0 Page 44

Relablty Relablty Relablty Estmato Of Stress- Stregth Relablty Model... Table 7 Maxmum lelhood estmate of survval probablty at varous tme pots T 0 50 75 00 0 St () 0.9659 0.787 4 0 0.060 0.076 Table 8 Hazard rate at varous tme pots T 0. 0. 0. 0.4 0.7 0.8 0.9 h(t) 867 50 58 0.486 0.455 0.45 0.988 0.740 07 0.89 Case () Stress has expoetal dstrbuto ad Stregth has mxture two of expoetal dstrbutos: 0.95 Fgure 0.9 0.85 0.8 0.75 0.7 5 5 0. 0. 0. 0.4 0.7 0.8 0.9 Stress Parameter 0.9 Fgure 0.85 0.8 0.75 0.7 5 5 0.45 0.4 0. 0. 0. 0.4 0.7 0.8 0.9 Stregth Parameter Case() Stress-Stregth has mxture two of expoetal dstrbutos: 0.9 Fgure 0.85 0.8 0.75 0.7 5 5 0.45 0.4 0. 0. 0. 0.4 0.7 0.8 0.9 Stress parameter Iss 50-005 November 0 Page 45

Hazard rate Relablty Estmato Of Stress- Stregth Relablty Model... 0.95 Fgure 4 0.9 0.85 0.8 0.75 0.7 5 5 0. 0. 0. 0.4 0.7 0.8 0.9 Stregth parameter Hazard rate fucto: 5 Fgure 5 5 0.45 0.4 0. 0. 0. 0.4 0.7 0.8 0.9 Lfe tme VI. CONCLUSION: The role of fte mxture of expoetal dstrbutos relablty aalyss s studed. We derve the relablty, whe the stregth X follows fte mxture of expoetal ad the stress Y taes expoetal ad fte mxture of expoetal. It has bee observed by the computatos ad graphs, case() f stress parameter creases the the value of relablty creases, f stregth parameter creases the the value of relablty decreases. Where as case() f stress parameter creases the the value of relablty decreases, f stregth parameter creases the the value of relablty creases. We developed estmates of parameters usg Maxmum lelhood estmato. The role of fte mxture of expoetal dstrbutos s llustrated usg a real lfe data o tme to death of two groups of leuaema patets. REFERENCES: [] Kapur, K.C ad Lamberso, L.R.: Relablty Egeerg Desg, Joh Wley ad Sos, Ic., U.K. (997). [] Sha, S.K. (986). Relablty ad Lfe Testg, Wley Easter Lmted, New Delh. [] Maya, T.Nar (007). O a fte mxture of Pareto ad beta dstrbutos. Ph.D Thess submtted to Coch Uversty of Scece ad Techology, February 007. [4] Fegal, P. ad Zele, M. (965). Estmato of survval probabltes wth cocomtat formato. Bometrcs, 86-88. Iss 50-005 November 0 Page 46