The current application of the Royston-Parmar model for prognostic modeling in health research: a scoping review

BackgroundPrognostic models incorporating survival analysis predict the risk (i.e., probability) of experiencing a future event over a specific time period. In 2002, Royston and Parmar described a type of flexible parametric survival model called the Royston-Parmar model in Statistics in Medicine, a model which fits a restricted cubic spline to flexibly model the baseline log cumulative hazard on the proportional hazards scale. This feature permits absolute measures of effect (e.g., hazard rates) to be estimated at all time points, an important feature when using the model. The Royston-Parmar model can also incorporate time-dependent effects and be used on different scales (e.g., proportional odds, probit). These features make the Royston-Parmar model attractive for prediction, yet their current uptake for prognostic modeling is unknown. Thus, the objectives were to conduct a scoping review of how the Royston-Parmar model has been applied to prognostic models in health research, to raise awareness of the model, to identify gaps in current reporting, and to offer model building considerations and reporting suggestions for other researchers.MethodsFive electronic databases and gray literature indexed in web sources from 2001 to 2016 were searched to identify articles for inclusion in the scoping review. Two reviewers independently screened 1429 articles, and after applying exclusion criteria through a two-step screening process, data from 12 studies were abstracted.ResultsSince 2001, only 12 studies were identified that used the Royston-Parmar model in some capacity for prognostic modeling, 10 of which used the model as the basis for their prognostic model. The restricted cubic spline varied across studies in the number of interior knots (range 1 to 6), and only three studies reported knot placement. Three studies provided details about the baseline function, with two studies using a figure and the third providing coefficients. However, no studies provided adequate information on their restricted cubic spline to permit others to validate or completely use the model.ConclusionsDespite the advantages of the Royston-Parmar model for prognostic models, they are not widely used in health research. Better reporting of details about the restricted cubic spline is needed, so the prognostic model can be used and validated by others.RegistrationThe protocol was registered with Open Science Framework (https://osf.io/r3232/).

[1]  Richard D Riley,et al.  Prediction of risk of recurrence of venous thromboembolism following treatment for a first unprovoked venous thromboembolism: systematic review, prognostic model and clinical decision rule, and economic evaluation. , 2016, Health technology assessment.

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

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

[4]  B. Møller,et al.  An empirical comparison of methods for predicting net survival. , 2016, Cancer epidemiology.

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

[6]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[7]  N. Obuchowski,et al.  Assessing the Performance of Prediction Models: A Framework for Traditional and Novel Measures , 2010, Epidemiology.

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

[9]  T. Peto,et al.  Predictors of First Recurrence of Clostridium difficile Infection: Implications for Initial Management , 2012, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[10]  Robert H. Halstead,et al.  Matrix Computations , 2011, Encyclopedia of Parallel Computing.

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

[12]  Gary S Collins,et al.  Transparent Reporting of a multivariable prediction model for Individual Prognosis Or Diagnosis (TRIPOD): Explanation and Elaboration , 2015, Annals of Internal Medicine.

[13]  Charles J. Stone,et al.  Additive Splines in Statistics , 2015 .

[14]  H. Arksey,et al.  Scoping studies: towards a methodological framework , 2005 .

[15]  Judith D. Goldberg,et al.  Applied Survival Analysis , 1999, Technometrics.

[16]  Christopher Dudley,et al.  Predicting patient survival after deceased donor kidney transplantation using flexible parametric modelling , 2016, BMC Nephrology.

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

[18]  E. Steyerberg,et al.  [Regression modeling strategies]. , 2011, Revista espanola de cardiologia.

[19]  Frank E. Harrell,et al.  Regression Modeling Strategies: With Applications to Linear Models, Logistic Regression, and Survival Analysis , 2001 .

[20]  K. Wilson,et al.  Predictive risk algorithms in a population setting: an overview , 2012, Journal of Epidemiology & Community Health.

[21]  Keith R Abrams,et al.  Flexible parametric joint modelling of longitudinal and survival data , 2012, Statistics in medicine.

[22]  Patricia McInerney,et al.  The Joanna Briggs Institute reviewers' manual 2015: methodology for JBI scoping reviews , 2015 .

[23]  F. Harrell,et al.  Regression modelling strategies for improved prognostic prediction. , 1984, Statistics in medicine.

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

[25]  C. Hermenegildo,et al.  Frailty and other geriatric conditions for risk stratification of older patients with acute coronary syndrome. , 2014, American heart journal.

[26]  D. Fuchs,et al.  Immunological markers of frailty predict outcomes beyond current risk scores in aortic stenosis following transcatheter aortic valve replacement: Role of neopterin and tryptophan , 2016 .

[27]  S. Berhane,et al.  Biomarker-based prognosis in hepatocellular carcinoma: validation and extension of the BALAD model , 2014, British Journal of Cancer.

[28]  M. Mansournia,et al.  Flexible parametric survival models built on age-specific antimüllerian hormone percentiles are better predictors of menopause , 2016, Menopause.

[29]  F E Harrell,et al.  The restricted cubic spline as baseline hazard in the proportional hazards model with step function time-dependent covariables. , 1995, Statistics in medicine.

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

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

[32]  Laurence L. George,et al.  The Statistical Analysis of Failure Time Data , 2003, Technometrics.

[33]  D. Mark,et al.  Clinical prediction models: are we building better mousetraps? , 2003, Journal of the American College of Cardiology.

[34]  M. Weinstock,et al.  Prognostic survival model for people diagnosed with invasive cutaneous melanoma , 2015, BMC Cancer.

[35]  R. D'Agostino,et al.  Validation of the Framingham coronary heart disease prediction scores: results of a multiple ethnic groups investigation. , 2001, JAMA.

[36]  David W. Hosmer,et al.  Applied Survival Analysis: Regression Modeling of Time-to-Event Data , 2008 .

[37]  K. Anderson,et al.  An updated coronary risk profile. A statement for health professionals. , 1991, Circulation.

[38]  M. Gail,et al.  Projecting individualized probabilities of developing breast cancer for white females who are being examined annually. , 1989, Journal of the National Cancer Institute.

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

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

[41]  Patrick Royston,et al.  Construction and validation of a prognostic model across several studies, with an application in superficial bladder cancer , 2004, Statistics in medicine.

[42]  L. Tao,et al.  [Transparent Reporting of a Multivariable Prediction Model for Individual Prognosis or Diagnosis]. , 2018, Zhonghua yi xue za zhi.

[43]  Frank E. Harrell,et al.  The restricted cubic spline hazard model , 1990 .

[44]  G. Collins,et al.  Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD Statement , 2015, BMC Medicine.

[45]  Keith R. Abrams,et al.  Joint modelling of longitudinal and survival data: incorporating delayed entry and an assessment of model misspecification , 2015, Statistics in medicine.

[46]  Patrick Royston Flexible alternatives to the Cox model, and more , 2001 .

[47]  Paul C Lambert,et al.  Flexible parametric modelling of cause-specific hazards to estimate cumulative incidence functions , 2013, BMC Medical Research Methodology.

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

[49]  Sally R. Hinchliffe,et al.  Extending the Flexible Parametric Survival Model for Competing Risks , 2013 .

[50]  Martha Sajatovic,et al.  Clinical Prediction Models , 2013 .

[51]  J. Aitken,et al.  Estimating the change in life expectancy after a diagnosis of cancer among the Australian population , 2015, BMJ Open.

[52]  J. Castillo,et al.  Population‐based prognostic factors for survival in patients with Burkitt lymphoma: An analysis from the Surveillance, Epidemiology, and End Results database , 2013, Cancer.

[53]  Mats Lambe,et al.  Estimating the loss in expectation of life due to cancer using flexible parametric survival models , 2013, Statistics in medicine.

[54]  Laura C Rosella,et al.  A population-based risk algorithm for the development of diabetes: development and validation of the Diabetes Population Risk Tool (DPoRT) , 2010, Journal of Epidemiology & Community Health.

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

[56]  J. Kalbfleisch,et al.  The Statistical Analysis of Failure Time Data: Kalbfleisch/The Statistical , 2002 .

[57]  G. Brier,et al.  External correspondence: Decompositions of the mean probability score , 1982 .

[58]  N. Breslow Covariance analysis of censored survival data. , 1974, Biometrics.