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.

Using Coxian Phase-Type Distributions to Identify Patient Characteristics for Duration of Stay in Hospital

Coxian phase-type distributions are a special type of Markov model that describes duration until an event occurs in terms of a process consisting of a sequence of latent phases. This paper considers the use of Coxian phase-type distributions for modelling patient duration of stay for the elderly in hospital and investigates the potential for using the resulting distribution as a classifying variable to identify common characteristics between different groups of patients according to their (anticipated) length of stay in hospital. The identification of common characteristics for patient length of stay groups would offer hospital managers and clinicians possible insights into the overall management and bed allocation of the hospital wards.

Simulating Coxian Phase-type Distributions for Patient Survival

Coxian phase-type distributions are a special type of Markov model that can be used to represent survival times in terms of phases through which an individual may progress until they eventually leave the system completely. Previous research has considered the Coxian phase-type distribution to be ideal in representing patient survival in hospital. However, problems exist in fitting the distributions. This paper investigates the problems that arise with the fitting process by simulating various Coxian phase-type models for the representation of patient survival and examining the estimated parameter values and eigenvalues obtained. The results indicate that numerical methods previously used for fitting the model parameters do not always converge. An alternative technique is therefore considered. All methods are influenced by the choice of initial parameter values. The investigation uses a data set of 1439 elderly patients and models their survival time, the length of time they spend in a UK hospital.

Investigating the efficiency of fitting Coxian phase-type distributions to health care data

Coxian phase-type distributions are becoming a popular means of representing survival times within a health care environment. They are favoured as they show a distribution as a system of phases and can allow for an easy visual representation of the rate of flow of patients through a system. Difficulties arise, however, in determining the parameter estimates of the Coxian phase-type distribution. This paper examines ways of making the fitting of the Coxian phase-type distribution less cumbersome by outlining different software packages and algorithms available to perform the fit and assessing their capabilities through a number of performance measures. The performance measures rate each of the methods and help in identifying the more efficient. Conclusions drawn from these performance measures suggest SAS to be the most robust package. It has a high rate of convergence in each of the four example model fits considered, short computational times, detailed output, convergence criteria options, along with a succinct ability to switch between different algorithms.

A Conditional Approach for Modelling Patient Re-admissions in Hospital using a mixture of Coxian phase-type distributions incorporating Bayes Theorem

The number of elderly patients requiring hospitalisation in Europe is rising. With a greater proportion of elderly people in the population comes a greater demand for health services and, in particular, hospital care. Thus, with a growing number of elderly patients requiring hospitalisation competing with non-elderly patients for a fixed (and in some cases, decreasing) number of hospital beds, this results in much longer waiting times for patients, often with a less satisfactory hospital experience. However, if a better understanding of the recurring nature of elderly patient movements between the community and hospital can be developed, then it may be possible for alternative provisions of care in the community to be put in place and thus prevent readmission to hospital. The research in this paper aims to model the multiple patient transitions between hospital and community by utilising a mixture of conditional Coxian phase-type distributions that incorporates Bayes' theorem. For the purpose of demonstration, the results of a simulation study are presented and the model is applied to hospital readmission data from the Lombardy region of Italy.

Actuarial applications of multivariate phase-type distributions : model calibration and credibility

Thèse numérisée par la Division de la gestion de documents et des archives de l'Université de Montréal

Experimenting with the Coxian Phase-Type Distribution to Uncover Suitable Fits

In the past few decades, Coxian phase-type distributions have become increasingly more popular as a means of representing survival times. In healthcare, they are considered suitable for modelling the length of stay of patients in hospital and more recently for modelling the patient waiting times in Accident and Emergency Departments. The Coxian phase-type distribution has not only been shown to provide a good representation of real survival data, but its interpretation seems reasonably initiative to the medical experts. The drawback, however, is fitting the distribution to the data. There have been many attempts at accurately estimating the Coxian phase-type parameters. This paper wishes to examine the most promising of the approaches reported in the literature to determine the most accurate. Three performance measures are introduced to assess the fitting process of the algorithms along with the likelihood values and AIC to examine the goodness of fit and complexity of the model. Previous research suggests that the fitting process is strongly influenced by the initial parameter estimates and the data itself being quite variable. To overcome this, one experiment in this research paper will use the same initial parameter values for each estimation and perform the fits on the data simulated from a Coxian phase-type distribution with known parameters.