A Tutorial on Multilevel Survival Analysis: Methods, Models and Applications. (24th March 2017)

• Record Type: Journal Article
• Authors: Austin, Peter C.
• Is Part Of: International statistical review. Volume 85:Number 2(2017)
• Page Start: 185
• Page End: 203
• Publication Date: 2017-03-24
• Subjects: Multilevel models -- hierarchical regression model -- survival analysis -- event history models -- Cox proportional hazards model -- clustered data -- health services research -- statistical software -- frailty models; Statistics -- Periodicals; Statistics -- Bibliography -- Periodicals; Statistics -- Bibliography; Statistics -- Periodicals; Statistique; Statistique -- Périodiques; Statistique -- Bibliographie -- Périodiques; Statistique -- Étude et enseignement -- Périodiques; Statistique -- Étude et enseignement -- Bibliographie -- Périodiques; Ressource Internet (Descripteur de forme); Périodique électronique (Descripteur de forme); 519.2
• DOI: 10.1111/insr.12214
• Languages: English
• ISSNs: 0306-7734

Abstract

Summary: Data that have a multilevel structure occur frequently across a range of disciplines, including epidemiology, health services research, public health, education and sociology. We describe three families of regression models for the analysis of multilevel survival data. First, Cox proportional hazards models with mixed effects incorporate cluster‐specific random effects that modify the baseline hazard function. Second, piecewise exponential survival models partition the duration of follow‐up into mutually exclusive intervals and fit a model that assumes that the hazard function is constant within each interval. This is equivalent to a Poisson regression model that incorporates the duration of exposure within each interval. By incorporating cluster‐specific random effects, generalised linear mixed models can be used to analyse these data. Third, after partitioning the duration of follow‐up into mutually exclusive intervals, one can use discrete time survival models that use a complementary log–log generalised linear model to model the occurrence of the outcome of interest within each interval. Random effects can be incorporated to account for within‐cluster homogeneity in outcomes. We illustrate the application of these methods using data consisting of patients hospitalised with a heart attack. We illustrate the application of these methods using three statistical programming languages (R, SAS and Stata).