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.

Bayesian Survival Analysis in STAN for Improved Measuring of Uncertainty in Parameter Estimates

Bayesian Survival Analysis Using the rstanarm R Package

A Bayesian Survival Model for Time-Varying Coefficients and Unobserved Heterogeneity

Fitting Bayesian item response models in Stata and Stan

Bayesian parametric models for survival prediction in medical applications

Bayesian Survival Model based on Moment Characterization

Maximum Likelihood Estimation of Flexible Survival Densities with Importance Sampling

Restricted Mean Survival Time Estimation Using Bayesian Nonparametric Dependent Mixture Models

Bayesian analyses of multiple random change points in survival models with applications to clinical trials

Flexible and efficient Bayesian pharmacometrics modeling using Stan and Torsten, Part I

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

High-dimensional Bayesian Parameter Estimation: Case Study for a Model of JAK2/STAT5 Signaling.

Dependent Dirichlet Processes for Analysis of a Generalized Shared Frailty Model

Dependent Dirichlet Processes for Analysis of aGeneralized Shared Frailty Model

Bayesian alternatives to null hypothesis significance testing in biomedical research: a non-technical introduction to Bayesian inference with JASP

A flexible parametric modelling framework for survival analysis

A two-stage approach for Bayesian joint models: reducing complexity while maintaining accuracy

The use of restricted mean survival time to estimate treatment effect under model misspecification, a simulation study

causalBETA: An R Package for Bayesian Semiparametric Causal Inference with Event-Time Outcomes

Semiparametric regression modelling of current status competing risks data: a Bayesian approach

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