DAGStat Event History Analysis.
|
|
- Albert Ward
- 6 years ago
- Views:
Transcription
1 DAGStat 2016 Event History Analysis 1 / 75
2 Schedule 9.00 Introduction Break Regression Models, Frailty and Multivariate Survival Lunch Time-Variation and Dynamic Covariates Break Competing Risks End Each session will consist of 45 minutes talks and 45 minutes computing exercises 2 / 75
3 Session 1: Introduction 1 Types of data 2 Survival analysis recap 3 Counting processes 4 Nelson-Aalen estimator 5 In R 3 / 75
4 Single event survival 1 : leukaemia data X X X X C X X C C X Years 1 Henderson et al, JASA, / 75
5 Competing risks: transplant data R C R D C R C C D R Years 5 / 75
6 Clustered survival 2 : retinipathy data X X X C X C X X X X C X C X C C C Months 2 Huster et al, Biometrics, / 75
7 Recurrent events I: ships data X X X X X X X X X X X X X C X X X X X C Age 7 / 75
8 Recurrent events II: more ships data X X X X X X X X X X X X X X X X C X X C X X X X X X X X C X C X X X X C Age 8 / 75
9 Complex event history 3 : diarrhoea data x xx xx O O x O x xx x x x x x x x O x x x x x x 3 Borgan et al, Scan J, Time 9 / 75
10 Single event survival definitions Time to event T (assumed continuous here) Cdf F (t) = P(T t), pdf f (t) = df (t)/dt Survival function S(t) = P(T > t) = 1 F (t) Hazard function Approximation Cumulative hazard α(t) = f (t) S(t) P (T [t, t + dt) T t) α(t)dt A(t) = t 0 α(u)du S(t) = exp{ A(t)} 10 / 75
11 Single event survival data basics Independent right-censoring Traditional notation: data (t i, δ i ) with δ i = 1 if event, δ i = 0 if censored Kaplan-Meier estimator 1 t < smallest observed failure time Ŝ(t) = ( ) i:t i t 1 d(t i ) n(t i ) otherwise (n(t i ), d(t i ) number at risk and number of events at t i ) Log-rank test (two groups) based on U = ( ) n 1 (t 1 ) d 1 (t i ) n 1 (t i ) + n 2 (t i ) {d 1(t 1 ) + d 2 (t i )} i Cox model α(t x) = α 0 (t)e βx 11 / 75
12 Counting processes N(t) = number of events that have occurred up to and including time t dn(t) = number of events in [t, t + dt) (0 or 1) History to t: F t History to just before t: F t Intensity P (dn(t) = 1 F t ) λ(t) = lim dt 0 dt Approximations P (dn(t) = 1 F t ) λ(t)dt E[dN(t) F t ] λ(t)dt Adding over lots of small intervals E[N(t)] = t 0 λ(u)du = Λ(t) 12 / 75
13 Martingales in under 140 characters M(t) = N(t) Λ(t) is a martingale E[M(t) F u ] = M(u) If H(t) is predictable W (t) = is a martingale Approximately t Var (W (t)) = 0 H(u)dM(u) t 0 H 2 (u)dλ(u) 13 / 75
14 Observed data and independent censoring At risk { 1 at risk just before t Y (t) = 0 otherwise α(t): underlying intensity (strictly α(t F t )) λ(t) = Y (t)α(t): intensity of observed counting process Cumulative versions A(t) = t 0 α(u)du Λ(t) = t 0 Independent censoring α(t + F t, Y (t)) = α(t + F t ) λ(u)du 14 / 75
15 Nelson-Aalen estimator Sample size n, no tied event times Now Informally N(t) = ˆα(t) = But for consistency n N i (t) Y (t) = i=1 n Y i (t) i=1 1 Y (t) dn(t) = 1, Y (t) > 0 0 dn(t) = 0 or Y (t) = 0 Â(t) = t where J(u) = I (Y (u) > 0) This is the Nelson-Aalen estimator 0 J(u) Y (u) dn(u) 15 / 75
16 Variance of Nelson-Aalen estimator Nelson-Aalen Â(t) = t We know M(t) = N(t) Λ(t) So And Estimated by Â(t) = t 0 0 J(u) Y (u) dn(u) J(u) t Y (u) dλ(u) + 0 ) t Var (Â(t) = 0 ) t Var ˆ (Â(t) = 0 J(u) Y (u) dm(u) J(u) Y 2 (u) dλ(u) J(u) Y 2 (u) dn(u) 16 / 75
17 In R Standard survival For counting processes Surv(time,status) Surv(time1,time2,status) (Intervals open on left, closed on right) Eg individal i = 1 has events at times 17 and 26 months and is right-censored at 36 months Consider as Old id New id time1 time2 status / 75
18 18 / 75
19 Session 2: Regression Models, Frailty and Multivariate Survival 1 Likelihood construction 2 Cox and Aalen models 3 Frailty for single-event survival 4 Frailty for clustered data 19 / 75
20 Likelihood for single event survival Observation t, δ Likelihood contribution f (t) δ = 1 L(t, δ) = S(t) δ = 0 Since α(t) = f (t)/s(t) L(t, δ) = α(t) δ S(t) = α(t) δ exp{ A(t)} where A(t) = t 0 α(u)du 20 / 75
21 Alternative derivation 0 Time X (t, δ = 1) Intervals I 1, I 2,..., I K of length dt, boundaries 0 = t 0 < t 1 < t 2... < t K = t L(t, δ = 1) P(No event in I 1 ) P(No event in I 2 past)... P(Event in I K past) = K 1 {1 α(t j )dt} α(t)dt α(t) exp{ A(t)} j=0 Similar argument if δ = 0 21 / 75
22 Likelihood for recurrent events 0 Time X X O Events (t 1, t 2 ) Censored t L(data) P(No event in I past) Empty intervals P(Event in I past) Occupied intervals α(t 1 )α(t 2 ) exp{ A(t)} 22 / 75
23 Likelihood for event history data In general L(data) = = K 1 k=0 K 1 k=0 P(data in[t k, t k + dt) F tk ) {P(events of interest in [t k, t k + dt) F tk ) P(other data in [t k, t k + dt) events of interest in [t k, t k + dt), F tk )} A partial likelihood is L(data) = K 1 k=0 P(events of interest in [t k, t k + dt) F tk ) 23 / 75
24 Regression models I: Cox Proportional Hazards α i (t x i ) = α 0 (t)e βx i Partial likelihood ( ) e βt x i L = i:event at t i j R i e βt x j where R i = R(t i ) = {k : Y k (t i ) = 1} is risk set Same works more generally, with α i (t F t ) = α(t x i ) 24 / 75
25 Regression models II: Aalen additive α i (t F t ) = β 0 (t) + β 1 (t)x 1 + β 2 (t)x 2... Inference on cumulative coefficients B j (t) = t 0 β j (u)du BASELINE X Day Day 25 / 75
26 Missing Data 26 / 75
27 Frailty Assume α(t x) = α 0 (t) exp{βx} α(t x obs, x miss ) = α 0 (t) exp{β obs x obs + β miss x miss } α(t x obs, W ) = α 0 (t) exp{β obs + W } α(t x obs, Z) = Zα 0 (t) exp{β obs } Z not observed Need to work with α(t x obs ) Assume E[Z] = 1 for identifiability 27 / 75
28 Frailty distributions Gamma. Can calculate S(t x) and α(t x) explicitly Log-Normal, so W = log Z Normal. Motivated by CLT for missing data Positive stable. Defined by Stays in PH family E[exp{ uz}] = exp{ u ν } 28 / 75
29 Frailty distributions Density Log normal Gamma Positive stable log(frailty) 29 / 75
30 Effect of frailty Survival No frailty Log normal frailty Gamma frailty Positive stable frailty Time 30 / 75
31 Leukaemia data: no frailty >coxph(surv(time,cens) age+male+wbc+dep,data=leukaemia) coef se(coef) z p age e+00 male e-01 wbc e-12 dep e / 75
32 Leukaemia data:frailty (n=1043) >coxph(surv(time,cens) age+male+wbc+dep +frailty(1:n),data=leukaemia) coef se(coef) se2 Chisq DF p age e+00 male e-01 wbc e-15 dep e-05 frailty(1:n) e+00 Variance of random effect= Warning message: / 75
33 Shared frailty for clustered data Eg pairs (T 1, T 2 ), (x 1, x 2 ) Hazards α 1 (t x 1, Z) = Zα 0 (t) exp{β 1 x 1 } α 2 (t x 2, Z) = Zα 0 (t) exp{β 2 x 2 } 33 / 75
34 34 / 75
35 Session 3: Time-Variation and Dynamic Covariates 1 Time-varying covariates 2 Time varying effects 3 Quick and dirty 4 Aalen revisited 5 Dynamic covariates 35 / 75
36 Time varying covariates Cox: α i (t x i ) = α 0 (t)e βx i Partial likelihood ( ) e βt x i L = i:event at t i j R i e βt x j Also works with Cox: α i (t x i ) = α 0 (t)e βx i (t) ( ) e βt x i (t i ) L = i:event at t i j R i e βt x j (t i ) 36 / 75
37 Ships data 1 Sales of 3908 ships 2 Fixed covariates: type (three), weight, speed 3 Time-varying covariates: owner number, price index 37 / 75
38 id type dwt speed owner index start stop cens / 75
39 fit1=coxph(surv(start,stop,cens) as.factor(type)*owner +dwt+speed+index,data=shipslong) coef se(coef) z as.factor(type) as.factor(type) owner dwt e e speed index as.factor(type)2:owner as.factor(type)3:owner / 75
40 Time varying effects Cox: α i (t x i ) = α 0 (t)e β(t)x i Various methods to estimate smooth β(t) Quick and dirty: changepoint α 0 (t)e β 1x t τ, α(t x) = α 0 (t)e β 2x t > τ. L(β 1, β 2 ) = i:t i τ ( e β 1x i j R i e β 1x j ) δi i:t i >τ ( e β 2x i j R i e β 2x j ) δi 40 / 75
41 Time-varying effects: quick and dirty 1 Split the time axis at a point τ 2 Fit a Cox model to times before τ. Get estimates ˆβ 1. 3 Fit a Cox model to times after τ. Get estimates ˆβ 2. 4 Do ˆβ 1 and ˆβ 2 seem to be very different? 5 Try various τ, compare (log) likelihoods 41 / 75
42 Leukaemia data >tau=1 >fit0=coxph(surv(time,cens) age+male+wbc+dep,data=leukaemia) >leuk1=leukaemia >leuk2=leukaemia >i1=leukaemia$time>tau >i2=leukaemia$time<=tau >leuk1$cens[i1]=0 >leuk2$cens[i2]=0 >fit1=coxph(surv(time,cens) age+male+wbc+dep,data=leuk1) >fit2=coxph(surv(time,cens) age+male+wbc+dep,data=leuk2) 42 / 75
43 Leukaemia data >fit0$loglik [1] > fit1$loglik [1] > fit2$loglik [1] > fit1$loglik[2]+fit2$loglik[2] [1] / 75
44 fit1 coef se z age male wbc dep fit2 coef se z age male wbc dep / 75
45 Aalen additive α i (t F t ) = β 0 (t) + β 1 (t)x i1 (t) + β 2 (t)x i2 (t)... Inference on cumulative coefficients B j (t) = t 0 β j (u)du Martingale theory for (cumulative) standard errors 45 / 75
46 Leukaemia data age male Time Time wbc dep Time Time 46 / 75
47 Recap Counting process N(t) History F t Intensity α(t F t ) At-risk Y (t) E[N(t)] = t 0 Y (u)α(u F u )du Static models α(t F t ) = α(t F 0 ) 47 / 75
48 Recap Counting process N(t) History F t Intensity α(t F t ) At-risk Y (t) E[N(t)] = t 0 Y (u)α(u F u )du Static models α(t F t ) = α(t F 0 ) Cox: α(t F 0 ) = α 0 (t)e βx Frailty α(t F 0, Z) = Zα 0 (t)e βx Aalen (constant): α(t F 0 ) = β 0 (t) + β 1 x 1 (t) +... Aalen (varying): α(t F 0 ) = β 0 (t) + β 1 (t)x 1 (t) +... Logistic: α(t F 0 ) = expit{β 0 (t) + β 1 (t)x 1 (t) +...} 47 / 75
49 Dynamic models E[N(t)] = t 0 Y (u)α(u F u )du Static models α(t F t ) = α(t F 0 ) Dynamic models incorporate F t Frailty: α(t F t ) = E[Z F t ] α 0 (t)e βx Dynamic covariate D t = g(f t ) Eg D t = Number of events before t Days at risk before t 48 / 75
50 Diarrhoea data Eg D t =previous episode rate (episodes/time) Example: test for effect of rain-affected accommodation (Wald) Model Rain-affected Previous episode rate No dynamic 3.70 Include D t / 75
51 Fixed X β XY (t) Y t 50 / 75
52 Fixed X β XY.D (t) Y t β XD (t) D t β DY (t) 51 / 75
53 Solution Assume Use D t = X γ t + Z t Ẑ t = D t ˆD t = D t X ˆγ t = D t X ( X T X ) 1 X T D t 52 / 75
54 X β XY (t) Y t Ẑ t β ZY (t) 53 / 75
55 Example: test for effect of rain-affected accommodation D t =previous episode rate (episodes/time) Model Rain-affected Previous episode rate No dynamic 3.70 Include D t Include Ẑ t / 75
56 55 / 75
57 Session 4: Competing Risks 1 Set-up 2 Cumulative incidence function 3 Cause-specific hazards 4 Subdistribution hazards 5 Words from the wise 56 / 75
58 Competing risks Events of more than one type Death from one of several causes Time to first event Latent failure time interpretation Independence not identifiable T = min{t 1, T 2,...} 57 / 75
59 Multistate interpretation C i (t)=state (0,1,2...) of person i at time t T = inf t>0 (C(t) 0) 58 / 75
60 Single-event survival recap S(t) = exp{ A(t)} Kaplan-Meier Ŝ(t) = i:t i t ( 1 d(t ) i) n(t i ) (n(t i ), d(t i ) number at risk and number of events at t i ) Nelson-Aalen Â(t) = i:t i t Assume independent right-censoring d(t i ) n(t i ) 59 / 75
61 Cumulative Incidence Function Two causes from now on T = time in state 0 C = C( ) = C(T ) CIF: F j (t) = P(T t, C = j) Marginal S(t) = 1 F 1 (t) F 2 (t) Estimated by usual Kaplan-Meier with event types pooled T is proper But lim S(t) = 0 t lim F j(t) = P(C = j) < 1 t 60 / 75
62 Naive Kaplan-Meier One cause at a time Treat other causes as censoring For cause j Ŝ j (t) = i:t i t ( 1 d ) j(t i ) n(t i ) (d j (t i ) number of events of type j at t i ) Does not in general estimate 1 F j (t) 61 / 75
63 Kidney transplant data Survival KM for tx failure KM for death Time (years) 62 / 75
64 Kidney transplant data Survival KM for tx failure KM for death Time (years) 63 / 75
65 Transition probabilities and cause-specific hazards P jk (s, t) = P(C(t) = k C(s) = j) P 00 (t) = S(t) α j (t) = lim dt 0 P 0j (t, t + dt)/dt Cumulative cause-specific hazard Can be estimated A j (t) = t 0 Â j (t) = i:t i t α j (u)du d j (t i ) n(t i ) And ˆF j (t) = i:t i t Ŝ(t i ) d j(t i ) n(t i ) 64 / 75
66 Kidney transplant data Survival KM for tx failure CIF for tx failure KM for death 1 CIF for death Time (years) 65 / 75
67 Kidney transplant data Survival KM for tx failure CIF for tx failure 1 KM for death CIF for death Time (years) 66 / 75
68 Categorical covariates CIF via cif=survfit(surv(time,cens,type="mstate") 1,data=kidney) As usual cif=survfit(surv(time,cens,type="mstate") capd,data=kidney) 67 / 75
69 Kidney transplant data: CIF Survival Failure, no CAPD Death,no CAPD Failure, CAPD Death, CAPD Time (years) 68 / 75
70 General covariates Cause-specific hazards can be modelled as usual α j (t x) = α 0j (t) exp{β j x} Estimation as usual 69 / 75
71 coxph(surv(time,censdeath) rage+capd+drmm+bmm,data=kidney) coef se z rage capd drmm bmm coxph(surv(time,censfail) rage+capd+drmm+bmm,data=kidney) coef se(coef) z rage capd drmm bmm / 75
72 Interpretation Depends on all causes F 1 (t) = = = t 0 t 0 t 0 S(u)α 1 (u x)du e { A 1(u x) A 2 (u x)} α 1 (u x)du e { A 01(u)e β 1 x A 02 (u)e β 2 x } α 01 (u)e β 1x du Better to calculate P 0j (0, t) for specific x Aalen-Johansen transition matrix P(t) = (I + da(t)) u t 71 / 75
73 State probabilities Recipient aged 60, no CAPD, drmm=1, bmm=1 P(t) OK Failed Dead Time t 72 / 75
74 Subdistribution hazards Standard survival α(t) = d log(s(t)) dt CIF F j (t) = P(T t, C = j) Define Interpretation α j (t) = d log (1 F j(t)) dt α j (t) P (C(t + dt) = j C(t) j) Compare with cause-specific α j (t) P (C(t + dt) = j C(t) = 0) Fine and Gray use proportional hazards for α j (t) Problem: estimation requires us to assume that individuals in state k j remain in the risk sets for transition to k 73 / 75
75 Other approaches P(T t, C = j) = P(T t C = j)p(c = j) P(T t, C = j) = P(C = j T t)p(t t) Parametric models 74 / 75
76 Final words Andersen & Keiding, SiM, Do not condition on the future 2 Do not regard individuals at risk after they have died 3 Stick to this world 75 / 75
In 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 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 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 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 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 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 informationLecture 7 Time-dependent Covariates in Cox Regression
Lecture 7 Time-dependent Covariates in Cox Regression So far, we ve been considering the following Cox PH model: λ(t Z) = λ 0 (t) exp(β Z) = λ 0 (t) exp( β j Z j ) where β j is the parameter for the the
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 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 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 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 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 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 informationSmoothing the Nelson-Aalen Estimtor Biostat 277 presentation Chi-hong Tseng
Smoothing the Nelson-Aalen Estimtor Biostat 277 presentation Chi-hong seng Reference: 1. Andersen, Borgan, Gill, and Keiding (1993). Statistical Model Based on Counting Processes, Springer-Verlag, p.229-255
More informationPhilosophy and Features of the mstate package
Introduction Mathematical theory Practice Discussion Philosophy and Features of the mstate package Liesbeth de Wreede, Hein Putter Department of Medical Statistics and Bioinformatics Leiden University
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 informationSurvival analysis in R
Survival analysis in R Niels Richard Hansen This note describes a few elementary aspects of practical analysis of survival data in R. For further information we refer to the book Introductory Statistics
More informationLecture 12. Multivariate Survival Data Statistics Survival Analysis. Presented March 8, 2016
Statistics 255 - Survival Analysis Presented March 8, 2016 Dan Gillen Department of Statistics University of California, Irvine 12.1 Examples Clustered or correlated survival times Disease onset in family
More informationADVANCED STATISTICAL ANALYSIS OF EPIDEMIOLOGICAL STUDIES. Cox s regression analysis Time dependent explanatory variables
ADVANCED STATISTICAL ANALYSIS OF EPIDEMIOLOGICAL STUDIES Cox s regression analysis Time dependent explanatory variables Henrik Ravn Bandim Health Project, Statens Serum Institut 4 November 2011 1 / 53
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 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 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 informationDefinitions and examples Simple estimation and testing Regression models Goodness of fit for the Cox model. Recap of Part 1. Per Kragh Andersen
Recap of Part 1 Per Kragh Andersen Section of Biostatistics, University of Copenhagen DSBS Course Survival Analysis in Clinical Trials January 2018 1 / 65 Overview Definitions and examples Simple estimation
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 informationUNIVERSITY OF CALIFORNIA, SAN DIEGO
UNIVERSITY OF CALIFORNIA, SAN DIEGO Estimation of the primary hazard ratio in the presence of a secondary covariate with non-proportional hazards An undergraduate honors thesis submitted to the Department
More informationAnalysis of competing risks data and simulation of data following predened subdistribution hazards
Analysis of competing risks data and simulation of data following predened subdistribution hazards Bernhard Haller Institut für Medizinische Statistik und Epidemiologie Technische Universität München 27.05.2013
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 informationAnalysis of Time-to-Event Data: Chapter 6 - Regression diagnostics
Analysis of Time-to-Event Data: Chapter 6 - Regression diagnostics Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term 2018/19 1/25 Residuals for the
More informationREGRESSION ANALYSIS FOR TIME-TO-EVENT DATA THE PROPORTIONAL HAZARDS (COX) MODEL ST520
REGRESSION ANALYSIS FOR TIME-TO-EVENT DATA THE PROPORTIONAL HAZARDS (COX) MODEL ST520 Department of Statistics North Carolina State University Presented by: Butch Tsiatis, Department of Statistics, NCSU
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 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 informationMultistate models in survival and event history analysis
Multistate models in survival and event history analysis Dorota M. Dabrowska UCLA November 8, 2011 Research supported by the grant R01 AI067943 from NIAID. The content is solely the responsibility of the
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 informationMulti-state Models: An Overview
Multi-state Models: An Overview Andrew Titman Lancaster University 14 April 2016 Overview Introduction to multi-state modelling Examples of applications Continuously observed processes Intermittently observed
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 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 informationOutline. Frailty modelling of Multivariate Survival Data. Clustered survival data. Clustered survival data
Outline Frailty modelling of Multivariate Survival Data Thomas Scheike ts@biostat.ku.dk Department of Biostatistics University of Copenhagen Marginal versus Frailty models. Two-stage frailty models: copula
More informationEstimating transition probabilities for the illness-death model The Aalen-Johansen estimator under violation of the Markov assumption Torunn Heggland
Estimating transition probabilities for the illness-death model The Aalen-Johansen estimator under violation of the Markov assumption Torunn Heggland Master s Thesis for the degree Modelling and Data Analysis
More informationSTAT331 Lebesgue-Stieltjes Integrals, Martingales, Counting Processes
STAT331 Lebesgue-Stieltjes Integrals, Martingales, Counting Processes This section introduces Lebesgue-Stieltjes integrals, and defines two important stochastic processes: a martingale process and a counting
More informationSurvival analysis in R
Survival analysis in R Niels Richard Hansen This note describes a few elementary aspects of practical analysis of survival data in R. For further information we refer to the book Introductory Statistics
More informationMulti-state models: prediction
Department of Medical Statistics and Bioinformatics Leiden University Medical Center Course on advanced survival analysis, Copenhagen Outline Prediction Theory Aalen-Johansen Computational aspects Applications
More informationMultistate models and recurrent event models
Multistate models Multistate models and recurrent event models Patrick Breheny December 10 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/22 Introduction Multistate models In this final lecture,
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 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 informationMultistate models and recurrent event models
and recurrent event models Patrick Breheny December 6 Patrick Breheny University of Iowa Survival Data Analysis (BIOS:7210) 1 / 22 Introduction In this final lecture, we will briefly look at two other
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 informationEmpirical Processes & Survival Analysis. The Functional Delta Method
STAT/BMI 741 University of Wisconsin-Madison Empirical Processes & Survival Analysis Lecture 3 The Functional Delta Method Lu Mao lmao@biostat.wisc.edu 3-1 Objectives By the end of this lecture, you will
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 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 informationGoodness-Of-Fit for Cox s Regression Model. Extensions of Cox s Regression Model. Survival Analysis Fall 2004, Copenhagen
Outline Cox s proportional hazards model. Goodness-of-fit tools More flexible models R-package timereg Forthcoming book, Martinussen and Scheike. 2/38 University of Copenhagen http://www.biostat.ku.dk
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 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 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 informationMultivariate Survival Data With Censoring.
1 Multivariate Survival Data With Censoring. Shulamith Gross and Catherine Huber-Carol Baruch College of the City University of New York, Dept of Statistics and CIS, Box 11-220, 1 Baruch way, 10010 NY.
More information9 Estimating the Underlying Survival Distribution for a
9 Estimating the Underlying Survival Distribution for a Proportional Hazards Model So far the focus has been on the regression parameters in the proportional hazards model. These parameters describe the
More informationFrailty Models and Copulas: Similarities and Differences
Frailty Models and Copulas: Similarities and Differences KLARA GOETHALS, PAUL JANSSEN & LUC DUCHATEAU Department of Physiology and Biometrics, Ghent University, Belgium; Center for Statistics, Hasselt
More informationFollow this and additional works at: Part of the Applied Mathematics Commons
Louisiana State University LSU Digital Commons LSU Master's Theses Graduate School 2004 Cox Regression Model Lindsay Sarah Smith Louisiana State University and Agricultural and Mechanical College, lsmit51@lsu.edu
More informationPart IV Extensions: Competing Risks Endpoints and Non-Parametric AUC(t) Estimation
Part IV Extensions: Competing Risks Endpoints and Non-Parametric AUC(t) Estimation Patrick J. Heagerty PhD Department of Biostatistics University of Washington 166 ISCB 2010 Session Four Outline Examples
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 informationSTAT331. Combining Martingales, Stochastic Integrals, and Applications to Logrank Test & Cox s Model
STAT331 Combining Martingales, Stochastic Integrals, and Applications to Logrank Test & Cox s Model Because of Theorem 2.5.1 in Fleming and Harrington, see Unit 11: For counting process martingales with
More informationCredit risk and survival analysis: Estimation of Conditional Cure Rate
Stochastics and Financial Mathematics Master Thesis Credit risk and survival analysis: Estimation of Conditional Cure Rate Author: Just Bajželj Examination date: August 30, 2018 Supervisor: dr. A.J. van
More informationLecture 7. Proportional Hazards Model - Handling Ties and Survival Estimation Statistics Survival Analysis. Presented February 4, 2016
Proportional Hazards Model - Handling Ties and Survival Estimation Statistics 255 - Survival Analysis Presented February 4, 2016 likelihood - Discrete Dan Gillen Department of Statistics University of
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 informationSurvival Analysis for Case-Cohort Studies
Survival Analysis for ase-ohort Studies Petr Klášterecký Dept. of Probability and Mathematical Statistics, Faculty of Mathematics and Physics, harles University, Prague, zech Republic e-mail: petr.klasterecky@matfyz.cz
More information( t) Cox regression part 2. Outline: Recapitulation. Estimation of cumulative hazards and survival probabilites. Ørnulf Borgan
Outline: Cox regression part 2 Ørnulf Borgan Department of Mathematics University of Oslo Recapitulation Estimation of cumulative hazards and survival probabilites Assumptions for Cox regression and check
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 informationLogistic regression model for survival time analysis using time-varying coefficients
Logistic regression model for survival time analysis using time-varying coefficients Accepted in American Journal of Mathematical and Management Sciences, 2016 Kenichi SATOH ksatoh@hiroshima-u.ac.jp Research
More information9. Estimating Survival Distribution for a PH Model
9. Estimating Survival Distribution for a PH Model Objective: Another Goal of the COX model Estimating the survival distribution for individuals with a certain combination of covariates. PH model assumption:
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 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 informationModelling geoadditive survival data
Modelling geoadditive survival data Thomas Kneib & Ludwig Fahrmeir Department of Statistics, Ludwig-Maximilians-University Munich 1. Leukemia survival data 2. Structured hazard regression 3. Mixed model
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 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 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 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 informationFaculty of Health Sciences. Cox regression. Torben Martinussen. Department of Biostatistics University of Copenhagen. 20. september 2012 Slide 1/51
Faculty of Health Sciences Cox regression Torben Martinussen Department of Biostatistics University of Copenhagen 2. september 212 Slide 1/51 Survival analysis Standard setup for right-censored survival
More informationDynamic Prediction of Disease Progression Using Longitudinal Biomarker Data
Dynamic Prediction of Disease Progression Using Longitudinal Biomarker Data Xuelin Huang Department of Biostatistics M. D. Anderson Cancer Center The University of Texas Joint Work with Jing Ning, Sangbum
More informationSTAT 526 Spring Final Exam. Thursday May 5, 2011
STAT 526 Spring 2011 Final Exam Thursday May 5, 2011 Time: 2 hours Name (please print): Show all your work and calculations. Partial credit will be given for work that is partially correct. Points will
More informationJOINT REGRESSION MODELING OF TWO CUMULATIVE INCIDENCE FUNCTIONS UNDER AN ADDITIVITY CONSTRAINT AND STATISTICAL ANALYSES OF PILL-MONITORING DATA
JOINT REGRESSION MODELING OF TWO CUMULATIVE INCIDENCE FUNCTIONS UNDER AN ADDITIVITY CONSTRAINT AND STATISTICAL ANALYSES OF PILL-MONITORING DATA by Martin P. Houze B. Sc. University of Lyon, 2000 M. A.
More informationBeyond GLM and likelihood
Stat 6620: Applied Linear Models Department of Statistics Western Michigan University Statistics curriculum Core knowledge (modeling and estimation) Math stat 1 (probability, distributions, convergence
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 informationUnderstanding product integration. A talk about teaching survival analysis.
Understanding product integration. A talk about teaching survival analysis. Jan Beyersmann, Arthur Allignol, Martin Schumacher. Freiburg, Germany DFG Research Unit FOR 534 jan@fdm.uni-freiburg.de It is
More informationA multi-state model for the prognosis of non-mild acute pancreatitis
A multi-state model for the prognosis of non-mild acute pancreatitis Lore Zumeta Olaskoaga 1, Felix Zubia Olaskoaga 2, Guadalupe Gómez Melis 1 1 Universitat Politècnica de Catalunya 2 Intensive Care Unit,
More informationResiduals and model diagnostics
Residuals and model diagnostics Patrick Breheny November 10 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/42 Introduction Residuals Many assumptions go into regression models, and the Cox proportional
More informationSurvival Analysis using Bivariate Archimedean Copulas. Krishnendu Chandra
Survival Analysis using Bivariate Archimedean Copulas Krishnendu Chandra Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy under the Executive Committee of the
More informationA note on the decomposition of number of life years lost according to causes of death
A note on the decomposition of number of life years lost according to causes of death Per Kragh Andersen Research Report 12/2 Department of Biostatistics University of Copenhagen A note on the decomposition
More informationPackage crrsc. R topics documented: February 19, 2015
Package crrsc February 19, 2015 Title Competing risks regression for Stratified and Clustered data Version 1.1 Author Bingqing Zhou and Aurelien Latouche Extension of cmprsk to Stratified and Clustered
More informationPart [1.0] Measures of Classification Accuracy for the Prediction of Survival Times
Part [1.0] Measures of Classification Accuracy for the Prediction of Survival Times Patrick J. Heagerty PhD Department of Biostatistics University of Washington 1 Biomarkers Review: Cox Regression Model
More informationLinear rank statistics
Linear rank statistics Comparison of two groups. Consider the failure time T ij of j-th subject in the i-th group for i = 1 or ; the first group is often called control, and the second treatment. Let n
More informationA Comparison of Different Approaches to Nonparametric Inference for Subdistributions
A Comparison of Different Approaches to Nonparametric Inference for Subdistributions Johannes Mertsching Johannes.Mertsching@gmail.com Master Thesis Supervision: dr. Ronald B. Geskus prof. Chris A.J. Klaassen
More informationOn the Breslow estimator
Lifetime Data Anal (27) 13:471 48 DOI 1.17/s1985-7-948-y On the Breslow estimator D. Y. Lin Received: 5 April 27 / Accepted: 16 July 27 / Published online: 2 September 27 Springer Science+Business Media,
More informationOutline. Frailty modelling of Multivariate Survival Data. Clustered survival data. Clustered survival data
Outline Frailty modelling of Multivariate Survival Data Thomas Scheike ts@biostat.ku.dk Department of Biostatistics University of Copenhagen Marginal versus Frailty models. Two-stage frailty models: copula
More informationDirect likelihood inference on the cause-specific cumulative incidence function: a flexible parametric regression modelling approach
Direct likelihood inference on the cause-specific cumulative incidence function: a flexible parametric regression modelling approach Sarwar I Mozumder 1, Mark J Rutherford 1, Paul C Lambert 1,2 1 Biostatistics
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 informationRobust estimates of state occupancy and transition probabilities for Non-Markov multi-state models
Robust estimates of state occupancy and transition probabilities for Non-Markov multi-state models 26 March 2014 Overview Continuously observed data Three-state illness-death General robust estimator Interval
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 informationLecture 9. Statistics Survival Analysis. Presented February 23, Dan Gillen Department of Statistics University of California, Irvine
Statistics 255 - Survival Analysis Presented February 23, 2016 Dan Gillen Department of Statistics University of California, Irvine 9.1 Survival analysis involves subjects moving through time Hazard may
More informationResampling methods for randomly censored survival data
Resampling methods for randomly censored survival data Markus Pauly Lehrstuhl für Mathematische Statistik und Wahrscheinlichkeitstheorie Mathematisches Institut Heinrich-Heine-Universität Düsseldorf Wien,
More informationContinuous case Discrete case General case. Hazard functions. Patrick Breheny. August 27. Patrick Breheny Survival Data Analysis (BIOS 7210) 1/21
Hazard functions Patrick Breheny August 27 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/21 Introduction Continuous case Let T be a nonnegative random variable representing the time to an event
More informationHypothesis Testing Based on the Maximum of Two Statistics from Weighted and Unweighted Estimating Equations
Hypothesis Testing Based on the Maximum of Two Statistics from Weighted and Unweighted Estimating Equations Takeshi Emura and Hisayuki Tsukuma Abstract For testing the regression parameter in multivariate
More information