A Tutorial on Multilevel Survival Analysis: Methods, Models and Applications

Peter C Austin, Peter C Austin

Abstract

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).

Keywords: Cox proportional hazards model; Multilevel models; clustered data; event history models; frailty models; health services research; hierarchical regression model; statistical software; survival analysis.

Figures

Figure 1
Figure 1
Variation in hospital-specific hazards and survival ( frailty model). SD, standard deviation. [Colour figure can be viewed at wileyonlinelibrary.com]
Figure 2
Figure 2
Variation in hazard functions across hospitals (piecewise exponential model). SD, standard deviation. [Colour figure can be viewed at wileyonlinelibrary.com]
Statistical software output 1
Statistical software output 1
SAS output for Cox frailty survival model (log-normal frailty distribution)
Statistical software output 2
Statistical software output 2
Stata output for PWE survival model
Statistical software output 3
Statistical software output 3
R output for discrete time mixed effects survival model
Statistical software output 4
Statistical software output 4
SAS output for discrete time mixed effects survival model with random intercept and random effect for cardiogenic shock

References

    1. Aalen OO, Borgan O, Gjessing HK. Survival and Event History Analysis. New York, NY: Springer; 2008.
    1. Aitkin M, Laird N, Francis B. A reanalysis of the Stanford heart transplant data. J Am Stat Assoc. 1983;78:264–274.
    1. Allison PD. Survival Analysis using SAS®: A Practical Guide. 2. Cary NC: SAS Institute; 2010.
    1. Altman DG, Andersen PK. Calculating the number needed to treat for trials where the outcome is time to an event. BMJ. 1999;319:1492–1495.
    1. Austin PC, Manca A, Zwarenstein M, Juurlink DN, Stanbrook MB. A substantial and confusing variation exists in handling of baseline covariates in randomized controlled trials: a review of trials published in leading medical journals. J Clin Epidemiol. 2010;63:142–153.
    1. Austin PC, Wagner P, Merlo J. The Median Hazard Ratio: A useful measure of variance and general contextual effects in multilevel survival analysis. Stat Med. 2017;36(6):928–938.
    1. Barber JS, Murphy SA, Axinn WG, Maples J. Discrete-time multilevel hazard analysis. Sociolo Method. 2000;30:201–235.
    1. Breslow N. Covariance analysis of censored survival data. Biometrics. 1974;30:89–99.
    1. Cox D, Oakes D. Analysis of Survival Data. London: Chapman & Hall; 1984.
    1. Crowther MJ, Look MP, Riley RD. Multilevel mixed effects parametric survival models using adaptive Gauss-Hermite quadrature with application to recurrent events and individual participant data meta-analysis. Stat Med. 2014;33:3844–3858.
    1. Crowther MJ, Riley RD, Staessen JA, Wang J, Gueyffier F, Lambert PC. Individual patient data meta-analysis of survival data using Poisson regression models. BMC Med Res Methodol. 2012;12:34.
    1. Duchateau L, Janssen P. The Frailty Model. New York, NY: Springer; 2008.
    1. Gail MH, Wieand S, Piantadosi S. Biased estimates of treatment effect in randomized experiments with nonlinear regressions and omitted covariates. Biometrika. 1984;7:431–444.
    1. Goldstein H. Multilevel Statistical Models. 4. West Sussex: John Wiley & Sons Ltd; 2011.
    1. Hougaard P. Analysis of Multivariate Survival Data. New York: Springer-Verlag; 2000.
    1. Hox JJ, Roberts JK. Handbook of Advanced Multilevel Analysis. New York: Routledge; 2011.
    1. Kalbfleisch JD, Prentice RL. The Statistical Analysis of Failure Time Data. 2. New York: John Wiley and Sons; 2002.
    1. Klein JP, Moeschberger ML. Survival Analysis: Techniques for Censored and Truncated Data. New York, NY: Springer-Verlag; 1997.
    1. Laird N, Olivier D. Covariance analysis of censored survival data using log-linear analysis techniques. J Am Stat Assoc. 1981;76:231–240.
    1. Lawless JF. Statistical Models and Methods for Lifetime Data. New York: John Wiley & Sons; 1982.
    1. Lin DY, Wei LJ. The robust inference for the proportional hazards model. J Am Stat Assoc. 1989;84:1074–1078.
    1. Mills M. Introducing Survival and Event History Analysis. Thousand Oaks, CA: Sage; 2011.
    1. Rabe-Hesketh S, Skrondal A. Multilevel and Longitudinal Modeling Using Stata, Volume 1: Continuous Responses. 3. Vol. 1. Stata Press; 2012a.
    1. Rabe-Hesketh S, Skrondal A. Multilevel and Longitudinal Modeling Using Stata, Volume 2: Categorical Responses, Counts, and Survival. 3. Vol. 2. Stata Press; 2012b.
    1. Raudenbush SW, Bryk AS. Hierarchical Linear Models: Applications and Data Analysis Methods. 2. Thousand Oaks: Sage Publications; 2002.
    1. Rodriguez G. Multilevel generalized linear models. In: de Leeuw J, Meijer E, editors. Handbook of Multilevel Analysis. New York: Springer; 2008. pp. 335–376.
    1. Rondeau V, Mazroui Y, Gonzalez JR. frailtypack: an R package for the analysis of correlated survival data with frailty models using penalized likelihood estimation or parametrical estimation. J Stat Softw. 2012;47(4):1–28.
    1. Singer JD, Willett JB. Applied Longitudinal Data Analysis. New York, NY: Oxford University Press; 2003.
    1. Snijders T, Bosker R. Multilevel Analysis: An Introduction to Basic and Advanced Multilevel Modeling. London: Sage Publications; 1999.
    1. Steele F. Multilevel discrete-time event history models with applications to the analysis of recurrent employment transitions. Australian & New Zealand, J Stat. 2011;53:1–26.
    1. Therneau TM, Grambsch PM. Modeling Survival Data: Extending the Cox Model. New York: Springer-Verlag; 2000.
    1. Tu JV, Austin P, Naylor CD. Temporal changes in the outcomes of acute myocardial infarction in Ontario, 1992–96. Can Med Assoc J. 1999;161:1257–1261.
    1. Tu JV, Austin PC, Walld R, Roos L, Agras J, McDonald KM. Development and validation of the Ontario acute myocardial infarction mortality prediction rules. J Am Coll Cardiol. 2001;37:992–997.
    1. Whitehead J. Fitting Cox’s regression model to survival data using GLIM. J R Stat Soc Ser C. 1980;29:268–275.
    1. Wienke A. Frailty Models in Survival Analysis. Boca Raton, FL: Chapman & Hall/CRC; 2011.

Source: PubMed

Sponsors and Collaborators

Medical Conditions

Drug Interventions

3
Subscribe