Bayesian Cox Regression
Bayesian Cox Proportional Hazards Regression · Also known as: Bayesian Cox PH model, Bayesian proportional hazards model, Bayesian survival regression, BCox
Bayesian Cox regression combines the Cox proportional hazards model for time-to-event data with Bayesian inference. Instead of point estimates, it produces full posterior distributions over the hazard ratios, naturally incorporating prior knowledge and providing coherent uncertainty quantification even with small samples or informative censoring.
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 Cox regression when you have time-to-event outcomes with censoring and need more than point estimates — especially when the sample is small or moderately sized, prior information from previous studies is available, complete uncertainty propagation is required (e.g., for decision analysis), or the goal is model comparison via information criteria. It is also preferable when classical maximum partial likelihood produces unstable estimates due to sparse data or many tied event times. Avoid it when the dataset is very large and computational cost is prohibitive, or when a classical Cox model already yields stable estimates and transparent frequentist inference is required by the research context.
Strengths & limitations
- Full posterior distributions over hazard ratios enable genuine probabilistic statements (e.g., P(HR > 1) = 0.95) rather than binary significance decisions.
- Prior distributions formally incorporate historical data or expert knowledge, reducing standard errors when sample sizes are small.
- Handles complex censoring, left-truncation, and irregular follow-up naturally within the probability model.
- Model comparison via WAIC or LOO-CV is straightforward and does not require nested-model assumptions.
- Posterior predictive survival curves with calibrated credible bands are directly available for individual prognosis.
- Computationally demanding: MCMC for the baseline hazard is far slower than classical partial likelihood maximization.
- Results are sensitive to the choice of prior on the baseline hazard, and misspecified priors can distort inference in small samples.
- Requires expertise in prior specification, MCMC diagnostics (Rhat, effective sample size), and model checking — more demanding than standard Cox regression.
Frequently asked
Do I need to specify a prior for the baseline hazard?
Yes, in a fully Bayesian Cox model the baseline hazard h0(t) must be assigned a prior. Common choices are a piecewise-constant hazard with gamma priors on each piece, or a gamma process prior. Software such as Stan or rstanarm handles this automatically with sensible defaults, but you should inspect the sensitivity of key posterior quantities to the prior specification.
How does this differ from the classical Cox model?
The classical Cox model uses partial likelihood to estimate β without modeling h0(t), giving point estimates and confidence intervals based on asymptotic theory. The Bayesian version places priors on both β and h0(t) and produces full posterior distributions, which are more informative about uncertainty and allow coherent probabilistic predictions.
Is the proportional hazards assumption still required?
Yes. The Bayesian framework does not relax the proportional hazards assumption. You must still check it using Schoenfeld residuals or time-varying coefficient tests; if it fails, consider a Bayesian accelerated failure time model or a model with time-by-covariate interactions.
Which software can I use for Bayesian Cox regression?
Stan (via the survival package or direct coding), rstanarm (survival extension), brms (censored-response families), and JAGS are all commonly used. R packages such as bayesm and BayesSurv also provide dedicated functionality.
How do I compare two Bayesian Cox models?
Use the widely applicable information criterion (WAIC) or leave-one-out cross-validation (LOO-CV) implemented in the loo R package. These penalize model complexity while rewarding predictive accuracy and do not require the models to be nested.
Sources
- Ibrahim, J. G., Chen, M.-H., & Sinha, D. (2001). Bayesian Survival Analysis. Springer. ISBN: 978-0387952772
- Cox, D. R. (1972). Regression models and life-tables. Journal of the Royal Statistical Society, Series B, 34(2), 187–220. link ↗
How to cite this page
ScholarGate. (2026, June 3). Bayesian Cox Proportional Hazards Regression. ScholarGate. https://scholargate.app/en/statistics/bayesian-cox-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 Generalized Linear ModelStatistics↔ compare
- Bayesian Mixed Effects ModelStatistics↔ compare
- Bayesian Survival regressionStatistics↔ compare
- Cox RegressionSurvival↔ compare
- Survival RegressionStatistics↔ compare
- Zero-inflated modelStatistics↔ compare