Survival Analysis: Weeks 2-3. Lu Tian and Richard Olshen Stanford University
|
|
- Beatrix Fisher
- 5 years ago
- Views:
Transcription
1 Survival Analysis: Weeks 2-3 Lu Tian and Richard Olshen Stanford University
2 2 Kaplan-Meier(KM) Estimator Nonparametric estimation of the survival function S(t) = pr(t > t) The nonparametric estimation is more robust and does not depend on any parametric assumption.
3 3 If there is no censoring S(t) can be consistently estimated by Ŝ(t) = n n i= I(T i > t). Ŝ is a discrete distribution with mass probability of n at observed times T,, T n Ŝ(t) is the nonparametric maximum likelihood estimator (NPMLE) for S(t).
4 4 CDF as NPMLE Assuming that F ( ) is discrete with mass probability at T < T 2 < < T n, where {T, T 2, } are observed times. Let f = pr(t = T ), f 2 = pr(t = T 2 ),. Objective: estimate f, f 2,. Method: maximize n i= f i subject to n i= f i =. Solution: ˆf = ˆf 2 = = ˆf n = n.
5 5 Data and Assumptions Data: {(U i, δ i ), i =,, n} where U i = min(t i, C i ) and δ i = I(T i C i ). Assumptions:. T,, T n i.i.d. F ( ) = S( ) 2. C,, C n i.i.d. G( ) 3. T i C i, i =,, n. Noninformative censoring!
6 6 If there is censoring Assuming that F ( ) is discrete with mass probability at v < v 2 <, where {v, v 2, } are observed times. Let f = pr(t = v ), f 2 = pr(t = v 2 ),. Objective: estimate f, f 2,.
7 7 Example Obs: 2, 2, 3 +, 5, 5 +, 7, 9, 6, 6, 8 +, where + means censored v = 2; v 2 = 3, v 3 = 5, v 4 = 7, v 5 = 9, v 6 = 6, v 7 = 8, v 8 = 8 + The likelihood function in terms of (f, f 2, ) : L(F ) = f 2 (f 3 +f 4 +f 5 +f 6 +f 7 +f 8 )f 3 (f 4 +f 5 +f 6 +f 7 +f 8 )f 4 f 5 f6 2 f 8, where f + f 2 + f 3 + f 4 + f 5 + f 6 + f 7 + f 8 =
8 8 Reparametrization tricks The discrete hazard function: h = pr(t = v ) and h j = pr(t = v j T > v j ), j > 2 For t [v j, v j+ ) S(t) = pr(t > t) = pr(t > v j ) = j ( h i ). i= For t = v j j f j = f(t) = pr(t = t) = h j ( h i ). i=
9 9 Reparametrization tricks The likelihood function in terms of (h, h 2, ) : L(F ) =h 2 {( h )( h 2 )} {( h )( h 2 )h 3 } {( h )( h 2 )( h 3 )} {( h )( h 2 )( h 3 )h 4 } {( h )( h 2 )( h 3 )( h 4 )h 5 } {( h )( h 2 )( h 3 )( h 4 )( h 5 )h 6 } 2 {( h )( h 2 )( h 3 )( h 4 )( h 5 )( h 6 )( h 7 )} =h 2 ( h ) 8 ( h 2 ) 8 h 3 ( h 3 ) 6 h 4 ( h 4 ) 4 h 5 ( h 5 ) 3 h 2 6( h 6 ) ( h 7 )
10 0 KM estimation The likelihood function L(F ) = j h d j j ( h j) Y (v j) d j where n d j = δ i I(U i = v j ) = #failures at v j Y (v j ) = i= n I(U i v j ) = # at risk at v j. i=
11 KM estimation ĥj = d j /Y (v j ) Ŝ(t) = t < v j i= ( ĥi), v j t < v j+ which is the Kaplan-Meier Estimator.
12 2 Example v j Y (v j ) d j ĥ j Ŝ(v j ) = j ( ĥi) = ˆP (T > v j ) i= / /7.69 (= ) 7 5 /5.55 (= ) 9 4 /4.4 (= ) /3.4 (=.4 3 )
13 3 KM estimation Suppose that v g denotes the largest v j for which Y (v j ) > 0.. if d g = Y (v j ), then Ŝ(t) = 0 for t v g 2. if d g < Y (v j ), then Ŝ(t) > 0 but not defined for t > v g. (Not identifiable beyond v g.) The survival distribution may not be estimable with right-censored data. Implicit extrapolation is sometimes used. The KM estimator can also be used to estimate the survival function for the censoring distribution.
14 4 KM estimation KM estimator is a special MLE Ŝ(t) = argmax F L(F ) where F is the CDF for all discrete random variables (nonparametric MLE).
15 5 Self-Consistency No censoring Ŝ(t) = n n i= I(T i > t) Right censoring: Ŝ(t) = n n i= E(I(T i > t) U i, δ i ). E(I(T i > t) U i, δ i = ) = I(U i > t) 2. E(I(T i > t) U i, δ i = 0) = S(t)/S(U i )I(t U i ) + I(U i > t) Self-consistency iteration: { Ŝ new (t) = n n i= The solution is still the KM estimator. I(U i > t) + ( δ i ) Ŝold(t) Ŝ old (U i ) I(U i t) }
16 6 Redistribution of Mass Step Step Step ( 8 70 ) 5 ( 8 70 ) 5 ( 8 70 ) 5 ( 8 70 ) 5 ( 8 70 ) Total Mass }{{}}{{}}{{} Assume this is somewhere > 8
17 7 Nelson-Aalen Estimator How to estimate the cumulative hazard function? Ĥ(t) = j ĥ i for v j t < v j+ i= H(t) = log{s(t)} log{ŝ(t)} = for v j t < v j+. j { log( ĥi)} i= j ĥ i = Ĥ(t) i=
18 8 Asymptotical properties of KM estimator As n, Ŝ(t) S(t) in probability. As n, n /2 {Ŝ(t) S(t)} converges to N(0, σ2 (t)) in distribution.
19 9 Asymptotical properties of KM estimator How to estimate the variance of Ŝ(t) ĥi is an estimated probability. The variance of ĥi can be approximated by ĥ i ( ĥi) Y (v i ) = d i(y (v i ) d i ) Y (v i ) 3 ĥi and ĥj are asymptotically independent (why?).
20 20 Asymptotical variance of KM estimator For v j t < v j+ (lnŝ(t) ) : Var j Var ( ) ln( ĥi) i= δ method = j ) Var (ĥi i= ( ĥi) 2 = j i= d i Y (v i )(Y (v i ) d i ).
21 2 Asymptotical variance of KM estimator Using the δ-method ) Var (Ŝ(t) (lnŝ(t) ) ) 2 Var (e lnŝ(t) (lnŝ(t) ) = Ŝ(t)2 Var Ŝ(t)2 j = ˆσ 2 (t) i= d i Y (v i )(Y (v i ) d i ) (v j t < v j+ ) Greenwood s formula
22 22 By-Product The by-product of the Greenwood s formula is the variance estimator for Nelson-Aalen Estimator: var(ĥ(t)) = j i= d i Y (v i )(Y (v i ) d i ), v j t < v j+
23 23 Confidence Interval Ŝ(t) ±.96ˆσ(t), drawbacks? By δ-method var(log( log(ŝ(t)))) = ˆσ 2 (t) (log(ŝ(t)))2ŝ(t)2 The confidence interval for Ŝ(t) [exp{ e log( log(ŝ(t))).96σ(t) ˆ log(ŝ(t))ŝ(t) }, exp{ e log( log(ŝ(t)))+.96ˆσ(t) log(ŝ(t))ŝ(t) } ]
24 24 Median survival time How to estimate the median survival time Solving Ŝ(ˆt M ) = /2, not always solvable! How to construct the CI for the median survival time? Suppose that. pr(ŝl(t) < S(t)) = pr(ŝu(t) > S(t)) = Ŝ L (ˆt ML ) = Ŝ U (ˆt MU ) = The confidence interval for t M is [ˆt ML, ˆt MU ].
25 25 CI for median survival time
26 26 Median survival time = pr(ŝl(t M ) < S(t M )) = pr(ŝl(t M ) < 0.5) = pr(ŝl(t M ) < ŜL(ˆt ML )) = pr(t M ˆt ML ) = pr(ŝu(t M ) > S(t M )) = pr(ŝu(t M ) > 0.5) = pr(ŝu(t M ) > ŜU(ˆt MU )) = pr(t M ˆt MU )
27 27 Restricted mean survival time The area under the survival curve is a nice summary for the curve The AUC µ = τ 0 S(t)dt : = τs(τ) + τ 0 µ = ts(t) tf(t)dt = µ can be estimated as τ 0 0 τ 0 + τ 0 tf(t)dt min(t, τ)f(t)dt = E{min(t, τ)} Ŝ(t)dt.
28 28 Restricted mean survival time The restricted mean survival time E{min(T, τ)} can also be estimated as ˆµ IP W = n n i= δ i + ( δ i )I(U i τ) T i τ, Ŝ C (T i τ) where ŜC( ) is a consistent estimator of the survival function of the censoring time C. (How?) Rational [ ] I(C i τ T i ) E T i τ T i Ŝ C (T i τ) (T i τ) P (C i τ T i T i ) S C (T i τ) = T i τ This type of estimator is called the inverse probability weighting estimator
29 29 Restricted mean survival time ˆµ and ˆµ IP W are equivalent! The variance of the estimated area under the survival curve is complicated (the derivation will be given later). n τ 0 { τ t S(u)du} 2 h(t)dt. P (U t)
30 30 Area under the KM curve
31 3 Model Checking Under the exponential distribution: log(s(t)) = λt Plot t vs. pattern. log(ŝ(t)) to visually check the expected linear Under the Weibull distribution log( log(s(t))) = p log(λ) + p log(t) Plot log(t) vs. log( log(ŝ(t))) to visually check the expected linear pattern.
32 32 Model Checking log( log(s(t))) log(t)
33 33 Coming Soon The distribution of Ŝ(t) as a process.
Survival 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 informationStatistics 262: Intermediate Biostatistics Non-parametric Survival Analysis
Statistics 262: Intermediate Biostatistics Non-parametric Survival Analysis Jonathan Taylor & Kristin Cobb Statistics 262: Intermediate Biostatistics p.1/?? Overview of today s class Kaplan-Meier Curve
More informationSTAT Sample Problem: General Asymptotic Results
STAT331 1-Sample Problem: General Asymptotic Results In this unit we will consider the 1-sample problem and prove the consistency and asymptotic normality of the Nelson-Aalen estimator of the cumulative
More informationEstimation for Modified Data
Definition. Estimation for Modified Data 1. Empirical distribution for complete individual data (section 11.) An observation X is truncated from below ( left truncated) at d if when it is at or below d
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 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 informationChapter 4 Fall Notations: t 1 < t 2 < < t D, D unique death times. d j = # deaths at t j = n. Y j = # at risk /alive at t j = n
Bios 323: Applied Survival Analysis Qingxia (Cindy) Chen Chapter 4 Fall 2012 4.2 Estimators of the survival and cumulative hazard functions for RC data Suppose X is a continuous random failure time with
More informationLecture 3. Truncation, length-bias and prevalence sampling
Lecture 3. Truncation, length-bias and prevalence sampling 3.1 Prevalent sampling Statistical techniques for truncated data have been integrated into survival analysis in last two decades. Truncation in
More informationExercises. (a) Prove that m(t) =
Exercises 1. Lack of memory. Verify that the exponential distribution has the lack of memory property, that is, if T is exponentially distributed with parameter λ > then so is T t given that T > t for
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 informationST745: Survival Analysis: Nonparametric methods
ST745: Survival Analysis: Nonparametric methods Eric B. Laber Department of Statistics, North Carolina State University February 5, 2015 The KM estimator is used ubiquitously in medical studies to estimate
More informationTwo-stage Adaptive Randomization for Delayed Response in Clinical Trials
Two-stage Adaptive Randomization for Delayed Response in Clinical Trials Guosheng Yin Department of Statistics and Actuarial Science The University of Hong Kong Joint work with J. Xu PSI and RSS Journal
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 informationLecture 6 PREDICTING SURVIVAL UNDER THE PH MODEL
Lecture 6 PREDICTING SURVIVAL UNDER THE PH MODEL The Cox PH model: λ(t Z) = λ 0 (t) exp(β Z). How do we estimate the survival probability, S z (t) = S(t Z) = P (T > t Z), for an individual with covariates
More informationPart III. Hypothesis Testing. III.1. Log-rank Test for Right-censored Failure Time Data
1 Part III. Hypothesis Testing III.1. Log-rank Test for Right-censored Failure Time Data Consider a survival study consisting of n independent subjects from p different populations with survival functions
More informationAnalysis of Time-to-Event Data: Chapter 2 - Nonparametric estimation of functions of survival time
Analysis of Time-to-Event Data: Chapter 2 - Nonparametric estimation of functions of survival time Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term
More informationConstrained estimation for binary and survival data
Constrained estimation for binary and survival data Jeremy M. G. Taylor Yong Seok Park John D. Kalbfleisch Biostatistics, University of Michigan May, 2010 () Constrained estimation May, 2010 1 / 43 Outline
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 informationSurvival Distributions, Hazard Functions, Cumulative Hazards
BIO 244: Unit 1 Survival Distributions, Hazard Functions, Cumulative Hazards 1.1 Definitions: The goals of this unit are to introduce notation, discuss ways of probabilistically describing the distribution
More informationOther Survival Models. (1) Non-PH models. We briefly discussed the non-proportional hazards (non-ph) model
Other Survival Models (1) Non-PH models We briefly discussed the non-proportional hazards (non-ph) model λ(t Z) = λ 0 (t) exp{β(t) Z}, where β(t) can be estimated by: piecewise constants (recall how);
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 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 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 informationTypical Survival Data Arising From a Clinical Trial. Censoring. The Survivor Function. Mathematical Definitions Introduction
Outline CHL 5225H Advanced Statistical Methods for Clinical Trials: Survival Analysis Prof. Kevin E. Thorpe Defining Survival Data Mathematical Definitions Non-parametric Estimates of Survival Comparing
More informationEfficient Semiparametric Estimators via Modified Profile Likelihood in Frailty & Accelerated-Failure Models
NIH Talk, September 03 Efficient Semiparametric Estimators via Modified Profile Likelihood in Frailty & Accelerated-Failure Models Eric Slud, Math Dept, Univ of Maryland Ongoing joint project with Ilia
More information1 Glivenko-Cantelli type theorems
STA79 Lecture Spring Semester Glivenko-Cantelli type theorems Given i.i.d. observations X,..., X n with unknown distribution function F (t, consider the empirical (sample CDF ˆF n (t = I [Xi t]. n Then
More informationChapter 17. Failure-Time Regression Analysis. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chapter 17 Failure-Time Regression Analysis William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on the authors
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 informationDAGStat Event History Analysis.
DAGStat 2016 Event History Analysis Robin.Henderson@ncl.ac.uk 1 / 75 Schedule 9.00 Introduction 10.30 Break 11.00 Regression Models, Frailty and Multivariate Survival 12.30 Lunch 13.30 Time-Variation and
More informationApplication of Time-to-Event Methods in the Assessment of Safety in Clinical Trials
Application of Time-to-Event Methods in the Assessment of Safety in Clinical Trials Progress, Updates, Problems William Jen Hoe Koh May 9, 2013 Overview Marginal vs Conditional What is TMLE? Key Estimation
More informationIn contrast, parametric techniques (fitting exponential or Weibull, for example) are more focussed, can handle general covariates, but require
Chapter 5 modelling Semi parametric We have considered parametric and nonparametric techniques for comparing survival distributions between different treatment groups. Nonparametric techniques, such as
More informationPENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA
PENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA Kasun Rathnayake ; A/Prof Jun Ma Department of Statistics Faculty of Science and Engineering Macquarie University
More information8. Parametric models in survival analysis General accelerated failure time models for parametric regression
8. Parametric models in survival analysis 8.1. General accelerated failure time models for parametric regression The accelerated failure time model Let T be the time to event and x be a vector of covariates.
More informationSurvival Analysis: Counting Process and Martingale. Lu Tian and Richard Olshen Stanford University
Survival Analysis: Counting Process and Martingale Lu Tian and Richard Olshen Stanford University 1 Lebesgue-Stieltjes Integrals G( ) is a right-continuous step function having jumps at x 1, x 2,.. b f(x)dg(x)
More informationST495: Survival Analysis: Maximum likelihood
ST495: Survival Analysis: Maximum likelihood Eric B. Laber Department of Statistics, North Carolina State University February 11, 2014 Everything is deception: seeking the minimum of illusion, keeping
More informationCIMAT Taller de Modelos de Capture y Recaptura Known Fate Survival Analysis
CIMAT Taller de Modelos de Capture y Recaptura 2010 Known Fate urvival Analysis B D BALANCE MODEL implest population model N = λ t+ 1 N t Deeper understanding of dynamics can be gained by identifying variation
More informationSurvival Analysis. STAT 526 Professor Olga Vitek
Survival Analysis STAT 526 Professor Olga Vitek May 4, 2011 9 Survival Data and Survival Functions Statistical analysis of time-to-event data Lifetime of machines and/or parts (called failure time analysis
More informationNonparametric Model Construction
Nonparametric Model Construction Chapters 4 and 12 Stat 477 - Loss Models Chapters 4 and 12 (Stat 477) Nonparametric Model Construction Brian Hartman - BYU 1 / 28 Types of data Types of data For non-life
More informationEconomics 508 Lecture 22 Duration Models
University of Illinois Fall 2012 Department of Economics Roger Koenker Economics 508 Lecture 22 Duration Models There is considerable interest, especially among labor-economists in models of duration.
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 informationPractice Exam 1. (A) (B) (C) (D) (E) You are given the following data on loss sizes:
Practice Exam 1 1. Losses for an insurance coverage have the following cumulative distribution function: F(0) = 0 F(1,000) = 0.2 F(5,000) = 0.4 F(10,000) = 0.9 F(100,000) = 1 with linear interpolation
More informationProduct-limit estimators of the survival function with left or right censored data
Product-limit estimators of the survival function with left or right censored data 1 CREST-ENSAI Campus de Ker-Lann Rue Blaise Pascal - BP 37203 35172 Bruz cedex, France (e-mail: patilea@ensai.fr) 2 Institut
More informationPackage Rsurrogate. October 20, 2016
Type Package Package Rsurrogate October 20, 2016 Title Robust Estimation of the Proportion of Treatment Effect Explained by Surrogate Marker Information Version 2.0 Date 2016-10-19 Author Layla Parast
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 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 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 informationLecture 2: Martingale theory for univariate survival analysis
Lecture 2: Martingale theory for univariate survival analysis In this lecture T is assumed to be a continuous failure time. A core question in this lecture is how to develop asymptotic properties when
More informationPart III Measures of Classification Accuracy for the Prediction of Survival Times
Part III Measures of Classification Accuracy for the Prediction of Survival Times Patrick J Heagerty PhD Department of Biostatistics University of Washington 102 ISCB 2010 Session Three Outline Examples
More informationβ j = coefficient of x j in the model; β = ( β1, β2,
Regression Modeling of Survival Time Data Why regression models? Groups similar except for the treatment under study use the nonparametric methods discussed earlier. Groups differ in variables (covariates)
More informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS3301 / MAS8311 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-10 1 13 The Cox proportional hazards model 13.1 Introduction In the
More informationA COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky
A COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky Empirical likelihood with right censored data were studied by Thomas and Grunkmier (1975), Li (1995),
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 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 informationECON 721: Lecture Notes on Duration Analysis. Petra E. Todd
ECON 721: Lecture Notes on Duration Analysis Petra E. Todd Fall, 213 2 Contents 1 Two state Model, possible non-stationary 1 1.1 Hazard function.......................... 1 1.2 Examples.............................
More informationDynamic Models Part 1
Dynamic Models Part 1 Christopher Taber University of Wisconsin December 5, 2016 Survival analysis This is especially useful for variables of interest measured in lengths of time: Length of life after
More informationStatistical Inference and Methods
Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 31st January 2006 Part VI Session 6: Filtering and Time to Event Data Session 6: Filtering and
More informationSTAT 6350 Analysis of Lifetime Data. Failure-time Regression Analysis
STAT 6350 Analysis of Lifetime Data Failure-time Regression Analysis Explanatory Variables for Failure Times Usually explanatory variables explain/predict why some units fail quickly and some units survive
More informationLet us use the term failure time to indicate the time of the event of interest in either a survival analysis or reliability analysis.
10.2 Product-Limit (Kaplan-Meier) Method Let us use the term failure time to indicate the time of the event of interest in either a survival analysis or reliability analysis. Let T be a continuous random
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 information18.465, further revised November 27, 2012 Survival analysis and the Kaplan Meier estimator
18.465, further revised November 27, 2012 Survival analysis and the Kaplan Meier estimator 1. Definitions Ordinarily, an unknown distribution function F is estimated by an empirical distribution function
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 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 informationEmpirical likelihood ratio with arbitrarily censored/truncated data by EM algorithm
Empirical likelihood ratio with arbitrarily censored/truncated data by EM algorithm Mai Zhou 1 University of Kentucky, Lexington, KY 40506 USA Summary. Empirical likelihood ratio method (Thomas and Grunkmier
More informationand Comparison with NPMLE
NONPARAMETRIC BAYES ESTIMATOR OF SURVIVAL FUNCTIONS FOR DOUBLY/INTERVAL CENSORED DATA and Comparison with NPMLE Mai Zhou Department of Statistics, University of Kentucky, Lexington, KY 40506 USA http://ms.uky.edu/
More informationA CHANGE-POINT MODEL FOR SURVIVAL DATA WITH LONG-TERM SURVIVORS
Statistica Sinica 19 (2009), 377-390 A CHANGE-POINT MODEL FOR SURVIVAL DATA WITH LONG-TERM SURVIVORS Xiaobing Zhao 1, Xianyi Wu 2 and Xian Zhou 3 1 Jiang Nan University, 2 East China Normal University
More informationEmpirical Likelihood in Survival Analysis
Empirical Likelihood in Survival Analysis Gang Li 1, Runze Li 2, and Mai Zhou 3 1 Department of Biostatistics, University of California, Los Angeles, CA 90095 vli@ucla.edu 2 Department of Statistics, The
More informationNew mixture models and algorithms in the mixtools package
New mixture models and algorithms in the mixtools package Didier Chauveau MAPMO - UMR 6628 - Université d Orléans 1ères Rencontres BoRdeaux Juillet 2012 Outline 1 Mixture models and EM algorithm-ology
More informationHacettepe Journal of Mathematics and Statistics Volume 45 (5) (2016), Abstract
Hacettepe Journal of Mathematics and Statistics Volume 45 (5) (2016), 1605 1620 Comparing of some estimation methods for parameters of the Marshall-Olkin generalized exponential distribution under progressive
More informationQuantile Regression for Residual Life and Empirical Likelihood
Quantile Regression for Residual Life and Empirical Likelihood Mai Zhou email: mai@ms.uky.edu Department of Statistics, University of Kentucky, Lexington, KY 40506-0027, USA Jong-Hyeon Jeong email: jeong@nsabp.pitt.edu
More informationSTAT Section 2.1: Basic Inference. Basic Definitions
STAT 518 --- Section 2.1: Basic Inference Basic Definitions Population: The collection of all the individuals of interest. This collection may be or even. Sample: A collection of elements of the population.
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 informationImproving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates
Improving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates Anastasios (Butch) Tsiatis Department of Statistics North Carolina State University http://www.stat.ncsu.edu/
More informationIntroduction to Statistical Analysis
Introduction to Statistical Analysis Changyu Shen Richard A. and Susan F. Smith Center for Outcomes Research in Cardiology Beth Israel Deaconess Medical Center Harvard Medical School Objectives Descriptive
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 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 informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 24 Paper 153 A Note on Empirical Likelihood Inference of Residual Life Regression Ying Qing Chen Yichuan
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 informationEfficiency of an Unbalanced Design in Collecting Time to Event Data with Interval Censoring
University of South Florida Scholar Commons Graduate Theses and Dissertations Graduate School 11-10-2016 Efficiency of an Unbalanced Design in Collecting Time to Event Data with Interval Censoring Peiyao
More informationInvestigation of goodness-of-fit test statistic distributions by random censored samples
d samples Investigation of goodness-of-fit test statistic distributions by random censored samples Novosibirsk State Technical University November 22, 2010 d samples Outline 1 Nonparametric goodness-of-fit
More informationMeei Pyng Ng 1 and Ray Watson 1
Aust N Z J Stat 444), 2002, 467 478 DEALING WITH TIES IN FAILURE TIME DATA Meei Pyng Ng 1 and Ray Watson 1 University of Melbourne Summary In dealing with ties in failure time data the mechanism by which
More informationLog-linearity for Cox s regression model. Thesis for the Degree Master of Science
Log-linearity for Cox s regression model Thesis for the Degree Master of Science Zaki Amini Master s Thesis, Spring 2015 i Abstract Cox s regression model is one of the most applied methods in medical
More informationBARTLETT IDENTITIES AND LARGE DEVIATIONS IN LIKELIHOOD THEORY 1. By Per Aslak Mykland University of Chicago
The Annals of Statistics 1999, Vol. 27, No. 3, 1105 1117 BARTLETT IDENTITIES AND LARGE DEVIATIONS IN LIKELIHOOD THEORY 1 By Per Aslak Mykland University of Chicago The connection between large and small
More informationImproved Simultaneous Estimation of Location and System Reliability via Shrinkage Ideas
University of South Carolina Scholar Commons Theses and Dissertations 2017 Improved Simultaneous Estimation of Location and System Reliability via Shrinkage Ideas Beidi Qiang University of South Carolina
More informationCox s proportional hazards model and Cox s partial likelihood
Cox s proportional hazards model and Cox s partial likelihood Rasmus Waagepetersen October 12, 2018 1 / 27 Non-parametric vs. parametric Suppose we want to estimate unknown function, e.g. survival function.
More informationEmpirical Likelihood
Empirical Likelihood Patrick Breheny September 20 Patrick Breheny STA 621: Nonparametric Statistics 1/15 Introduction Empirical likelihood We will discuss one final approach to constructing confidence
More informationAnalysis of Time-to-Event Data: Chapter 4 - Parametric regression models
Analysis of Time-to-Event Data: Chapter 4 - Parametric regression models Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term 2018/19 1/25 Right censored
More informationOptimum Hybrid Censoring Scheme using Cost Function Approach
Optimum Hybrid Censoring Scheme using Cost Function Approach Ritwik Bhattacharya 1, Biswabrata Pradhan 1, Anup Dewanji 2 1 SQC and OR Unit, Indian Statistical Institute, 203, B. T. Road, Kolkata, PIN-
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 informationPhD course in Advanced survival analysis. One-sample tests. Properties. Idea: (ABGK, sect. V.1.1) Counting process N(t)
PhD course in Advanced survival analysis. (ABGK, sect. V.1.1) One-sample tests. Counting process N(t) Non-parametric hypothesis tests. Parametric models. Intensity process λ(t) = α(t)y (t) satisfying Aalen
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 informationFundamentals of Reliability Engineering and Applications
Fundamentals of Reliability Engineering and Applications E. A. Elsayed elsayed@rci.rutgers.edu Rutgers University Quality Control & Reliability Engineering (QCRE) IIE February 21, 2012 1 Outline Part 1.
More informationChapter 15. System Reliability Concepts and Methods. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chapter 15 System Reliability Concepts and Methods William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on
More informationYou know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?
You know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?) I m not goin stop (What?) I m goin work harder (What?) Sir David
More informationTHESIS for the degree of MASTER OF SCIENCE. Modelling and Data Analysis
PROPERTIES OF ESTIMATORS FOR RELATIVE RISKS FROM NESTED CASE-CONTROL STUDIES WITH MULTIPLE OUTCOMES (COMPETING RISKS) by NATHALIE C. STØER THESIS for the degree of MASTER OF SCIENCE Modelling and Data
More informationAFT Models and Empirical Likelihood
AFT Models and Empirical Likelihood Mai Zhou Department of Statistics, University of Kentucky Collaborators: Gang Li (UCLA); A. Bathke; M. Kim (Kentucky) Accelerated Failure Time (AFT) models: Y = log(t
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 informationDuration Analysis. Joan Llull
Duration Analysis Joan Llull Panel Data and Duration Models Barcelona GSE joan.llull [at] movebarcelona [dot] eu Introduction Duration Analysis 2 Duration analysis Duration data: how long has an individual
More informationMaximum likelihood estimation of a log-concave density based on censored data
Maximum likelihood estimation of a log-concave density based on censored data Dominic Schuhmacher Institute of Mathematical Statistics and Actuarial Science University of Bern Joint work with Lutz Dümbgen
More informationConsistency of the Semi-parametric MLE under the Cox Model with Linearly Time-dependent Covariates and Interval-censored Data
Journal of Advanced Statistics, Vol. 4, No. 1, March 2019 https://dx.doi.org/10.22606/jas.2019.41001 1 Consistency of the Semi-parametric MLE under the Cox Model with Linearly Time-dependent Covariates
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 information