Capturing simple and complex time-dependent effects using flexible parametric survival models: A simulation study

Abstract Non-proportional hazards are common within time-to-event data and can be modeled using restricted cubic splines in flexible parametric survival models. This simulation study assesses the ability of these models in capturing non-proportional hazards, and the ability of the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) in selecting degrees of freedom. The simulation results for scenarios with differing complexities showed little bias in the survival and hazard functions for simple scenarios; bias increased in complex scenarios when fewer degrees of freedom were modeled. Neither AIC nor BIC consistently performed better and both generally selected models with little bias.

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

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

[3]  Michael J Crowther,et al.  A general framework for parametric survival analysis , 2014, Statistics in medicine.

[4]  Adrian E. Raftery,et al.  Bayesian Model Averaging: A Tutorial , 2016 .

[5]  Adrian E. Raftery,et al.  Bayesian Model Averaging , 1998 .

[6]  C. Z. Wei On Predictive Least Squares Principles , 1992 .

[7]  P. Lambert,et al.  The loss in expectation of life after colon cancer: a population-based study , 2015, BMC Cancer.

[8]  R. Simon,et al.  Flexible regression models with cubic splines. , 1989, Statistics in medicine.

[9]  Miguel A Hernán,et al.  The hazards of hazard ratios. , 2010, Epidemiology.

[10]  Adrian E. Raftery,et al.  Bayesian model averaging: a tutorial (with comments by M. Clyde, David Draper and E. I. George, and a rejoinder by the authors , 1999 .

[11]  Paul W Dickman,et al.  Estimating and modelling cure in population-based cancer studies within the framework of flexible parametric survival models , 2011, BMC medical research methodology.

[12]  V. Carstairs,et al.  Deprivation and health in Scotland. , 1990, Health bulletin.

[13]  P. Dieppe,et al.  Failure rates of stemmed metal-on-metal hip replacements: analysis of data from the National Joint Registry of England and Wales , 2012, The Lancet.

[14]  Paul C. Lambert,et al.  Strcs: A Command for Fitting Flexible Parametric Survival Models on the Log-hazard Scale , 2016 .

[15]  P. Lambert,et al.  Continued improvement in survival of acute myeloid leukemia patients: an application of the loss in expectation of life , 2016, Blood Cancer Journal.

[16]  D. Madigan,et al.  Correction to: ``Bayesian model averaging: a tutorial'' [Statist. Sci. 14 (1999), no. 4, 382--417; MR 2001a:62033] , 2000 .

[17]  Ralf Bender,et al.  Generating survival times to simulate Cox proportional hazards models , 2005, Statistics in medicine.

[18]  R. A. Leibler,et al.  On Information and Sufficiency , 1951 .

[19]  G. Kitagawa,et al.  Bootstrapping Log Likelihood and EIC, an Extension of AIC , 1997 .

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

[21]  David R. Anderson,et al.  Model selection and multimodel inference : a practical information-theoretic approach , 2003 .

[22]  H. Akaike,et al.  Information Theory and an Extension of the Maximum Likelihood Principle , 1973 .

[23]  P. Lambert,et al.  Temporal trends in the proportion cured for cancer of the colon and rectum: A population‐based study using data from the Finnish Cancer Registry , 2007, International journal of cancer.

[24]  P. Royston,et al.  A New Proposal for Multivariable Modelling of Time‐Varying Effects in Survival Data Based on Fractional Polynomial Time‐Transformation , 2007, Biometrical journal. Biometrische Zeitschrift.

[25]  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.

[26]  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.

[27]  P. Lambert,et al.  Estimating the cure proportion of malignant melanoma, an alternative approach to assess long term survival: a population-based study. , 2014, Cancer epidemiology.

[28]  P. Carbone Breast cancer adjuvant therapy , 1990, Cancer.

[29]  Paul C. Lambert,et al.  The use of restricted cubic splines to approximate complex hazard functions in the analysis of time-to-event data: a simulation study , 2015 .

[30]  Henry E. Kyburg,et al.  Comparison of Approaches , 1974 .

[31]  Paulo César Emiliano,et al.  Information criteria: How do they behave in different models? , 2014, Comput. Stat. Data Anal..

[32]  K. Czene,et al.  Prognosis of patients with breast cancer: causes of death and effects of time since diagnosis, age, and tumor characteristics. , 2011, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[33]  Paul C. Lambert,et al.  stgenreg: A Stata Package for General Parametric Survival Analysis , 2013 .

[34]  R. Tibshirani,et al.  Varying‐Coefficient Models , 1993 .

[35]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[36]  Robert Gray,et al.  Flexible Methods for Analyzing Survival Data Using Splines, with Applications to Breast Cancer Prognosis , 1992 .

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

[38]  Aileen Clarke,et al.  Setting benchmark revision rates for total hip replacement: analysis of registry evidence , 2015, BMJ : British Medical Journal.

[39]  G. Schwarz Estimating the Dimension of a Model , 1978 .

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

[41]  Paul C Lambert,et al.  Flexible parametric models for relative survival, with application in coronary heart disease , 2007, Statistics in medicine.

[42]  Jana Vogel,et al.  Model Selection And Multimodel Inference , 2016 .

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

[44]  Chong Gu Smoothing Spline Anova Models , 2002 .

[45]  M. Abrahamowicz,et al.  Joint estimation of time‐dependent and non‐linear effects of continuous covariates on survival , 2007, Statistics in medicine.

[46]  Y. Pawitan In all likelihood : statistical modelling and inference using likelihood , 2002 .

[47]  Willy Wynant,et al.  Flexible estimation of survival curves conditional on non‐linear and time‐dependent predictor effects , 2016, Statistics in medicine.

[48]  P. Lambert,et al.  Robustness of individual and marginal model-based estimates: A sensitivity analysis of flexible parametric models , 2019, Cancer Epidemiology.

[49]  P. Pynsent,et al.  Mortality and implant revision rates of hip arthroplasty in patients with osteoarthritis: registry based cohort study , 2012, BMJ : British Medical Journal.

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

[51]  Michal Abrahamowicz,et al.  Time-Dependent Hazard Ratio: Modeling and Hypothesis Testing with Application in Lupus Nephritis , 1996 .

[52]  David R. Anderson,et al.  Multimodel Inference , 2004 .

[53]  Jong-Hyeon Jeong,et al.  Breast cancer adjuvant therapy: time to consider its time-dependent effects. , 2011, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[54]  P. Grambsch,et al.  Proportional hazards tests and diagnostics based on weighted residuals , 1994 .

[55]  David R. Cox,et al.  Regression models and life tables (with discussion , 1972 .

[56]  G. Kitagawa,et al.  Generalised information criteria in model selection , 1996 .