A generalized Brown-Proschan model for preventive and corrective maintenance. Laurent DOYEN.
|
|
- Donna Reeves
- 5 years ago
- Views:
Transcription
1 A generalized Brown-Proschan model for preventive and corrective maintenance Jean Kuntzmann Laboratory Grenoble University France 1 of 28
2 I. Introduction The dependability of complex repairable systems depends strongly on the efficiency of maintenance actions. Corrective Maintenance (CM): After a failure, put the system into a state in which it can perform its function again. Preventive Maintenance (PM): When the system is operating, intends to slow down the wear process. Basic maintenance effect assumptions: As Bad As Old (ABAO): restores the system in the same state it was just before maintenance. As Good As New (AGAN): restores the system as if it were new. Reality is between these two extreme cases: Imperfect maintenance. 2 of 28
3 Brown-Proschan model [83], CM is : AGAN with probability p, ABAO with probability 1 p. repair effects (AGAN or ABAO) are mutually independent and independent of already observed failure times. CM effects can be characterized with random variables: { 1 if the i th CM is AGAN B i = 0 if the i th CM is ABAO Statistical studies when the {B i } i 1 are known: Whitaker-Samaniego[89], Hollander-Presnell-Sethuraman [92], Kvam-Singh- Whitaker [02], Bathe-Franz [96], Agustin-Pena[99],... 3 of 28
4 Practical purposes : {B i } i 0 are hidden variables. p = 0 : ABAO p characterizes maintenance efficiency: p = 1 : AGAN 0 < p < 1 : imperfect Joint assessment of CM efficiency and intrinsic wear-out: Lim [98]; Lim-Lie[00]; Lim-Lu-Park[98]; Langseth-Lindqvist [04]. Presentation aim: Generalize the BP model to PM effects Assess PM efficiency and intrinsic wear-out when PM effects are unknown Compute reliability indicators. 4 of 28
5 II/ A maintenance data set PM and CM times of a subsystem of a fossil-fired thermal plant of EDF: Cumulative number of failures 5 of 28
6 III. Notations 0 PM PM c CM CM PM CM T m t C 1 C 2 C 3 C 4 u 1 = 0 u 2 = 0 u 3 = 1 u 4 = 0 τ 1 τ 2 τ 3 N t T 1 T 2 T 3 6 of 28
7 IV. Model assumptions λ(t) the failure rate of the new unmaintained system. Λ(t) = t 0 λ(s) ds, Maintenance durations are not taken into account, PM are done at deterministic times, all the PM times are known CM are done at random times, CM are observed over [c, T ] (c < T ) CM effects are ABAO, PM effects follow a Brown-Proschan model, i.e. PM effects are : mutually independent and independent of previous failure times, AGAN with probability p, ABAO with probability 1 p. B i = { 1 if the i th PM is AGAN 0 if the i th PM is ABAO P (B i = 1 T i, B i 1 ) = P (B i = 1) = p 7 of 28
8 V. Joint assessment of intrinsic wear-out and PM efficiency a/ Maximum likelihood method Likelihood associated to a single observation of the failure process over [c, t]: L t (θ) = f T nt N t =n t (t nt ) P (N t = n t ) Maximum Likelihood Estimator (MLE): ˆθ = arg max L T (θ) = arg max log(l T (θ)) θ θ Notations: L c (θ) = 1 L τ,t (θ) = f T nt nτ,n t τ=n t n τ (t nτ +1 τ,..., t nt τ) the likelihood associated to the system new in τ and observed over[c τ, t]. 8 of 28
9 Let us denote : D m t = (1 p) m t m c τ m t j t λ(t j τ m ) Two equivalent equations for likelihood recursive computation: L t (θ) = m t m=0 p 1l {m 0} D m t forward computation algorithm. L t (θ) = e (Λ(t τ m) Λ((c τ m ) τ m )) L c τm (θ) p 1l {τ t} Dτ 0 m e (Λ(τ m c) Λ(c)) L τm,t(θ) τ {τ 1,...,τ mt,t} backward computation algorithm. 9 of 28
10 Proof of the equation used for the forward computation: {B m = 1, {B j = 0} m<j mt }{{} τ m is the last AGAN PM since t } 0 m mt is a partition of the probability space. Given B m = 1, the CM process can be divided into 2 independent processes. L t (θ) = [ m t p 1l {m 0} (1 p) m t m }{{} m=0 P (B m =1,{B j =0} m<j mt ) c τ m t j t λ(t j τ m ) e (Λ(t τ m) Λ((c τ m ) τ m )) }{{} Given B m =1,{B j =0} m<j mt, CM times after τ m follow a NHPP initialized in τ m ] L c τm (θ) }{{} Given B m =1, maintenances before τ m follow a BP PM-ABAO CM model 10 of 28
11 b/ MLE combined with moment estimation : λ(t) = αβt β 1 Power Law Process: E[N T ] = α(t β c β ) α β = N T /(T β c β ). BP model, E[N t ] = α S t, where S t = m τ τ {τ m c+1,...,τ mt,t} k=0 p 1l {k 0} (1 p) m τ k ((τ τ k ) β ((τ mτ c) τ k ) β ) MLE combined with moment estimation of α: α β,p = N T S T E[ α β,p ] = α ( β, p) = arg max (β,p) log(l T ( α β,p, β, p)) and α = α β, p Dimension reduction of the likelihood maximization space. 11 of 28
12 c/ Individual PM efficiency assessment p characterized the average global PM efficiency. The m th PM effect can be characterized by: π θ m = P (B m = 1 N T = n T, T nt = t nt ) which verifies: π θ m = p L c τ m (θ) L τm,t (θ) L T (θ) It can naturally be estimated by: π ˆθ m L τm (ˆθ): intermediate computing values of the forward algorithm L τm,t (ˆθ): intermediate computing values of the backward algorithm 12 of 28
13 d/ Expectation-Maximization algorithm Complete likelihood: L c t (θ) = f T nt N t =n t,b nt =b nt (t nt ) P (N t = n t B nt = b nt )P (B nt = b nt ) EM algorithm: Expectation (E) step: Compute Q(θ θ k ) = E θ k [log(l c t (θ)) N] Maximization (M) step: θ k+1 = arg max θ θ k k local likelihood maxima Q(θ θ k ) 13 of 28
14 Complete likelihood: L c t (θ) = N t n=1 λ(a tn ) e t c λ(a s ) ds m t m=1 p B m (1 p) 1 B m where A s is the virtual age, ie the time elapsed since the last perfect PM m s ]τ m, τ m+1 ] [c, T ] λ(a s ) = [λ(s τ m i )] 1l {Bτm i } i=0 where Bτ i m = at time τ m, the last AGAN PM time is τ m i 14 of 28
15 Complete likelihood: E θ [log(l c t (θ)) N]=Q(p θ) {}}{ L c t (θ) = N t n=1 λ(a tn ) e t c λ(a s ) ds }{{} E θ [log(. ) N]=Q 2 (λ θ) m t m=1 p B m (1 p) 1 B m }{{} E θ [log(. ) N]=Q 1 (p θ) where A s is the virtual age, ie the time elapsed since the last perfect PM m s ]τ m, τ m+1 ] [c, T ] λ(a s ) = [λ(s τ m i )] 1l {Bτm i } i=0 where Bτ i m = at time τ m, the last AGAN PM time is τ m i 15 of 28
16 Complete likelihood: E θ [log(l c t (θ)) N]=Q(p θ) {}}{ L c t (θ) = N t n=1 λ(a tn ) e t c λ(a s ) ds }{{} E θ [log(. ) N]=Q 2 (λ θ) m t m=1 p B m (1 p) 1 B m }{{} E θ [log(. ) N]=Q 1 (p θ) where A s is the virtual age, ie the time elapsed since the last perfect PM m s ]τ m, τ m+1 ] [c, T ] λ(a s ) = [λ(s τ m i )] 1l {Bτm i } i=0 where Bτ i m = at time τ m, the last AGAN PM time is τ m i Q 2 function of πm,i [1l θ = Eθ {B i τm } N T = n T, T nt = t nt [ ] Q 1 function of πm,0 θ = πθ m = E θ B m N T = n T, T nt = t nt ] 16 of 28
17 E step: Compute for 0 m m T and 0 i m, π 0,0 (θ) = 1 π θ m,i = Eθ [1l {B i τm } N T = n T, T nt = t nt ] for m {1,..., m T }, π θ m,0 = πθ m for 0 m < m T and 0 i m, π θ m+1,i+1 = πθ m,i p1+1l {m>i} (1 p) i (c τ m i )<t j τ m+1 λ(t j τ m i ) e (Λ((c τ m+1) τ m i ) Λ((c τ m i ) τ m i )) L c τ m i (θ)l τm+1,t (θ) L T (θ) 17 of 28
18 M step: λ(t) = αβt β 1 p k+1 = β k+1 = arg max β [ mt ] m=1 πθ k m,0 /m T τ {τ m c+1,...,τ mt,t } with S k (β) = [ n T (log(n T β) log(s k (β))) + (β 1) m τ i=0 τ {τ m c+1,...,τ mt,t } π θ k m τ,m τ i m τ i=0 (c τ mτ )<t j τ π θ k m τ,m τ i ] log(t j τ i ) ((τ τ i ) β ((c τ mτ ) τ i) β ) α k+1 = n T /S k (β k+1 ) 18 of 28
19 VI. Reliability indicators Failure intensity: λ t = lim t 0 1 t P (N t+ t N t T Nt, N t ) λ t = m t m=0 Dm C Kt e (Λ(t τ m) Λ((c τ m ) τ m ) L c τm (θ) m t m=0 Dm C Kt λ(t τ m ) e (Λ(t τ m) Λ((c τ m ) τ m ) L c τm (θ) Cumulative failure intensity: Λ t = t 0 λ s ds K t Λ t = log L Ck (θ) log L t (θ) L C + (θ) L k 1 C + (θ) Kt k=1 19 of 28
20 Reliability: P (T NT +1 > s T NT, N T ) = exp( (Λ s Λ T )) Expected cumulative number of failures: E[N s T NT, N T ] = N T + + m s i=m T +1 [ mt (Λ(s τ i ) + E[N τ i T NT, N T ] N T )p(1 p) m s i ] (Λ(s τ i ) Λ(T τ i ))(1 p) m s m T πm θ T,m T i i=0 20 of 28
21 VII. Application The BP PM-ABAO CM model is implemented in MARS (Maintenance Assessment of Repairable Systems): a free software developped by LJK (Grenoble university) and EDF (French electricity utility). 21 of 28
22 MLE combined with moment estimation: α = , β = 1.95, p = 0.602, log(lt ( α, β, p)) = MLE computed with EM algorithm or direct maximization: ˆα = , ˆβ = 1.94, ˆp = 0.614, log(lt (ˆα, ˆβ, ˆp)) = π ˆθ m : Wearing out state at time c failures surge ABAO ABAO AGAN AGAN AGAN 22 of 28
23 Quality of estimation: θ = (1.9e 6, 1.9, 0.61) } NEB: Normalized Empirical Bias estimated over NESD: Norm. Emp. Standard Deviation simulations θ : MLE estimator when the B i are known Initial intensity: λ(t) = αβt β 1 or λ(t) = β/η(t/η) β 1 EM is more robust than direct likelihood maximization NEB NESD α η β p α η β p θ 590% 36% 8.2% -9.4% 3300% 130% 35% 50% ˆθ 580% 37% 9.2% -12% 3300% 130% 36% 51% θ 18000% -24% -12% 0% 5000% 35% 22% 25% 23 of 28
24 Quality of individual PM efficiency estimation: Empirical mean Empirical standard deviation m th PM π θ m B m p B m π θ m B m p B m e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e-1 Mean -7.5 e e e e-1 Individual PM efficiency estimation is relevent after c 24 of 28
25 Failure intensity for θ = ˆθ Cumulative failure intensity and number of failures for θ = ˆθ of 28
26 Forecast reliability for θ = ˆθ Forecast cumulative number of failures for θ = ˆθ x x of 28
27 VIII. Prospects Maintenance times optimization Consider a more general distribution over ], 1] for the B i. Develop confidence intervals and tests for the BP model of 28
28 References Agustin MZN and Peˆna EA. Order statistics properties, random generation, and goodness-of-fit testing for a minimal repair model. Journal of American Statistical Association 1999; 94: Bathe F and Franz J. Modeling of repairable systems with various degrees of repair. Metrika 1996; 43: Brown M and Proschan F. Imperfect repair. Journal of Applied Probability 1983; 20: Kumar U and Klefsjö B. Reliability analysis of hydraulic systems of LHD machines using the power law process model. Reliability Engineering and System Safety 1992; 35: Kvam PH, Singh H and Whitaker LR. Estimating distributions with increasing failure rate in an imperfect repair model. Lifetime Data Analysis 2002; 8: Hollander M, Presnell P and Sethuraman J. Nonparametric methods for imperfect repair models. The Annals of Statistics 1992; 20: Langseth H and Lindqvist BH. A maintenance model for components exposed to several failure mechanisms and imperfect repair. In Mathematical and Statistical Methods in Reliability, B.H. Lindqvist, K.L. Doksum (eds). World Scientific: 2004; Lim JH, Lu KL, Park DH. Bayesian imperfect repair model. Communications in Statistics - Theory and Methods 1998; 27: Lim TJ. Estimating system reliability with fully masked data under Brown-Proschan imperfect repair model. Reliability Engineering and System Safety 1998; 59: Lim TJ and Lie CH. Analysis of system reliability with dependent repair modes. IEEE Transactions on reliability 2000; 49: Whitaker LR and Samaniego FJ. Estimating the reliability of systems subject to imperfect repair. American Statistical Association 1989; 84: Journal of 28 of 28
Time-varying failure rate for system reliability analysis in large-scale railway risk assessment simulation
Time-varying failure rate for system reliability analysis in large-scale railway risk assessment simulation H. Zhang, E. Cutright & T. Giras Center of Rail Safety-Critical Excellence, University of Virginia,
More informationCHAPTER 1 A MAINTENANCE MODEL FOR COMPONENTS EXPOSED TO SEVERAL FAILURE MECHANISMS AND IMPERFECT REPAIR
CHAPTER 1 A MAINTENANCE MODEL FOR COMPONENTS EXPOSED TO SEVERAL FAILURE MECHANISMS AND IMPERFECT REPAIR Helge Langseth and Bo Henry Lindqvist Department of Mathematical Sciences Norwegian University of
More informationCHAPTER 1. Stochastic modelling of the effect of preventive and corrective. maintenance on repairable systems reliability
CHAPTER 1 Stochastic modelling of the effect of preventive and corrective maintenance on repairable systems reliability Laurent Doyen and Olivier Gaudoin Institut National Polytechnique de Grenoble Laboratoire
More informationXLVII SIMPÓSIO BRASILEIRO DE PESQUISA OPERACIONAL
Static and Dynamics Maintenance Policies for Repairable Systems under Imperfect Repair Models Maria Luíza Guerra de Toledo Escola Nacional de Ciências Estatísticas - IBGE Rua André Cavalcanti, 106 - Santa
More informationOptimal maintenance for minimal and imperfect repair models
Optimal maintenance for minimal and imperfect repair models Gustavo L. Gilardoni Universidade de Brasília AMMSI Workshop Grenoble, January 2016 Belo Horizonte Enrico Colosimo Marta Freitas Rodrigo Padilha
More informationBasic math for biology
Basic math for biology Lei Li Florida State University, Feb 6, 2002 The EM algorithm: setup Parametric models: {P θ }. Data: full data (Y, X); partial data Y. Missing data: X. Likelihood and maximum likelihood
More informationDetermination of the optimal periodic maintenance policy under imperfect repair assumption
Universidade Federal de Minas Gerais Escola de Engenharia Departamento de Engenharia de Produção Programa de Pós-Graduação em Engenharia de Produção Determination of the optimal periodic maintenance policy
More informationStatic and Dynamic Maintenance Policies for Repairable Systems under Imperfect Repair Models
Static and Dynamic Maintenance Policies for Repairable Systems under Imperfect Repair Models Gustavo L. Gilardoni (UnB) Joint work with Maria Luiza G. de Toledo (ENCE-IBGE), Enrico Colosimo (UFMG) e Marta
More informationIntroduction to repairable systems STK4400 Spring 2011
Introduction to repairable systems STK4400 Spring 2011 Bo Lindqvist http://www.math.ntnu.no/ bo/ bo@math.ntnu.no Bo Lindqvist Introduction to repairable systems Definition of repairable system Ascher and
More informationA nonparametric monotone maximum likelihood estimator of time trend for repairable systems data
A nonparametric monotone maximum likelihood estimator of time trend for repairable systems data Knut Heggland 1 and Bo H Lindqvist Department of Mathematical Sciences Norwegian University of Science and
More informationA novel repair model for imperfect maintenance
IMA Journal of Management Mathematics (6) 7, 35 43 doi:.93/imaman/dpi36 Advance Access publication on July 4, 5 A novel repair model for imperfect maintenance SHAOMIN WU AND DEREK CLEMENTS-CROOME School
More informationThe information complexity of sequential resource allocation
The information complexity of sequential resource allocation Emilie Kaufmann, joint work with Olivier Cappé, Aurélien Garivier and Shivaram Kalyanakrishan SMILE Seminar, ENS, June 8th, 205 Sequential allocation
More informationNon-Parametric Tests for Imperfect Repair
Non-Parametric Tests for Imperfect Repair GUSTAVO GILARDONI Departamento de Estatística, UnB Prédio CIC/EST, Campus Darcy Ribeiro Brasília - DF, 70910-900 gilardon@gmail.com ENRICO COLOSIMO Departamento
More informationExamination paper for TMA4275 Lifetime Analysis
Department of Mathematical Sciences Examination paper for TMA4275 Lifetime Analysis Academic contact during examination: Ioannis Vardaxis Phone: 95 36 00 26 Examination date: Saturday May 30 2015 Examination
More informationA Model of Human Capital Accumulation and Occupational Choices. A simplified version of Keane and Wolpin (JPE, 1997)
A Model of Human Capital Accumulation and Occupational Choices A simplified version of Keane and Wolpin (JPE, 1997) We have here three, mutually exclusive decisions in each period: 1. Attend school. 2.
More informationReliability Growth in JMP 10
Reliability Growth in JMP 10 Presented at Discovery Summit 2012 September 13, 2012 Marie Gaudard and Leo Wright Purpose of Talk The goal of this talk is to provide a brief introduction to: The area of
More informationModeling of Dependence Between Critical Failure and Preventive Maintenance: The Repair Alert Model
Modeling of Dependence Between Critical Failure and Preventive Maintenance: The Repair Alert Model Bo Lindqvist Helge Langseth Bård Støve bo@math.ntnu.no Department of Mathematical Sciences Norwegian University
More informationBAYESIAN MODELING OF DYNAMIC SOFTWARE GROWTH CURVE MODELS
BAYESIAN MODELING OF DYNAMIC SOFTWARE GROWTH CURVE MODELS Zhaohui Liu, Nalini Ravishanker, University of Connecticut Bonnie K. Ray, IBM Watson Research Center Department of Mathematical Sciences, IBM Watson
More informationLECTURE 5 NOTES. n t. t Γ(a)Γ(b) pt+a 1 (1 p) n t+b 1. The marginal density of t is. Γ(t + a)γ(n t + b) Γ(n + a + b)
LECTURE 5 NOTES 1. Bayesian point estimators. In the conventional (frequentist) approach to statistical inference, the parameter θ Θ is considered a fixed quantity. In the Bayesian approach, it is considered
More informatione 4β e 4β + e β ˆβ =0.765
SIMPLE EXAMPLE COX-REGRESSION i Y i x i δ i 1 5 12 0 2 10 10 1 3 40 3 0 4 80 5 0 5 120 3 1 6 400 4 1 7 600 1 0 Model: z(t x) =z 0 (t) exp{βx} Partial likelihood: L(β) = e 10β e 10β + e 3β + e 5β + e 3β
More informationThe Expectation-Maximization Algorithm
1/29 EM & Latent Variable Models Gaussian Mixture Models EM Theory The Expectation-Maximization Algorithm Mihaela van der Schaar Department of Engineering Science University of Oxford MLE for Latent Variable
More informationReview of Maximum Likelihood Estimators
Libby MacKinnon CSE 527 notes Lecture 7, October 7, 2007 MLE and EM Review of Maximum Likelihood Estimators MLE is one of many approaches to parameter estimation. The likelihood of independent observations
More informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
More informationOptimal goodness-of-fit tests for recurrent event data
Lifetime Data Anal (211) 17:49 432 DOI 1.17/s1985-11-9193-1 Optimal goodness-of-fit tests for recurrent event data Russell S. Stocker IV Akim Adekpedjou Received: 13 October 29 / Accepted: 14 February
More informationEstimation of reliability parameters from Experimental data (Parte 2) Prof. Enrico Zio
Estimation of reliability parameters from Experimental data (Parte 2) This lecture Life test (t 1,t 2,...,t n ) Estimate θ of f T t θ For example: λ of f T (t)= λe - λt Classical approach (frequentist
More informationSoftware Reliability Growth Modelling using a Weighted Laplace Test Statistic
Software Reliability Growth Modelling using a Weighted Laplace Test Statistic Yan Luo Torsten Bergander A. Ben Hamza Concordia Institute for Information Systems Engineering Concordia University, Montréal,
More informationA STUDY OF ASYMPTOTIC AVAILABILITY MODELING FOR A FAILURE AND A REPAIR RATES FOLLOWING A WEIBULL DISTRIBUTION
A STUDY OF ASYMPTOTIC AVAILABILITY MODELING FOR A FAILURE AND A REPAIR RATES FOLLOWING A WEIBULL DISTRIBUTION Salem Bahri a, Fethi Ghribi b, Habib Ben Bacha a,c a Electro Mechanical systems laboratory
More informationChapter 9 Part II Maintainability
Chapter 9 Part II Maintainability 9.4 System Repair Time 9.5 Reliability Under Preventive Maintenance 9.6 State-Dependent Systems with Repair C. Ebeling, Intro to Reliability & Maintainability Chapter
More informationSolutions For Stochastic Process Final Exam
Solutions For Stochastic Process Final Exam (a) λ BMW = 20 0% = 2 X BMW Poisson(2) Let N t be the number of BMWs which have passes during [0, t] Then the probability in question is P (N ) = P (N = 0) =
More informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationParameters Estimation for a Linear Exponential Distribution Based on Grouped Data
International Mathematical Forum, 3, 2008, no. 33, 1643-1654 Parameters Estimation for a Linear Exponential Distribution Based on Grouped Data A. Al-khedhairi Department of Statistics and O.R. Faculty
More informationCHAPTER 1 COMPETING RISK MODELING IN RELIABILITY TIM BEDFORD
CHAPTER 1 COMPETING RISK MODELING IN RELIABILITY TIM BEDFORD Department of Management Science Strathclyde University Glasgow, UK E-mail: tim.bedford@strath.ac.uk This paper gives a review of some work
More informationA STATISTICAL TEST FOR MONOTONIC AND NON-MONOTONIC TREND IN REPAIRABLE SYSTEMS
A STATISTICAL TEST FOR MONOTONIC AND NON-MONOTONIC TREND IN REPAIRABLE SYSTEMS Jan Terje Kvaløy Department of Mathematics and Science, Stavanger University College, P.O. Box 2557 Ullandhaug, N-491 Stavanger,
More informationÚstav teorie informace a automatizace. Academy of Sciences of the Czech Republic Institute of Information Theory and Automation RESEARCH REPORT
Akademie věd České republiky Ústav teorie informace a automatizace Academy of Sciences of the Czech Republic Institute of Information Theory and Automation RESEARCH REPORT Jaroslav Ševčík and Petr Volf
More informationMATHEMATICAL AND STATISTICAL METHODS IN RELIABILITY. Bo H. Lindqvist Kjell A. Doksum Editors
MATHEMATICAL AND STATISTICAL METHODS IN RELIABILITY Bo H. Lindqvist Kjell A. Doksum Editors 1 CONTENTS PREFACE CONTRIBUTORS xi xv I RELIABILITY THEORY IN THE PAST AND PRESENT CENTURIES 1 1 MATHEMATICAL
More informationQualifying Exam in Probability and Statistics. https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf
Part : Sample Problems for the Elementary Section of Qualifying Exam in Probability and Statistics https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf Part 2: Sample Problems for the Advanced Section
More informationMaintenance of Repairable Systems
1 Maintenance of Repairable Systems Bo Henry Lindqvist 1.1 Introduction A commonly used definition of a repairable system (Ascher and Feingold [3]) states that this is a system which, after failing to
More informationSection 7.5: Nonstationary Poisson Processes
Section 7.5: Nonstationary Poisson Processes Discrete-Event Simulation: A First Course c 2006 Pearson Ed., Inc. 0-13-142917-5 Discrete-Event Simulation: A First Course Section 7.5: Nonstationary Poisson
More informationImprecise Software Reliability
Munich, September 2009 Denitions and assumptions of the software reliability models Well-known models Standard denitions of software reliability A fault in a software is an incorrect step, process, or
More informationModelling and Analysis of Recurrent Event Data
Modelling and Analysis of Recurrent Event Data Edsel A. Peña Department of Statistics University of South Carolina Research support from NIH, NSF, and USC/MUSC Collaborative Grants Joint work with Prof.
More informationIntroduction to queuing theory
Introduction to queuing theory Claude Rigault ENST claude.rigault@enst.fr Introduction to Queuing theory 1 Outline The problem The number of clients in a system The client process Delay processes Loss
More informationChap 4. Software Reliability
Chap 4. Software Reliability 4.2 Reliability Growth 1. Introduction 2. Reliability Growth Models 3. The Basic Execution Model 4. Calendar Time Computation 5. Reliability Demonstration Testing 1. Introduction
More informationSequence labeling. Taking collective a set of interrelated instances x 1,, x T and jointly labeling them
HMM, MEMM and CRF 40-957 Special opics in Artificial Intelligence: Probabilistic Graphical Models Sharif University of echnology Soleymani Spring 2014 Sequence labeling aking collective a set of interrelated
More informationBayesian Nonparametric Inference Methods for Mean Residual Life Functions
Bayesian Nonparametric Inference Methods for Mean Residual Life Functions Valerie Poynor Department of Applied Mathematics and Statistics, University of California, Santa Cruz April 28, 212 1/3 Outline
More informationProblem Selected Scores
Statistics Ph.D. Qualifying Exam: Part II November 20, 2010 Student Name: 1. Answer 8 out of 12 problems. Mark the problems you selected in the following table. Problem 1 2 3 4 5 6 7 8 9 10 11 12 Selected
More informationDegradation data analysis for samples under unequal operating conditions: a case study on train wheels
Degradation data analysis for samples under unequal operating conditions: a case study on train wheels Marta AFONSO FREITAS UNIVERSIDADE FEDERAL DE MINAS GERAIS DEPARTMENT OF PRODUCTION ENGINEERING marta.afreitas@gmail.com
More informationCHAPTER 1 STATISTICAL MODELLING AND INFERENCE FOR COMPONENT FAILURE TIMES UNDER PREVENTIVE MAINTENANCE AND INDEPENDENT CENSORING
CHAPTER 1 STATISTICAL MODELLING AND INFERENCE FOR COMPONENT FAILURE TIMES UNDER PREVENTIVE MAINTENANCE AND INDEPENDENT CENSORING Bo Henry Lindqvist and Helge Langseth Department of Mathematical Sciences
More informationSpring 2017 Econ 574 Roger Koenker. Lecture 14 GEE-GMM
University of Illinois Department of Economics Spring 2017 Econ 574 Roger Koenker Lecture 14 GEE-GMM Throughout the course we have emphasized methods of estimation and inference based on the principle
More informationEmpirical Bayes Estimation of Reliability
Empirical Bayes Estimation of Reliability Introduction Assessment of the reliability of various types of equipment relies on statistical inference about characteristics of reliability such as reliability
More informationLecture 25: Review. Statistics 104. April 23, Colin Rundel
Lecture 25: Review Statistics 104 Colin Rundel April 23, 2012 Joint CDF F (x, y) = P [X x, Y y] = P [(X, Y ) lies south-west of the point (x, y)] Y (x,y) X Statistics 104 (Colin Rundel) Lecture 25 April
More informationA Higher-Order Interactive Hidden Markov Model and Its Applications Wai-Ki Ching Department of Mathematics The University of Hong Kong
A Higher-Order Interactive Hidden Markov Model and Its Applications Wai-Ki Ching Department of Mathematics The University of Hong Kong Abstract: In this talk, a higher-order Interactive Hidden Markov Model
More informationSTAT 512 sp 2018 Summary Sheet
STAT 5 sp 08 Summary Sheet Karl B. Gregory Spring 08. Transformations of a random variable Let X be a rv with support X and let g be a function mapping X to Y with inverse mapping g (A = {x X : g(x A}
More informationReliability Analysis of Tampered Failure Rate Load-Sharing k-out-of-n:g Systems
Reliability Analysis of Tampered Failure Rate Load-Sharing k-out-of-n:g Systems Suprasad V. Amari Relex Software Corporation 540 Pellis Road Greensburg, PA 15601 USA Krishna B. Misra RAMS Consultants 71
More informationON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION
ON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION ZHENLINYANGandRONNIET.C.LEE Department of Statistics and Applied Probability, National University of Singapore, 3 Science Drive 2, Singapore
More informationRepairable Systems Reliability Trend Tests and Evaluation
Repairable Systems Reliability Trend Tests and Evaluation Peng Wang, Ph.D., United Technologies Research Center David W. Coit, Ph.D., Rutgers University Keywords: repairable system, reliability trend test,
More informationAn Integral Measure of Aging/Rejuvenation for Repairable and Non-repairable Systems
An Integral Measure of Aging/Rejuvenation for Repairable and Non-repairable Systems M.P. Kaminskiy and V.V. Krivtsov Abstract This paper introduces a simple index that helps to assess the degree of aging
More informationSTAT215: Solutions for Homework 2
STAT25: Solutions for Homework 2 Due: Wednesday, Feb 4. (0 pt) Suppose we take one observation, X, from the discrete distribution, x 2 0 2 Pr(X x θ) ( θ)/4 θ/2 /2 (3 θ)/2 θ/4, 0 θ Find an unbiased estimator
More informationDynamic Asset Allocation - Identifying Regime Shifts in Financial Time Series to Build Robust Portfolios
Downloaded from orbit.dtu.dk on: Jan 22, 2019 Dynamic Asset Allocation - Identifying Regime Shifts in Financial Time Series to Build Robust Portfolios Nystrup, Peter Publication date: 2018 Document Version
More informationLecture 22 Survival Analysis: An Introduction
University of Illinois Department of Economics Spring 2017 Econ 574 Roger Koenker Lecture 22 Survival Analysis: An Introduction There is considerable interest among economists in models of durations, which
More informationToday. Calculus. Linear Regression. Lagrange Multipliers
Today Calculus Lagrange Multipliers Linear Regression 1 Optimization with constraints What if I want to constrain the parameters of the model. The mean is less than 10 Find the best likelihood, subject
More informationA hidden semi-markov model for the occurrences of water pipe bursts
A hidden semi-markov model for the occurrences of water pipe bursts T. Economou 1, T.C. Bailey 1 and Z. Kapelan 1 1 School of Engineering, Computer Science and Mathematics, University of Exeter, Harrison
More informationVariable inspection plans for continuous populations with unknown short tail distributions
Variable inspection plans for continuous populations with unknown short tail distributions Wolfgang Kössler Abstract The ordinary variable inspection plans are sensitive to deviations from the normality
More informationConfidence Intervals for Reliability Growth Models with Small Sample Sizes. Reliability growth models, Order statistics, Confidence intervals
Confidence Intervals for Reliability Growth Models with Small Sample Sizes John Quigley, Lesley Walls University of Strathclyde, Glasgow, Scotland Key Words Reliability growth models, Order statistics,
More informationTwo optimization problems in a stochastic bandit model
Two optimization problems in a stochastic bandit model Emilie Kaufmann joint work with Olivier Cappé, Aurélien Garivier and Shivaram Kalyanakrishnan Journées MAS 204, Toulouse Outline From stochastic optimization
More informationReview. December 4 th, Review
December 4 th, 2017 Att. Final exam: Course evaluation Friday, 12/14/2018, 10:30am 12:30pm Gore Hall 115 Overview Week 2 Week 4 Week 7 Week 10 Week 12 Chapter 6: Statistics and Sampling Distributions Chapter
More informationBandit models: a tutorial
Gdt COS, December 3rd, 2015 Multi-Armed Bandit model: general setting K arms: for a {1,..., K}, (X a,t ) t N is a stochastic process. (unknown distributions) Bandit game: a each round t, an agent chooses
More informationExact Linear Likelihood Inference for Laplace
Exact Linear Likelihood Inference for Laplace Prof. N. Balakrishnan McMaster University, Hamilton, Canada bala@mcmaster.ca p. 1/52 Pierre-Simon Laplace 1749 1827 p. 2/52 Laplace s Biography Born: On March
More informationECE 275A Homework 7 Solutions
ECE 275A Homework 7 Solutions Solutions 1. For the same specification as in Homework Problem 6.11 we want to determine an estimator for θ using the Method of Moments (MOM). In general, the MOM estimator
More informationQualifying Exam in Probability and Statistics. https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf
Part 1: Sample Problems for the Elementary Section of Qualifying Exam in Probability and Statistics https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf Part 2: Sample Problems for the Advanced Section
More informationTail bound inequalities and empirical likelihood for the mean
Tail bound inequalities and empirical likelihood for the mean Sandra Vucane 1 1 University of Latvia, Riga 29 th of September, 2011 Sandra Vucane (LU) Tail bound inequalities and EL for the mean 29.09.2011
More informationSome Applications of Hawkes Processes for Order Book Modelling
Some Applications of Hawkes Processes for Order Book Modelling Ioane Muni Toke Ecole Centrale Paris Chair of Quantitative Finance First Workshop on Quantitative Finance and Economics International Christian
More informationSTATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL. A Thesis. Presented to the. Faculty of. San Diego State University
STATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL A Thesis Presented to the Faculty of San Diego State University In Partial Fulfillment of the Requirements for the Degree
More informationPROBABILITY DISTRIBUTIONS
Review of PROBABILITY DISTRIBUTIONS Hideaki Shimazaki, Ph.D. http://goo.gl/visng Poisson process 1 Probability distribution Probability that a (continuous) random variable X is in (x,x+dx). ( ) P x < X
More informationIntroduction to Machine Learning
Introduction to Machine Learning Linear Regression Varun Chandola Computer Science & Engineering State University of New York at Buffalo Buffalo, NY, USA chandola@buffalo.edu Chandola@UB CSE 474/574 1
More informationStatistical Inference Using Progressively Type-II Censored Data with Random Scheme
International Mathematical Forum, 3, 28, no. 35, 1713-1725 Statistical Inference Using Progressively Type-II Censored Data with Random Scheme Ammar M. Sarhan 1 and A. Abuammoh Department of Statistics
More informationBayesian Modeling of Accelerated Life Tests with Random Effects
Bayesian Modeling of Accelerated Life Tests with Random Effects Ramón V. León Avery J. Ashby Jayanth Thyagarajan Joint Statistical Meeting August, 00 Toronto, Canada Abstract We show how to use Bayesian
More informationthe Presence of k Outliers
RESEARCH ARTICLE OPEN ACCESS On the Estimation of the Presence of k Outliers for Weibull Distribution in Amal S. Hassan 1, Elsayed A. Elsherpieny 2, and Rania M. Shalaby 3 1 (Department of Mathematical
More informationIntroduction to Reliability Theory (part 2)
Introduction to Reliability Theory (part 2) Frank Coolen UTOPIAE Training School II, Durham University 3 July 2018 (UTOPIAE) Introduction to Reliability Theory 1 / 21 Outline Statistical issues Software
More informationAn Analysis of Record Statistics based on an Exponentiated Gumbel Model
Communications for Statistical Applications and Methods 2013, Vol. 20, No. 5, 405 416 DOI: http://dx.doi.org/10.5351/csam.2013.20.5.405 An Analysis of Record Statistics based on an Exponentiated Gumbel
More informationBayesian and Frequentist Methods in Bandit Models
Bayesian and Frequentist Methods in Bandit Models Emilie Kaufmann, Telecom ParisTech Bayes In Paris, ENSAE, October 24th, 2013 Emilie Kaufmann (Telecom ParisTech) Bayesian and Frequentist Bandits BIP,
More informationStatistical Analysis of Competing Risks With Missing Causes of Failure
Proceedings 59th ISI World Statistics Congress, 25-3 August 213, Hong Kong (Session STS9) p.1223 Statistical Analysis of Competing Risks With Missing Causes of Failure Isha Dewan 1,3 and Uttara V. Naik-Nimbalkar
More informationSurvival Analysis: Weeks 2-3. Lu Tian and Richard Olshen Stanford University
Survival Analysis: Weeks 2-3 Lu Tian and Richard Olshen Stanford University 2 Kaplan-Meier(KM) Estimator Nonparametric estimation of the survival function S(t) = pr(t > t) The nonparametric estimation
More informationWeibull Reliability Analysis
Weibull Reliability Analysis = http://www.rt.cs.boeing.com/mea/stat/reliability.html Fritz Scholz (425-865-3623, 7L-22) Boeing Phantom Works Mathematics &Computing Technology Weibull Reliability Analysis
More informationFor more information about how to cite these materials visit
Author(s): Kerby Shedden, Ph.D., 2010 License: Unless otherwise noted, this material is made available under the terms of the Creative Commons Attribution Share Alike 3.0 License: http://creativecommons.org/licenses/by-sa/3.0/
More informationDesign of Optimal Bayesian Reliability Test Plans for a Series System
Volume 109 No 9 2016, 125 133 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://wwwijpameu ijpameu Design of Optimal Bayesian Reliability Test Plans for a Series System P
More informationStatistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach
Statistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach Jae-Kwang Kim Department of Statistics, Iowa State University Outline 1 Introduction 2 Observed likelihood 3 Mean Score
More informationBayesian Methods. David S. Rosenberg. New York University. March 20, 2018
Bayesian Methods David S. Rosenberg New York University March 20, 2018 David S. Rosenberg (New York University) DS-GA 1003 / CSCI-GA 2567 March 20, 2018 1 / 38 Contents 1 Classical Statistics 2 Bayesian
More informationABC methods for phase-type distributions with applications in insurance risk problems
ABC methods for phase-type with applications problems Concepcion Ausin, Department of Statistics, Universidad Carlos III de Madrid Joint work with: Pedro Galeano, Universidad Carlos III de Madrid Simon
More informationBIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation
BIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation Yujin Chung November 29th, 2016 Fall 2016 Yujin Chung Lec13: MLE Fall 2016 1/24 Previous Parametric tests Mean comparisons (normality assumption)
More informationCS Lecture 18. Expectation Maximization
CS 6347 Lecture 18 Expectation Maximization Unobserved Variables Latent or hidden variables in the model are never observed We may or may not be interested in their values, but their existence is crucial
More informationLecture 5 Models and methods for recurrent event data
Lecture 5 Models and methods for recurrent event data Recurrent and multiple events are commonly encountered in longitudinal studies. In this chapter we consider ordered recurrent and multiple events.
More informationProbabilistic Time Series Classification
Probabilistic Time Series Classification Y. Cem Sübakan Boğaziçi University 25.06.2013 Y. Cem Sübakan (Boğaziçi University) M.Sc. Thesis Defense 25.06.2013 1 / 54 Problem Statement The goal is to assign
More informationPattern Recognition and Machine Learning. Bishop Chapter 9: Mixture Models and EM
Pattern Recognition and Machine Learning Chapter 9: Mixture Models and EM Thomas Mensink Jakob Verbeek October 11, 27 Le Menu 9.1 K-means clustering Getting the idea with a simple example 9.2 Mixtures
More informationHypothesis testing: theory and methods
Statistical Methods Warsaw School of Economics November 3, 2017 Statistical hypothesis is the name of any conjecture about unknown parameters of a population distribution. The hypothesis should be verifiable
More informationA Multi-level Trend-Renewal Process for Modeling Systems with Recurrence Data
Statistics Preprints Statistics 4-2015 A Multi-level Trend-Renewal Process for Modeling Systems with Recurrence Data Zhibing Xu Virginia Tech Yili Hong Virginia Tech William Q. Meeker Iowa State University,
More informationKernel density estimation with doubly truncated data
Kernel density estimation with doubly truncated data Jacobo de Uña-Álvarez jacobo@uvigo.es http://webs.uvigo.es/jacobo (joint with Carla Moreira) Department of Statistics and OR University of Vigo Spain
More informationModeling and Analysis of Recurrent Event Data
Edsel A. Peña (pena@stat.sc.edu) Department of Statistics University of South Carolina Columbia, SC 29208 New Jersey Institute of Technology Conference May 20, 2012 Historical Perspective: Random Censorship
More informationHigh Dimensional Empirical Likelihood for Generalized Estimating Equations with Dependent Data
High Dimensional Empirical Likelihood for Generalized Estimating Equations with Dependent Data Song Xi CHEN Guanghua School of Management and Center for Statistical Science, Peking University Department
More informationSoftware Reliability & Testing
Repairable systems Repairable system A reparable system is obtained by glueing individual non-repairable systems each around a single failure To describe this gluing process we need to review the concept
More informationAN INTEGRAL MEASURE OF AGING/REJUVENATION FOR REPAIRABLE AND NON REPAIRABLE SYSTEMS
R&RAA # 1 (Vol.1) 8, March AN INEGRAL MEASURE OF AGING/REJUVENAION FOR REPAIRABLE AND NON REPAIRABLE SYSEMS M.P. Kaminskiy and V.V. Krivtsov Abstract his paper introduces a simple index that helps to assess
More information