Outline. Frailty modelling of Multivariate Survival Data. Clustered survival data. Clustered survival data
|
|
- Jesse Mitchell
- 5 years ago
- Views:
Transcription
1 Outline Frailty modelling of Multivariate Survival Data Thomas Scheike Department of Biostatistics University of Copenhagen Marginal versus Frailty models. Two-stage frailty models: copula modelling. Example and Exercises. November 212 Clustered survival data For k = 1,..., K, i = 1,..., n, let T ik and C ik be the failure and censoring times for the ith individual in the kth cluster and let X ik (t) be a p-vector of covariates. Put T k = ( T 1k,..., T nk ), C k = (C 1k,..., C nk ), X k (t) = (X 1k (t),..., X nk (t)). We assume that ( T k, C k, X k ( )), k = 1,..., K are independent and identically distributed variables and these variables follow the model described in the following. The right-censored failure time is denoted T ik = T ik C ik and as usual we let Y ik (t) = 1(T ik t) Clustered survival data Two basic approaches : 1) Frailty models: given random effects "Z" the survival times are independent with hazards λ (t)z exp(x T i β) 2) Marginal approach, given covariates the marginal intensity is (need intensity to construct likelihood of data) Y i (t)λ (t) exp(x T i β) This gives population regression parameters, but one can also characterize the dependende using copula (or underlying frailty distribution). Assume for example that given frailty Z, the intensities are Zh(Λ, exp(xi T β)) for a specific choice of h. More on this later.
2 Frailty models Characterizing dependence ( T k, C k, X k ( ), V k ), k = 1,, K are assumed to be i.i.d. variables. Censoring, conditional on V k and covariates, is assumed to be independent and noninformative on V k. When the frailty variable V k is gamma-distributed with mean one and variance 1 we get the. The shared frailty model is specified solely with respect to the (unobserved) filtration H t = k Hk t, where H k t = σ{n ik (s), Y ik (s), X ik (s), V k : i = 1, n, s t}, and often it is assumed to be a proportional hazards model We shall talk about estimation for the random effects models, but first some general remarks on dependence measures. Given two (T 1, T 2 ) survival times, we can compute Pearson correlation, Spearman correlation, and the Kendall s tau. Different ways of interpreting the degree of association. Spearman correlation depends on aspects of the bivariate distribution, but not the marginals. The same property is kept by Kendall s tau. These measures are rank-based. λ H ik (t) = Y ik(t)v k λ (t) exp(xik T (t)β). (1) Kendall s tau Dependence Kendall s tau is defined given an i.i.d. copy, and the degree of difference between such two i.i.d. copies E(I ((T 1 T 1 )(T 2 T 2 ) > )) E(I ((T 1 T 1 )(T 2 T 2 ) < )) which is just the concordance probability - discordance probability Given a bivariate distrubtion written on copula-form F (s, t) = C(F 1 (s), F 2 (t)) where C is the copula specifying bivariate distributions on [, 1] 2, and F j is the jth marginal. Kendall s tau is computed as τ = C(u, v)c(du, dv) Y ik (t)v k λ (t) exp(xik T (t)β). (2) with V k i.i.d. Γ with mean 1, and variance 1/, then τ = in the sense that with fixed covariates we see this Kendall s tau.
3 The frailty model When V k i.i.d. Γ with mean 1, and variance 1/, then the observed intensity given covariates become N k. Y ik (t) i exp(x ik T (t)β)λ (T ik ) λ (t) exp(xik T (t)β). (3) And the marginal intensity becomes: N ik (t ) Y ik (t) exp(xik T (t)β)λ (T ik ) λ (t) exp(xik T (t)β). (4) What is the observed hazard, and how is it different from marginal intensities? As a part of the next exercise we shall also derive these expressions. Frailty models Estimation : NPMLE for Frailty model (V k dλ (s) exp(xki T β)) N ki (s) exp ( k i s t Th V k exp(x T ki β)dλ (s) (5) and with V k Gamma distributed with mean 1, and variance 1/, giving the additional term f (V k ) (6) k For fixed write up the EM-algorithm for this model. How do we get standard errors for this model. To compute E-step we need E(Vk Data), and this conditional mean (and distribution) can be found from full likelihood. Recall that a Gamma distribution with parameters λ, α (Γ(λ, α)) has mean λα and variance λα 2 and density (Γ(λ)α λ ) 1 x λ 1 exp( x/α). Write up full EM-algorithm for the model where is not known. Frailty models, Attenuation > fit<-coxph(surv(time,status)~adult*trt+frailty(id), data=diabetes) > fit Call: coxph(formula = Surv(time, status) ~ adult * trt + frailty(id), data = diabetes) coef se(coef) se2 Chisq DF p adult trt frailty(id) adult:trt Iterations: 6 outer, 25 Newton-Raphson Variance of random effect=.926 I-likelihood = -847 Fit the model yoursel, and try to simplify the regression part, is the frailty variance significant, and what is it. What is the related Kendall s tau? Compare with the marginal esimates. How do the regression effects differ? If dim(x)=1 and gamma frailty with variance and mean 1 then the relative risk is Y (t)λ (t) exp (X T 1 β) 1 + Λ (t) ( exp (β) 1 + Λ (t) 1 + exp (β)λ (t) ), (7) where Λ (t) = t λ (s) ds. exp (β) at time t = and tends to 1 as t >. V positive stable with Laplace transform ϕ (t) = exp ( t ), < 1, then the marginal intensity is Y (t)λ (t)λ (t) 1 exp (X T β),
4 Frailty models, Attenuation Frailty models Assume that V is a positive stochastic variable with Laplace transform ϕ (t), that the covariate is one-dimensional, and that β >. The relative risk in the marginal model is exp (β) (D log ϕ )(exp ((X + 1)β)Λ (t)) (D log ϕ )(exp (X β)λ (t)) = exp (β)k(t), and we see that k(t) 1 if and only if (D log ϕ )(exp ((X + 1)β)Λ (t)) (D log ϕ )(exp (X β)λ (t)). The latter inequality holds if log (ϕ ) is convex, which is the case since D 2 log (ϕ )(t) = E(V 2 h(t, V )) E(Vh(t, V )) 2 with h(t, V ) = exp ( tv )/E(exp ( tv )). Two sets of frailty models quite similar: Frailty model Frailty model has subject specific regression effects, that relies on choice of frailty (ofcourse). Frailty parameter model can in principle be identified solely on univariate data because marginal depends on. Frailty parameter primarily identified from correlation in practice. Frailty parameter has Kendall s tau interpretation. Frailty model more easy to extend to various other settings. Two-stage model, Later today Marginals fixed and regression parameters gives population effects. Frailty accounts solely for dependence. Frailty parameter has Kendall s tau interpretation. Two-stage modelling Two-stage modelling To be explicit, assume that the marginal and conditional intensities are λ F ik ik (t) = λ ik(t), λ H ik (t) = V kλ ik (t), where we assume that λ ik (t) is predictable with respect to the marginal filtration. One may show that the relationship between the above two intensities is Establish the connection between the the two sets of intensities: Hint : compute the marginal survival for both models! λ ik (t) = Y ik (t)( λ ik (t))(d log ϕ )( t λ ik (t) = Y ik(t)( λ ik (t)) exp ( t λ ik (s) ds), t λ ik (s) ds)(dϕ 1 )(exp ( λ ik (s) ds)),
5 Marginal models Right-censored failure times T ik = T ik C ik, Y ik (t) = 1(T ik t), N ik (t) = 1(T ik t, T ik = T ik ) Marginal (intensity) model is a model given Cox model: F ik t = σ{n ik (s), Y ik (s), X ik (s) : s t}, (8) λ F ik ik (t) = Y ik(t)λ (t) exp (Xik T (t)β). (9) It is important to note that (9) is not the intensity with respect to the observed filtration F t = Ft k, (1) k where F k t = σ{n ik (s), Y ik (s), X ik (s) : i = 1, n, s t} is the information generated by observing all the individuals in the kth cluster. Now, with respect to (unobserved) filtration where is H t = k H k t, (11) H k t = σ{n ik (s), Y ik (s), X ik (s), V k : i = 1, n, s t}, λ H ik (t) = V kλ ik (t,, λ ( )), (12) referred to as the, and so that the marginal intensities are on Cox form λ F ik ik (t) = Y ik(t)λ (t) exp (Xik T (t)β), (13) Characterizing dependence Copula models P( T 1 > t 1,..., T n > t n ) = C (S 1 (t 1 ),..., S n (t n )), where S j, j = 1,..., n, denotes the marginal survivor functions. All multivariate distribution has this form. Archimedean copula model family, C (u 1,..., u n ) = ϕ (ϕ 1 (u 1) + + ϕ 1 (u n)) for some non-negative convex decreasing function ϕ with ϕ () = 1. Then λ ik is λ ik (t,, λ ( )) = Y ik (t)λ (t) exp(xik T (t)β) t exp( 1 exp(xik T (s)βt )λ (s)ds), It is of interest to find the intensities with respect to the observed filtration F t given in (1). It can be shown that these are where λ F ik (t) = Y ik(t)λ (t) exp(x T ik (t)β)f ik(t), (14) f ik (t) = ( + N k (t ) n t f k (t) = 1 + (exp( 1 j=1 )(exp( 1 t λ (s) exp(x T f k (t) ik (s)β)ds) ), Y jk (s)λ (s) exp(xjk T (s)β)ds) 1).
6 The observed (partial) log-likelihood function is K k=1 log(1 + N k(t ) )dn k (t) + K [ + N k (τ) ] n log(1 + 1 k=1 i=1 K n k=1 i=1 log(y ik (t) λ ik (t))dn ik (t) Y ik (t) λ ik (t)dt), (15) Terms depending on in (15) gives 1 K ( K k=1 log(1 + 1 N k (t ))dn k (t) + K k=1 i=1 n 1 N ik (τ)h ik K ) ( + N k (τ)) log(r k ()), (16) k=1 where λ ik (t) = λ (t)e X T ik (t)β exp( 1 t e X T ik (s)β λ (s)ds). H ik = Y ik (t)e X T ik (t)β dλ (t), R k () = 1 + n (exp( 1 H ik ) 1). i=1 Now, by replacing H ik with Ĥ ik = Y ik (t) exp(x T ik (t) ˆβ I )d ˆΛ I (t) in (16), we obtain the pseudo log likelihood for, and maximizing this function in gives the two-stage estimator of. Under some regularity conditions, [?] showed consistency and asymptotic normality of this estimator. This is equivalent to fitting directly the MLE for the Clayton-Oakes copula form C (u 1,..., u n ) = ϕ (ϕ 1 (u 1) + + ϕ 1 (u n)) Derive this Copula form for the Clayton-Oakes Gamma frailty model with Cox marginals. Write up the likelihood for censored bivariate survival data based on this represenation.
7 Two-stage modelling Two-stage modelling data(diabetes) # Marginal Cox model with treat as covariate marg <- cox.aalen(surv(time,status) ~prop(treat)+cluster(id),data=diabetes) fit<-two.stage(marg,data=diabetes,theta=1) summary(fit) Dependence parameter for Clayton-Oakes-Glidden model Variance SE z P-val Kendall s tau intercept Marginal Cox-Aalen model fit Proportional Cox terms : Coef. SE Robust SE D2log(L)^-1 z P-val prop(treat) e-6 # Stratification after adult theta.des<-model.matrix(~-1+factor(adult),diabetes); fit.s2<-two.stage(marg,data=diabetes,nit=4,theta=1,theta.des=theta.des) summary(fit.s2) Dependence parameter for Clayton-Oakes-Glidden model Variance SE z P-val Kendall s tau factor(adult) factor(adult) Marginal Cox-Aalen model fit Proportional Cox terms : Coef. SE Robust SE D2log(L)^-1 z P-val prop(treat) e-6 Two-stage modelling Likelihoods # test for same variance among the two strata theta.des<-model.matrix(~factor(adult),diabetes); fit.s3<-two.stage(marg,data=diabetes,nit=4,theta=1,theta.des=theta.des) summary(fit.s3) Dependence parameter for Clayton-Oakes-Glidden model Variance SE z P-val Kendall s tau (Intercept) factor(adult) Marginal Cox-Aalen model fit Proportional Cox terms : Coef. SE Robust SE D2log(L)^-1 z P-val prop(treat) e-6 Generally given af bivariate survival distribution G(x, y) and censored bivariate survival data (X i, δ i, Y i, ν i ), the likelihood is written as G(X i, Y i ) (1 δ i )(1 ν i ) (D x G)(X i, Y i ) δ i (1 ν i ) (17) (D y G)(X i, Y i ) (1 δ i )ν i (D X D y G)(X i, Y i ) δ i ν i (18) # to fit model without covariates, beta.fixed=1, but still need prop term! fit<-two.stage(marg,data=diabetes,theta=.95,detail=,beta.fixed=1) summary(fit) Dependence parameter for Clayton-Oakes-Glidden model Variance SE z P-val Kendall s tau intercept
8 Goodness of fit of Gamma Frailty Different frailt distribution have different consequences for the type of dependence. GOF One idea (Shih, 1998, Shih and Louis 1995) Compare likelihoods from different models. Score process. Under Gamma-distribution we earlier computed E(Z k History) = γ k (t) = N k i exp(x T ik (t)β)λ (T ik ) Idea is to look at γ k (t) that has mean 1 under the gamma distribution. Let W (t) = K 1/2 k (γ k (t) 1) that is asymptotically Gaussian with a variance we can also estiamte under the null. This can be resampled and a supremum test can be constructed. Frailty models Considering the twin.csv data of menarche ages for pairs of twins. Look at the two.stage help pages to see different examples of its use. First estimate the marginal effect of cohort, and assess wheter or not the zygosity affects the marginal models. Fit a frailty model for the overall data, to asses the effect of cohort. Use the phmm and coxph program and compare the estimates of the cumulative baseline. Compare also the baseline with that from the marginal model for mono and dyzygotic twins. Fit separate models for monozygotic and dizygotic twins. Estimate also the marginal baselines using this model, and compare them. Is there a genetic effect on the timing of menarche? Report a Kendall s tau for this data for mono- and dizygotic twins and overall. Is the dependence different for mono- versus dizygotic twins. How can you make survival predictions for the monozygotic and dyzygotic twins? Consider the (simulated) twin data of menarche age for twins. Fit marginal models for the mono- and the dizygotic twins. Report a Kendall s tau. How does this Kendall s tau compare with that from the standard frailty model? Validate the fit of the model, by validating the marginal model. What about the frailty part? Compare formally the dependence between mono- and dyzygotic twins.
Outline. 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 informationGoodness-Of-Fit for Cox s Regression Model. Extensions of Cox s Regression Model. Survival Analysis Fall 2004, Copenhagen
Outline Cox s proportional hazards model. Goodness-of-fit tools More flexible models R-package timereg Forthcoming book, Martinussen and Scheike. 2/38 University of Copenhagen http://www.biostat.ku.dk
More 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 informationCox s proportional hazards/regression model - model assessment
Cox s proportional hazards/regression model - model assessment Rasmus Waagepetersen September 27, 2017 Topics: Plots based on estimated cumulative hazards Cox-Snell residuals: overall check of fit Martingale
More 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 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 informationFaculty of Health Sciences. Cox regression. Torben Martinussen. Department of Biostatistics University of Copenhagen. 20. september 2012 Slide 1/51
Faculty of Health Sciences Cox regression Torben Martinussen Department of Biostatistics University of Copenhagen 2. september 212 Slide 1/51 Survival analysis Standard setup for right-censored survival
More informationSSUI: Presentation Hints 2 My Perspective Software Examples Reliability Areas that need work
SSUI: Presentation Hints 1 Comparing Marginal and Random Eects (Frailty) Models Terry M. Therneau Mayo Clinic April 1998 SSUI: Presentation Hints 2 My Perspective Software Examples Reliability Areas that
More 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 informationLecture 12. Multivariate Survival Data Statistics Survival Analysis. Presented March 8, 2016
Statistics 255 - Survival Analysis Presented March 8, 2016 Dan Gillen Department of Statistics University of California, Irvine 12.1 Examples Clustered or correlated survival times Disease onset in family
More 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 informationProportional hazards model for matched failure time data
Mathematical Statistics Stockholm University Proportional hazards model for matched failure time data Johan Zetterqvist Examensarbete 2013:1 Postal address: Mathematical Statistics Dept. of Mathematics
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 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 informationMultivariate Survival Data With Censoring.
1 Multivariate Survival Data With Censoring. Shulamith Gross and Catherine Huber-Carol Baruch College of the City University of New York, Dept of Statistics and CIS, Box 11-220, 1 Baruch way, 10010 NY.
More 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 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 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 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 informationTime-dependent coefficients
Time-dependent coefficients Patrick Breheny December 1 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/20 Introduction As we discussed previously, stratification allows one to handle variables that
More informationDependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline.
Practitioner Course: Portfolio Optimization September 10, 2008 Before we define dependence, it is useful to define Random variables X and Y are independent iff For all x, y. In particular, F (X,Y ) (x,
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 informationPairwise dependence diagnostics for clustered failure-time data
Biometrika Advance Access published May 13, 27 Biometrika (27), pp. 1 15 27 Biometrika Trust Printed in Great Britain doi:1.193/biomet/asm24 Pairwise dependence diagnostics for clustered failure-time data
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 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 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 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 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 informationModelling Dependence with Copulas and Applications to Risk Management. Filip Lindskog, RiskLab, ETH Zürich
Modelling Dependence with Copulas and Applications to Risk Management Filip Lindskog, RiskLab, ETH Zürich 02-07-2000 Home page: http://www.math.ethz.ch/ lindskog E-mail: lindskog@math.ethz.ch RiskLab:
More informationLecture 8 Stat D. Gillen
Statistics 255 - Survival Analysis Presented February 23, 2016 Dan Gillen Department of Statistics University of California, Irvine 8.1 Example of two ways to stratify Suppose a confounder C has 3 levels
More 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 informationEstimation and Goodness of Fit for Multivariate Survival Models Based on Copulas
Estimation and Goodness of Fit for Multivariate Survival Models Based on Copulas by Yildiz Elif Yilmaz A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the
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 informationASYMPTOTIC PROPERTIES AND EMPIRICAL EVALUATION OF THE NPMLE IN THE PROPORTIONAL HAZARDS MIXED-EFFECTS MODEL
Statistica Sinica 19 (2009), 997-1011 ASYMPTOTIC PROPERTIES AND EMPIRICAL EVALUATION OF THE NPMLE IN THE PROPORTIONAL HAZARDS MIXED-EFFECTS MODEL Anthony Gamst, Michael Donohue and Ronghui Xu University
More informationFinancial Econometrics and Volatility Models Copulas
Financial Econometrics and Volatility Models Copulas Eric Zivot Updated: May 10, 2010 Reading MFTS, chapter 19 FMUND, chapters 6 and 7 Introduction Capturing co-movement between financial asset returns
More informationLecture 5 Models and methods for recurrent event data
Lecture 5 Models and methods for recurrent event data Recurrent and multiple events are commonly encountered in longitudinal studies. In this chapter we consider ordered recurrent and multiple events.
More informationSurvival analysis in R
Survival analysis in R Niels Richard Hansen This note describes a few elementary aspects of practical analysis of survival data in R. For further information we refer to the book Introductory Statistics
More 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 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 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 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 informationA Goodness-of-fit Test for Semi-parametric Copula Models of Right-Censored Bivariate Survival Times
A Goodness-of-fit Test for Semi-parametric Copula Models of Right-Censored Bivariate Survival Times by Moyan Mei B.Sc. (Honors), Dalhousie University, 2014 Project Submitted in Partial Fulfillment of the
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 informationPhD course in Advanced survival analysis. One-sample tests. Properties. Idea: (ABGK, sect. V.1.1) Counting process N(t)
PhD course in Advanced survival analysis. (ABGK, sect. V.1.1) One-sample tests. Counting process N(t) Non-parametric hypothesis tests. Parametric models. Intensity process λ(t) = α(t)y (t) satisfying Aalen
More informationSurvival Analysis using Bivariate Archimedean Copulas. Krishnendu Chandra
Survival Analysis using Bivariate Archimedean Copulas Krishnendu Chandra Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy under the Executive Committee of the
More informationModelling Dependent Credit Risks
Modelling Dependent Credit Risks Filip Lindskog, RiskLab, ETH Zürich 30 November 2000 Home page:http://www.math.ethz.ch/ lindskog E-mail:lindskog@math.ethz.ch RiskLab:http://www.risklab.ch Modelling Dependent
More informationA Measure of Association for Bivariate Frailty Distributions
journal of multivariate analysis 56, 6074 (996) article no. 0004 A Measure of Association for Bivariate Frailty Distributions Amita K. Manatunga Emory University and David Oakes University of Rochester
More informationSAMPLE SIZE ESTIMATION FOR SURVIVAL OUTCOMES IN CLUSTER-RANDOMIZED STUDIES WITH SMALL CLUSTER SIZES BIOMETRICS (JUNE 2000)
SAMPLE SIZE ESTIMATION FOR SURVIVAL OUTCOMES IN CLUSTER-RANDOMIZED STUDIES WITH SMALL CLUSTER SIZES BIOMETRICS (JUNE 2000) AMITA K. MANATUNGA THE ROLLINS SCHOOL OF PUBLIC HEALTH OF EMORY UNIVERSITY SHANDE
More informationMaximum likelihood estimation for Cox s regression model under nested case-control sampling
Biostatistics (2004), 5, 2,pp. 193 206 Printed in Great Britain Maximum likelihood estimation for Cox s regression model under nested case-control sampling THOMAS H. SCHEIKE Department of Biostatistics,
More informationNonparametric estimation of linear functionals of a multivariate distribution under multivariate censoring with applications.
Chapter 1 Nonparametric estimation of linear functionals of a multivariate distribution under multivariate censoring with applications. 1.1. Introduction The study of dependent durations arises in many
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
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 informationComposite likelihood and two-stage estimation in family studies
Biostatistics (2004), 5, 1,pp. 15 30 Printed in Great Britain Composite likelihood and two-stage estimation in family studies ELISABETH WREFORD ANDERSEN The Danish Epidemiology Science Centre, Statens
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 informationCopula modeling for discrete data
Copula modeling for discrete data Christian Genest & Johanna G. Nešlehová in collaboration with Bruno Rémillard McGill University and HEC Montréal ROBUST, September 11, 2016 Main question Suppose (X 1,
More informationSurvival analysis in R
Survival analysis in R Niels Richard Hansen This note describes a few elementary aspects of practical analysis of survival data in R. For further information we refer to the book Introductory Statistics
More informationDependence. MFM Practitioner Module: Risk & Asset Allocation. John Dodson. September 11, Dependence. John Dodson. Outline.
MFM Practitioner Module: Risk & Asset Allocation September 11, 2013 Before we define dependence, it is useful to define Random variables X and Y are independent iff For all x, y. In particular, F (X,Y
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 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 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 informationST745: Survival Analysis: Cox-PH!
ST745: Survival Analysis: Cox-PH! Eric B. Laber Department of Statistics, North Carolina State University April 20, 2015 Rien n est plus dangereux qu une idee, quand on n a qu une idee. (Nothing is more
More informationCox s proportional hazards model and Cox s partial likelihood
Cox s proportional hazards model and Cox s partial likelihood Rasmus Waagepetersen October 12, 2018 1 / 27 Non-parametric vs. parametric Suppose we want to estimate unknown function, e.g. survival function.
More 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 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 informationRelative-risk regression and model diagnostics. 16 November, 2015
Relative-risk regression and model diagnostics 16 November, 2015 Relative risk regression More general multiplicative intensity model: Intensity for individual i at time t is i(t) =Y i (t)r(x i, ; t) 0
More information(a) (3 points) Construct a 95% confidence interval for β 2 in Equation 1.
Problem 1 (21 points) An economist runs the regression y i = β 0 + x 1i β 1 + x 2i β 2 + x 3i β 3 + ε i (1) The results are summarized in the following table: Equation 1. Variable Coefficient Std. Error
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 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 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 information1 Introduction. 2 Residuals in PH model
Supplementary Material for Diagnostic Plotting Methods for Proportional Hazards Models With Time-dependent Covariates or Time-varying Regression Coefficients BY QIQING YU, JUNYI DONG Department of Mathematical
More informationA copula goodness-of-t approach. conditional probability integral transform. Daniel Berg 1 Henrik Bakken 2
based on the conditional probability integral transform Daniel Berg 1 Henrik Bakken 2 1 Norwegian Computing Center (NR) & University of Oslo (UiO) 2 Norwegian University of Science and Technology (NTNU)
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 informationCox regression: Estimation
Cox regression: Estimation Patrick Breheny October 27 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/19 Introduction The Cox Partial Likelihood In our last lecture, we introduced the Cox partial
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 informationFrailty Modeling for Spatially Correlated Survival Data, with Application to Infant Mortality in Minnesota By: Sudipto Banerjee, Mela. P.
Frailty Modeling for Spatially Correlated Survival Data, with Application to Infant Mortality in Minnesota By: Sudipto Banerjee, Melanie M. Wall, Bradley P. Carlin November 24, 2014 Outlines of the talk
More informationReduced-rank hazard regression
Chapter 2 Reduced-rank hazard regression Abstract The Cox proportional hazards model is the most common method to analyze survival data. However, the proportional hazards assumption might not hold. The
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 informationMultilevel Models in Matrix Form. Lecture 7 July 27, 2011 Advanced Multivariate Statistical Methods ICPSR Summer Session #2
Multilevel Models in Matrix Form Lecture 7 July 27, 2011 Advanced Multivariate Statistical Methods ICPSR Summer Session #2 Today s Lecture Linear models from a matrix perspective An example of how to do
More informationOn a connection between the Bradley-Terry model and the Cox proportional hazards model
On a connection between the Bradley-Terry model and the Cox proportional hazards model Yuhua Su and Mai Zhou Department of Statistics University of Kentucky Lexington, KY 40506-0027, U.S.A. SUMMARY This
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 informationOutline. Cox's regression model Goodness-of-t methods. Cox's proportional hazards model: Survival analysis
Outline Cox's regression model Goodness-of-t methods Giuliana Cortese gico@biostat.ku.dk Advanced Survival Analysis 21 Copenhagen Cox's proportional hazards model: assumptions Goodness-of-t for Cox's model:
More informationOn the Breslow estimator
Lifetime Data Anal (27) 13:471 48 DOI 1.17/s1985-7-948-y On the Breslow estimator D. Y. Lin Received: 5 April 27 / Accepted: 16 July 27 / Published online: 2 September 27 Springer Science+Business Media,
More 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 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 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 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 informationfrailtyem: An R Package for Estimating Semiparametric Shared Frailty Models
frailtyem: An R Package for Estimating Semiparametric Shared Frailty Models Theodor Adrian Balan Leiden University Medical Center Hein Putter Leiden University Medical Center Abstract When analyzing correlated
More informationTwo-level lognormal frailty model and competing risks model with missing cause of failure
University of Iowa Iowa Research Online Theses and Dissertations Spring 2012 Two-level lognormal frailty model and competing risks model with missing cause of failure Xiongwen Tang University of Iowa Copyright
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 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 informationEnsemble estimation and variable selection with semiparametric regression models
Ensemble estimation and variable selection with semiparametric regression models Sunyoung Shin Department of Mathematical Sciences University of Texas at Dallas Joint work with Jason Fine, Yufeng Liu,
More informationSurvival Analysis: Counting Process and Martingale. Lu Tian and Richard Olshen Stanford University
Survival Analysis: Counting Process and Martingale Lu Tian and Richard Olshen Stanford University 1 Lebesgue-Stieltjes Integrals G( ) is a right-continuous step function having jumps at x 1, x 2,.. b f(x)dg(x)
More informationWeb-based Supplementary Material for A Two-Part Joint. Model for the Analysis of Survival and Longitudinal Binary. Data with excess Zeros
Web-based Supplementary Material for A Two-Part Joint Model for the Analysis of Survival and Longitudinal Binary Data with excess Zeros Dimitris Rizopoulos, 1 Geert Verbeke, 1 Emmanuel Lesaffre 1 and Yves
More informationSTATISTICAL ANALYSIS OF MULTIVARIATE INTERVAL-CENSORED FAILURE TIME DATA
STATISTICAL ANALYSIS OF MULTIVARIATE INTERVAL-CENSORED FAILURE TIME DATA A Dissertation Presented to the Faculty of the Graduate School University of Missouri-Columbia In Partial Fulfillment Of the Requirements
More informationSEMIPARAMETRIC REGRESSION WITH TIME-DEPENDENT COEFFICIENTS FOR FAILURE TIME DATA ANALYSIS
Statistica Sinica 2 (21), 853-869 SEMIPARAMETRIC REGRESSION WITH TIME-DEPENDENT COEFFICIENTS FOR FAILURE TIME DATA ANALYSIS Zhangsheng Yu and Xihong Lin Indiana University and Harvard School of Public
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 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 informationTied survival times; estimation of survival probabilities
Tied survival times; estimation of survival probabilities Patrick Breheny November 5 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/22 Introduction Tied survival times Introduction Breslow approximation
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 informationLecture 2: Martingale theory for univariate survival analysis
Lecture 2: Martingale theory for univariate survival analysis In this lecture T is assumed to be a continuous failure time. A core question in this lecture is how to develop asymptotic properties when
More information