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.

Between-within effects in survival models with cross-classified clustering: Application to the evaluation of the effectiveness of medical devices

A review and empirical comparison of causal inference methods for clustered observational data with application to the evaluation of the effectiveness of medical devices

Overview of model validation for survival regression model with competing risks using melanoma study data

Using cross-validation to evaluate predictive accuracy of survival risk classifiers based on high-dimensional data

Weibull Regression and Machine Learning Survival Models: Methodology, Comparison, and Application to Biomedical Data Related to Cardiac Surgery

A stratification approach using logit-based models for confounder adjustment in the study of continuous outcomes

Targeted Estimation of Heterogeneous Treatment Effect in Observational Survival Analysis

A flexible approach for causal inference with multiple treatments and clustered survival outcomes

Application of Cox Model to predict the survival of patients with Chronic Heart Failure: A latent class regression approach

Nonparametric Causal Survival Analysis with Clustered Interference

MODELLING COUNTY LEVEL BREAST CANCER SURVIVAL DATA USING A COVARIATE-ADJUSTED FRAILTY PROPORTIONAL HAZARDS MODEL

Causal effect estimation in survival analysis with high dimensional confounders

Bayesian parametric models for survival prediction in medical applications

Modeling county level breast cancer survival data using a covariate-adjusted frailty proportional hazards model

Hierarchical Bayesian modeling of heterogeneous outcome variance in cluster randomized trials

Interpretable (not just posthoc-explainable) heterogeneous survivor bias-corrected treatment effects for assignment of postdischarge interventions to prevent readmissions

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

Joint Modeling of Longitudinal and Survival Data with Censored Single-index Varying Coefficient Models

SEMIPARAMETRIC LATENT-CLASS MODELS FOR MULTIVARIATE LONGITUDINAL AND SURVIVAL DATA

Evaluating analytic models for individually randomized group treatment trials with complex clustering in nested and crossed designs