Modeling and Analysis of Recurrent Event Data
|
|
- Griselda Hudson
- 6 years ago
- Views:
Transcription
1 Edsel A. Peña Department of Statistics University of South Carolina Columbia, SC New Jersey Institute of Technology Conference May 20, 2012
2 Historical Perspective: Random Censorship Model (RCM) Random Observables: T 1, T 2,...,T n IID F C 1, C 2,...,C n IID G F and G not related {T i } {C i } (Z 1, δ 1 ), (Z 2, δ 2 ),...,(Z n, δ n ) Z i = T i C i and δ i = I {T i C i } Goal: To make inference on the distribution F or the hazard Λ = df/ F, or functionals (mean, median, etc.).
3 Nonparametric Inference F F = space of continuous distributions Λ = log F; F = 1 F = exp( Λ); λ = Λ N(s) = i I {Z i s; δ i = 1} and Y (s) = i I {Z i s} NAE: ˆΛ = dn Y and PLE: ˆ F = [ 1 dn Y ] Properties of ˆΛ and ˆF well-known, e.g., biased; consistent; and when normalized, weakly convergent to Gaussian processes. Avar(ˆ F(t)) = 1 n F(t) 2 t 0 df(s) F(s) 2Ḡ(s)
4 An Informative RCM Koziol-Green (KG) Model (1976): β 0, Ḡ(t) = F(t) β Lehmann-type alternatives Proportional hazards: Z i = min(t i, C i ) F β+1 Λ G = βλ F δ i = I {T i C i } BER(1/(β + 1)) An important characterizing property: Z i δ i
5 MMC Data: First Event Times and Tau s Migratory Motor Complex (MMC) data set from Aalen and Husebye, Stat Med, MMC Data Set Unit Number Calendar Time
6 MMC Data: Hazard Estimates KG model holds if and only if Λ τ Λ T1 MMC Data Set Cumulative Hazard Estimate First Event Time and Tau
7 KG Model s Utility Chen, Hollander and Langberg (JASA, 1982): exploited independence between Z i and δ i to obtain the exact bias, variance, and MSE functions of PLE. Comparisons with asymptotic results. Cheng and Lin (1987): exploited semiparametric nature of KG model to obtain a more efficient estimator of F compared to the PLE. Hollander and Peña (1988): obtained better confidence bands for F under KG model. Csorgo & Faraway (1998): KG model not practically viable, but serves as a mathematical specimen (a lá yeast) for examining exact properties of procedures and for providing pinpoint assessment of efficiency losses/gains.
8 Recurrent Events: Some Examples admission to hospital due to chronic disease tumor re-occurrence migraine attacks alcohol or drug (eg cocaine) addiction machine failure or discovery of a bug in a software commission of a criminal act by a delinquent minor! major disagreements between a couple non-life insurance claim drop of 200 points in DJIA during trading day publication of a research paper by a professor
9 Full MMC Data: With Recurrences n = 19 subjects; event = end of migratory motor complex cycle; random length of monitoring period per subject. MMC Data Set Unit Number Calendar Time
10 Data Accrual: One Subject
11 Some Aspects in Recurrent Data random monitoring length (τ). random # of events (K) and sum-quota constraint: k K K = max k : T j τ with T j τ < j=1 j=1 Basic Observable: (K, τ, T 1, T 2,...,T K, τ S K ) always a right-censored observation. dependent and informative censoring. K+1 effects of covariates, frailties, interventions after each event, and accumulation of events. j=1 T j
12 Simplest Model: One Subject T 1, T 2,... IID F: (renewal model) perfect interventions after each event τ G F and G not related no covariates (X) no frailties (Z) F could be parametric or nonparametric Peña, Strawderman and Hollander (JASA, 01): nonparametric estimation of F.
13 Nonparametric Estimation of F Y (t) = N(t) = n K i I {T ij t} i=1 j=1 n K i I {T ij t} + I {τ i S iki t} i=1 j=1 t dn(w) GNAE : Λ(t) = 0 Y (w) t [ GPLE : F(t) = 1 dn(w) ] Y (w) 0
14 Main Asymptotic Result k kth Convolution: F (k) (t) = Pr T j t Renewal Function: ρ(t) = ν(t) = 1 Ḡ(t) σ 2 (t) = F(t) 2 t 0 t j=1 F (k) (t) k=1 ρ(w t)dg(w) df(w) F(w) 2Ḡ(w)[1 + ν(w)] Theorem (JASA, 01): n( F(t) F(t)) GP(0, σ 2 (t))
15 Extending KG Model: Recurrent Setting Wanted: a tractable model with monitoring time informative about F. Potential to refine analysis of efficiency gains/losses. Idea: Why not simply generalize the KG model for the RCM. Generalized KG Model (GKG) for Recurrent Events: β > 0, Ḡ(t) = F(t) β with β unknown, and F the common inter-event time distribution function. Remark: τ may also represent system failure/death, while the recurrent event could be shocks to the system. Remark: Association (within unit) could be modeled through a frailty.
16 Estimation Issues and Some Questions How to semiparametrically estimate β, Λ, and F? Parametric estimation in Adekpedjou, Peña, and Quiton (2010, JSPI). How much efficiency loss is incurred when the informative monitoring model structure is ignored? How much penalty is incurred with Single-event analysis relative to Recurrent-event analysis? In particular, what is the efficiency loss for estimating F when using the nonparametric estimator in PSH (2001) relative to the semiparametric estimator that exploits the informative monitoring structure?
17 Basic Processes S ij = j k=1 T ik N i (s) = I {S ij s} j=1 Y i (s) = I {τ i s} R i (s) = s S in = backward recurrence time i (s ) A i (s) = s 0 Y i (v)λ[r i (v)]dv N τ i (s) = I {τ i s} Y τ i (s) = I {τ i s}
18 Transformed Processes Z i (s, t) = I {R i (s) t} N i (s, t) = Y i (s, t) = A i (s, t) = s 0 N i (s ) j=1 s 0 N Z i (v, t)n i (s)) i (dv) = I {T ij t} j=1 I {T ij t} + I {(s τ i ) S in i (s ) t} t Z i (v, t)a i (dv) = Y i (s, w)λ(w)dw 0
19 Aggregated Processes N(s, t) = Y (s, t) = A(s, t) = n N i (s, t) i=1 n Y i (s, t) i=1 n A i (s, t) i=1 N τ (s) = Y τ (s) = n Ni τ (s) i=1 n i=1 Y τ i (s)
20 First, Assume β Known Via Method-of-Moments Approach, estimator of Λ: ˆΛ(s, t β) = t 0 { N(s, dw) + N τ } (dw) Y (s, w) + βy τ (w) Using product-integral representation of F in terms of Λ, estimator of F: ˆ F(s, t β) = t w=0 { 1 N(s, dw) + } Nτ (dw) Y (s, w) + βy τ (w)
21 Estimating β: Profile Likelihood MLE Profile Likelihood: L P (s ; β) = β Nτ (s ) {[ n s { i=1 Estimator of β: [ s v=0 { v=0 1 Y (s, v) + βy τ (v) 1 Y (s, v) + βy τ (v) ˆβ = arg max β L P(s ; β) } N τ i ( v) ] } Ni (s, v) ]}
22 Is It A Legitimate Profile Likelihood? Let C be the space of all cumulative hazard functions defined on [0, ). Theorem Let L(s ; Λ, β) be the full likelihood function given by L = { n s i=1 v=0 exp { } [ ] N Y i (v) i (v)λ[r i (v)]dv [Y τ i (v)βλ(v)dv] Nτ i (v) [ n s i=1 0 ]} s Y i (v)λ[r i (v)]dv + Yi τ (v)βλ(v)dv. 0 then L P (s ; β) = max Λ C L(s ; Λ, β).
23 Estimators of Λ and F Estimator of Λ: ˆΛ(s, t) = ˆΛ(s, t ˆβ) = Estimator of F: ˆ F(s, t) = ˆ F(s, t ˆβ) = t 0 { } N(s, dw) + N τ (dw) Y (s, w) + ˆβY τ (w) { } t 1 N(s, dw) + N τ (dw) Y (s, w) + ˆβY τ (w) w=0
24 Illustrative Data (n = 30): GKG[Wei(2,.1), β =.2] Recurrent Event Data Subject Calendar Time
25 Estimates of β and F ˆβ =.2331 Estimates of SF for Illustrative Data (Blue=GKG; Red=PSH; Green=TRUE) Estimate Survivor Probability Time
26 Properties of Estimators G s (w) = G(w)I {w < s} + I {w s} E{Y 1 (s, t)} y(s, t) = F(t)Ḡ s (t) + F(t) E{Y τ 1 (t)} y τ (t) = F(t) β t ρ(w t)dg s (w) True Values = (F 0, Λ 0, β 0 ) y 0 (s, t) = y(s, t; Λ 0, β 0 ) y τ 0 (s) = y τ (s; Λ 0, β 0 )
27 Existence, Consistency, Normality Theorem There is a sequence of ˆβ that is consistent, and ˆΛ(s, ) and ˆ F(s, ) are both uniformly strongly consistent. Theorem As n, we have n(ˆβ β 0 ) N(0, [I P (s ; Λ 0, β 0 )] 1 ) with I P (s ; Λ 0, β 0 ) = 1 β 0 s 0 y τ 0 (v)y 0(s, v) y 0 (s, v) + β 0 y τ 0 (v)λ 0(v)dv.
28 Weak Convergence of ˆΛ(s, ) Theorem As n, { n[ˆλ(s, t) Λ 0 (t)] : t [0, t ]} converges weakly to a zero-mean Gaussian process with variance function σ 2ˆΛ (s, t) = [ s 0 [ t 0 t 0 Λ 0 (dv) y 0 (s, v) + β 0 y τ 0 (v)+ y 0 (s, v)y τ 0 (v) β 0 [y 0 (s, v) + β 0 y τ 0 (v)]λ 0(dv) y τ 0 (v) y 0 (s, v) + β 0 y τ 0 (v)λ 0(dv) ] 2. ] 1 Remark: The last product term is the effect of estimating β. It inflates the asymptotic variance.
29 Weak Convergence of ˆ F(s, ) and F(s, ) Corollary As n, { n[ˆ F(s, t) F 0 (t)] : t [0, t ]} converges weakly to a zero-mean Gaussian process whose variance function is σ 2ˆF (s, t) = F 0 (t) 2 σ 2ˆΛ (s, t) F 0 (t) 2 σ 2 Λ(s, t). Recall/Compare! Theorem (PSH, 2001) As n, { n[ F(s, t) F 0 (t)] : t [0, t ]} converges weakly to a zero-mean Gaussian process whose variance function is σ 2 F(s, t) = F 0 (t) 2 t 0 Λ 0 (dv) y 0 (s, v).
30 Asymptotic Relative Efficiency: β 0 Known If we know β 0 : ARE{ F(s, t) : ˆ F(s, t β 0 )} = { t } 1 Λ 0 (dw) 0 y 0 (s, w) { t } Λ 0 (dw) y 0 (s, w) + β 0 y0 τ(w) 0 Clearly, less than or equal to unity, as is to be expected.
31 Case of Exponential F: β 0 Known Theorem If F 0 (t) = exp{ θ 0 t} for t 0 and s, then Also, t 0, ARE{ F(, t) : ˆ F(, t β 0 )} = { } 1 1 du F 0 (t) (1 + β 0 )u 2+β 0 { } 1 du (1 + β 0 )u 2+β 0 + β 2 0 u 1+β. 0 F 0 (t) ARE{ F(, t) : ˆ F(, t; β 0 )} 1 + β β 0 + β0 2.
32 ARE-Plots; β 0 {.1,.3,.5,.7,.9, 1.0, 1.5} Known; F = Exponential Efficiency(PSH:GKG0) F(t)
33 Case of β 0 Unknown As to be expected, if β 0 is known, then the estimator exploiting the GKG structure is more efficient. Question: Does this dominance hold true still if β 0 is now estimated? Theorem Under the GKG model, for all ( F 0, β 0 ) with β 0 > 0, F(s, t) is asymptotically dominated by ˆ F(s, t) in the sense that ARE( F(s, t) : ˆ F(s, t)) 1. Proof. Neat application of Cauchy-Schwartz Inequality.
34 ARE-Plots; β 0 {.1,.3,.5,.7,.9, 1.0, 1.5} Unknown; F = Exponential Efficiency(PSH:GKG) F(t)
35 Assessing Asymptotic Approximations under Exponential n = 50 n = 100 β 0 Mean SE ASE Mean SE ASE Mean: mean of the simulated values of ˆβ. SE: standard error of the simulated values of ˆβ. ASE: asymptotic standard error obtained from earlier formula replications were performed.
36 Simulations under Weibull F: ˆβ Results n = 30,α = 2 n = 50,α = 2 n = 30,α =.9 n = 50,α =.9 β 0 Mean SE Mean SE Mean SE Mean SE
37 Simulated Rel Eff of F : ˆ F under a Weibull F with α = 2 Simulated Efficiencies of PLEPSH relative to PLEGKG Effi(PSH:GKG) F(t)
38 Simulated Rel Eff of F : ˆ F under a Weibull F with α = 0.9 Simulated Efficiencies of PLEPSH relative to PLEGKG Effi(PSH:GKG) F(t)
39 Mis-Specified Model: F 0 is Exp(1) and G 0 is Uniform[0, θ] Simulated Efficiencies of PLEPSH relative to PLEGKG Effi(PSH:GKG) Theta = 1 Theta = 2 Theta = 4 Theta = F(t)
40 Mis-Specified Model: F 0 is Weibull(2, 1) and G 0 is Exp(θ) Simulated Efficiencies of PLEPSH relative to PLEGKG Effi(PSH:GKG) Theta = 1.2 Theta = 1 Theta = 0.5 Theta = F(t)
41 Briefly: On More General Recurrent Event Models Must be able to model the following: comparing treatment groups. effect of covariates or regressor variables. effect of interventions after each event. effect of accumulating event occurrences. modeling association of inter-event times within a unit or subject.
42 Recall: Data Accrual for One Subject
43 Bladder Cancer Data Set: Groups [control (red), thiotepa (blue)]; Covariates; in WLW (89) Subject Calendar Time
44 A General Class of Dynamic Recurrent Event Models P{dN i (s) = 1 F s, Z i } = Z i λ 0 [E i (s)]ρ[(s, N i (s ));α]ψ[x i(s)β]ds λ 0 ( ): baseline hazard rate function; nonparametric. E i ( ): predictable effective (real) age process; effect of interventions. Ei (s) = s: minimal intervention. Ei (s) = s S in i (s ): perfect intervention. ρ(, ; α): known function; effect of event occurrences. ψ( ): known link function; effect of covariate X i. Z i : unobserved frailty variable; distribution H( ; η). Model Parameters: (λ 0 ( ), α, β, η).
45 Some Inferential Aspects for General Model Peña and Hollander (2004): proposed and discussed model, especially its generality. Peña, Slate and Gonzalez (2007): semiparametric estimation of model parameters. Stocker and Peña (2007): parametric estimation of model parameters. Peña (2012?): asymptotic properties of the estimators of semiparametric estimators of model parameters. Other inference methods: goodness-of-fit; group comparisons; model validation; estimating parameters in effective age process; Bayesian Approach. Challenge: need to gather effective age or real age of subjects/units/components in real studies.
46 Some Concluding Thoughts Revisited RCM with one event and Koziol-Green (KG) model. Revisited nonparametric estimation of inter-event distribution with recurrent event data. Extended KG model to recurrent event settings. Semiparametric estimation under this GKG recurrent event model. Examined properties of estimators under GKG structure. Efficiency losses of fully nonparametric estimator relative to GKG-based estimator. How about a Bayesian Approach to Inference? Touched on general dynamic models for recurrent event data.
47 Acknowledgements Work on GKG model joint with Akim Adekpedjou (Missouri University of Science and Technology). Research support from the National Science Foundation and National Institutes of Health. Thanks to Professor Sundar Subramanian and the organizers for inviting me to this conference.
Semiparametric Estimation for a Generalized KG Model with R. Model with Recurrent Event Data
Semiparametric Estimation for a Generalized KG Model with Recurrent Event Data Edsel A. Peña (pena@stat.sc.edu) Department of Statistics University of South Carolina Columbia, SC 29208 MMR Conference Beijing,
More informationModelling and Analysis of Recurrent Event Data
Modelling and Analysis of Recurrent Event Data Edsel A. Peña Department of Statistics University of South Carolina Research support from NIH, NSF, and USC/MUSC Collaborative Grants Joint work with Prof.
More informationRecurrent Event Data: Models, Analysis, Efficiency
Recurrent Event Data: Models, Analysis, Efficiency Edsel A. Peña Department of Statistics University of South Carolina Columbia, SC 29208 Talk at Brown University March 10, 2008 Recurrent Event Data:Models,
More informationTime-to-Event Modeling and Statistical Analysis
PeñaTalkGT) p. An Appetizer The important thing is not to stop questioning. Curiosity has its own reason for existing. One cannot help but be in awe when he contemplates the mysteries of eternity, of life,
More informationModeling and Non- and Semi-Parametric Inference with Recurrent Event Data
Modeling and Non- and Semi-Parametric Inference with Recurrent Event Data Edsel A. Peña Department of Statistics University of South Carolina Columbia, SC 29208 Talk at Columbia University November 19,
More informationDynamic Models in Reliability and Survival Analysis
Dynamic Models in Reliability and Survival Analysis Edsel A. Peña Department of Statistics University of South Carolina (USC) Joint work with E. Slate (MUSC) MMR 2004, June 21-24, 2004 pena@stat.sc.edu
More informationIntroduction to repairable systems STK4400 Spring 2011
Introduction to repairable systems STK4400 Spring 2011 Bo Lindqvist http://www.math.ntnu.no/ bo/ bo@math.ntnu.no Bo Lindqvist Introduction to repairable systems Definition of repairable system Ascher and
More informatione 4β e 4β + e β ˆβ =0.765
SIMPLE EXAMPLE COX-REGRESSION i Y i x i δ i 1 5 12 0 2 10 10 1 3 40 3 0 4 80 5 0 5 120 3 1 6 400 4 1 7 600 1 0 Model: z(t x) =z 0 (t) exp{βx} Partial likelihood: L(β) = e 10β e 10β + e 3β + e 5β + e 3β
More 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 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 informationEfficient Semiparametric Estimators via Modified Profile Likelihood in Frailty & Accelerated-Failure Models
NIH Talk, September 03 Efficient Semiparametric Estimators via Modified Profile Likelihood in Frailty & Accelerated-Failure Models Eric Slud, Math Dept, Univ of Maryland Ongoing joint project with Ilia
More information1 Glivenko-Cantelli type theorems
STA79 Lecture Spring Semester Glivenko-Cantelli type theorems Given i.i.d. observations X,..., X n with unknown distribution function F (t, consider the empirical (sample CDF ˆF n (t = I [Xi t]. n Then
More informationEstimating Load-Sharing Properties in a Dynamic Reliability System. Paul Kvam, Georgia Tech Edsel A. Peña, University of South Carolina
Estimating Load-Sharing Properties in a Dynamic Reliability System Paul Kvam, Georgia Tech Edsel A. Peña, University of South Carolina Modeling Dependence Between Components Most reliability methods are
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 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 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 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 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 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 informationNon- and Semi-parametric Bayesian Inference with Recurrent Events and Coherent Systems Data
University of South Carolina Scholar Commons Theses and Dissertations 8-9-2014 Non- and Semi-parametric Bayesian Inference with Recurrent Events and Coherent Systems Data A. K. M. Fazlur Rahman University
More informationFrailty Models and Copulas: Similarities and Differences
Frailty Models and Copulas: Similarities and Differences KLARA GOETHALS, PAUL JANSSEN & LUC DUCHATEAU Department of Physiology and Biometrics, Ghent University, Belgium; Center for Statistics, Hasselt
More informationSurvival Times (in months) Survival Times (in months) Relative Frequency. Relative Frequency
Smooth Goodness-of-Fit Tests in Hazard-Based Models by Edsel A. Pe~na Department of Statistics University of South Carolina at Columbia E-Mail: pena@stat.sc.edu Univ. of Virginia Colloquium Talk Department
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 informationGoodness-of-Fit Tests With Right-Censored Data by Edsel A. Pe~na Department of Statistics University of South Carolina Colloquium Talk August 31, 2 Research supported by an NIH Grant 1 1. Practical Problem
More informationSurvival Times (in months) Survival Times (in months) Relative Frequency. Relative Frequency
Smooth Goodness-of-Fit Tests in Hazard-Based Models by Edsel A. Pe~na Department of Statistics University of South Carolina at Columbia E-Mail: pena@stat.sc.edu Univ. of Georgia Colloquium Talk Department
More informationSTAT Sample Problem: General Asymptotic Results
STAT331 1-Sample Problem: General Asymptotic Results In this unit we will consider the 1-sample problem and prove the consistency and asymptotic normality of the Nelson-Aalen estimator of the cumulative
More informationST495: Survival Analysis: Maximum likelihood
ST495: Survival Analysis: Maximum likelihood Eric B. Laber Department of Statistics, North Carolina State University February 11, 2014 Everything is deception: seeking the minimum of illusion, keeping
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 01: Introduction and Overview
Introduction to Empirical Processes and Semiparametric Inference Lecture 01: Introduction and Overview Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations
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 informationConcepts and Tests for Trend in Recurrent Event Processes
JIRSS (2013) Vol. 12, No. 1, pp 35-69 Concepts and Tests for Trend in Recurrent Event Processes R. J. Cook, J. F. Lawless Department of Statistics and Actuarial Science, University of Waterloo, Ontario,
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 informationA Regression Model For Recurrent Events With Distribution Free Correlation Structure
A Regression Model For Recurrent Events With Distribution Free Correlation Structure J. Pénichoux(1), A. Latouche(2), T. Moreau(1) (1) INSERM U780 (2) Université de Versailles, EA2506 ISCB - 2009 - Prague
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 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 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 informationSurvival Analysis: Weeks 2-3. Lu Tian and Richard Olshen Stanford University
Survival Analysis: Weeks 2-3 Lu Tian and Richard Olshen Stanford University 2 Kaplan-Meier(KM) Estimator Nonparametric estimation of the survival function S(t) = pr(t > t) The nonparametric estimation
More informationlog T = β T Z + ɛ Zi Z(u; β) } dn i (ue βzi ) = 0,
Accelerated failure time model: log T = β T Z + ɛ β estimation: solve where S n ( β) = n i=1 { Zi Z(u; β) } dn i (ue βzi ) = 0, Z(u; β) = j Z j Y j (ue βz j) j Y j (ue βz j) How do we show the asymptotics
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 informationOptimal goodness-of-fit tests for recurrent event data
Lifetime Data Anal (211) 17:49 432 DOI 1.17/s1985-11-9193-1 Optimal goodness-of-fit tests for recurrent event data Russell S. Stocker IV Akim Adekpedjou Received: 13 October 29 / Accepted: 14 February
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 informationFrom semi- to non-parametric inference in general time scale models
From semi- to non-parametric inference in general time scale models Thierry DUCHESNE duchesne@matulavalca Département de mathématiques et de statistique Université Laval Québec, Québec, Canada Research
More informationCTDL-Positive Stable Frailty Model
CTDL-Positive Stable Frailty Model M. Blagojevic 1, G. MacKenzie 2 1 Department of Mathematics, Keele University, Staffordshire ST5 5BG,UK and 2 Centre of Biostatistics, University of Limerick, Ireland
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 25: Semiparametric Models
Introduction to Empirical Processes and Semiparametric Inference Lecture 25: Semiparametric Models Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations
More informationSemiparametric maximum likelihood estimation in normal transformation models for bivariate survival data
Biometrika (28), 95, 4,pp. 947 96 C 28 Biometrika Trust Printed in Great Britain doi: 1.193/biomet/asn49 Semiparametric maximum likelihood estimation in normal transformation models for bivariate survival
More informationBeyond GLM and likelihood
Stat 6620: Applied Linear Models Department of Statistics Western Michigan University Statistics curriculum Core knowledge (modeling and estimation) Math stat 1 (probability, distributions, convergence
More 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 informationSTAT 331. Accelerated Failure Time Models. Previously, we have focused on multiplicative intensity models, where
STAT 331 Accelerated Failure Time Models Previously, we have focused on multiplicative intensity models, where h t z) = h 0 t) g z). These can also be expressed as H t z) = H 0 t) g z) or S t z) = e Ht
More informationYou know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?
You know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?) I m not goin stop (What?) I m goin work harder (What?) Sir David
More informationHarvard University. Harvard University Biostatistics Working Paper Series
Harvard University Harvard University Biostatistics Working Paper Series Year 2008 Paper 85 Semiparametric Maximum Likelihood Estimation in Normal Transformation Models for Bivariate Survival Data Yi Li
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 informationOutline. Frailty modelling of Multivariate Survival Data. Clustered survival data. Clustered survival data
Outline Frailty modelling of Multivariate Survival Data Thomas Scheike ts@biostat.ku.dk Department of Biostatistics University of Copenhagen Marginal versus Frailty models. Two-stage frailty models: copula
More informationFRAILTY MODELS FOR MODELLING HETEROGENEITY
FRAILTY MODELS FOR MODELLING HETEROGENEITY By ULVIYA ABDULKARIMOVA, B.Sc. A Thesis Submitted to the School of Graduate Studies in Partial Fulfillment of the Requirements for the Degree Master of Science
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 informationReview. DS GA 1002 Statistical and Mathematical Models. Carlos Fernandez-Granda
Review DS GA 1002 Statistical and Mathematical Models http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall16 Carlos Fernandez-Granda Probability and statistics Probability: Framework for dealing with
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 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 informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS330 / MAS83 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-0 8 Parametric models 8. Introduction In the last few sections (the KM
More informationTests of independence for censored bivariate failure time data
Tests of independence for censored bivariate failure time data Abstract Bivariate failure time data is widely used in survival analysis, for example, in twins study. This article presents a class of χ
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 informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 02: Overview Continued
Introduction to Empirical Processes and Semiparametric Inference Lecture 02: Overview Continued Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations Research
More informationOn the generalized maximum likelihood estimator of survival function under Koziol Green model
On the generalized maximum likelihood estimator of survival function under Koziol Green model By: Haimeng Zhang, M. Bhaskara Rao, Rupa C. Mitra Zhang, H., Rao, M.B., and Mitra, R.C. (2006). On the generalized
More informationAnalysis of Time-to-Event Data: Chapter 4 - Parametric regression models
Analysis of Time-to-Event Data: Chapter 4 - Parametric regression models Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term 2018/19 1/25 Right censored
More informationGOODNESS-OF-FIT TESTS FOR ARCHIMEDEAN COPULA MODELS
Statistica Sinica 20 (2010), 441-453 GOODNESS-OF-FIT TESTS FOR ARCHIMEDEAN COPULA MODELS Antai Wang Georgetown University Medical Center Abstract: In this paper, we propose two tests for parametric models
More information,..., θ(2),..., θ(n)
Likelihoods for Multivariate Binary Data Log-Linear Model We have 2 n 1 distinct probabilities, but we wish to consider formulations that allow more parsimonious descriptions as a function of covariates.
More informationEfficiency of Profile/Partial Likelihood in the Cox Model
Efficiency of Profile/Partial Likelihood in the Cox Model Yuichi Hirose School of Mathematics, Statistics and Operations Research, Victoria University of Wellington, New Zealand Summary. This paper shows
More informationOutline. Frailty modelling of Multivariate Survival Data. Clustered survival data. Clustered survival data
Outline Frailty modelling of Multivariate Survival Data Thomas Scheike ts@biostat.ku.dk Department of Biostatistics University of Copenhagen Marginal versus Frailty models. Two-stage frailty models: copula
More informationA generalized Brown-Proschan model for preventive and corrective maintenance. Laurent DOYEN.
A generalized Brown-Proschan model for preventive and corrective maintenance laurent.doyen@iut2.upmf-grenoble.fr Jean Kuntzmann Laboratory Grenoble University France 1 of 28 I. Introduction The dependability
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 informationCan we do statistical inference in a non-asymptotic way? 1
Can we do statistical inference in a non-asymptotic way? 1 Guang Cheng 2 Statistics@Purdue www.science.purdue.edu/bigdata/ ONR Review Meeting@Duke Oct 11, 2017 1 Acknowledge NSF, ONR and Simons Foundation.
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 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 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 informationFigure 36: Respiratory infection versus time for the first 49 children.
y BINARY DATA MODELS We devote an entire chapter to binary data since such data are challenging, both in terms of modeling the dependence, and parameter interpretation. We again consider mixed effects
More informationBayesian estimation of the discrepancy with misspecified parametric models
Bayesian estimation of the discrepancy with misspecified parametric models Pierpaolo De Blasi University of Torino & Collegio Carlo Alberto Bayesian Nonparametrics workshop ICERM, 17-21 September 2012
More informationHacettepe Journal of Mathematics and Statistics Volume 45 (5) (2016), Abstract
Hacettepe Journal of Mathematics and Statistics Volume 45 (5) (2016), 1605 1620 Comparing of some estimation methods for parameters of the Marshall-Olkin generalized exponential distribution under progressive
More 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 informationQualifying Exam in Probability and Statistics. https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf
Part : Sample Problems for the Elementary Section of Qualifying Exam in Probability and Statistics https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf Part 2: Sample Problems for the Advanced Section
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 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 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 informationStatistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation
Statistics - Lecture One Charlotte Wickham wickham@stat.berkeley.edu http://www.stat.berkeley.edu/~wickham/ Outline 1. Basic ideas about estimation 2. Method of Moments 3. Maximum Likelihood 4. Confidence
More informationCHAPTER 1 A MAINTENANCE MODEL FOR COMPONENTS EXPOSED TO SEVERAL FAILURE MECHANISMS AND IMPERFECT REPAIR
CHAPTER 1 A MAINTENANCE MODEL FOR COMPONENTS EXPOSED TO SEVERAL FAILURE MECHANISMS AND IMPERFECT REPAIR Helge Langseth and Bo Henry Lindqvist Department of Mathematical Sciences Norwegian University of
More informationPropensity Score Weighting with Multilevel Data
Propensity Score Weighting with Multilevel Data Fan Li Department of Statistical Science Duke University October 25, 2012 Joint work with Alan Zaslavsky and Mary Beth Landrum Introduction In comparative
More information11 Survival Analysis and Empirical Likelihood
11 Survival Analysis and Empirical Likelihood The first paper of empirical likelihood is actually about confidence intervals with the Kaplan-Meier estimator (Thomas and Grunkmeier 1979), i.e. deals with
More informationStock Sampling with Interval-Censored Elapsed Duration: A Monte Carlo Analysis
Stock Sampling with Interval-Censored Elapsed Duration: A Monte Carlo Analysis Michael P. Babington and Javier Cano-Urbina August 31, 2018 Abstract Duration data obtained from a given stock of individuals
More informationUsing Estimating Equations for Spatially Correlated A
Using Estimating Equations for Spatially Correlated Areal Data December 8, 2009 Introduction GEEs Spatial Estimating Equations Implementation Simulation Conclusion Typical Problem Assess the relationship
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 informationON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION
ON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION ZHENLINYANGandRONNIET.C.LEE Department of Statistics and Applied Probability, National University of Singapore, 3 Science Drive 2, Singapore
More informationTABLE OF CONTENTS CHAPTER 1 COMBINATORIAL PROBABILITY 1
TABLE OF CONTENTS CHAPTER 1 COMBINATORIAL PROBABILITY 1 1.1 The Probability Model...1 1.2 Finite Discrete Models with Equally Likely Outcomes...5 1.2.1 Tree Diagrams...6 1.2.2 The Multiplication Principle...8
More informationThe consequences of misspecifying the random effects distribution when fitting generalized linear mixed models
The consequences of misspecifying the random effects distribution when fitting generalized linear mixed models John M. Neuhaus Charles E. McCulloch Division of Biostatistics University of California, San
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 informationModelling the risk process
Modelling the risk process Krzysztof Burnecki Hugo Steinhaus Center Wroc law University of Technology www.im.pwr.wroc.pl/ hugo Modelling the risk process 1 Risk process If (Ω, F, P) is a probability space
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 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 informationPrevious lecture. P-value based combination. Fixed vs random effects models. Meta vs. pooled- analysis. New random effects testing.
Previous lecture P-value based combination. Fixed vs random effects models. Meta vs. pooled- analysis. New random effects testing. Interaction Outline: Definition of interaction Additive versus multiplicative
More informationGOODNESS-OF-FIT TEST FOR RANDOMLY CENSORED DATA BASED ON MAXIMUM CORRELATION. Ewa Strzalkowska-Kominiak and Aurea Grané (1)
Working Paper 4-2 Statistics and Econometrics Series (4) July 24 Departamento de Estadística Universidad Carlos III de Madrid Calle Madrid, 26 2893 Getafe (Spain) Fax (34) 9 624-98-49 GOODNESS-OF-FIT TEST
More informationStatistical Analysis of Competing Risks With Missing Causes of Failure
Proceedings 59th ISI World Statistics Congress, 25-3 August 213, Hong Kong (Session STS9) p.1223 Statistical Analysis of Competing Risks With Missing Causes of Failure Isha Dewan 1,3 and Uttara V. Naik-Nimbalkar
More informationSession 9: Introduction to Sieve Analysis of Pathogen Sequences, for Assessing How VE Depends on Pathogen Genomics Part I
Session 9: Introduction to Sieve Analysis of Pathogen Sequences, for Assessing How VE Depends on Pathogen Genomics Part I Peter B Gilbert Vaccine and Infectious Disease Division, Fred Hutchinson Cancer
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 informationQuantile Regression for Residual Life and Empirical Likelihood
Quantile Regression for Residual Life and Empirical Likelihood Mai Zhou email: mai@ms.uky.edu Department of Statistics, University of Kentucky, Lexington, KY 40506-0027, USA Jong-Hyeon Jeong email: jeong@nsabp.pitt.edu
More information