flexsurv: A Platform for Parametric Survival Modeling in R.

flexsurv is an R package for fully-parametric modeling of survival data. Any parametric time-to-event distribution may be fitted if the user supplies a probability density or hazard function, and ideally also their cumulative versions. Standard survival distributions are built in, including the three and four-parameter generalized gamma and F distributions. Any parameter of any distribution can be modeled as a linear or log-linear function of covariates. The package also includes the spline model of Royston and Parmar (2002), in which both baseline survival and covariate effects can be arbitrarily flexible parametric functions of time. The main model-fitting function, flexsurvreg, uses the familiar syntax of survreg from the standard survival package (Therneau 2016). Censoring or left-truncation are specified in 'Surv' objects. The models are fitted by maximizing the full log-likelihood, and estimates and confidence intervals for any function of the model parameters can be printed or plotted. flexsurv also provides functions for fitting and predicting from fully-parametric multi-state models, and connects with the mstate package (de Wreede, Fiocco, and Putter 2011). This article explains the methods and design principles of the package, giving several worked examples of its use.

[1]  P. Hougaard,et al.  Frailty models for survival data , 1995, Lifetime data analysis.

[2]  Patrick Royston,et al.  Flexible Parametric Alternatives to the Cox Model, and more , 2001 .

[3]  Tatiana Benaglia,et al.  Survival extrapolation in the presence of cause specific hazards , 2015, Statistics in medicine.

[4]  Competing Risks and Multi-State Models : Overview and Case Studies , 2018 .

[5]  N. Keiding,et al.  Multi-state models for event history analysis , 2002, Statistical methods in medical research.

[6]  Micha Mandel,et al.  Simulation-Based Confidence Intervals for Functions With Complicated Derivatives , 2013 .

[7]  Saralees Nadarajah,et al.  A new R package for actuarial survival models , 2013, Comput. Stat..

[8]  Richard J. Cook,et al.  Statistical Issues in Modeling Chronic Disease in Cohort Studies , 2014 .

[9]  P. Grambsch,et al.  A Package for Survival Analysis in S , 1994 .

[10]  P. Royston,et al.  Building multivariable prognostic and diagnostic models: transformation of the predictors by using fractional polynomials , 1999 .

[11]  Michael J. Crowther,et al.  Multilevel mixed effects parametric survival analysis , 2013, 1709.06633.

[12]  Nicholas R Latimer,et al.  Survival Analysis for Economic Evaluations Alongside Clinical Trials—Extrapolation with Patient-Level Data , 2013, Medical decision making : an international journal of the Society for Medical Decision Making.

[13]  E. Stacy.  A Generalization of the Gamma Distribution , 1962 .

[14]  P. Grambsch.  Survival and Event History Analysis: A Process Point of View by AALEN, O. O., BORGAN, O., and GJESSING, H. K. , 2009 .

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

[16]  R. Rigby,et al.  Generalized additive models for location, scale and shape , 2005 .

[17]  David J. Lunn,et al.  Survival extrapolation using the poly-Weibull model , 2015, Statistical methods in medical research.

[18]  H Putter,et al.  Tutorial in biostatistics: competing risks and multi‐state models , 2007, Statistics in medicine.

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

[20]  Andrew C Titman,et al.  Flexible Nonhomogeneous Markov Models for Panel Observed Data , 2011, Biometrics.

[21]  J. Klein,et al.  Statistical Models Based On Counting Process , 1994 .

[22]  Nancy Reid,et al.  On “A conversation with Sir David Cox” , 1994, Issue 5.2, Spring 2023.

[23]  Christopher H. Jackson,et al.  Multi-State Models for Panel Data: The msm Package for R , 2011 .

[24]  Hein Putter,et al.  Reduced-rank proportional hazards regression and simulation-based prediction for multi-state models. , 2008, Statistics in medicine.

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

[26]  H Putter,et al.  Reduced rank proportional hazards model for competing risks. , 2006, Biostatistics.

[27]  Patrick Royston,et al.  Flexible Parametric Alternatives to the Cox Model: Update , 2004 .

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

[29]  C. Wild,et al.  Vector Generalized Additive Models , 1996 .

[30]  D Heng,et al.  Bronchiolitis obliterans syndrome: incidence, natural history, prognosis, and risk factors. , 1998, The Journal of heart and lung transplantation : the official publication of the International Society for Heart Transplantation.

[31]  H. Müller,et al.  Hazard rate estimation under random censoring with varying kernels and bandwidths. , 1994, Biometrics.

[32]  J. Kalbfleisch,et al.  The Analysis of Panel Data under a Markov Assumption , 1985 .

[33]  R. Prentice.  A LOG GAMMA MODEL AND ITS MAXIMUM LIKELIHOOD ESTIMATION , 1974 .

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

[35]  Hein Putter,et al.  mstate: An R Package for the Analysis of Competing Risks and Multi-State Models , 2011 .

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

[37]  P. Royston,et al.  Regression using fractional polynomials of continuous covariates: parsimonious parametric modelling. , 1994 .

[38]  Ross L. Prentice,et al.  Discrimination among some parametric models , 1975 .

[39]  Christopher Cox,et al.  The generalized F distribution: An umbrella for parametric survival analysis , 2008, Statistics in medicine.

[40]  F Louzada-Neto,et al.  Polyhazard Models for Lifetime Data , 1999, Biometrics.

[41]  Hein Putter,et al.  The mstate package for estimation and prediction in non- and semi-parametric multi-state and competing risks models , 2010, Comput. Methods Programs Biomed..

[42]  Francesca Ieva,et al.  Multi-state modelling of repeated hospitalisation and death in patients with heart failure: The use of large administrative databases in clinical epidemiology , 2017, Statistical methods in medical research.

[43]  P. Royston,et al.  Selection of important variables and determination of functional form for continuous predictors in multivariable model building , 2007, Statistics in medicine.

[44]  Karline Soetaert,et al.  Solving Differential Equations in R: Package deSolve , 2010 .