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.

Bayes Factors for Choosing among Six Common Survival Models.

spBayesSurv: Fitting Bayesian Spatial Survival Models Using R

Parametric Regression Model for Survival Data: Weibull Regression Model As an Example

Semi-parametric Regression Model for Survival Data: Graphical Visualization with R.

Bayesian Variable Selection In Semiparametric Proportional Hazards Model For High Dimensional Survival Data

PROPORTIONAL HAZARDS MODELS WITH CONTINUOUS MARKS

Bayesian parametric models for survival prediction in medical applications

A scalable and flexible Cox proportional hazards model for high-dimensional survival prediction and functional selection

Bayesian Semiparametric Joint Model of Multivariate Longitudinal and Survival Data with Dependent Censoring.

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

Efficient Estimation for the Proportional Hazards Model with Competing Risks and Current Status Data

Penalized variable selection in multi-parameter regression survival modeling

Modified Cox Models: A Simulation Study on Different Survival Distributions, Censoring Rates, and Sample Sizes

Bayesian transformation cure frailty models with multivariate failure time data

FastSurvival: Hidden Computational Blessings in Training Cox Proportional Hazards Models

Predicting Survival Probabilities with Semiparametric Transformation Models

A flexible parametric modelling framework for survival analysis

Additive Hazards Model For Competing Risks Analysis Of The Case-Cohort Design

Optimal subsampling for parametric accelerated failure time models with massive survival data

Bayesian pairwise meta‐analysis of time‐to‐event outcomes in the presence of non‐proportional hazards: A simulation study of flexible parametric, piecewise exponential and fractional polynomial models

A Proportional Hazards Model for the Subdistribution of a Competing Risk