Lecture 9. Statistics Survival Analysis. Presented February 23, Dan Gillen Department of Statistics University of California, Irvine
|
|
- Emery Walters
- 5 years ago
- Views:
Transcription
1 Statistics Survival Analysis Presented February 23, 2016 Dan Gillen Department of Statistics University of California, Irvine 9.1
2 Survival analysis involves subjects moving through time Hazard may change over time Relative hazard may change over time Factors (covariates) determining relative hazard may change over time Why should we only examine covariates that are fixed at their baseline values? 9.2
3 Examples 1. In a study of Medicare patients, interest is on the effects of heart failure on mortality (age at death) We wish to allow relative hazard to change as subjects experience sentinel events such as: ischemia arrhythmia pump failure Health status of subjects is dynamic 2. In HIV studies of the time from infection (sero-conversion) to AIDS, the risk of AIDS-onset varies dramatically with the subject s CD-4 T-cell count, which is a smooth function of time 3. The relative hazard between high and low risk AML (leukemia) patients may decline with time 9.3
4 These effects can be expressed with time-varying covariates: 1. x i (t) = ischemia i (t) indicator for ith subject for ischemia hx before age t 2. x i (t) = CD4 i (t) ith HIV patient s CD4 count at time t since HIV-seroconversion 3. x i1 = I(risk i = high) = baseline high risk indicator, and x i2 (t) = x i1 log(t) is the interaction of risk status with time Conceptually, these effects are easy to include in a proportional hazards. Why? 9.4
5 The partial likelihood for the proportional hazards with time-fixed covariates is L P = failure times j ( ) exp(β T x (j) ) i R (j) exp(β T x i ) L P compares the covariate x (j) of the failed subject to the covariates x i in the risk set R (j) at time t (j) But, there is absolutely no reason that the covariate x i for the ith subject has to be the same for all times t (j). 9.5
6 The partial likelihood for the proportional hazards with time-varying covariates is L P = failure times j ( ) exp{β T x (j) (t (j) )} i R (j) exp{β T x i (t (j) )} where x i (t (j) ) is the value of the ith subject s covariate at time t (j) L P now compares the covariate x (j) (t (j) ) at time t (j) of the failed subject to the covariates x i (t (j) ) at time t (j) in the risk set R (j) 9.6
7 To formalize, consider the data: (y 1, δ 1, x 1 (t)),, (y n, δ n, x n (t)), where y i = observed follow-up time δ i = failure (vs. censoring) indicator x i (t) = vector covariate function of time t Note that x i (t) is set of values that (may) vary with t. R (j) is the set of subjects at risk for failure at t (j) Q: who is in the risk set? A: the covariate values at risk at t (j) are {x i (t (j) ) : i R (j) } 9.7
8 The partial likelihood is constructed by comparing the risk given x (j) (t (j) ) to the risk given all other x i (t (j) )s for i R (j) : L (j) = Pr{subject with x (j) (t (j) ) fails at t (j) } Pr{ some subject in R (j) failed at t (j) } and L P = L (j) failure times j The is λ i {t x i ( )} = λ 0 (t) exp{β T x i (t)} 9.8
9 Note: The cumulative hazard function for the ith subject is then Λ i {t x i ( )} = = t 0 t 0 λ i {s x i ( )} ds λ 0 (s) exp{β T x i (s)} ds which no longer factors into Λ 0 (t) and exp(β T x i ) The baseline cumulative hazard Λ 0 (t) now refers to a subject with x i (t) = 0 for all t λ0 (t), Λ 0 (t), S 0 (t) may no longer have a clear interpretation Recovering Λi {t x i ( )} is more difficult than before, and perhaps less meaningful 9.9
10 Consider the bone-marrow transplant data (Section 1.3, K & M): time origin is time of transplant failure event is relapse or death Intermediate events that could effect the hazard for death or relapse include: development of acute graft-versus-host disease (agvhd) development of chronic graft-versus-host disease (cgvhd) return of platelet count to self-sustaining level (platelet recovery) 9.10
11 We are interested here in the effect of platelet recovery (iplate), a binary variable taking value 0 at t = 0 and, for some subjects, taking value 1 at some time tplate. Lets look at the data > bmt <- read.csv( " STAT295/Data/bmt.csv" ) > bmt <- bmt[ order(bmt$tnodis), ] > bmt <- bmt[,-1] > bmt <- cbind( 1:dim(bmt)[1], bmt ) > names(bmt)[1] <- "id" 9.11
12 ### Center age of patient and donor for later interpretation > bmt$agep.c = bmt$agep - 28 > bmt$aged.c = bmt$aged - 28 > bmt[1:25,c("id","tnodis","inodis", "iplate", "tplate")] id tnodis inodis iplate tplate
13 Data analysis set-up For a subject who never experiences platelet recovery, iplate= 0, for all time t up to the end of observation tnodis. This is like one observation with one exit time. For a subject experiencing platelet recovery, think of the data as two observations: first, we have an observation with iplate= 0, entering at time t = 0 and being censored at time t =tplate. iplate= 0 for this observation then, we have an observation with iplate= 1, entering the just just after time t =tplate, and exiting at time t =tnodis. iplate= 1 for this observation 9.13
14 Data analysis set-up ### Create a duplicate record for each subject > bmt.tvc <- bmt[rep(1:dim(bmt)[1],each=2),] > id.1 <-!duplicated( bmt.tvc$id ) > id.2 <- duplicated( bmt.tvc$id ) ### Deal with first record, set start, stop, event indicator, ### and covariate value > bmt.tvc$start[ id.1 bmt.tvc$iplate==0 ] <- 0 > bmt.tvc$stop[ id.1 & bmt.tvc$iplate==1 ] <- bmt.tvc[ id.1 & bmt.tvc$iplate==1, ]$tplate > bmt.tvc$inodis[ id.1 & bmt.tvc$iplate==1 ] <- 0 > bmt.tvc$iplate[ id.1 ] <- 0 ### Deal with second record, set start and stop > bmt.tvc$start[ id.2 & bmt.tvc$iplate==1 ] <- bmt.tvc[ id.2 & bmt.tvc$iplate==1, ]$tplate > bmt.tvc$stop[ id.2 ] <- bmt.tvc[ id.2, ]$tnodis 9.14
15 Data analysis set-up ### Remove records with missing stop value ### (these did not change iplate status) > bmt.tvc <- bmt.tvc[!is.na(bmt.tvc$stop), ] > bmt.tvc[1:20,c("id","tnodis", "start", "stop", "inodis", "iplate")] id tnodis start stop inodis iplate
16 Data analysis set-up Interpretation of what we have so far: id indicates multiple observations on a single subject start is the entry time for each observation stop is the exit time for each observation the exit time for the previous observation is used as the entry time for the following observations To create the Surv() response we can now specify Surv(start, stop, event) ### Look at the Surv() object > Surv( bmt.tvc$start, bmt.tvc$stop, bmt.tvc$inodis) [1] ( 0, 1 ] ( 0, 2 ] ( 0, 10 ] ( 0, 16 ] ( 0, 16+] ( 16, 32 ] ( 0, 35 ] 9.16
17 then proceeds by calling coxph() as before, with the expanded" Surv() object ### Fit Cox with time-varying covariate iplate > fit <- coxph( Surv( start, stop, inodis) ~ fab + agep.c*aged.c + factor(g) + imtx + iplate, data=bmt.tvc ) > summary(fit) coef exp(coef) se(coef) z Pr(> z ) fab ** agep.c aged.c factor(g) ** factor(g) imtx iplate ** agep.c:aged.c ** exp(coef) exp(-coef) lower.95 upper.95 fab agep.c aged.c factor(g) factor(g) imtx iplate agep.c:aged.c
18 From the Wald test on the previous page, we obtain a p-value of , indicating that platelet recovery is significantly associated with increased survival. We can also perform a likelihood ratio test... ### Fit Cox without time-varying covariate for LRT > fit.red <- coxph( Surv( start, stop, inodis) ~ fab + agep.c*aged.c + factor(g) + imtx, data=bmt.tvc ) > anova( fit.red, fit ) Analysis of Deviance Table Cox : response is Surv(start, stop, inodis) Model 1: ~ fab + agep.c * aged.c + factor(g) + imtx Model 2: ~ fab + agep.c * aged.c + factor(g) + imtx + iplate loglik Chisq Df P(> Chi ) * 9.18
19 Interpretation: Adjusting for FAB, group, ages of patient and donor and MTX assignment, a randomly sampled subject having had a platelet recovery by any given time t has about 2/5 the risk of relapse / death than a (similar) subject at the same time t who has not yet experienced platelet recovery. Conclusion: Platelet recovery appears to be an important indicator for successful recovery, or at least, for delayed relapse. 9.19
20 Recall we had a with several covariates, and we were considering the effect of MTX (0 or 1), indicating receipt of a GVH prophylactic We had concluded via a graphical investigation: The hazards for the two MTX groups do not appear to be proportional Can we / test this conclusion? 9.20
21 Our original was λ i (t x i ) = λ 0 (t) exp(η i ) where η i = β T x i x ij = ages of patient and donor, fab, group, and x i7 = imtx i β 7 is the estimated log hazard ratio at all times t for imtx= 1 versus imtx= 0, holding other x s constant 9.21
22 Consider the alternative λ i (t imtx i = 1)/λ i (t imtx i = 0) = t β 8 exp(β 7 ) This expresses a departure from the proportional hazards assumption: β 8 = 0 t β 8 = 1 t so λ i(t imtx = 1) λ i (t imtx = 0) = exp(β 7) Testing H 0 : β 8 = 0 will provide a test of the proportional hazards assumption 9.22
23 The is now and x i7 = imtx i x i8 = imtx i log(t) ie. a covariate by (log) time interaction... Q: How do we test this in R? A: Create a data set with one observation for each person at each observed failure time: entry: just after last failure time t(j 1) exit: observed failure time t(j) Note: Fitting such a is computationally intensive and can be difficult for large datasets (with many failures)
24 Data analysis set-up Here is one way we can obtain the necessary dataset... ## ##### Create dataset to look at an interaction with time u.evtimes <- unique( bmt$tnodis[ bmt$inodis==1 ] ) num.event <- length( u.evtimes ) bmt.texpand <- bmt[, c("id", "tnodis", "inodis", "fab", "agep.c", "aged.c", "g", "imtx") ] ## ##### Replicate each record for the number of observed failures bmt.texpand <- bmt.texpand[ rep(bmt.texpand$id,each=num.event), ] bmt.texpand$start <- rep( c(0,u.evtimes[1:(num.event-1)]), sum(!duplicated(bmt.texpand$id)) ) bmt.texpand$stop <- rep( u.evtimes, sum(!duplicated(bmt.texpand$id)) ) ## ##### Remove unnecessary rows and create event indicators ## bmt.texpand <- bmt.texpand[ bmt.texpand$tnodis > bmt.texpand$start, ] bmt.texpand <- bmt.texpand[ dim(bmt.texpand)[1]:1, ] bmt.texpand$stop <- ifelse(!duplicated(bmt.texpand$id), bmt.texpand$tnodis, bmt.texpand$stop) bmt.texpand$inodis <- ifelse(!duplicated(bmt.texpand$id), bmt.texpand$inodis,0 ) bmt.texpand <- bmt.texpand[ dim(bmt.texpand)[1]:1, ] 9.24
25 Now, our scientific question is whether the effect of imtx as measured by the hazard ratio changes with respect to time This can be parameterized and tested using a multiplicative interaction ### Test whether effect of imtx varies with time > fit <- coxph( Surv( start, stop, inodis) ~ fab + agep.c*aged.c + factor(g) + imtx + imtx:log(stop), data=bmt.texpand ) > summary(fit) coef exp(coef) se(coef) z Pr(> z ) fab ** agep.c aged.c factor(g) ** factor(g) imtx * agep.c:aged.c ** imtx:log(stop) * 9.25
26 ### Test whether effect of imtx varies with time > fit <- coxph( Surv( start, stop, inodis) ~ fab + agep.c*aged.c + factor(g) + imtx + imtx:log(stop), data=bmt.texpand ) > summary(fit) exp(coef) exp(-coef) lower.95 upper.95 fab agep.c aged.c factor(g) factor(g) imtx agep.c:aged.c imtx:log(stop)
27 The Wald test indicates a time-varying effect with imtx (z-statistic of and p-value of ) We can also look at a likelihood ratio test... ### Fit Cox without interaction for LRT Analysis of Deviance Table Cox : response is Surv(start, stop, inodis) Model 1: ~ fab + agep.c * aged.c + factor(g) + imtx Model 2: ~ fab + agep.c * aged.c + factor(g) + imtx + imtx:log(stop) loglik Chisq Df P(> Chi ) * --- Signif. codes: 0 Ô***Õ Ô**Õ 0.01 Ô*Õ 0.05 Ô.Õ 0.1 Ô Õ
28 Interpretation: When log(t) = 0, i.e. at t = 1, the estimated relative hazard for an imtx= 1 subject, compared to an imtx= 0 subject is 16.5, CI [1.7,157] At 1 year (t = 365), the estimated log-relative hazard for an imtx= 1 subject, compared to aimtx= 0 subject is log(16.48) + log(.607) log(365) = 0.143, giving a relative hazard estimate of exp( 0.143) =
29 The lincontr.coxph() makes this easy... Conclusion: The relative risk of relapse or death for MTX patients is very high early in the study, but drops off later in the study. The hazards are not proportional (X 2 = 5.98 on 1 df) ### Estimated effect of imtx at 1 year (365 days) > lincontr.coxph( fit, contr.names=c("imtx","imtx:log(stop)"), contr.coef=c(1,log(365)) ) Test of H_0: exp( 1*imtx *imtx:log(stop) ) = 1 : exp( Est ) se.est zstat pval ci95.lo ci95.hi
30 Modeling with time-varying covariates With time-varying covariates, it is much more difficult to interpret the in terms of survival time baseline survival function may not have a real meaning you are always safe to maintain interpretation at the hazard level Keep your eyes on the risk sets! Time-varying covariates methods can be used to test the proportional hazards assumptions 9.30
31 Modeling with time-varying covariates Model-building strategy: Follow similar strategy to that for time-fixed covariates Treat true time-varying covariates on equal footing with other covariates Treat interactions of time-fixed covariates with log(t) (used to test proportional hazards assumption) on equal footing with other interactions: Examine their effects later in building strategy Go for parsimony, if possible 9.31
32 Modeling with time-varying covariates Need to be careful about the functional form of t in the : Good ideal to categorize time. For example, consider the hazard ratio associated with your predictor of interest over 0-6 months, 6-12 months, etc. (See Homework 4) Can also consider Schoenfeld resdiduals...coming soon 9.32
Lecture 8 Stat D. Gillen
Statistics 255 - Survival Analysis Presented February 23, 2016 Dan Gillen Department of Statistics University of California, Irvine 8.1 Example of two ways to stratify Suppose a confounder C has 3 levels
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 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 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 11. Interval Censored and. Discrete-Time Data. Statistics Survival Analysis. Presented March 3, 2016
Statistics 255 - Survival Analysis Presented March 3, 2016 Motivating Dan Gillen Department of Statistics University of California, Irvine 11.1 First question: Are the data truly discrete? : Number of
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 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 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 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 informationChapter 7 Fall Chapter 7 Hypothesis testing Hypotheses of interest: (A) 1-sample
Bios 323: Applied Survival Analysis Qingxia (Cindy) Chen Chapter 7 Fall 2012 Chapter 7 Hypothesis testing Hypotheses of interest: (A) 1-sample H 0 : S(t) = S 0 (t), where S 0 ( ) is known survival function,
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 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 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 informationThe coxvc_1-1-1 package
Appendix A The coxvc_1-1-1 package A.1 Introduction The coxvc_1-1-1 package is a set of functions for survival analysis that run under R2.1.1 [81]. This package contains a set of routines to fit Cox models
More informationTime-dependent coefficients
Time-dependent coefficients Patrick Breheny December 1 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/20 Introduction As we discussed previously, stratification allows one to handle variables that
More informationCox s proportional hazards/regression model - model assessment
Cox s proportional hazards/regression model - model assessment Rasmus Waagepetersen September 27, 2017 Topics: Plots based on estimated cumulative hazards Cox-Snell residuals: overall check of fit Martingale
More informationLecture 10. Diagnostics. Statistics Survival Analysis. Presented March 1, 2016
Statistics 255 - Survival Analysis Presented March 1, 2016 Dan Gillen Department of Statistics University of California, Irvine 10.1 Are model assumptions correct? Is the proportional hazards assumption
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 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 informationPackage threg. August 10, 2015
Package threg August 10, 2015 Title Threshold Regression Version 1.0.3 Date 2015-08-10 Author Tao Xiao Maintainer Tao Xiao Depends R (>= 2.10), survival, Formula Fit a threshold regression
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 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 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 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 informationRelative-risk regression and model diagnostics. 16 November, 2015
Relative-risk regression and model diagnostics 16 November, 2015 Relative risk regression More general multiplicative intensity model: Intensity for individual i at time t is i(t) =Y i (t)r(x i, ; t) 0
More informationFitting Cox Regression Models
Department of Psychology and Human Development Vanderbilt University GCM, 2010 1 Introduction 2 3 4 Introduction The Partial Likelihood Method Implications and Consequences of the Cox Approach 5 Introduction
More informationTime-dependent covariates
Time-dependent covariates Rasmus Waagepetersen November 5, 2018 1 / 10 Time-dependent covariates Our excursion into the realm of counting process and martingales showed that it poses no problems to introduce
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 informationChapter 4 Regression Models
23.August 2010 Chapter 4 Regression Models The target variable T denotes failure time We let x = (x (1),..., x (m) ) represent a vector of available covariates. Also called regression variables, regressors,
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 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 informationStat 642, Lecture notes for 04/12/05 96
Stat 642, Lecture notes for 04/12/05 96 Hosmer-Lemeshow Statistic The Hosmer-Lemeshow Statistic is another measure of lack of fit. Hosmer and Lemeshow recommend partitioning the observations into 10 equal
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 informationEstimating Causal Effects of Organ Transplantation Treatment Regimes
Estimating Causal Effects of Organ Transplantation Treatment Regimes David M. Vock, Jeffrey A. Verdoliva Boatman Division of Biostatistics University of Minnesota July 31, 2018 1 / 27 Hot off the Press
More informationFrailty Modeling for clustered survival data: a simulation study
Frailty Modeling for clustered survival data: a simulation study IAA Oslo 2015 Souad ROMDHANE LaREMFiQ - IHEC University of Sousse (Tunisia) souad_romdhane@yahoo.fr Lotfi BELKACEM LaREMFiQ - IHEC University
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 informationConsider Table 1 (Note connection to start-stop process).
Discrete-Time Data and Models Discretized duration data are still duration data! Consider Table 1 (Note connection to start-stop process). Table 1: Example of Discrete-Time Event History Data Case Event
More information[Part 2] Model Development for the Prediction of Survival Times using Longitudinal Measurements
[Part 2] Model Development for the Prediction of Survival Times using Longitudinal Measurements Aasthaa Bansal PhD Pharmaceutical Outcomes Research & Policy Program University of Washington 69 Biomarkers
More informationMultivariable Fractional Polynomials
Multivariable Fractional Polynomials Axel Benner September 7, 2015 Contents 1 Introduction 1 2 Inventory of functions 1 3 Usage in R 2 3.1 Model selection........................................ 3 4 Example
More information5. Parametric Regression Model
5. Parametric Regression Model The Accelerated Failure Time (AFT) Model Denote by S (t) and S 2 (t) the survival functions of two populations. The AFT model says that there is a constant c > 0 such that
More informationSurvival Analysis. 732G34 Statistisk analys av komplexa data. Krzysztof Bartoszek
Survival Analysis 732G34 Statistisk analys av komplexa data Krzysztof Bartoszek (krzysztof.bartoszek@liu.se) 10, 11 I 2018 Department of Computer and Information Science Linköping University Survival analysis
More informationSTAT 526 Spring Midterm 1. Wednesday February 2, 2011
STAT 526 Spring 2011 Midterm 1 Wednesday February 2, 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
More informationLecture 14: Introduction to Poisson Regression
Lecture 14: Introduction to Poisson Regression Ani Manichaikul amanicha@jhsph.edu 8 May 2007 1 / 52 Overview Modelling counts Contingency tables Poisson regression models 2 / 52 Modelling counts I Why
More informationModelling counts. Lecture 14: Introduction to Poisson Regression. Overview
Modelling counts I Lecture 14: Introduction to Poisson Regression Ani Manichaikul amanicha@jhsph.edu Why count data? Number of traffic accidents per day Mortality counts in a given neighborhood, per week
More informationCURE FRACTION MODELS USING MIXTURE AND NON-MIXTURE MODELS. 1. Introduction
ØÑ Å Ø Ñ Ø Ð ÈÙ Ð Ø ÓÒ DOI: 10.2478/v10127-012-0001-4 Tatra Mt. Math. Publ. 51 (2012), 1 9 CURE FRACTION MODELS USING MIXTURE AND NON-MIXTURE MODELS Jorge A. Achcar Emílio A. Coelho-Barros Josmar Mazucheli
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 informationStatistical Methods III Statistics 212. Problem Set 2 - Answer Key
Statistical Methods III Statistics 212 Problem Set 2 - Answer Key 1. (Analysis to be turned in and discussed on Tuesday, April 24th) The data for this problem are taken from long-term followup of 1423
More information11 November 2011 Department of Biostatistics, University of Copengen. 9:15 10:00 Recap of case-control studies. Frequency-matched studies.
Matched and nested case-control studies Bendix Carstensen Steno Diabetes Center, Gentofte, Denmark http://staff.pubhealth.ku.dk/~bxc/ Department of Biostatistics, University of Copengen 11 November 2011
More informationPower and Sample Size Calculations with the Additive Hazards Model
Journal of Data Science 10(2012), 143-155 Power and Sample Size Calculations with the Additive Hazards Model Ling Chen, Chengjie Xiong, J. Philip Miller and Feng Gao Washington University School of Medicine
More informationMultivariable Fractional Polynomials
Multivariable Fractional Polynomials Axel Benner May 17, 2007 Contents 1 Introduction 1 2 Inventory of functions 1 3 Usage in R 2 3.1 Model selection........................................ 3 4 Example
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 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 informationCase-control studies
Matched and nested case-control studies Bendix Carstensen Steno Diabetes Center, Gentofte, Denmark b@bxc.dk http://bendixcarstensen.com Department of Biostatistics, University of Copenhagen, 8 November
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 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 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 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 informationStatistics in medicine
Statistics in medicine Lecture 4: and multivariable regression Fatma Shebl, MD, MS, MPH, PhD Assistant Professor Chronic Disease Epidemiology Department Yale School of Public Health Fatma.shebl@yale.edu
More informationThe nltm Package. July 24, 2006
The nltm Package July 24, 2006 Version 1.2 Date 2006-07-17 Title Non-linear Transformation Models Author Gilda Garibotti, Alexander Tsodikov Maintainer Gilda Garibotti Depends
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 informationLongitudinal + Reliability = Joint Modeling
Longitudinal + Reliability = Joint Modeling Carles Serrat Institute of Statistics and Mathematics Applied to Building CYTED-HAROSA International Workshop November 21-22, 2013 Barcelona Mainly from Rizopoulos,
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 informationMultivariate Survival Analysis
Multivariate Survival Analysis Previously we have assumed that either (X i, δ i ) or (X i, δ i, Z i ), i = 1,..., n, are i.i.d.. This may not always be the case. Multivariate survival data can arise in
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 informationFaculty of Health Sciences. Regression models. Counts, Poisson regression, Lene Theil Skovgaard. Dept. of Biostatistics
Faculty of Health Sciences Regression models Counts, Poisson regression, 27-5-2013 Lene Theil Skovgaard Dept. of Biostatistics 1 / 36 Count outcome PKA & LTS, Sect. 7.2 Poisson regression The Binomial
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 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 information7/28/15. Review Homework. Overview. Lecture 6: Logistic Regression Analysis
Lecture 6: Logistic Regression Analysis Christopher S. Hollenbeak, PhD Jane R. Schubart, PhD The Outcomes Research Toolbox Review Homework 2 Overview Logistic regression model conceptually Logistic regression
More informationLogistic Regression. Interpretation of linear regression. Other types of outcomes. 0-1 response variable: Wound infection. Usual linear regression
Logistic Regression Usual linear regression (repetition) y i = b 0 + b 1 x 1i + b 2 x 2i + e i, e i N(0,σ 2 ) or: y i N(b 0 + b 1 x 1i + b 2 x 2i,σ 2 ) Example (DGA, p. 336): E(PEmax) = 47.355 + 1.024
More informationMultinomial Logistic Regression Models
Stat 544, Lecture 19 1 Multinomial Logistic Regression Models Polytomous responses. Logistic regression can be extended to handle responses that are polytomous, i.e. taking r>2 categories. (Note: The word
More informationSection IX. Introduction to Logistic Regression for binary outcomes. Poisson regression
Section IX Introduction to Logistic Regression for binary outcomes Poisson regression 0 Sec 9 - Logistic regression In linear regression, we studied models where Y is a continuous variable. What about
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 informationChapter 7: Hypothesis testing
Chapter 7: Hypothesis testing Hypothesis testing is typically done based on the cumulative hazard function. Here we ll use the Nelson-Aalen estimate of the cumulative hazard. The survival function is used
More informationExtensions of Cox Model for Non-Proportional Hazards Purpose
PhUSE 2013 Paper SP07 Extensions of Cox Model for Non-Proportional Hazards Purpose Jadwiga Borucka, PAREXEL, Warsaw, Poland ABSTRACT Cox proportional hazard model is one of the most common methods used
More informationProportional hazards regression
Proportional hazards regression Patrick Breheny October 8 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/28 Introduction The model Solving for the MLE Inference Today we will begin discussing regression
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 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 informationIndications and outcomes after UD HSCT
Indications and outcomes after UD HSCT Jakob R Passweg 1 Impact on Outcome: Patient Age, Disease Severity I II III Title: Optimization of Therapy for Severe AplasticAnemia Based on Clinical, Biological
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 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 informationLongitudinal Modeling with Logistic Regression
Newsom 1 Longitudinal Modeling with Logistic Regression Longitudinal designs involve repeated measurements of the same individuals over time There are two general classes of analyses that correspond to
More informationLecture 12: Effect modification, and confounding in logistic regression
Lecture 12: Effect modification, and confounding in logistic regression Ani Manichaikul amanicha@jhsph.edu 4 May 2007 Today Categorical predictor create dummy variables just like for linear regression
More informationA fast routine for fitting Cox models with time varying effects
Chapter 3 A fast routine for fitting Cox models with time varying effects Abstract The S-plus and R statistical packages have implemented a counting process setup to estimate Cox models with time varying
More informationST 732, Midterm Solutions Spring 2019
ST 732, Midterm Solutions Spring 2019 Please sign the following pledge certifying that the work on this test is your own: I have neither given nor received aid on this test. Signature: Printed Name: There
More informationST745: Survival Analysis: Cox-PH!
ST745: Survival Analysis: Cox-PH! Eric B. Laber Department of Statistics, North Carolina State University April 20, 2015 Rien n est plus dangereux qu une idee, quand on n a qu une idee. (Nothing is more
More informationLecture 1. Introduction Statistics Statistical Methods II. Presented January 8, 2018
Introduction Statistics 211 - Statistical Methods II Presented January 8, 2018 linear models Dan Gillen Department of Statistics University of California, Irvine 1.1 Logistics and Contact Information Lectures:
More informationSimple techniques for comparing survival functions with interval-censored data
Simple techniques for comparing survival functions with interval-censored data Jinheum Kim, joint with Chung Mo Nam jinhkim@suwon.ac.kr Department of Applied Statistics University of Suwon Comparing survival
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 informationSSUI: Presentation Hints 2 My Perspective Software Examples Reliability Areas that need work
SSUI: Presentation Hints 1 Comparing Marginal and Random Eects (Frailty) Models Terry M. Therneau Mayo Clinic April 1998 SSUI: Presentation Hints 2 My Perspective Software Examples Reliability Areas that
More informationExtensions of Cox Model for Non-Proportional Hazards Purpose
PhUSE Annual Conference 2013 Paper SP07 Extensions of Cox Model for Non-Proportional Hazards Purpose Author: Jadwiga Borucka PAREXEL, Warsaw, Poland Brussels 13 th - 16 th October 2013 Presentation Plan
More informationSurvival models and health sequences
Survival models and health sequences Walter Dempsey University of Michigan July 27, 2015 Survival Data Problem Description Survival data is commonplace in medical studies, consisting of failure time information
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 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 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 informationSimple logistic regression
Simple logistic regression Biometry 755 Spring 2009 Simple logistic regression p. 1/47 Model assumptions 1. The observed data are independent realizations of a binary response variable Y that follows a
More informationECLT 5810 Linear Regression and Logistic Regression for Classification. Prof. Wai Lam
ECLT 5810 Linear Regression and Logistic Regression for Classification Prof. Wai Lam Linear Regression Models Least Squares Input vectors is an attribute / feature / predictor (independent variable) The
More informationSingle-level Models for Binary Responses
Single-level Models for Binary Responses Distribution of Binary Data y i response for individual i (i = 1,..., n), coded 0 or 1 Denote by r the number in the sample with y = 1 Mean and variance E(y) =
More informationfor Time-to-event Data Mei-Ling Ting Lee University of Maryland, College Park
Threshold Regression for Time-to-event Data Mei-Ling Ting Lee University of Maryland, College Park MLTLEE@UMD.EDU Outline The proportional hazards (PH) model is widely used in analyzing time-to-event data.
More informationDescription Syntax for predict Menu for predict Options for predict Remarks and examples Methods and formulas References Also see
Title stata.com stcrreg postestimation Postestimation tools for stcrreg Description Syntax for predict Menu for predict Options for predict Remarks and examples Methods and formulas References Also see
More information