Multilevel mixed effects parametric survival analysis

With the release of Stata 14 came the mestreg command to fit multilevel mixed effects parametric survival models, assuming normally distributed random effects, estimated with maximum likelihood utilising Gaussian quadrature. In this article, I present the user written stmixed command, which serves as both an alternative and a complimentary program for the fitting of multilevel parametric survival models, to mestreg. The key extensions include incorporation of the flexible parametric Royston-Parmar survival model, and the ability to fit multilevel relative survival models. The methods are illustrated with a commonly used dataset of patients with kidney disease suffering recurrent infections, and a simulated example, illustrating a simple approach to simulating clustered survival data using survsim (Crowther and Lambert 2012, 2013).

[1]  Paul C. Lambert,et al.  Further Development of Flexible Parametric Models for Survival Analysis , 2009 .

[2]  Richard D Riley,et al.  Individual patient data meta-analysis of survival data using Poisson regression models , 2012, BMC Medical Research Methodology.

[3]  Sam Harper,et al.  Use of relative and absolute effect measures in reporting health inequalities: structured review , 2012, BMJ : British Medical Journal.

[4]  R. Gutierrez.  Parametric Frailty and Shared Frailty Survival Models , 2002 .

[5]  C. Mcgilchrist,et al.  Regression with frailty in survival analysis. , 1991, Biometrics.

[6]  Paula R Williamson,et al.  Investigating heterogeneity in an individual patient data meta‐analysis of time to event outcomes , 2005, Statistics in medicine.

[7]  T. Mok,et al.  Gefitinib or carboplatin-paclitaxel in pulmonary adenocarcinoma. , 2009, The New England journal of medicine.

[8]  D.,et al.  Regression Models and Life-Tables , 2022 .

[9]  Mark C Simmonds,et al.  Meta-analysis of individual patient data from randomized trials: a review of methods used in practice , 2005, Clinical trials.

[10]  Catherine Legrand,et al.  Baseline and treatment effect heterogeneity for survival times between centers using a random effects accelerated failure time model with flexible error distribution , 2007, Statistics in medicine.

[11]  Benjamin Djulbegovic,et al.  A Flexible Alternative to the Cox Proportional Hazards Model for Assessing the Prognostic Accuracy of Hospice Patient Survival , 2012, PloS one.

[12]  Paul C. Lambert,et al.  Flexible Parametric Survival Analysis Using Stata: Beyond the Cox Model , 2011 .

[13]  Enzo Coviello,et al.  Estimating and Modeling Relative Survival , 2015 .

[14]  Paul C. Lambert,et al.  Simulating Complex Survival Data , 2012 .

[15]  S. Rabe-Hesketh,et al.  Reliable Estimation of Generalized Linear Mixed Models using Adaptive Quadrature , 2002 .

[16]  P. Janssen,et al.  A Bayesian approach to jointly estimate centre and treatment by centre heterogeneity in a proportional hazards model , 2005, Statistics in medicine.

[17]  P. Royston,et al.  Flexible parametric proportional‐hazards and proportional‐odds models for censored survival data, with application to prognostic modelling and estimation of treatment effects , 2002, Statistics in medicine.

[18]  D. Bates,et al.  Approximations to the Log-Likelihood Function in the Nonlinear Mixed-Effects Model , 1995 .

[19]  Richard D Riley,et al.  Multilevel mixed effects parametric survival models using adaptive Gauss–Hermite quadrature with application to recurrent events and individual participant data meta‐analysis , 2014, Statistics in medicine.

[20]  Henrik Møller,et al.  Quantifying differences in breast cancer survival between England and Norway. , 2011, Cancer epidemiology.

[21]  Paul W Dickman,et al.  Regression models for relative survival , 2004, Statistics in medicine.

[22]  Michael J Crowther,et al.  Simulating biologically plausible complex survival data , 2013, Statistics in medicine.

[23]  Lei Liu,et al.  The use of Gaussian quadrature for estimation in frailty proportional hazards models , 2008, Statistics in medicine.

[24]  T. Yamaguchi,et al.  Proportional hazards models with random effects to examine centre effects in multicentre cancer clinical trials , 2002, Statistical methods in medical research.

[25]  Ben Goldacre,et al.  Survival Analysis of Irish Amyotrophic Lateral Sclerosis Patients Diagnosed from 1995–2010 , 2013, PloS one.