Abstract:Background: Evidence-based treatment decisions in medicine are made founded on population-level evidence obtained during randomized clinical trials. In an era of personalized medicine, these decisions should be based on the predicted benefit of a treatment on a patient-level. Survival prediction models play a central role as they incorporate the time-to-event and censoring. In medical applications uncertainty is critical especially when treatments differ in their side effect profiles or costs. Additionally, models must be adapted to local populations without diminishing performance and often without the original training data available due to privacy concern. Both points are supported by Bayesian models-yet they are rarely used. The aim of this work is to evaluate Bayesian parametric survival models on public datasets including cardiology, infectious diseases, and oncology. Materials and methods: Bayesian parametric survival models based on the Exponential and Weibull distribution were implemented as a Python package. A linear combination and a neural network were used for predicting the parameters of the distributions. A superiority design was used to assess whether Bayesian models are better than commonly used models such as Cox Proportional Hazards, Random Survival Forest, and Neural Network-based Cox Proportional Hazards. In a secondary analysis, overfitting was compared between these models. An equivalence design was used to assess whether the prediction performance of Bayesian models after model updating using Bayes rule is equivalent to retraining on the full dataset. Results: In this study, we found that Bayesian parametric survival models perform as good as state-of-the art models while requiring less hyperparameters to be tuned and providing a measure of the uncertainty of the predictions. In addition, these models were less prone to overfitting. Furthermore, we show that updating these models using Bayes rule yields equivalent performance compared to models trained on combined original and new datasets. Conclusions: Bayesian parametric survival models are non-inferior to conventional survival models while requiring less hyperparameter tuning, being less prone to overfitting, and allowing model updating using Bayes rule. Further, the Bayesian models provide a measure of the uncertainty on the statistical inference, and, in particular, on the prediction.

Flexible Bayesian excess hazard models using low-rank thin plate splines

A Unifying Framework for Flexible Excess Hazard Modeling with Applications in Cancer Epidemiology

Flexible Bayesian survival modeling with semiparametric time-dependent and shape-restricted covariate effects

Flexible estimation in cure survival models using Bayesian P-splines

On a general structure for hazard-based regression models: An application to population-based cancer research

Extended excess hazard models for spatially dependent survival data

Flexible joint model for time-to-event and non-Gaussian longitudinal outcomes

A flexible parametric modelling framework for survival analysis

On models for the estimation of the excess mortality hazard in case of insufficiently stratified life tables

Functional proportional hazards mixture cure model and its application to modelling the association between cancer mortality and physical activity in NHANES 2003-2006

Bayesian regularization for flexible baseline hazard functions in Cox survival models

Individual frailty excess hazard models in cancer epidemiology

Penalized Bayesian forward continuation ratio model with application to high-dimensional data with discrete survival outcomes

A multi-state model incorporating estimation of excess hazards and multiple time scales

Bayesian approach for flexible modeling of semicompeting risks data

Nonlinear association structures in flexible Bayesian additive joint models

A Bayesian survival treed hazards model using latent Gaussian processes

A Proportional Hazards Model for the Subdistribution of a Competing Risk

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

Bayesian parametric models for survival prediction in medical applications