Modeling of Dependence Between Critical Failure and Preventive Maintenance: The Repair Alert Model
|
|
- Lauren Cunningham
- 5 years ago
- Views:
Transcription
1 Modeling of Dependence Between Critical Failure and Preventive Maintenance: The Repair Alert Model Bo Lindqvist Helge Langseth Bård Støve Department of Mathematical Sciences Norwegian University of Science and Technology Workshop Delft June 2003 p.1/20
2 Purpose of talk Explicit modeling of maintenance in repairable systems: Several failure mechanisms Imperfect repair Degraded failures Preventive Maintenance (PM) (scheduled/unscheduled) Several failure mechanisms Periodically tested components Workshop Delft June 2003 p.2/20
3 Simultaneous Modeling of Time to Failure and PM 0 X (failure time) Competing risk problem: Z (PM) X is the (potential) failure time of an item Z is the time of a (potential) PM action before time X Possible outcomes: Failure at X, no PM, Z is not observed PM at Z (< X), while X is not observed Typical observations: N independent pairs {min(x, Z), I(Z < X)} are observed, represented as x 1,..., x m ; z 1,..., z n where N = m + n Workshop Delft June 2003 p.3/20
4 Cooke s Random Signs Censoring 0 X (failure time) Definition: Motivation: Z (PM) The event {Z < X} is independent of X Suppose the item emits some warning of deterioration, prior to failure. If warning signal is observed, then the item will be preventively maintained and the PM variable Z is observed. If the event of observing the signal is independent of the item s age, then random signs censoring is appropriate. Workshop Delft June 2003 p.4/20
5 Definition: Motivation: Repair Alert Model The pair (X, Z) satisfies the requirements of the Repair Alert Model provided the following two conditions both hold. 1. Random signs censoring, that is {Z < X} is independent of X 2. There exists an increasing function G(x) such that for all x > 0, P (Z z Z < X, X = x) = G(z) 0 z x G(x), The function G(z) is called the cumulative repair alert function. Its derivative g(z) is called the repair alert function. G(x) G(z) Z = z X = x Workshop Delft June 2003 p.5/20
6 Some Notation F X (x) = P (X x) (Distribution function) F Z (z) = P (Z z) FX (x) = P (X x, X < Z) (Subdistribution function) FZ (z) = P (Z z, Z < X) F X (x) = P (X x X < Z) (Conditional subdistribution function) F Z (z) = P (Z z Z < X) Only the subdistribution functions are in general identifiable from competing risk data Under Random Signs Censoring we have by definition F X (x) = P (X x X < Z) = P (X x) = F X (x) making F X (x) identifiable. Workshop Delft June 2003 p.6/20
7 Identification of G( ) in Repair Alert Model Recall assumptions: {Z < X} is independent of X P (Z z Z < X, X = x) = G(z) G(x), 0 z x From this: F Z (z) = P (Z z Z < X) = = 0 0 P (Z z Z < X, X = x)p (x X x + dx Z < X) min( G(z) G(x), 1) f X(x)dx = F X (z) + G(z) z f X (x) G(x) dx Differentiate: f Z (z) = f X (z) + g(z) z f X (x) G(x) dx G(z) f X(z) G(z) = g(z) z f X (x) G(x) dx. Combine: F Z (z) = F X (z) + G(z) f Z (z) g(z) Workshop Delft June 2003 p.7/20
8 Identification of G( ) in Repair Alert Model Basic formula: F Z (z) = F X (z) + G(z) f Z (z) g(z) Rearranging: g(z) G(z) = fz (z) F Z (z) F X (z). Integrating from fixed point a > 0, assuming G(a) > 0: w a g(z) G(z) dz = w a f Z (z) F Z (z) F X (z) dz G(w) G(a) dy y = FZ (w) F Z (a) dy 1 y F X ( F Z (y)). G(w) FZ (w) G(a) = exp{ F Z (a) dy }. 1 y F X ( F Z (y)) Hence, G(z) is identifiable from data, modulo a constant. Workshop Delft June 2003 p.8/20
9 Nonparametric Estimation of Cumulative Repair Function Let x 1,..., x m and z 1,..., z n be the observed X and Z. F X (t) is estimated by (1) ˆF X (t) = i m for x i t < x i+1, i = 0, 1,..., m. With t = 1 F Z (y) we get (2) ˆF X ( F 1 Z (y)) = i m for F Z (x i ) y < F Z (x i+1 ), i = 0, 1,..., m. Thus (3) FZ (x j ) F Z (x 1 ) y ˆF X ( dy F 1 Z (y)) = j 1 i=1 FZ (x i+1 ) F Z (x i ) dy y i/m = j 1 i=1 ln F Z (x i+1 ) i/m F Z (x i ) i/m. yielding the estimator (4) Ĝ(x j ) Ĝ(x 1 ) = j 1 i=1 ˆ F Z (x i+1 ) i/m. ˆ F Z (x i ) i/m The estimator for F Z (t) is Workshop Delft June 2003 p.9/20
10 Simulated example: Repair alert model Failure time X is exponentially distributed, failure rate λ = 1 q = P (Z < X) = 0.5 Cumulative repair function is G(x) = x, i.e. Z given Z < X, X = x is uniform on (0, x) m + n = Cumulative Repair Alert Function 6 Log Log Cumulative Repair Alert Function m+n= m+n= q= q=0.5 log(g) X exponentially distributed X independent of Z<X Z Z<X,X=x uniform on (0,x) 2 X exponentially distributed X independent of Z<X Z Z<X,X=x uniform on (0,x) Time log(time) Workshop Delft June 2003 p.10/20
11 Parametric statistical inference for repair alert model Suppose we observe N independent copies of the pair This gives data x 1,..., x m ; z 1,..., z n, with obvious meaning, where n + m = N. {min(x, Z), I(Z < X)} We are interested in estimating Density of X, f X (x) (for example exponential, Weibull, etc.) Repair probability q Repair function g(x) LIKELIHOOD FUNCTION Likelihood contribution from an observation is under the repair alert model: f(x i, X < Z) = (1 q)f X (x i ) for x i f(z i, Z < X) = q f Z (z i ) for z i = qg(z i ) z i (f X (x)/g(x))dx Workshop Delft June 2003 p.11/20
12 An Exponential-Power Repair Alert Model f X (x) = λe λx, q = P (Z < X), G(x) = x β Likelihood contributions: Complete log-likelihood: f(x i, X < Z) = (1 q)f X (x i ) = (1 q)λe λx i f(z i, Z < X) = qg(z i ) = qβz β 1 i z i z i = qλβ(λz i ) β 1 (f X (x)/g(x))dx λe λx x β dx λz i w β e w dw l(λ, β, q) = m ln(1 q) + n ln q + (n + m) ln λ + n ln β λ n (β 1) ln(λz i ) + i=1 n i=1 m x i + i=1 ln( w β e w dw). λz i Workshop Delft June 2003 p.12/20
13 Simulated parametric example N = m + n = 100 Parameter True value Estimate Lower bound Upper bound λ β q Maximum likelihood estimates and approximate 95% confidence intervals for simulated data The 95% confidence region for λ (horizontal axis) and β from the simulated data Workshop Delft June 2003 p.13/20
14 The EM-algorithm general description Y is the observed data; X is a piece of unknown data; θ is the parameter of interest; and l C (θ; Y, X ) is the hypothetical complete-data log-likelihood, defined for all θ Ω. Starting with an initial parameter value θ (0) Ω, the EM algorithm repeats the following two steps until convergence. E-step: Compute l (j) (θ) = E X Y,θ (j 1) [l C (θ; Y, X )], where the expectation is taken with respect to the conditional distribution of the missing data X given the observed data Y, and the current numerical value θ (j 1) is used in evaluating the expected value. M-step: Find θ (j) Ω that maximizes l (j) (θ). Workshop Delft June 2003 p.14/20
15 The EM-algorithm in our model Likelihood contributions to complete likelihood are now simplified: Resulting iterative algorithm: f(x i, X < Z) = (1 q)f X (x i ) = (1 q)λe λx i f(x i, z i, Z < X) = qλe λx i βzβ 1 i x β i ˆq = n N ˆλ j+1 = N m i=1 x i + (1/ˆλ j ) n i=1 ˆλj zi ˆλj zi w ( ˆβ j 1) e w dw w ˆβ j e w dw ˆβ j+1 = n ˆλj z ln(w)w ˆβ j e w dw i i=1 w ˆβ j e w dw ˆλj zi n ln(ˆλ j z i ) Workshop Delft June 2003 p.15/20
16 Example: VHF data (Mendenhall and Hader, 1958) Times to failure for ARC-1 VHF communication transmitter-receivers of a single commercial airline. X = time to confirmed failure, m = 218 Z = time to unconfirmed failure (censoring), n = X Z Workshop Delft June 2003 p.16/20
17 Example: VHF data (Mendenhall and Hader, 1958) Left: Estimated subsurvival functions for X, thin line, and Z, thick line Right: Estimated conditional subsurvival functions for X, thin line, Z, thick line, estimated Φ(t) = P (Z < X X > t, Z > t), and the estimated CSF corresponding to an independent exponential model Null hypothesis of independent exponentials is accepted at p-value.15 Workshop Delft June 2003 p.17/20
18 Estimated failure rate for X, VHF data Assuming Failure rate for X Repair Alert, X exponentially distributed X, Z independent exponentials Bounds for the failure rate of X. Thick line: accounting for sampling fluctuations. Thin line: without sampling fluctuations Workshop Delft June 2003 p.18/20
19 VHF data: Nonparametric Repair Alert Model The estimated ln G(x j ) G(x 1 ) plotted against ln(x j) Workshop Delft June 2003 p.19/20
20 VHF data: Parametric Repair Alert Model MODEL Failure time X is exponentially distributed, failure rate λ Cumulative repair function is G(x) = x β Parameter Estimate Lower bound Upper bound λ β q estimates and approximate 95% confidence Maximum likelihood intervals for parametric repair alert model Workshop Delft June 2003 p.20/20
CHAPTER 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 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 informationCompeting risks for repairable systems: A data study
Competing risks for repairable systems: A data study Helge Langseth and Bo Henry Lindqvist Department of Mathematical Sciences, Norwegian University of Science and Technology, Norway Abstract We analyze
More informationLikelihood Construction, Inference for Parametric Survival Distributions
Week 1 Likelihood Construction, Inference for Parametric Survival Distributions In this section we obtain the likelihood function for noninformatively rightcensored survival data and indicate how to make
More informationMultistate Modeling and Applications
Multistate Modeling and Applications Yang Yang Department of Statistics University of Michigan, Ann Arbor IBM Research Graduate Student Workshop: Statistics for a Smarter Planet Yang Yang (UM, Ann Arbor)
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 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 informationExponential Family and Maximum Likelihood, Gaussian Mixture Models and the EM Algorithm. by Korbinian Schwinger
Exponential Family and Maximum Likelihood, Gaussian Mixture Models and the EM Algorithm by Korbinian Schwinger Overview Exponential Family Maximum Likelihood The EM Algorithm Gaussian Mixture Models Exponential
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 informationCensoring and Truncation - Highlighting the Differences
Censoring and Truncation - Highlighting the Differences Micha Mandel The Hebrew University of Jerusalem, Jerusalem, Israel, 91905 July 9, 2007 Micha Mandel is a Lecturer, Department of Statistics, The
More informationFailure modeling and maintenance optimization for a railway line
Failure modeling and maintenance optimization for a railway line Per Hokstad *, Helge Langseth SINTEF Technology and Society, Department of Safety and Reliability, N-7465 Trondheim, Norway Bo H. Lindqvist,
More informationStatistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation
Statistics - Lecture One Charlotte Wickham wickham@stat.berkeley.edu http://www.stat.berkeley.edu/~wickham/ Outline 1. Basic ideas about estimation 2. Method of Moments 3. Maximum Likelihood 4. Confidence
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 20, 2012 Introduction Introduction The world of the model-builder
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 informationEXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY
EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 00 MODULE : Statistical Inference Time Allowed: Three Hours Candidates should answer FIVE questions. All questions carry equal marks. The
More information1 Degree distributions and data
1 Degree distributions and data A great deal of effort is often spent trying to identify what functional form best describes the degree distribution of a network, particularly the upper tail of that distribution.
More informationf X (x) = λe λx, , x 0, k 0, λ > 0 Γ (k) f X (u)f X (z u)du
11 COLLECTED PROBLEMS Do the following problems for coursework 1. Problems 11.4 and 11.5 constitute one exercise leading you through the basic ruin arguments. 2. Problems 11.1 through to 11.13 but excluding
More informationAnalysis of Progressive Type-II Censoring. in the Weibull Model for Competing Risks Data. with Binomial Removals
Applied Mathematical Sciences, Vol. 5, 2011, no. 22, 1073-1087 Analysis of Progressive Type-II Censoring in the Weibull Model for Competing Risks Data with Binomial Removals Reza Hashemi and Leila Amiri
More informationAnalysis of Type-II Progressively Hybrid Censored Data
Analysis of Type-II Progressively Hybrid Censored Data Debasis Kundu & Avijit Joarder Abstract The mixture of Type-I and Type-II censoring schemes, called the hybrid censoring scheme is quite common in
More informationMathematical statistics
October 1 st, 2018 Lecture 11: Sufficient statistic Where are we? Week 1 Week 2 Week 4 Week 7 Week 10 Week 14 Probability reviews Chapter 6: Statistics and Sampling Distributions Chapter 7: Point Estimation
More informationStatistics for scientists and engineers
Statistics for scientists and engineers February 0, 006 Contents Introduction. Motivation - why study statistics?................................... Examples..................................................3
More informationMTH739U/P: Topics in Scientific Computing Autumn 2016 Week 6
MTH739U/P: Topics in Scientific Computing Autumn 16 Week 6 4.5 Generic algorithms for non-uniform variates We have seen that sampling from a uniform distribution in [, 1] is a relatively straightforward
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 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 information11 Survival Analysis and Empirical Likelihood
11 Survival Analysis and Empirical Likelihood The first paper of empirical likelihood is actually about confidence intervals with the Kaplan-Meier estimator (Thomas and Grunkmeier 1979), i.e. deals with
More informationParameter Estimation of Power Lomax Distribution Based on Type-II Progressively Hybrid Censoring Scheme
Applied Mathematical Sciences, Vol. 12, 2018, no. 18, 879-891 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.8691 Parameter Estimation of Power Lomax Distribution Based on Type-II Progressively
More informationDepartment of Statistics, School of Mathematical Sciences, Ferdowsi University of Mashhad, Iran.
JIRSS (2012) Vol. 11, No. 2, pp 191-202 A Goodness of Fit Test For Exponentiality Based on Lin-Wong Information M. Abbasnejad, N. R. Arghami, M. Tavakoli Department of Statistics, School of Mathematical
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 informationTime-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 informationSurvival Analysis I (CHL5209H)
Survival Analysis Dalla Lana School of Public Health University of Toronto olli.saarela@utoronto.ca January 7, 2015 31-1 Literature Clayton D & Hills M (1993): Statistical Models in Epidemiology. Not really
More informationOn the Comparison of Fisher Information of the Weibull and GE Distributions
On the Comparison of Fisher Information of the Weibull and GE Distributions Rameshwar D. Gupta Debasis Kundu Abstract In this paper we consider the Fisher information matrices of the generalized exponential
More informationSTAT 6350 Analysis of Lifetime Data. Probability Plotting
STAT 6350 Analysis of Lifetime Data Probability Plotting Purpose of Probability Plots Probability plots are an important tool for analyzing data and have been particular popular in the analysis of life
More informationTMA 4275 Lifetime Analysis June 2004 Solution
TMA 4275 Lifetime Analysis June 2004 Solution Problem 1 a) Observation of the outcome is censored, if the time of the outcome is not known exactly and only the last time when it was observed being intact,
More informationA Bivariate Weibull Regression Model
c Heldermann Verlag Economic Quality Control ISSN 0940-5151 Vol 20 (2005), No. 1, 1 A Bivariate Weibull Regression Model David D. Hanagal Abstract: In this paper, we propose a new bivariate Weibull regression
More informationDUBLIN CITY UNIVERSITY
DUBLIN CITY UNIVERSITY SAMPLE EXAMINATIONS 2017/2018 MODULE: QUALIFICATIONS: Simulation for Finance MS455 B.Sc. Actuarial Mathematics ACM B.Sc. Financial Mathematics FIM YEAR OF STUDY: 4 EXAMINERS: Mr
More informatione x3 dx dy. 0 y x 2, 0 x 1.
Problem 1. Evaluate by changing the order of integration y e x3 dx dy. Solution:We change the order of integration over the region y x 1. We find and x e x3 dy dx = y x, x 1. x e x3 dx = 1 x=1 3 ex3 x=
More informationHybrid Censoring Scheme: An Introduction
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline 1 2 3 4 5 Outline 1 2 3 4 5 What is? Lifetime data analysis is used to analyze data in which the time
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 informationChing-Han Hsu, BMES, National Tsing Hua University c 2015 by Ching-Han Hsu, Ph.D., BMIR Lab. = a + b 2. b a. x a b a = 12
Lecture 5 Continuous Random Variables BMIR Lecture Series in Probability and Statistics Ching-Han Hsu, BMES, National Tsing Hua University c 215 by Ching-Han Hsu, Ph.D., BMIR Lab 5.1 1 Uniform Distribution
More informationHybrid Censoring; An Introduction 2
Hybrid Censoring; An Introduction 2 Debasis Kundu Department of Mathematics & Statistics Indian Institute of Technology Kanpur 23-rd November, 2010 2 This is a joint work with N. Balakrishnan Debasis Kundu
More informationPump failure data. Pump Failures Time
Outline 1. Poisson distribution 2. Tests of hypothesis for a single Poisson mean 3. Comparing multiple Poisson means 4. Likelihood equivalence with exponential model Pump failure data Pump 1 2 3 4 5 Failures
More informationSpring 2012 Math 541B Exam 1
Spring 2012 Math 541B Exam 1 1. A sample of size n is drawn without replacement from an urn containing N balls, m of which are red and N m are black; the balls are otherwise indistinguishable. Let X denote
More informationNotes largely based on Statistical Methods for Reliability Data by W.Q. Meeker and L. A. Escobar, Wiley, 1998 and on their class notes.
Unit 2: Models, Censoring, and Likelihood for Failure-Time Data Notes largely based on Statistical Methods for Reliability Data by W.Q. Meeker and L. A. Escobar, Wiley, 1998 and on their class notes. Ramón
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 informationIntegration by Substitution
November 22, 2013 Introduction 7x 2 cos(3x 3 )dx =? 2xe x2 +5 dx =? Chain rule The chain rule: d dx (f (g(x))) = f (g(x)) g (x). Use the chain rule to find f (x) and then write the corresponding anti-differentiation
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER.
Two hours MATH38181 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER EXTREME VALUES AND FINANCIAL RISK Examiner: Answer any FOUR
More informationMaximum Smoothed Likelihood for Multivariate Nonparametric Mixtures
Maximum Smoothed Likelihood for Multivariate Nonparametric Mixtures David Hunter Pennsylvania State University, USA Joint work with: Tom Hettmansperger, Hoben Thomas, Didier Chauveau, Pierre Vandekerkhove,
More informationChapter 4: CONTINUOUS RANDOM VARIABLES AND PROBABILITY DISTRIBUTIONS
Chapter 4: CONTINUOUS RANDOM VARIABLES AND PROBABILITY DISTRIBUTIONS Part 4: Gamma Distribution Weibull Distribution Lognormal Distribution Sections 4-9 through 4-11 Another exponential distribution example
More informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS330 / MAS83 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-0 8 Parametric models 8. Introduction In the last few sections (the KM
More informationMultivariable Calculus and Matrix Algebra-Summer 2017
Multivariable Calculus and Matrix Algebra-Summer 017 Homework 4 Solutions Note that the solutions below are for the latest version of the problems posted. For those of you who worked on an earlier version
More informationReliability Engineering I
Happiness is taking the reliability final exam. Reliability Engineering I ENM/MSC 565 Review for the Final Exam Vital Statistics What R&M concepts covered in the course When Monday April 29 from 4:30 6:00
More informationSPRING 2007 EXAM C SOLUTIONS
SPRING 007 EXAM C SOLUTIONS Question #1 The data are already shifted (have had the policy limit and the deductible of 50 applied). The two 350 payments are censored. Thus the likelihood function is L =
More informationRandom variables. DS GA 1002 Probability and Statistics for Data Science.
Random variables DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Motivation Random variables model numerical quantities
More informationMath Review Sheet, Fall 2008
1 Descriptive Statistics Math 3070-5 Review Sheet, Fall 2008 First we need to know about the relationship among Population Samples Objects The distribution of the population can be given in one of the
More information1 Inverse Transform Method and some alternative algorithms
Copyright c 2016 by Karl Sigman 1 Inverse Transform Method and some alternative algorithms Assuming our computer can hand us, upon demand, iid copies of rvs that are uniformly distributed on (0, 1), it
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 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 informationStatistical inference for Markov deterioration models of bridge conditions in the Netherlands
Statistical inference for Markov deterioration models of bridge conditions in the Netherlands M.J.Kallen & J.M. van Noortwijk HKV Consultants, Lelystad, and Delft University of Technology, Delft, Netherlands
More informationMATH Midterm 1 Sample 4
1. (15 marks) (a) (4 marks) Given the function: f(x, y) = arcsin x 2 + y, find its first order partial derivatives at the point (0, 3). Simplify your answers. Solution: Compute the first order partial
More informationAccelerated life testing in the presence of dependent competing causes of failure
isid/ms/25/5 April 21, 25 http://www.isid.ac.in/ statmath/eprints Accelerated life testing in the presence of dependent competing causes of failure Isha Dewan S. B. Kulathinal Indian Statistical Institute,
More informationProbability Theory and Statistics. Peter Jochumzen
Probability Theory and Statistics Peter Jochumzen April 18, 2016 Contents 1 Probability Theory And Statistics 3 1.1 Experiment, Outcome and Event................................ 3 1.2 Probability............................................
More informationBMIR Lecture Series on Probability and Statistics Fall, 2015 Uniform Distribution
Lecture #5 BMIR Lecture Series on Probability and Statistics Fall, 2015 Department of Biomedical Engineering and Environmental Sciences National Tsing Hua University s 5.1 Definition ( ) A continuous random
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 informationSurvival Analysis. Lu Tian and Richard Olshen Stanford University
1 Survival Analysis Lu Tian and Richard Olshen Stanford University 2 Survival Time/ Failure Time/Event Time We will introduce various statistical methods for analyzing survival outcomes What is the survival
More informationA Quantile Implementation of the EM Algorithm and Applications to Parameter Estimation with Interval Data
A Quantile Implementation of the EM Algorithm and Applications to Parameter Estimation with Interval Data Chanseok Park The author is with the Department of Mathematical Sciences, Clemson University, Clemson,
More informationST5212: Survival Analysis
ST51: Survival Analysis 8/9: Semester II Tutorial 1. A model for lifetimes, with a bathtub-shaped hazard rate, is the exponential power distribution with survival fumction S(x) =exp{1 exp[(λx) α ]}. (a)
More informationThis paper is not to be removed from the Examination Halls
~~ST104B ZA d0 This paper is not to be removed from the Examination Halls UNIVERSITY OF LONDON ST104B ZB BSc degrees and Diplomas for Graduates in Economics, Management, Finance and the Social Sciences,
More informationSurvival Regression Models
Survival Regression Models David M. Rocke May 18, 2017 David M. Rocke Survival Regression Models May 18, 2017 1 / 32 Background on the Proportional Hazards Model The exponential distribution has constant
More informationSurvival Analysis. Stat 526. April 13, 2018
Survival Analysis Stat 526 April 13, 2018 1 Functions of Survival Time Let T be the survival time for a subject Then P [T < 0] = 0 and T is a continuous random variable The Survival function is defined
More informationEstimation Under Multivariate Inverse Weibull Distribution
Global Journal of Pure and Applied Mathematics. ISSN 097-768 Volume, Number 8 (07), pp. 4-4 Research India Publications http://www.ripublication.com Estimation Under Multivariate Inverse Weibull Distribution
More informationDepartment of Statistical Science FIRST YEAR EXAM - SPRING 2017
Department of Statistical Science Duke University FIRST YEAR EXAM - SPRING 017 Monday May 8th 017, 9:00 AM 1:00 PM NOTES: PLEASE READ CAREFULLY BEFORE BEGINNING EXAM! 1. Do not write solutions on the exam;
More informationDr. Maddah ENMG 617 EM Statistics 10/15/12. Nonparametric Statistics (2) (Goodness of fit tests)
Dr. Maddah ENMG 617 EM Statistics 10/15/12 Nonparametric Statistics (2) (Goodness of fit tests) Introduction Probability models used in decision making (Operations Research) and other fields require fitting
More informationA generalized Brown-Proschan model for preventive and corrective maintenance. Laurent DOYEN.
A generalized Brown-Proschan model for preventive and corrective maintenance laurent.doyen@iut2.upmf-grenoble.fr Jean Kuntzmann Laboratory Grenoble University France 1 of 28 I. Introduction The dependability
More informationSemiparametric Regression
Semiparametric Regression Patrick Breheny October 22 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Introduction Over the past few weeks, we ve introduced a variety of regression models under
More informationMaster s Written Examination
Master s Written Examination Option: Statistics and Probability Spring 05 Full points may be obtained for correct answers to eight questions Each numbered question (which may have several parts) is worth
More informationA Recursive Formula for the Kaplan-Meier Estimator with Mean Constraints
Noname manuscript No. (will be inserted by the editor) A Recursive Formula for the Kaplan-Meier Estimator with Mean Constraints Mai Zhou Yifan Yang Received: date / Accepted: date Abstract In this note
More informationSOLUTION FOR HOMEWORK 8, STAT 4372
SOLUTION FOR HOMEWORK 8, STAT 4372 Welcome to your 8th homework. Here you have an opportunity to solve classical estimation problems which are the must to solve on the exam due to their simplicity. 1.
More informationStatistics 135 Fall 2008 Final Exam
Name: SID: Statistics 135 Fall 2008 Final Exam Show your work. The number of points each question is worth is shown at the beginning of the question. There are 10 problems. 1. [2] The normal equations
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More informationReliability Analysis of Repairable Systems with Dependent Component Failures. under Partially Perfect Repair
Reliability Analysis of Repairable Systems with Dependent Component Failures under Partially Perfect Repair Qingyu Yang & Nailong Zhang Department of Industrial and System Engineering Wayne State University
More informationApplied Statistics and Probability for Engineers. Sixth Edition. Chapter 4 Continuous Random Variables and Probability Distributions.
Applied Statistics and Probability for Engineers Sixth Edition Douglas C. Montgomery George C. Runger Chapter 4 Continuous Random Variables and Probability Distributions 4 Continuous CHAPTER OUTLINE Random
More informationChapter 4 Continuous Random Variables and Probability Distributions
Applied Statistics and Probability for Engineers Sixth Edition Douglas C. Montgomery George C. Runger Chapter 4 Continuous Random Variables and Probability Distributions 4 Continuous CHAPTER OUTLINE 4-1
More informationFormal groups. Peter Bruin 2 March 2006
Formal groups Peter Bruin 2 March 2006 0. Introduction The topic of formal groups becomes important when we want to deal with reduction of elliptic curves. Let R be a discrete valuation ring with field
More informationSurvival Analysis Math 434 Fall 2011
Survival Analysis Math 434 Fall 2011 Part IV: Chap. 8,9.2,9.3,11: Semiparametric Proportional Hazards Regression Jimin Ding Math Dept. www.math.wustl.edu/ jmding/math434/fall09/index.html Basic Model Setup
More informationConstant Stress Partially Accelerated Life Test Design for Inverted Weibull Distribution with Type-I Censoring
Algorithms Research 013, (): 43-49 DOI: 10.593/j.algorithms.01300.0 Constant Stress Partially Accelerated Life Test Design for Mustafa Kamal *, Shazia Zarrin, Arif-Ul-Islam Department of Statistics & Operations
More informationA Sequential Bayesian Approach with Applications to Circadian Rhythm Microarray Gene Expression Data
A Sequential Bayesian Approach with Applications to Circadian Rhythm Microarray Gene Expression Data Faming Liang, Chuanhai Liu, and Naisyin Wang Texas A&M University Multiple Hypothesis Testing Introduction
More informationSOLUTIONS TO MATH68181 EXTREME VALUES AND FINANCIAL RISK EXAM
SOLUTIONS TO MATH68181 EXTREME VALUES AND FINANCIAL RISK EXAM Solutions to Question A1 a) The marginal cdfs of F X,Y (x, y) = [1 + exp( x) + exp( y) + (1 α) exp( x y)] 1 are F X (x) = F X,Y (x, ) = [1
More informationContinuous Distributions
Chapter 5 Continuous Distributions 5.1 Density and Distribution Functions In many situations random variables can take any value on the real line or in a certain subset of the real line. For concrete examples,
More informationSTAT331. Cox s Proportional Hazards Model
STAT331 Cox s Proportional Hazards Model In this unit we introduce Cox s proportional hazards (Cox s PH) model, give a heuristic development of the partial likelihood function, and discuss adaptations
More informationMathematics Qualifying Examination January 2015 STAT Mathematical Statistics
Mathematics Qualifying Examination January 2015 STAT 52800 - Mathematical Statistics NOTE: Answer all questions completely and justify your derivations and steps. A calculator and statistical tables (normal,
More informationSTAB57: Quiz-1 Tutorial 1 (Show your work clearly) 1. random variable X has a continuous distribution for which the p.d.f.
STAB57: Quiz-1 Tutorial 1 1. random variable X has a continuous distribution for which the p.d.f. is as follows: { kx 2.5 0 < x < 1 f(x) = 0 otherwise where k > 0 is a constant. (a) (4 points) Determine
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 informationEXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS
EXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS Course code: 8D Date: Thursday April 8 th, Time: 4h 7h Place: AUD 3 Read this first! Write your name and student identification number on each paper
More informationFULL LIKELIHOOD INFERENCES IN THE COX MODEL
October 20, 2007 FULL LIKELIHOOD INFERENCES IN THE COX MODEL BY JIAN-JIAN REN 1 AND MAI ZHOU 2 University of Central Florida and University of Kentucky Abstract We use the empirical likelihood approach
More informationExercises and Answers to Chapter 1
Exercises and Answers to Chapter The continuous type of random variable X has the following density function: a x, if < x < a, f (x), otherwise. Answer the following questions. () Find a. () Obtain mean
More informationPairwise rank based likelihood for estimating the relationship between two homogeneous populations and their mixture proportion
Pairwise rank based likelihood for estimating the relationship between two homogeneous populations and their mixture proportion Glenn Heller and Jing Qin Department of Epidemiology and Biostatistics Memorial
More informationStochastic Calculus. Kevin Sinclair. August 2, 2016
Stochastic Calculus Kevin Sinclair August, 16 1 Background Suppose we have a Brownian motion W. This is a process, and the value of W at a particular time T (which we write W T ) is a normally distributed
More informationYou may use your calculator and a single page of notes.
LAST NAME (Please Print): KEY FIRST NAME (Please Print): HONOR PLEDGE (Please Sign): Statistics 111 Midterm 2 This is a closed book exam. You may use your calculator and a single page of notes. The room
More informationEstimation of Parameters of the Weibull Distribution Based on Progressively Censored Data
International Mathematical Forum, 2, 2007, no. 41, 2031-2043 Estimation of Parameters of the Weibull Distribution Based on Progressively Censored Data K. S. Sultan 1 Department of Statistics Operations
More informationEconomics 520. Lecture Note 19: Hypothesis Testing via the Neyman-Pearson Lemma CB 8.1,
Economics 520 Lecture Note 9: Hypothesis Testing via the Neyman-Pearson Lemma CB 8., 8.3.-8.3.3 Uniformly Most Powerful Tests and the Neyman-Pearson Lemma Let s return to the hypothesis testing problem
More information