Open Access
2010 Bayesian adaptive B-spline estimation in proportional hazards frailty models
Emmanuel Sharef, Robert L. Strawderman, David Ruppert, Mark Cowen, Lakshmi Halasyamani
Electron. J. Statist. 4: 606-642 (2010). DOI: 10.1214/10-EJS566


Frailty models derived from the proportional hazards regression model are frequently used to analyze clustered right-censored survival data. We propose a semiparametric Bayesian methodology for this purpose, modeling both the unknown baseline hazard and density of the random effects using mixtures of B-splines. The posterior distributions for all regression coefficients and spline parameters are obtained using Markov Chain Monte Carlo (MCMC). The methodology permits the use of weighted mixtures of parametric and nonparametric components in modeling the hazard function and frailty distribution; in addition, the spline knots may also be selected adaptively using reversible-jump MCMC. Simulations indicate that the method produces smooth and accurate posterior hazard and frailty density estimates. The Bayesian approach not only produces point estimators that outperform existing approaches in certain circumstances, but also offers a wealth of information about the parameters of interest in the form of MCMC samples from the joint posterior probability distribution. We illustrate the adaptability of the method with data from a study of congestive heart failure.


Download Citation

Emmanuel Sharef. Robert L. Strawderman. David Ruppert. Mark Cowen. Lakshmi Halasyamani. "Bayesian adaptive B-spline estimation in proportional hazards frailty models." Electron. J. Statist. 4 606 - 642, 2010.


Published: 2010
First available in Project Euclid: 6 July 2010

zbMATH: 1329.62406
MathSciNet: MR2660535
Digital Object Identifier: 10.1214/10-EJS566

Keywords: frailty distribution , Hazard regression , heart failure , knot selection , random effect density , re-hospitalization , reversible-jump MCMC , Survival analysis

Rights: Copyright © 2010 The Institute of Mathematical Statistics and the Bernoulli Society

Back to Top