English

Full Bayesian inference with hazard mixture models

Statistics Theory 2016-05-04 v4 Methodology Statistics Theory

Abstract

Bayesian nonparametric inferential procedures based on Markov chain Monte Carlo marginal methods typically yield point estimates in the form of posterior expectations. Though very useful and easy to implement in a variety of statistical problems, these methods may suffer from some limitations if used to estimate non-linear functionals of the posterior distribution. The main goal of the present paper is to develop a novel methodology that extends a well-established marginal procedure designed for hazard mixture models, in order to draw approximate inference on survival functions that is not limited to the posterior mean but includes, as remarkable examples, credible intervals and median survival time. Our approach relies on a characterization of the posterior moments that, in turn, is used to approximate the posterior distribution by means of a technique based on Jacobi polynomials. The inferential performance of our methodology is analyzed by means of an extensive study of simulated data and real data consisting of leukemia remission times. Although tailored to the survival analysis context, the procedure we introduce can be adapted to a range of other models for which moments of the posterior distribution can be estimated.

Keywords

Cite

@article{arxiv.1405.6628,
  title  = {Full Bayesian inference with hazard mixture models},
  author = {Julyan Arbel and Antonio Lijoi and Bernardo Nipoti},
  journal= {arXiv preprint arXiv:1405.6628},
  year   = {2016}
}

Comments

27 pages, 5 figures

R2 v1 2026-06-22T04:23:27.756Z