SEMIPARAMETRIC APPROACHES TO INFERENCE IN JOINT MODELS FOR LONGITUDINAL AND TIME-TO-EVENT DATA
|
|
- Dana Joseph
- 5 years ago
- Views:
Transcription
1 SEMIPARAMETRIC APPROACHES TO INFERENCE IN JOINT MODELS FOR LONGITUDINAL AND TIME-TO-EVENT DATA Marie Davidian and Anastasios A. Tsiatis davidian/ tsiatis/ In part joint work with Xiao Song, Department of Biostatistics, University of Washington Semparametric Inference in Joint Models 1
2 OUTLINE 1. Introduction and example 2. Joint models 3. Conditional score approach 4. Semiparametric likelihood approach 5. Simulation evidence 6. Example, revisited 7. Discussion Semparametric Inference in Joint Models 2
3 1. INTRODUCTION AND EXAMPLE Longitudinal studies in medical research: Time-to-event, e.g. death, progression to AIDS, etc. (possibly censored) Intermittent longitudinal measurements (time-dependent) Time-independent covariates, e.g. treatment group, age at baseline, gender, etc Questions of interest: Interrelationships between these variables Semparametric Inference in Joint Models 3
4 AIDS Clinical Trials Group (ACTG) 175: Compare 4 antiretroviral regimens in asymptomatic HIV-infected patients 2467 subjects, Primary objective: Compare on basis of time to AIDS or death CD4 counts approximately every 12 weeks Subsequent objective 1: Characterize within-subject patterns of CD4 change Subsequent objective 2: Characterize relationship between features of CD4 profiles and survival Semparametric Inference in Joint Models 4
5 CD4 profiles: 10 randomly chosen subjects log CD week Semparametric Inference in Joint Models 5
6 Ideally: To address these objectives, would like to have for all subjects i T i = time to event X i (u) = longitudinal trajectory for all times u 0 Z i = vector of time-independent covariates Semparametric Inference in Joint Models 6
7 Objective 1: Characterize within-subject patterns of CD4 change Estimate aspects of average longitudinal trajectories; e.g. E{X i (u)}, var{x i (u)} Trajectory of the average person Difficulties/complications: Data are collected intermittently Possibly measured with error Informative censoring: no measurement available after death! [E.g., Wu and Carroll (1988), Hogan and Laird (1997)] Observe only W i (t ij ), j =1,...,m i,t ij T i W i (t ij )=X i (t ij )+e(t ij ) Semparametric Inference in Joint Models 7
8 Objective 2: Characterize relationship between features of CD4 profiles and survival Establish relationships between T i and X i (u), u 0, and Z i History X H i (t) ={X i(u),u t} Proportional hazards model: λ i (u) = lim du 0 du 1 pr{u T i <u+ du T i u, Xi H (u),z i } = λ 0 (u)exp{γx i (u)+η T Z i } Interest in estimating γ, η, (q 1). Semiparametric model (baseline hazard λ 0 (u) isunspecified) Focus here: Objective 2 Semparametric Inference in Joint Models 8
9 Popular focus: Surrogate marker problem Can we replace the time-to-event endpoint by CD4 in assessing treatment efficacy? Roughly speaking (Prentice conditions) Treatment is effective on time-to-event Treatment has effect on marker; e.g. CD4 more on treatment than placebo Effect of treatment should manifest through the marker; i.e., risk of event given specific marker trajectory should be of treatment For example, λ 0 (u)exp{γx i (u)+η T Z i }, If γ<0andη =0? X i (u) isasurrogate marker Semparametric Inference in Joint Models 9
10 Complication 1: Time-to-event may be censored Underlying censoring variable C i (T i,c i ) = (survival,censoring) time for i Observe V i =min(t i,c i ), i = I(T i C i ) Cox (1972, 1975) partial likelihood: maximize [ n exp{γx i (V i )+η T Z i } n j=1 exp{γx j(v i )+η T Z j }I(V j V i ) i=1 ] i Semparametric Inference in Joint Models 10
11 Complication 2: Presumes X i (u),i=1,...,n is available at all failure times X i (u) only available intermittently Subject to error and intra-subject variation Naive approach: Extrapolate X i (u) at failure times u using available longitudinal data, substitute in partial likelihood E.g. use most recent observed, error-prone covariate value (in place of unobserved, true covariate) Last Value Carried Forward (LVCF) Such substitution biased estimation (Prentice, 1982) Semparametric Inference in Joint Models 11
12 2. JOINT MODELS Popular approach: Joint longitudinal data-survival model True longitudinal data process: mixed effects model (+ stochastic process) Normal measurement error Proportional hazards model depending on true longitudinal data process Semparametric Inference in Joint Models 12
13 Longitudinal data process: Popular models X i (u) =α 0i + α 1i u, α i =(α 0i,α 1i ) T X i (u) =α 0i + α 1i u + α 2i u α pi u p X i (u) =α 0i + α 1i u + U i (u) U i (u) stationary Gaussian process (Henderson, Diggle, Zeger, et al.) U i (u) integrated Ornstein-Uhlenbeck (IOU) process (Taylor, DeGruttola, et al.) Introduce within-subject autocorrelation, allow time-varying trend Semparametric Inference in Joint Models 13
14 Ingredients of joint model: Two linked model components Longitudinal data: At times t i = {t ij, j =1,...,m i } W i (t ij )=X i (t ij )+e(t ij ) e(t ij ) N(0,σ 2 ) of α i, U i (u) α i usually assumed normal (more later... ) Time-to-event data: Proportional hazards model λ i (u) = lim du 0 du 1 pr{u T i <u+ du T i u, Xi H (u),z i } = λ 0 (u)exp{γx i (u)+η T Z i } Dependence on α i through X i (u) Semparametric Inference in Joint Models 14
15 Philosophical interlude: From biological point of view, form of X i (u) may be dictated by beliefs in underlying biological mechanisms E.g., X i (u) =α 0i + α 1i u smooth trend is predominant feature associated with prognosis E.g., X i (u) =α 0i + α 1i u + U i (u) local good or bad periods predominant feature associated with prognosis From empirical point of view, represent hazard in terms of relevant features of longitudinal process Higher-order polynomials, splines versus stochastic process Ease of implementation From practical point of view intermittency motivates assumption that within-subject autocorrelation negligible U i (u) absorbed intoe i (u) Semparametric Inference in Joint Models 15
16 Observed data: (V i, i,w i,z i,t i ), i =1,...,m V i =min(t i,c i ), i = I(T i C i ),W i = {W (t i1 ),...,W(t imi )} T Implementation: Assuming normal α i Two-stage ( Regression calibration ): (i) EBLUPs from normal mixed model fit at each survival time, (ii) insert in Cox partial likelihood reduces but doesn t eliminate bias (Pawitan & Self, 1993; Tsiatis et al., 1995; Dafni & Tsiatis, 1998; Lavalley and DeGruttola, 1996; Bycott & Taylor, 1998) Likelihood: DeGruttola & Tu (1994), Wulfsohn & Tsiatis (1997), Henderson et al. (2000) Bayesian: Faucett & Thomas (1996), Xu and Zeger (2001), Wang and Taylor (2001), Brown & Ibrahim (2003) Issues: Sensitivity to assumptions Computational complexity Semparametric Inference in Joint Models 16
17 ACTG 175: Assuming X i (u) = log CD4, X i (u) =α 0i + α 1i u Raw individual OLS estimates intercept slope Semparametric Inference in Joint Models 17
18 3. CONDITIONAL SCORE APPROACH Our objective: X i (u) =f(u) T α i ; e.g., X i (u) =α 0i + α 1i u Simple method (computationally and conceptually straightforward) that yields consistent, asymptotically normal estimator for γ, η... And that furthermore makes no distributional assumption about the underlying random effects α i Semiparametric model/method (λ 0 (u) and distribution of α i unspecified) Our approach: Exploit the conditional score idea of Stefanski & Carroll (1987) for GLIMs Condition away the random effects α i Semparametric Inference in Joint Models 18
19 Assume: σ 2 known for now Notation: ˆX i (u) OLS estimator for X i (u) based on t i (u) =(t ij <u) (requires 2 measurements prior to u) dn i (u) =I(u V i <u+ du, i =1,t i2 u) puts point mass at u for observed death time if after 2nd measurement Y i (u) =I(V i u, t i2 u) =at risk with 2 measurements at u Semparametric Inference in Joint Models 19
20 Assume: Hazard relationship satisfies λ i (u) = lim i <u+ du T i u, C i u, t i (u), W i (u),α i,z i ) du 0 = lim i <u+ du T i u, C i u, t i (u), ẽ i (u),α i,z i ) du 0 = lim i <u+ du T i u, α i,z i ) du 0 = λ 0 (u)exp{γx i (u)+η T Z i }. W i (u) ={W i (t ij ): t ij <u}, ẽ i (u) ={e i (t ij ): t ij <u} Also assume: Distribution of e i (t ij ) given measurement at t ij,at risk at t ij, measurement history prior to t ij, α i, Z i is N (0,σ 2 ) May be shown with additional assumptions on censoring, timing of measurements ˆX i (u) Y i (u) =1,t i (u),α i,z i N{X i (u),σ 2 θ i (u)} σ 2 θ i (u) =prediction variance [depends on t i (u)] Semparametric Inference in Joint Models 20
21 At any time u: Conditional on being at risk at u pr{dn i (u) =r, ˆX i (u) =x Y i (u) =1,α i,z i,t i (u)} =pr{dn i (u) =r Y i (u) =1, ˆX i (u) =x, α i,z i,t i (u)} pr{ ˆX i (u) =x Y i (u) =1,α i,z i,t i (u)} 1st piece: Bernoulli[λ 0 (u)du exp{γx i (u)+η T Z i }] 2nd piece: N{X i (u),σ 2 θ i (u)} At time u, likelihood {dn i (u), ˆX i (u)} Y i (u) =1,α i,z i,t i (u): To order du [ { γσ 2 θ i (u)dn i (u)+ =exp X i (u) ˆX }] i (u) σ 2 θ i (u) { {λ 0(u)exp(η T Z i )du} dn i(u) exp ˆX } i 2(u)+X2 i (u) {2πσ 2 θ i (u)} 1/2 2σ 2 θ i (u) Semparametric Inference in Joint Models 21
22 Thus: Sufficient statistic for α i (conditional on Y i (u) =1) S i (u, γ, σ 2 )=γσ 2 θ i (u)dn i (u)+ ˆX i (u) Suggestion: Conditioning on S i (u, γ, σ 2 ) would remove the dependence of the conditional distribution on (the nuisance parameter ) α i Form estimating equations in same spirit as the partial likelihood score equations Semparametric Inference in Joint Models 22
23 Usual partial likelihood equations: X i (u) known Intensity lim du 0 du 1 pr{dn i (u) =1 X i (u),z i,y i (u)} = λ 0 (u)exp{γx i (u)+η T Z i }Y i (u) =λ 0 (u)e 0i (u, γ, η) Estimating equations n {X i (u),zi T } T {dn i (u) E 0i (u, γ, η, σ 2 )λ 0 (u)du} =0 i=1 n {dn i (u) E 0i (u, γ, η, σ 2 )λ 0 (u)du} =0 i=1 E 0 (u, γ, η) = n j=1 E 0j(u, γ, η), ˆλ 0 (u)du = n j=1 dn j(u)/e 0 (u, γ, η) Substitute ˆλ 0 (u)du, rearrange, solve in (γ,η) n [ {X i (u),zi T } T E ] 1(u, γ, η) dn i (u) =0 E 0 (u, γ, η) i=1 E 1j (u, γ, η) ={X i (u),zi T } T E 0j (u, γ, η), E 1 (u, γ, η) = n j=1 E 1j(u, γ, η) Semparametric Inference in Joint Models 23
24 Conditional score estimating equations: By analogy Conditional intensity canshow lim du 0 du 1 pr{dn i (u) =1 S i (u, γ, σ 2 ),Z i,t i (u),y i (u)} = λ 0 (u)exp{γs i (u, γ, σ 2 ) γ 2 σ 2 θ i (u)/2+η T Z i }Y i (u) = λ 0 (u)e 0i(u, γ, η, σ 2 ) Estimating equations n {S i (u, γ, σ 2 ),Zi T } T {dn i (u) E0i(u, γ, η, σ 2 )λ 0 (u)du} =0 i=1 n {dn i (u) E0i(u, γ, η, σ 2 )λ 0 (u)du} =0 i=1 E 0 (u, γ, η, σ 2 )= n j=1 E 0j(u, γ, η, σ 2 ), ˆλ 0 (u)du = n j=1 dn j(u)/e 0 (u, γ, η, σ 2 ) Semparametric Inference in Joint Models 24
25 Substitute ˆλ 0 (u)du, rearrange solve in (γ,η) n [ {S i (u, γ, σ 2 ),Zi T } T E 1(u, γ, η, σ 2 ] ) E0 (u, γ, η, dn i (u) =0 σ2 ) i=1 E1j(u, γ, η, σ 2 )={S j (u, γ, σ 2 ),Zj T } T E0j(u, γ, η, σ 2 ), E1 (u, γ, η, σ 2 )= n j=1 E 1j(u, γ, η, σ 2 ) Remarks: Easy to compute Can substitute ˆσ 2 based on OLS residuals Consistent, asymptotically normal, SEs via sandwich Reduces to usual partial likelihood equations when σ 2 =0[so ˆX i (u) =X i (u)] Tsiatis and Davidian (2001, Biometrika), Song, Davidian, and Tsiatis (2002, Biostatistics) Semparametric Inference in Joint Models 25
26 4. SEMIPARAMETRIC LIKELIHOOD APPROACH Potential drawbacks of conditional score approach: Possible inefficiency Includes any distribution for α i Alternative approach: Likelihood formulation Take X i (u) =f T (u)α i as before Make realistic but not overly restrictive assumption on α i Assume: α i have density p(α i )inaclassof smooth densities H (Gallant and Nychka, 1987) Densities in H do not have jumps, kinks or oscillations but may be skewed, multi-modal, fat- or thin-tailed relative to the normal (and the normal is H) α i = g(µ, Z i )+Rz i, R lower triangular, z i has density h H Semparametric Inference in Joint Models 26
27 Representation of h H: h(z) =P (z)ϕ 2 q (z)+small lower bound for tail behavior P (z) is an infinite-dimensional polynomial ϕ q (z) isq-variate standard normal density Practical approximation: Truncate h K (z) =PK(z)ϕ 2 q (z) P K (z) iskth order polynomial; e.g., for K =2 P K (z) =a 00 + a 10 z 1 + a 01 z 2 + a 20 z1 2 + a 02 z2 2 + a 11 z 1 z 2 Vector of coefficients a must satisfy h K (z) dz =1 K =0isstandard normal α i N{g(µ, Z i ),RR T } SemiNonParametric (SNP) Semparametric Inference in Joint Models 27
28 n i=1 As before: Proportional hazard model assumption, normal errors same as for conditional score approach Implementation issues: Polar coordinate reparameterization in terms of vector of angles φ to enforce h K (z) dz = 1, numerical stability Parameters of interest: Ω = {λ 0 ( ),γ,η,σ 2,µ,R,φ} K controls degree of flexibility, must be chosen somehow Need a likelihood function Likelihood: Under certain assumptions on nature of censoring, timing of measurements, likelihood as function of Ω [ λ 0 (V i )exp{γx i (V i )+η T Z i }] i exp [ [ 1 (2πσ 2 ) exp m i/2 m i j=1 Vi 0 ] λ 0 (u)exp{γx i (u)+η T Z i } du ] {W i (t ij ) X i (t ij )} 2 h 2σ 2 K (z i ) dz i Semparametric Inference in Joint Models 28
29 Implementation: EM algorithm E-step: Intractable integration by quadrature M-step: Maximization in (µ, R, φ) and(γ,η,σ 2,λ 0 ) separates; one-step NR update for (γ,η) SEs, CIs by profile likelihood Choice of K: Adaptive based on inspection of information criteria If l K ( Ω) is maximized loglikelihood for fixed K, N = total number of observations, Ω (p 1), minimize { l K ( Ω) + pc(n)}/n AIC, C(N) = 1; BIC, C(N) = log N/2; Hannan-Quinn (HQ), C(N) = log log N AIC prefers larger models, BIC smaller, HQ intermediate Song, Davidian, and Tsiatis (2002), Biometrics Semparametric Inference in Joint Models 29
30 5. SIMULATION EVIDENCE Simple situation: η =0 X i (u) =α 0i + α 1i u E(α i )=(4.173, ) T, γ = 1.0, λ 0 (u) =1, C i exp(1/110), additional censoring at 80 weeks, σ 2 =0.6 Nominal t ij =(0, 2, 4, 8, 16, 24, 32, 40, 48, 56, 64, 72, 80), 10% missing Case 1: α i bivariate mixture of normals Case 2: α i bivariate normal Semparametric Inference in Joint Models 30
31 Methods: I, Ideal, X i (u) known for all u LV, LVCF CS, Conditional score, σ 2 estimated K = 0, likelihood with normal α i SNP, K chosen by HQ Next slide: 200 Monte Carlo replications, n = 200, σ 2 =0.60 MC SD, Average estimated SE, CP of 95% Wald CI Semparametric Inference in Joint Models 31
32 I LV CS K =0 SNP Case 1: Mixture Scenario γ SD SE CP Case 2: Normal Scenario γ SD SE CP Semparametric Inference in Joint Models 32
33 Observations: Conditional score is inefficient (but much faster... ) Interesting robustness to specification of distribution of α i Not shown: Likelihood inference on aspects of α i compromised by incorrect distributional assumption Advantage of likelihood inference: Insight on distribution of α i model refinements Semparametric Inference in Joint Models 33
34 αi αi αi αi1 density density IWSM2003 True and MC average density estimates: Semparametric Inference in Joint Models 34
35 Intercept marginal Average and raw estimates: density α i0 density α i0 Semparametric Inference in Joint Models 35
36 6. EXAMPLE, REVISITED Model: W (t ij )=X i (t ij )+e(t ij ), X i (u) =α 0i + α 1i u, e(t ij ) N(0,σ 2 ) α i = µ 0 (1 Z i )+µ 1 Z i + Rz i, Z i = I(Trt=ZDV) λ i (u) =λ 0 (u)exp{γx i (u)+ηz i } 2467 subjects Conditional Score: Estimate (γ,η,σ 2 ) with no assumptions on α i (so no model for α i ) Likelihood: Assume h Hand estimate {γ,η,σ 2,µ,R,λ 0 (u)} with K =0, 1, 2, 3, 4; all criteria chose K =3or4 Semparametric Inference in Joint Models 36
37 CS K =0 K =2 K =3 K =4 γ (0.207) (0.091) (0.092) (0.092) (0.092) η (0.264) (0.132) (0.132) (0.132) (0.131) loglike AIC HQ BIC Semparametric Inference in Joint Models 37
38 αi αi density 800 IWSM2003 Estimated joint density, K =4: Semparametric Inference in Joint Models 38
39 Estimated slope marginal, K =0, 2, 3, 4: density α i1 Semparametric Inference in Joint Models 39
40 7. DISCUSSION Naive methods yield biased inferences Regression calibration may also yield bias Likelihood or methods like conditional score are required, yield unbiased inferences Conditional score is easy to compute, readily extends to multiple longitudinal processes, and does not require as restrictive assumptions on censoring and timing of measurements But is inefficient, does not accommodate additional stochastic process, only permits inference on hazard parameters Full likelihood approaches are computationally intensive; robustness to violation of assumptions still unclear Philosophical issue: Modeling the longitudinal process X i (u) Semparametric Inference in Joint Models 40
Semiparametric Mixed Effects Models with Flexible Random Effects Distribution
Semiparametric Mixed Effects Models with Flexible Random Effects Distribution Marie Davidian North Carolina State University davidian@stat.ncsu.edu www.stat.ncsu.edu/ davidian Joint work with A. Tsiatis,
More informationJOINT MODELING OF LONGITUDINAL AND TIME-TO-EVENT DATA: AN OVERVIEW
Statistica Sinica 14(2004), 809-834 JOINT MODELING OF LONGITUDINAL AND TIME-TO-EVENT DATA: AN OVERVIEW Anastasios A. Tsiatis and Marie Davidian North Carolina State University Abstract: A common objective
More informationJoint Modeling of Longitudinal and Time-to-Event Data: An Overview
Joint Modeling of Longitudinal and Time-to-Event Data: An Overview Anastasios A. Tsiatis and Marie Davidian Department of Statistics, North Carolina State University Box 8203, Raleigh, North Carolina 27695-8203,
More informationSome New Methods for Latent Variable Models and Survival Analysis. Latent-Model Robustness in Structural Measurement Error Models.
Some New Methods for Latent Variable Models and Survival Analysis Marie Davidian Department of Statistics North Carolina State University 1. Introduction Outline 3. Empirically checking latent-model robustness
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 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 informationReview Article Analysis of Longitudinal and Survival Data: Joint Modeling, Inference Methods, and Issues
Journal of Probability and Statistics Volume 2012, Article ID 640153, 17 pages doi:10.1155/2012/640153 Review Article Analysis of Longitudinal and Survival Data: Joint Modeling, Inference Methods, and
More informationModelling Survival Events with Longitudinal Data Measured with Error
Modelling Survival Events with Longitudinal Data Measured with Error Hongsheng Dai, Jianxin Pan & Yanchun Bao First version: 14 December 29 Research Report No. 16, 29, Probability and Statistics Group
More informationA TWO-STAGE LINEAR MIXED-EFFECTS/COX MODEL FOR LONGITUDINAL DATA WITH MEASUREMENT ERROR AND SURVIVAL
A TWO-STAGE LINEAR MIXED-EFFECTS/COX MODEL FOR LONGITUDINAL DATA WITH MEASUREMENT ERROR AND SURVIVAL Christopher H. Morrell, Loyola College in Maryland, and Larry J. Brant, NIA Christopher H. Morrell,
More informationJoint modelling of longitudinal measurements and event time data
Biostatistics (2000), 1, 4,pp. 465 480 Printed in Great Britain Joint modelling of longitudinal measurements and event time data ROBIN HENDERSON, PETER DIGGLE, ANGELA DOBSON Medical Statistics Unit, Lancaster
More information2 X. Song et al. 1. Introduction In longitudinal studies, information is collected on a time-to-event (e.g. `survival') and timedependent and time-ind
Biostatistics (2001), 1, 1, pp. 1 32 Printed ingreat Britain An estimator for the proportional hazards model with multiple longitudinal covariates measured with error By XIAO SONG, MARIE DAVIDIAN AND ANASTASIOS
More informationImproving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates
Improving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates Anastasios (Butch) Tsiatis Department of Statistics North Carolina State University http://www.stat.ncsu.edu/
More informationBayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features. Yangxin Huang
Bayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features Yangxin Huang Department of Epidemiology and Biostatistics, COPH, USF, Tampa, FL yhuang@health.usf.edu January
More informationLatent-model Robustness in Joint Models for a Primary Endpoint and a Longitudinal Process
Latent-model Robustness in Joint Models for a Primary Endpoint and a Longitudinal Process Xianzheng Huang, 1, * Leonard A. Stefanski, 2 and Marie Davidian 2 1 Department of Statistics, University of South
More informationMixed model analysis of censored longitudinal data with flexible random-effects density
Biostatistics (2012), 13, 1, pp. 61 73 doi:10.1093/biostatistics/kxr026 Advance Access publication on September 13, 2011 Mixed model analysis of censored longitudinal data with flexible random-effects
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 informationParametric Joint Modelling for Longitudinal and Survival Data
Parametric Joint Modelling for Longitudinal and Survival Data Paraskevi Pericleous Doctor of Philosophy School of Computing Sciences University of East Anglia July 2016 c This copy of the thesis has been
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 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 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 informationEstimation of Conditional Kendall s Tau for Bivariate Interval Censored Data
Communications for Statistical Applications and Methods 2015, Vol. 22, No. 6, 599 604 DOI: http://dx.doi.org/10.5351/csam.2015.22.6.599 Print ISSN 2287-7843 / Online ISSN 2383-4757 Estimation of Conditional
More informationOn Estimating the Relationship between Longitudinal Measurements and Time-to-Event Data Using a Simple Two-Stage Procedure
Biometrics DOI: 10.1111/j.1541-0420.2009.01324.x On Estimating the Relationship between Longitudinal Measurements and Time-to-Event Data Using a Simple Two-Stage Procedure Paul S. Albert 1, and Joanna
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 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 informationEstimating the Mean Response of Treatment Duration Regimes in an Observational Study. Anastasios A. Tsiatis.
Estimating the Mean Response of Treatment Duration Regimes in an Observational Study Anastasios A. Tsiatis http://www.stat.ncsu.edu/ tsiatis/ Introduction to Dynamic Treatment Regimes 1 Outline Description
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 informationSemiparametric Regression
Semiparametric Regression Patrick Breheny October 22 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Introduction Over the past few weeks, we ve introduced a variety of regression models under
More informationConditional Estimation for Generalized Linear Models When Covariates Are Subject-specific Parameters in a Mixed Model for Longitudinal Measurements
Conditional Estimation for Generalized Linear Models When Covariates Are Subject-specific Parameters in a Mixed Model for Longitudinal Measurements Erning Li, Daowen Zhang, and Marie Davidian Department
More informationPersonalized Treatment Selection Based on Randomized Clinical Trials. Tianxi Cai Department of Biostatistics Harvard School of Public Health
Personalized Treatment Selection Based on Randomized Clinical Trials Tianxi Cai Department of Biostatistics Harvard School of Public Health Outline Motivation A systematic approach to separating subpopulations
More informationDiscussion of Missing Data Methods in Longitudinal Studies: A Review by Ibrahim and Molenberghs
Discussion of Missing Data Methods in Longitudinal Studies: A Review by Ibrahim and Molenberghs Michael J. Daniels and Chenguang Wang Jan. 18, 2009 First, we would like to thank Joe and Geert for a carefully
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 informationA Sampling of IMPACT Research:
A Sampling of IMPACT Research: Methods for Analysis with Dropout and Identifying Optimal Treatment Regimes Marie Davidian Department of Statistics North Carolina State University http://www.stat.ncsu.edu/
More informationStatistical Methods for Joint Analysis of Survival Time and Longitudinal Data
Statistical Methods for Joint Analysis of Survival Time and Longitudinal Data Jaeun Choi A dissertation sumitted to the faculty of the University of North Carolina at Chapel Hill in partial fulfillment
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 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 informationFUNCTIONAL PRINCIPAL COMPONENT ANALYSIS FOR LONGITUDINAL AND SURVIVAL DATA
Statistica Sinica 17(27), 965-983 FUNCTIONAL PRINCIPAL COMPONENT ANALYSIS FOR LONGITUDINAL AND SURVIVAL DATA Fang Yao University of Toronto Abstract: This paper proposes a nonparametric approach for jointly
More informationA Bayesian Nonparametric Approach to Monotone Missing Data in Longitudinal Studies with Informative Missingness
A Bayesian Nonparametric Approach to Monotone Missing Data in Longitudinal Studies with Informative Missingness A. Linero and M. Daniels UF, UT-Austin SRC 2014, Galveston, TX 1 Background 2 Working model
More informationJoint Modeling of Survival and Longitudinal Data: Likelihood Approach Revisited
Biometrics 62, 1037 1043 December 2006 DOI: 10.1111/j.1541-0420.2006.00570.x Joint Modeling of Survival and Longitudinal Data: Likelihood Approach Revisited Fushing Hsieh, 1 Yi-Kuan Tseng, 2 and Jane-Ling
More informationSemiparametric Generalized Linear Models
Semiparametric Generalized Linear Models North American Stata Users Group Meeting Chicago, Illinois Paul Rathouz Department of Health Studies University of Chicago prathouz@uchicago.edu Liping Gao MS Student
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 informationApproximation of Survival Function by Taylor Series for General Partly Interval Censored Data
Malaysian Journal of Mathematical Sciences 11(3): 33 315 (217) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal Approximation of Survival Function by Taylor
More informationIntegrated Likelihood Estimation in Semiparametric Regression Models. Thomas A. Severini Department of Statistics Northwestern University
Integrated Likelihood Estimation in Semiparametric Regression Models Thomas A. Severini Department of Statistics Northwestern University Joint work with Heping He, University of York Introduction Let Y
More informationBayesian Semiparametric Methods for Joint Modeling Longitudinal and Survival Data
Bayesian Semiparametric Methods for Joint Modeling Longitudinal and Survival Data Tim Hanson Division of Biostatistics University of Minnesota and Adam Branscum University of Kentucky Departments of Biostatistics,
More information6 Pattern Mixture Models
6 Pattern Mixture Models A common theme underlying the methods we have discussed so far is that interest focuses on making inference on parameters in a parametric or semiparametric model for the full data
More informationWeb-based Supplementary Materials for A Robust Method for Estimating. Optimal Treatment Regimes
Biometrics 000, 000 000 DOI: 000 000 0000 Web-based Supplementary Materials for A Robust Method for Estimating Optimal Treatment Regimes Baqun Zhang, Anastasios A. Tsiatis, Eric B. Laber, and Marie Davidian
More informationSemiparametric Models for Joint Analysis of Longitudinal Data and Counting Processes
Semiparametric Models for Joint Analysis of Longitudinal Data and Counting Processes by Se Hee Kim A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial
More informationSeparate and Joint Modeling of Longitudinal and Event Time Data Using Standard Computer Packages
Separate and Joint Modeling of Longitudinal and Event Time Data Using Standard Computer Packages Xu GUO and Bradley P. CARLIN Many clinical trials and other medical and reliability studies generate both
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
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 informationOptimal Treatment Regimes for Survival Endpoints from a Classification Perspective. Anastasios (Butch) Tsiatis and Xiaofei Bai
Optimal Treatment Regimes for Survival Endpoints from a Classification Perspective Anastasios (Butch) Tsiatis and Xiaofei Bai Department of Statistics North Carolina State University 1/35 Optimal Treatment
More informationNon-maximum likelihood estimation and statistical inference for linear and nonlinear mixed models
Optimum Design for Mixed Effects Non-Linear and generalized Linear Models Cambridge, August 9-12, 2011 Non-maximum likelihood estimation and statistical inference for linear and nonlinear mixed models
More informationGENERALIZED LINEAR MIXED MODELS AND MEASUREMENT ERROR. Raymond J. Carroll: Texas A&M University
GENERALIZED LINEAR MIXED MODELS AND MEASUREMENT ERROR Raymond J. Carroll: Texas A&M University Naisyin Wang: Xihong Lin: Roberto Gutierrez: Texas A&M University University of Michigan Southern Methodist
More informationOn Measurement Error Problems with Predictors Derived from Stationary Stochastic Processes and Application to Cocaine Dependence Treatment Data
On Measurement Error Problems with Predictors Derived from Stationary Stochastic Processes and Application to Cocaine Dependence Treatment Data Yehua Li Department of Statistics University of Georgia Yongtao
More informationJoint Modelling of Accelerated Failure Time and Longitudinal Data
Joint Modelling of Accelerated Failure Time and Longitudinal Data YI-KUAN TSENG, FUSHING HSIEH and JANE-LING WANG yktseng@wald.ucdavis.edu fushing@wald.ucdavis.edu wang@wald.ucdavis.edu Department of Statistics,
More informationIntegrated likelihoods in survival models for highlystratified
Working Paper Series, N. 1, January 2014 Integrated likelihoods in survival models for highlystratified censored data Giuliana Cortese Department of Statistical Sciences University of Padua Italy Nicola
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 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 informationOn the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models
On the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models Thomas Kneib Department of Mathematics Carl von Ossietzky University Oldenburg Sonja Greven Department of
More informationEstimation with clustered censored survival data with missing covariates in the marginal Cox model
Estimation with clustered censored survival data with missing covariates in the marginal Cox model Michael Parzen 1, Stuart Lipsitz 2,Amy Herring 3, and Joseph G. Ibrahim 3 (1) Emory University (2) Harvard
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 informationBiost 518 Applied Biostatistics II. Purpose of Statistics. First Stage of Scientific Investigation. Further Stages of Scientific Investigation
Biost 58 Applied Biostatistics II Scott S. Emerson, M.D., Ph.D. Professor of Biostatistics University of Washington Lecture 5: Review Purpose of Statistics Statistics is about science (Science in the broadest
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 informationChapter 2 Inference on Mean Residual Life-Overview
Chapter 2 Inference on Mean Residual Life-Overview Statistical inference based on the remaining lifetimes would be intuitively more appealing than the popular hazard function defined as the risk of immediate
More informationUSING TRAJECTORIES FROM A BIVARIATE GROWTH CURVE OF COVARIATES IN A COX MODEL ANALYSIS
USING TRAJECTORIES FROM A BIVARIATE GROWTH CURVE OF COVARIATES IN A COX MODEL ANALYSIS by Qianyu Dang B.S., Fudan University, 1995 M.S., Kansas State University, 2000 Submitted to the Graduate Faculty
More informationConstrained Maximum Likelihood Estimation for Model Calibration Using Summary-level Information from External Big Data Sources
Constrained Maximum Likelihood Estimation for Model Calibration Using Summary-level Information from External Big Data Sources Yi-Hau Chen Institute of Statistical Science, Academia Sinica Joint with Nilanjan
More informationA penalized likelihood approach to joint modeling of longitudinal measurements and time-to-event data
Statistics and Its Interface Volume 1 2008 33 45 A penalized likelihood approach to joint modeling of longitudinal measurements and time-to-event data Wen Ye, Xihong Lin and Jeremy M. G. Taylor Recently
More informationLikelihood-based inference with missing data under missing-at-random
Likelihood-based inference with missing data under missing-at-random Jae-kwang Kim Joint work with Shu Yang Department of Statistics, Iowa State University May 4, 014 Outline 1. Introduction. Parametric
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 informationGroup Sequential Tests for Delayed Responses. Christopher Jennison. Lisa Hampson. Workshop on Special Topics on Sequential Methodology
Group Sequential Tests for Delayed Responses Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj Lisa Hampson Department of Mathematics and Statistics,
More informationLatent Variable Model for Weight Gain Prevention Data with Informative Intermittent Missingness
Journal of Modern Applied Statistical Methods Volume 15 Issue 2 Article 36 11-1-2016 Latent Variable Model for Weight Gain Prevention Data with Informative Intermittent Missingness Li Qin Yale University,
More informationPairwise rank based likelihood for estimating the relationship between two homogeneous populations and their mixture proportion
Pairwise rank based likelihood for estimating the relationship between two homogeneous populations and their mixture proportion Glenn Heller and Jing Qin Department of Epidemiology and Biostatistics Memorial
More 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 informationOn the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models
On the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models Thomas Kneib Department of Mathematics Carl von Ossietzky University Oldenburg Sonja Greven Department of
More informationSurvival Prediction Under Dependent Censoring: A Copula-based Approach
Survival Prediction Under Dependent Censoring: A Copula-based Approach Yi-Hau Chen Institute of Statistical Science, Academia Sinica 2013 AMMS, National Sun Yat-Sen University December 7 2013 Joint work
More informationStatistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach
Statistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach Jae-Kwang Kim Department of Statistics, Iowa State University Outline 1 Introduction 2 Observed likelihood 3 Mean Score
More informationStatistical Methods for Alzheimer s Disease Studies
Statistical Methods for Alzheimer s Disease Studies Rebecca A. Betensky, Ph.D. Department of Biostatistics, Harvard T.H. Chan School of Public Health July 19, 2016 1/37 OUTLINE 1 Statistical collaborations
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 informationContents. Part I: Fundamentals of Bayesian Inference 1
Contents Preface xiii Part I: Fundamentals of Bayesian Inference 1 1 Probability and inference 3 1.1 The three steps of Bayesian data analysis 3 1.2 General notation for statistical inference 4 1.3 Bayesian
More informationPartly Conditional Survival Models for Longitudinal Data
UW Biostatistics Working Paper Series 12-19-23 Partly Conditional Survival Models for Longitudinal Data Yingye Zheng Fred Hutchinson Cancer Research Center, yzheng@fhcrc.org Patrick Heagerty 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 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 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 informationAnalysing repeated measurements whilst accounting for derivative tracking, varying within-subject variance and autocorrelation: the xtiou command
Analysing repeated measurements whilst accounting for derivative tracking, varying within-subject variance and autocorrelation: the xtiou command R.A. Hughes* 1, M.G. Kenward 2, J.A.C. Sterne 1, K. Tilling
More informationStat 579: Generalized Linear Models and Extensions
Stat 579: Generalized Linear Models and Extensions Linear Mixed Models for Longitudinal Data Yan Lu April, 2018, week 15 1 / 38 Data structure t1 t2 tn i 1st subject y 11 y 12 y 1n1 Experimental 2nd subject
More information18 Bivariate normal distribution I
8 Bivariate normal distribution I 8 Example Imagine firing arrows at a target Hopefully they will fall close to the target centre As we fire more arrows we find a high density near the centre and fewer
More informationEstimating subgroup specific treatment effects via concave fusion
Estimating subgroup specific treatment effects via concave fusion Jian Huang University of Iowa April 6, 2016 Outline 1 Motivation and the problem 2 The proposed model and approach Concave pairwise fusion
More informationBayesian Quantile Regression for Longitudinal Studies with Nonignorable Missing Data
Biometrics 66, 105 114 March 2010 DOI: 10.1111/j.1541-0420.2009.01269.x Bayesian Quantile Regression for Longitudinal Studies with Nonignorable Missing Data Ying Yuan and Guosheng Yin Department of Biostatistics,
More informationModels for Multivariate Panel Count Data
Semiparametric Models for Multivariate Panel Count Data KyungMann Kim University of Wisconsin-Madison kmkim@biostat.wisc.edu 2 April 2015 Outline 1 Introduction 2 3 4 Panel Count Data Motivation Previous
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 informationFractional Hot Deck Imputation for Robust Inference Under Item Nonresponse in Survey Sampling
Fractional Hot Deck Imputation for Robust Inference Under Item Nonresponse in Survey Sampling Jae-Kwang Kim 1 Iowa State University June 26, 2013 1 Joint work with Shu Yang Introduction 1 Introduction
More informationGeneralized Linear. Mixed Models. Methods and Applications. Modern Concepts, Walter W. Stroup. Texts in Statistical Science.
Texts in Statistical Science Generalized Linear Mixed Models Modern Concepts, Methods and Applications Walter W. Stroup CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint
More informationIntegrated approaches for analysis of cluster randomised trials
Integrated approaches for analysis of cluster randomised trials Invited Session 4.1 - Recent developments in CRTs Joint work with L. Turner, F. Li, J. Gallis and D. Murray Mélanie PRAGUE - SCT 2017 - Liverpool
More informationPQL Estimation Biases in Generalized Linear Mixed Models
PQL Estimation Biases in Generalized Linear Mixed Models Woncheol Jang Johan Lim March 18, 2006 Abstract The penalized quasi-likelihood (PQL) approach is the most common estimation procedure for the generalized
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 informationGeneralized Linear Models with Functional Predictors
Generalized Linear Models with Functional Predictors GARETH M. JAMES Marshall School of Business, University of Southern California Abstract In this paper we present a technique for extending generalized
More informationCensoring mechanisms
Censoring mechanisms Patrick Breheny September 3 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Fixed vs. random censoring In the previous lecture, we derived the contribution to the likelihood
More informationDESIGN CONSIDERATIONS FOR COMPLEX SURVIVAL MODELS
DESIGN CONSIDERATIONS FOR COMPLEX SURVIVAL MODELS Liddy Miaoli Chen A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial fulfillment of the requirements
More informationA Study of Joinpoint Models for Longitudinal Data
UNLV Theses, Dissertations, Professional Papers, and Capstones 8-1-2014 A Study of Joinpoint Models for Longitudinal Data Libo Zhou University of Nevada, Las Vegas, zhoul@unlv.nevada.edu Follow this and
More informationApproximate Bayesian Computation
Approximate Bayesian Computation Michael Gutmann https://sites.google.com/site/michaelgutmann University of Helsinki and Aalto University 1st December 2015 Content Two parts: 1. The basics of approximate
More informationStatistical Distribution Assumptions of General Linear Models
Statistical Distribution Assumptions of General Linear Models Applied Multilevel Models for Cross Sectional Data Lecture 4 ICPSR Summer Workshop University of Colorado Boulder Lecture 4: Statistical Distributions
More information