A fully Bayesian approach to inference for Coxian phase-type distributions with covariate dependent mean
Phase-type distributions represent the time to absorption for a finite state Markov chain in continuous time, generalising the exponential distribution and providing a flexible and useful modelling tool. We present a new reversible jump Markov chain Monte Carlo scheme for performing a fully Bayesian analysis of the popular Coxian subclass of phase-type models; the convenient Coxian representation involves fewer parameters than a more general phase-type model. The key novelty of our approach is that we model covariate dependence in the mean whilst using the Coxian phase-type model as a very general residual distribution. Such incorporation of covariates into the model has not previously been attempted in the Bayesian literature. A further novelty is that we also propose a reversible jump scheme for investigating structural changes to the model brought about by the introduction of Erlang phases. Our approach addresses more questions of inference than previous Bayesian treatments of this model and is automatic in nature. We analyse an example dataset comprising lengths of hospital stays of a sample of patients collected from two Australian hospitals to produce a model for a patient's expected length of stay which incorporates the effects of several covariates. This leads to interesting conclusions about what contributes to length of hospital stay with implications for hospital planning. We compare our results with an alternative classical analysis of these data.
References listed on IDEAS
Please report citation or reference errors to , or , if you are the registered author of the cited work, log in to your RePEc Author Service profile, click on "citations" and make appropriate adjustments.:
- Ajay Jasra & David A. Stephens & Christopher C. Holmes, 2007. "Population-Based Reversible Jump Markov Chain Monte Carlo," Biometrika, Biometrika Trust, vol. 94(4), pages 787-807.
- H. Xie & T. J. Chaussalet & P. H. Millard, 2005. "A continuous time Markov model for the length of stay of elderly people in institutional long-term care," Journal of the Royal Statistical Society Series A, Royal Statistical Society, vol. 168(1), pages 51-61.
- C. P. Robert & T. Rydén & D. M. Titterington, 2000. "Bayesian inference in hidden Markov models through the reversible jump Markov chain Monte Carlo method," Journal of the Royal Statistical Society Series B, Royal Statistical Society, vol. 62(1), pages 57-75.
When requesting a correction, please mention this item's handle: RePEc:eee:csdana:v:53:y:2009:i:12:p:4311-4321. See general information about how to correct material in RePEc.
For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: (Shamier, Wendy)
If references are entirely missing, you can add them using this form.