for Time-to-event Data Mei-Ling Ting Lee University of Maryland, College Park
|
|
- Dorothy Craig
- 5 years ago
- Views:
Transcription
1 Threshold Regression for Time-to-event Data Mei-Ling Ting Lee University of Maryland, College Park
2 Outline The proportional hazards (PH) model is widely used in analyzing time-to-event data. It requires, however, a strong PH assumption. Threshold regression (TR) is an alternative model without the PH assumption. Brief introduction of the TR model. Extensions and Examples Connections between TR and PH models
3 Hazard Function: The instantaneous failure rate Pr( t T < t +Δ t T t ) ht () = lim Δ t 0 Δt t ht () = f () t St ()
4 The Proportional Hazards (PH) Model For any covariates x 1 and x 2, the PH model assumes that the hazard functions are proportional to a baseline hazard h 0 (t). h t x h t x ( 1) = 0( )exp( 1β ) h ( t x 2 ) = h 0 ( t ) exp( x 2 β ) h ( t x 1) e x p ( x 1β ) = = h ( t x 2 ) exp( x 2 β ) exp[( x 1 x 2) ] β The hazard ratio of the PH model has a simple e a a d a o o e ode as a s p e form that does not depend on time t.
5 Hazard rate is really an elusive concept, especially when one tries to interpret its shape considered as a function of time. Aalen, Gjessing (2001), Understanding the shape of the hazard rate: a process point of view, Statistical Science, 16, S i l d Aalen, Borgan, and Gjessing (2008). Survival and Event History Analysis: a Process Point of view. Springer.
6 Consequences of the PH Assumption As a consequence of the proportional hazards assumption, the survivorship functions of the two groups must be parallel.
7 Example 1 Time to infection of kidney dialysis i patients t with different catheterization procedures (Nahman et al 1992, Klein & Moesberger 2003) Surgical group: 43 patients utilized a surgically placed catheter Percutaneous group: 76 patients utilized a percutaneous placement of their catheter Th t ti i d fi d b th ti t The event time is defined by the time to cutaneous exit-site infection.
8 Example 1 Kaplan-Meier Estimate t versus PH Cox Model
9 Example 1 Generalized Gamma Threshold Regression 9
10
11 Threshold Regression Model can estimate hazard ratios at different time points timevalue=2 timevalue=20 xi: threg,lny0(i.group) mu(i.group) failure(infection) /* */ time(time) hr(group) timevalue(2) graph(hz) xi: threg,lny0(i.group) mu(i.group) failure(infection) /* */ time(time) hr(group) timevalue(20) graph(hz) Hazard Rate Hazard Rate 0 0 h(t) h(t) Calendar Time Calendar Time reference _Igroup_2==1 reference _Igroup_2==1 Proportional Hazard Cox Model can only estimate t a constant t hazard ratio stset time, failure(infection) xi: stcox i.group
12 Example 2 Time to progression of disease for patients with limited Stage II or IIIA ovarian carcinoma (Fleming et al., 1980; Hess, 1995) Low-grade tumor group: 15 patients with low-grade or well-differentiated tumors High-grade tumor group: 20 patients with high-grade grade or undifferentiated tumors 12
13 Example 2 Kaplan-Meier Estimates 13
14 Example 2 Generalized Gamma Threshold Regression Threshold 14
15 Example 3 Time to death for patients with bile duct cancer (Fleming et al., 1980; Lawless, 2003) Control group: 25 patients without treatment Treatment group: 22 patients with a combined treatment of radiation and the drug 5-fluorouracil (5-FU) 15
16 Example 3 Kaplan-Meier Estimates 16
17 Example 3 Generalized Gamma Threshold Regression 17
18 Threshold Regression is based on the First-hitting-time (FHT) Model Modeling event times by a stochastic process reaching a boundary
19 Threshold Regression is based on First-hitting-time time Methodology : Example: Engineering Equipment fails when its cumulative wear first reaches a failure threshold. Example: Health and Medicine A patient dies of heart disease when condition of the heart deteriorates to a critical state. Example: Social Sciences Divorce occurs when tensions in a marriage reach a breaking point. A review article: Modeling Event Times by Stochastic Process Reaching a Boundary (Lee & Whitmore 2006, Statistical Sciences)
20 y 0 Process Sample path Y(t) 0 0 time t S First hitting time S of a fixed boundary at level zero for a stochastic process of interest Y(t)
21 Process Y(t) y 0 Sample paths Y(t) 0 Time t S L Two sample paths of a stochastic process of interest: (1) One path experiences failure at first hitting time S (2) One path is surviving at end of follow up at time L
22 First Hitting Time (FHT) Models {Y(t)} : the stochastic process of interest Β : the threshold or boundary set First hitting time S defined by S = inf { t : Y(t) B }
23 Examples of First-hitting-time (FHT) Models: Wiener diffusion to a fixed boundary Progress of multiple l myeloma until death Renewal process to a fixed count of renewal events Time to the n th epileptic seizure Semi-Markov process to an absorbing state Multi-state model for disease with death as an absorbing state
24 First Hitting Time (FHT) Models The first hitting time (FHT) model describes many time-to-event applications The stochastic process of interest {Y(t)} may represent the latent (unobservable) health status of a subject. The threshold constitutes the critical level of the process that triggers the failure event (e.g., symptomatic cancer, death). The event occurs when health status {Y(t)} first decreases to the zero threshold.
25 Parameters for the FHT Model Model parameters for a latent process Y(t) : Process parameters: θ = (μ, σ 2 ), where μ is the mean drift and σ 2 is the variance Baseline level of process: Y(0) = y 0 Because Y(t) () is latent, we set σ 2 = 1.
26 Threshold Regression: FHT Model with Regression Structure Including covariates to FHT model parameters θ = (μ, σ 2 ) using link functions, where σ 2 = 1. Mean drift parameters: μ = g 1 ( x 1, x 2, ) Baseline level: y(0) = g 2 (x 1, x 2, ) A variety of link functions g have been considered: Linear combinations of covariates X1,, Xp Polynomial combinations of X1,, Xp Semi-parametric regression splines Penalized regression splines Random effects and Bayesian models 26
27 Likelihood Inference for the FHT Model The likelihood contribution of each sample subject is as follows. If the subject fails at S=s: f (s y 0, μ) = Pr [ first-hitting-time in (s, s+ds) ] If the subject survives beyond time L: 1- F (L y 0,μ) = Pr [ no first-hitting-time before L ]
28 n { } i i i i ln L, x = d ln f t, x + 1 d ln F t, x. where ( θ ) ( θ ) ( ) ( θ ) d i i= 1 is the failure indicator for subject t is a censored survival time t = s if subject i fails i i i i ( ) f and F denote the FHT p.d.f and complementary c.d.f.
29 Extensions of Threshold Regression (TR) 1. Process Y(t): Wiener process, gamma process, etc 2. Boundary: straight lines or curves 3. Time scale: calendar time, running or analytical time
30 Example 4 Comparing TR and PH for two possible scenarios Scenario 1: A favorable surgical intervention might double the health level of a patient (relative to control) at the outset of the study but leaves the rate of decline of health unchanged on both study arms. Scenario 2: Drug treatment for chronic disease may slow the rate of decline in health (relative to control) but leaves the initial health level unchanged on both study arms. 30
31 Scenario 1 Scenario Process Level y 1-1 Treatment t Process Level 1 y Treatment t -1 Control Control Time Time
32 Example 4 Consider a simulated randomized clinical trial with two study arms and n=200 subjects on each arm. (control: arm=0; treatment: arm=1). The TR model is a Wiener diffusion model. Regression links for the TR model ln (y 0 )=α α 0 + α 1 arm μ = β 0 + β 1 arm In Scenario 1: α 1 =ln(2) = , β 1 =0 In Scenario 2: α 1 = ln (1) = 0, β 1 =
33 Example 4 (simulation) The TR model: ln (y 0 ) = α 0 + α 1 arm μ = β 0 + β 1 arm The mean survival time on any arm is given by E(S)= y 0 / μ Both scenarios have the same mean of 1 year survival on the control arm and 2 years on the treatment arm. 33
34 A Comparison of TR and PH Regression Results
35 Example 4: Modeling Results for Scenario 1 In scenario 1, the PH assumption is untenable. The analyst would then need to proceed to some other model. If the TR model is considered, the table shows a significant positive treatment effect for ln y 0 but an insignificant effect for μ (as expected). The significant coefficient for α 1 translates into an estimated t exp(0.7752)= multiplier li for initial health level y 0 (note that the true value is 2, i.e. the treatment arm is twice as large at y 0 )
36 Example 4: Modeling Results for Scenario 2 In scenario 2, the PH assumption is tenable. The analyst using PH regression would discover a significant hazard ratio of , indicating that the treatment is about 56% of the control risk. However, the question of how the treatment is acting on the subject is not known from the PH model. If the TR model is considered, as expected, the table shows a significant ifi positive treatment t t effect for μ but an insignificant effect for ln y 0 (). The significant ifi coefficient i for β 1 translates t into an estimated multiplier in the value of μ (the treatment is about half of the control). 36
37 2 2 Scenario Control hazard Treatment hazard Log-hazard ratio Time (years) Time (years) Scenario Control hazard Log-hazard ratio Treatment hazard Time (years) Time (years)
38 Connection of the TR model with the PH model Most survival distributions are hitting time distributions for stochastic processes (the basic TR model) Families of PH functions can be generated by varying time scales or boundaries of a TR model Th f il f PH f ti b d d The same family of PH functions can be produced by different TR mode
39 Connection of the TR model with PH model Families of PH functions can be generated by 1. Varying process Y and time scales of a TR model 2. Varying process Y and boundaries of a TR model. Hence, the same family of PH functions can be produced by different TR models. For details, see M-LT Lee, Whitmore GA. Proportional hazards and threshold regression: their theoretical and practical connections. Lifetime Data Analysis (2010). 39
40 Benefits of the TR Model Probing the causal forces behind a hazard function is always a worthwhile endeavor. A TR model represents more fundamental knowledge about the underlying science than its corresponding PH model. These insights would deepen the investigator s understanding of the case application.
41 Whitmore, G. A. (1986). First-passage-time models for duration data: Regression structures and competing risks. The Statistician Aalen O.O. and Gjessing H.K. (2001). Understanding the shape of the hazard rate: a process point of view. Statistical Science, 16: Lee, M.-L. T. and G. A. Whitmore (2006). Threshold regression for survival analysis: modeling event times by a stochastic process reaching a boundary. Statistical Sciences. Aalen O.O., Borgon O, and Gjessing H.K (2008). Survival and Event History Analysis: A process Point of View. Springer. 41
42 Yu Z, Tu W, Lee M-LT (2009). A Semiparametric Threshold Regression Analysis of Sexually Transmitted Infections in Adolescent Women, Statistics in Medicine. Pennell M, GA Whitmore, Lee M-LT (2009). Bayesian random effects threshold regression with application to survival data with nonproportional hazards, Biostatistics. Lee M-LT, Whitmore GA (2010). Proportional hazards and threshold regression: their theoretical and practical connections. Lifetime Data Analysis. Lee M-LT Whitmore GA Rosner B (2010) Threshold Lee M LT, Whitmore GA, Rosner B. (2010) Threshold regression for survival data with time-varying covariates. Statistics in Medicine.
43 Nested boundaries Z=0 Z=1 Z=2 Process Level y Representative sample path for z=2 Time
44 .4 Z=0 Z=1 Z=2.2.3 Hazard h(t) Time x t
45 Boundary b(t z) y(t z) Process Level y 0 0 s Time t
46 y(t z) Process Level y Overshoot y(s z)-b(s z) Boundary b(s z) Boundary b(t z) Time s t
47 Aalen O.O. and Gjessing H.K. (2001). Understanding the shape of the hazard rate: a process point of view. Statistical Science, 16: Lee, M.-L. T. and G. A. Whitmore (2006). Threshold regression for survival analysis: modeling event times by a stochastic process reaching a boundary. Statistical Sciences. Aalen O.O., Borgon O, and Gjessing H.K (2008). Survival and Event History Analysis: A process Point of View. Springer. Yu Z, Tu W, Lee M-LT (2009). A Semiparametric Threshold h Regression Analysis of Sexually Transmitted Infections in Adolescent Women, Statistics in Medicine (in press). Pennell M, GA Whitmore, Lee M-LT (2009). Bayesian random effects threshold regression with application to survival data with nonproportional hazards, Biostatistics, (In Press).
48
49 The Connection of Threshold Regression Models With Cox Proportional Hazard Models Mei-Ling Ting Lee, University of Maryland Alex Whitmore, McGill University Lifetime Data Analysis 2010
50 Relation of Cox Proportional Hazards and Threshold Regression The class of survival functions that are firsthitting-time time (FHT) distributions includes most distributions of practical interest Almost any family of proportional hazard (PH) functions can be generated from an FHT context by an appropriate choice of process time or boundary for the stochastic process defining the FHT.
51 Cox regression can be considered d as a special case of threshold regression (TR) for the class of ffhtdi distributions tib ti that thave the proportional hazard property Suitable methods allow TR to replicate Cox PH regression with either fixed or time-varying i covariates
52 Generating PH Functions in an FHT Context Three building blocks for threshold regression 1. A stochastic process that describes the evolution of the subject's underlying health state 2. A critical level, condition or boundary (the threshold) that triggers the event of interest when it is reached by the process for the first time (the FHT) 3. A time scale on which the process unfolds.
53 Generating PH Functions in an FHT Context The three building blocks of TR give many combinations of conditions that can generate the same family of Cox proportional hazard functions: We now show two methods of construction of a PH family from an FHT context.
54 1. Generating PH by Varying Process Time Let denote a process time function, defined as a function of calendar time t, conditional on covariate z. Let the survival function for FHT context be denoted by: We have the following two correspondences: Varying process time:
55 Generating PH by Varying Process Time Equating the two right-hand hand sides and solving for gives process time functions that generate a family of PH functions.
56 Examples: Poisson process: A Poisson process with unit intensity has the following exponential survival time distribution: Solving for the process time function gives the following. Observe that the process time is proportional to the cumulative baseline hazard function H 0 (t).
57 Examples: Brownian motion process: A Brownian motion process with a unit variance, starting at x 0, has the following survival time distribution: where Φ denotes the standard normal c.d.f. Solving for the process time function gives:
58 2. Generating PH by Varying Boundaries A family of PH functions can be generated by varying the boundary function. Example: Brownian motion process The following identity links the boundary b(t z) to the c.d.f. of survival time F(t z) for a given covariate vector z in standard Brownian motion starting at zero:
59 Generating PH by Varying Boundaries Replacing F(t z) ( by the c.d.f. for each member of a PH family and solving for boundary function b(t z) ( )gives a set of nested boundaries that generate FHTs with proportional hazard functions.
60 Example: The next figure shows a nested set of three boundaries for standard Brownian motion that correspond to three levels of a fixed covariate z. Covariate z might represent, for example, level of occupational exposure to a carcinogen The figure shows an illustrative sample path for a case where exposure level z=2 generates a survival time S.
61 10 0 b0 b1 b2 Nested boundaries z=0 z=1 z=2-10 Simulated path x(t) for z= x S Three boundaries in Brownian motion that correspond to three levels of covariate z
62 Example: The next figure shows the hazard functions for survival time corresponding to the preceding three boundaries for Brownian motion. The boundaries and process have been chosen p so the hazard functions are, in fact, proportional for the three levels of fixed covariate z.
63 haz0 haz1 haz2 h(t) ().3.2 z=1 z=0.1 z= x Hazard functions generated by the FHT model. The functions are proportional p for the three levels of covariate z.
64 The threg package in R
65 Install threg package in R Open R in Windows From menu Packages > > Install package(s) from local zip files..., select the file "threg_1.0.1.zip 1 zip" and click "Open. Type the following command in R prompt to load this package: > library(threg) Now you can use the threg package in R.
66 An example The following example can be found in the help document of the threg package, by typing the following command in R prompt: >?threg
67 The myeloma dataset *The myeloma dataset is a simulated dataset with 49 observations. time fail age gender treat Variable Name survival time Description censoring indicator variable 0 indicates a right censored time and 1 indicates an observed event time. age of the corresponding subject. It is a continuous variable. gender of the corresponding subject. It is a categorical variable. treatment type of the corresponding subject. It is a categorical variable with 4 levels.
68 Load the myeloma dataset load the data "myeloma in the package: > data("myeloma") Transform the "gender" and "treat" variables into factor variables "f.gender gender and "ftreat" f.treat. > f.gender< factor(myeloma$gender) > f.treat< factor(myeloma$treat) f bind the two new factor variables "f.gender" and "f.treat" to the dataset "myeloma". > myeloma< cbind(myeloma,f.gender,f.treat) fgenderftreat)
69 Fit a threshold regression model by using threg function Fit a threshold h regression model, dlby using age and f.gender to predict parameter lny0 in the threshold regression model, and by using age, g, f.gender and f.treat to predict parameter mu in the threshold regression model: > threg(data=myeloma, t + lny0=age+f.gender, + mu=age+f.gender+f.treat, gender+f + time=time, + failure=fail + )
70 Output (regression coefficient estimates, log likelihood value and AIC value) $coef.estest Coef. Est. Std. Err. z P> z lny0: (Intercept) lny0: age lny0: f.gender(1) mu: (Intercept) mu: age mu: f.gender(1) mu: f.treat(1) mu: f.treat(2) mu: f.treat(4) $log.likelihood [1] $AIC [1]
71 Hazard ratio calculation Fit the same threshold h regression model dlas previous, and calculate the hazard ratios for the factor variable "f.treat", at a fixed time point 3, for a specified scenario: age=50, gender=0: > threg(data=myeloma, + lny0=age+f.gender, + mu=age+f.gender+f.treat, + time=time, + failure=fail, + hr=f.treat, + timevalue=3, + scenario=age(50)+f.gender(0) + )
72 Output (hazard ratios) (some output tomitted) $scenario [1] "f.gender = 0 age=50 at time = 3" $hazard.ratio Haz. Ratio f.treat(1) f.treat(2) f.treat(4)
73 Plot of hazard functions Generate the overlaid lidestimated t hazard dfunctions corresponding to different levels of "treat", for a specified scenario: age=50,,gender=0: > threg(data=myeloma, + lny0=age+f.gender, + mu=age+f.gender+f.treat, + time=time, + failure=fail, + hr=f.treat, + timevalue=3, + scenario=age(50)+f.gender(0), + graph="hz" + )
74 Output (estimated hazard functions )
75 Plot of survival functions Generate the overlaid lidestimated t survival lfunctions corresponding to different levels of "treat", for a specified scenario: age=50,,gender=0: > threg(data=myeloma, + lny0=age+f.gender, + mu=age+f.gender+f.treat, + time=time, + failure=fail, + hr=f.treat, + timevalue=3, + scenario=age(50)+f.gender(0), + graph= sv" + )
76 Output (estimated survival functions )
77 Plot of probability density functions Generate the overlaid estimated t probability bilit density functions corresponding to different levels of "treat", for a specified scenario: age=50,,gender=0: > threg(data=myeloma, + lny0=age+f.gender, + mu=age+f.gender+f.treat, + time=time, + failure=fail, + hr=f.treat, + timevalue=3, + scenario=age(50)+f.gender(0), + graph= ds" + )
78 Output (estimated probability density functions )
Package 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 informationBayesian Analysis for Markers and Degradation
Bayesian Analysis for Markers and Degradation Mei-Ling Ting Lee 12, Maria Shubina 1, and Alan Zaslavsky 13 1 Biostatistics Department, Harvard School of Public Health, Boston, USA 2 Channing Laboratory,
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 informationRegression Modeling of Time to Event Data Using the Ornstein-Uhlenbeck Process
Regression Modeling of Time to Event Data Using the Ornstein-Uhlenbeck Process Dissertation Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the Graduate School
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. 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 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 general mixed model approach for spatio-temporal regression data
A general mixed model approach for spatio-temporal regression data Thomas Kneib, Ludwig Fahrmeir & Stefan Lang Department of Statistics, Ludwig-Maximilians-University Munich 1. Spatio-temporal regression
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 informationIndividualized Treatment Effects with Censored Data via Nonparametric Accelerated Failure Time Models
Individualized Treatment Effects with Censored Data via Nonparametric Accelerated Failure Time Models Nicholas C. Henderson Thomas A. Louis Gary Rosner Ravi Varadhan Johns Hopkins University July 31, 2018
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 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 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 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 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 information1 The problem of survival analysis
1 The problem of survival analysis Survival analysis concerns analyzing the time to the occurrence of an event. For instance, we have a dataset in which the times are 1, 5, 9, 20, and 22. Perhaps those
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 informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
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 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 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 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 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 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 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 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 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 informationA Bayesian Nonparametric Approach to Causal Inference for Semi-competing risks
A Bayesian Nonparametric Approach to Causal Inference for Semi-competing risks Y. Xu, D. Scharfstein, P. Mueller, M. Daniels Johns Hopkins, Johns Hopkins, UT-Austin, UF JSM 2018, Vancouver 1 What are semi-competing
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 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 informationJoint Modeling of Longitudinal Item Response Data and Survival
Joint Modeling of Longitudinal Item Response Data and Survival Jean-Paul Fox University of Twente Department of Research Methodology, Measurement and Data Analysis Faculty of Behavioural Sciences Enschede,
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 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 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 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 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 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 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 informationInstantaneous geometric rates via Generalized Linear Models
The Stata Journal (yyyy) vv, Number ii, pp. 1 13 Instantaneous geometric rates via Generalized Linear Models Andrea Discacciati Karolinska Institutet Stockholm, Sweden andrea.discacciati@ki.se Matteo Bottai
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 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 informationLOGISTIC REGRESSION Joseph M. Hilbe
LOGISTIC REGRESSION Joseph M. Hilbe Arizona State University Logistic regression is the most common method used to model binary response data. When the response is binary, it typically takes the form of
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 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 informationPASS Sample Size Software. Poisson Regression
Chapter 870 Introduction Poisson regression is used when the dependent variable is a count. Following the results of Signorini (99), this procedure calculates power and sample size for testing the hypothesis
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 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 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 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 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 informationExam C Solutions Spring 2005
Exam C Solutions Spring 005 Question # The CDF is F( x) = 4 ( + x) Observation (x) F(x) compare to: Maximum difference 0. 0.58 0, 0. 0.58 0.7 0.880 0., 0.4 0.680 0.9 0.93 0.4, 0.6 0.53. 0.949 0.6, 0.8
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 informationSurvival Analysis. Lu Tian and Richard Olshen Stanford University
1 Survival Analysis Lu Tian and Richard Olshen Stanford University 2 Survival Time/ Failure Time/Event Time We will introduce various statistical methods for analyzing survival outcomes What is the survival
More 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 informationVariable Selection in Competing Risks Using the L1-Penalized Cox Model
Virginia Commonwealth University VCU Scholars Compass Theses and Dissertations Graduate School 2008 Variable Selection in Competing Risks Using the L1-Penalized Cox Model XiangRong Kong Virginia Commonwealth
More informationPrerequisite: STATS 7 or STATS 8 or AP90 or (STATS 120A and STATS 120B and STATS 120C). AP90 with a minimum score of 3
University of California, Irvine 2017-2018 1 Statistics (STATS) Courses STATS 5. Seminar in Data Science. 1 Unit. An introduction to the field of Data Science; intended for entering freshman and transfers.
More informationAn Overview of Methods for Applying Semi-Markov Processes in Biostatistics.
An Overview of Methods for Applying Semi-Markov Processes in Biostatistics. Charles J. Mode Department of Mathematics and Computer Science Drexel University Philadelphia, PA 19104 Overview of Topics. I.
More informationPractical considerations for survival models
Including historical data in the analysis of clinical trials using the modified power prior Practical considerations for survival models David Dejardin 1 2, Joost van Rosmalen 3 and Emmanuel Lesaffre 1
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 informationCategorical and Zero Inflated Growth Models
Categorical and Zero Inflated Growth Models Alan C. Acock* Summer, 2009 *Alan C. Acock, Department of Human Development and Family Sciences, Oregon State University, Corvallis OR 97331 (alan.acock@oregonstate.edu).
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 informationBayesian Nonparametric Inference Methods for Mean Residual Life Functions
Bayesian Nonparametric Inference Methods for Mean Residual Life Functions Valerie Poynor Department of Applied Mathematics and Statistics, University of California, Santa Cruz April 28, 212 1/3 Outline
More informationAdaptive Prediction of Event Times in Clinical Trials
Adaptive Prediction of Event Times in Clinical Trials Yu Lan Southern Methodist University Advisor: Daniel F. Heitjan May 8, 2017 Yu Lan (SMU) May 8, 2017 1 / 19 Clinical Trial Prediction Event-based trials:
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 informationBIOS 312: Precision of Statistical Inference
and Power/Sample Size and Standard Errors BIOS 312: of Statistical Inference Chris Slaughter Department of Biostatistics, Vanderbilt University School of Medicine January 3, 2013 Outline Overview and Power/Sample
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 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 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 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 informationPackage SimSCRPiecewise
Package SimSCRPiecewise July 27, 2016 Type Package Title 'Simulates Univariate and Semi-Competing Risks Data Given Covariates and Piecewise Exponential Baseline Hazards' Version 0.1.1 Author Andrew G Chapple
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 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 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 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 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 informationInstantaneous geometric rates via generalized linear models
Instantaneous geometric rates via generalized linear models Andrea Discacciati Matteo Bottai Unit of Biostatistics Karolinska Institutet Stockholm, Sweden andrea.discacciati@ki.se 1 September 2017 Outline
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 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 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 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 informationDoes Cox analysis of a randomized survival study yield a causal treatment effect?
Published in final edited form as: Lifetime Data Analysis (2015), 21(4): 579 593 DOI: 10.1007/s10985-015-9335-y Does Cox analysis of a randomized survival study yield a causal treatment effect? Odd O.
More informationRegularization in Cox Frailty Models
Regularization in Cox Frailty Models Andreas Groll 1, Trevor Hastie 2, Gerhard Tutz 3 1 Ludwig-Maximilians-Universität Munich, Department of Mathematics, Theresienstraße 39, 80333 Munich, Germany 2 University
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 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 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 informationBIAS OF MAXIMUM-LIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY
BIAS OF MAXIMUM-LIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY Ingo Langner 1, Ralf Bender 2, Rebecca Lenz-Tönjes 1, Helmut Küchenhoff 2, Maria Blettner 2 1
More informationISQS 5349 Spring 2013 Final Exam
ISQS 5349 Spring 2013 Final Exam Name: General Instructions: Closed books, notes, no electronic devices. Points (out of 200) are in parentheses. Put written answers on separate paper; multiple choices
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 informationMultilevel Statistical Models: 3 rd edition, 2003 Contents
Multilevel Statistical Models: 3 rd edition, 2003 Contents Preface Acknowledgements Notation Two and three level models. A general classification notation and diagram Glossary Chapter 1 An introduction
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 informationParametric Maximum Likelihood Estimation of Cure Fraction Using Interval-Censored Data
Columbia International Publishing Journal of Advanced Computing (2013) 1: 43-58 doi:107726/jac20131004 Research Article Parametric Maximum Likelihood Estimation of Cure Fraction Using Interval-Censored
More informationEvaluating the value of structural heath monitoring with longitudinal performance indicators and hazard functions using Bayesian dynamic predictions
Evaluating the value of structural heath monitoring with longitudinal performance indicators and hazard functions using Bayesian dynamic predictions C. Xing, R. Caspeele, L. Taerwe Ghent University, Department
More informationPubH 7470: STATISTICS FOR TRANSLATIONAL & CLINICAL RESEARCH
PubH 7470: STATISTICS FOR TRANSLATIONAL & CLINICAL RESEARCH The First Step: SAMPLE SIZE DETERMINATION THE ULTIMATE GOAL The most important, ultimate step of any of clinical research is to do draw inferences;
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 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 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 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 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 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 information