Bayesian Survival Regression
Also known as: Bayesian time-to-event regression, Bayesian parametric survival model, Bayesian survival analysis, Bayesian accelerated failure time model
Bayesian Survival Regression combines parametric or semiparametric survival models — such as Weibull, log-normal, or Cox proportional hazards — with Bayesian inference. Instead of point estimates, it produces full posterior distributions for regression coefficients and the baseline hazard, naturally handling censored observations and incorporating prior knowledge about event times or covariate effects.
Read the full method
Sign in with a free account to read this section.
Method map
The neighbourhood of related methods — select a node to explore.
When to use it
Use Bayesian survival regression when your outcome is time to an event (death, relapse, failure, default) and observations may be censored. It is particularly advantageous when prior information from historical trials or expert knowledge is available, when the sample is small and frequentist uncertainty intervals are unreliable, or when you need full predictive distributions for individual patients or items. Avoid it when you have no defensible prior (use Cox regression instead), when computational resources are severely limited, or when a simple Kaplan-Meier or log-rank comparison suffices because no covariate adjustment is needed.
Strengths & limitations
- Produces full posterior distributions for all parameters, giving genuine probabilistic statements about hazard ratios and survival probabilities.
- Naturally incorporates prior information from historical data or expert knowledge through the prior distribution.
- Handles censored and incomplete observations directly within the likelihood without ad hoc adjustments.
- Predictive inference for new subjects propagates all parameter uncertainty, yielding calibrated credible intervals.
- Supports complex hierarchical or multilevel structures (e.g., patients nested in clinics) within a unified Bayesian framework.
- Flexible choice of baseline hazard family (parametric: Weibull, log-normal; semiparametric: gamma process prior on cumulative hazard).
- Computationally expensive: MCMC sampling can take hours or days on large datasets or complex models.
- Results depend on the prior specification; a poorly chosen prior can distort inference, especially with small samples.
- Requires careful convergence diagnostics (R-hat, effective sample size) and expertise to identify MCMC failures.
- Interpreting and communicating posterior distributions is more involved than reporting a hazard ratio and p-value.
- Proportional-hazards assumption still applies unless time-varying coefficients or flexible hazard functions are explicitly modelled.
Frequently asked
How does Bayesian survival regression differ from the classical Cox model?
The classical Cox model estimates a single set of coefficient values by maximising the partial likelihood. Bayesian survival regression places prior distributions on the coefficients and baseline hazard, then uses MCMC to obtain a full posterior distribution. This yields credible intervals, predictive distributions, and the ability to incorporate prior knowledge — but at higher computational cost.
What prior should I use for the regression coefficients?
A weakly informative normal prior (mean 0, standard deviation 1–2 on the log-hazard-ratio scale) is a safe default. It regularises estimates without strongly constraining them. If you have historical data supporting specific effect sizes, a more informative normal or t prior centred on those values can improve estimation. Avoid completely flat priors on unbounded parameters.
Can I use this method with interval-censored or left-censored data?
Yes. The Bayesian likelihood framework extends naturally to interval censoring (event occurred in a known interval) and left censoring (event occurred before observation began) by adjusting the likelihood contribution for each observation accordingly. Stan and JAGS both support these censoring types via data augmentation or truncated distributions.
How do I check whether the proportional-hazards assumption holds?
Even within a Bayesian framework, you can compute Schoenfeld residuals from posterior draws, plot log-log survival curves, or extend the model with time-varying coefficients. If the assumption is violated, a Bayesian accelerated failure time (AFT) model or a flexible baseline hazard (e.g., piecewise exponential with many segments) may be more appropriate.
Which software can fit Bayesian survival regression?
Stan (via the brms or rstanarm R packages) and JAGS are the most widely used. brms provides a formula interface similar to lm/lmer and supports Weibull, log-normal, and Cox-type models with Bayesian estimation. Python users can use PyMC or Stan's cmdstanpy interface.
Sources
- Ibrahim, J. G., Chen, M.-H., & Sinha, D. (2001). Bayesian Survival Analysis. Springer. ISBN: 978-0387952772
- Gelman, A., Carlin, J. B., Stern, H. S., Dunson, D. B., Vehtari, A., & Rubin, D. B. (2013). Bayesian Data Analysis (3rd ed.). CRC Press. ISBN: 978-1439840955
How to cite this page
ScholarGate. (2026, June 3). Bayesian Survival Regression. ScholarGate. https://scholargate.app/en/statistics/bayesian-survival-regression
Which method?
Set this method beside its closest kin and read them side by side — the library lays the books on the table; the choice is yours.
- Bayesian Cox RegressionStatistics↔ compare
- Bayesian Generalized Linear ModelStatistics↔ compare
- Bayesian Mixed Effects ModelStatistics↔ compare
- Cox RegressionSurvival↔ compare
- Survival RegressionStatistics↔ compare