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.

[1]  Juan Eloy Ruiz-Castro,et al.  Modelling a reliability system governed by discrete phase-type distributions , 2008, Reliab. Eng. Syst. Saf..

[2]  L. Schmickler.  Meda: mixed erlang distributions as phase-type representations of empirical distribution functions. , 1992 .

[3]  W. R. Schucany,et al.  Minimum Distance and Robust Estimation , 1980 .

[4]  D. Cox.  A use of complex probabilities in the theory of stochastic processes , 1955, Mathematical Proceedings of the Cambridge Philosophical Society.

[5]  Marcel F. Neuts,et al.  Matrix-geometric solutions in stochastic models - an algorithmic approach , 1982 .

[6]  H. Akaike.  A new look at the statistical model identification , 1974 .

[7]  Malcolm J. Faddy,et al.  On inferring the number of phases in a coxian phase-type distribution , 1998 .

[8]  Hassan S. Bakouch,et al.  Probability, Markov chains, queues, and simulation , 2011 .

[9]  M. Faddy,et al.  A structured compartmental model for drug kinetics. , 1993, Biometrics.

[10]  Sally McClean,et al.  Analysing data on lengths of stay of hospital patients using phase‐type distributions , 1999 .

[11]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[12]  Adele H. Marshall,et al.  Simulating Coxian phase-type distributions for patient survival , 2009, Int. Trans. Oper. Res..

[13]  Sally McClean,et al.  Conditional phase-type distributions for modelling patient length of stay in hospital , 2003 .

[14]  M. J. Faddy,et al.  PENALISED MAXIMUM LIKELIHOOD ESTIMATION OF THE PARAMETERS IN A COXIAN PHASE-TYPE DISTRIBUTION , 2002 .

[15]  Adele H Marshall,et al.  Using Coxian Phase-Type Distributions to Identify Patient Characteristics for Duration of Stay in Hospital , 2004, Health care management science.

[16]  J. E. Ruiz-Castro,et al.  A Multiple‐Absorbent Markov Process in Survival Studies: Application to Breast Cancer , 2003 .

[17]  Reinhard Augustin,et al.  Characteristics of the COX-distribution , 1982, PERV.

[18]  Christophe Dutang,et al.  actuar: An R Package for Actuarial Science , 2008 .

[19]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[20]  M. A. Johnson,et al.  Selecting Parameters of Phase Distributions: Combining Nonlinear Programming, Heuristics, and Erlang Distributions , 1993, INFORMS J. Comput..

[21]  M. Faddy.  Examples of fitting structured phase–type distributions , 1994 .

[22]  H. D. Miller,et al.  The Theory Of Stochastic Processes , 1977, The Mathematical Gazette.

[23]  Ren Asmussen,et al.  Fitting Phase-type Distributions via the EM Algorithm , 1996 .

[24]  Alma Riska,et al.  Efficient fitting of long-tailed data sets into hyperexponential distributions , 2002, Global Telecommunications Conference, 2002. GLOBECOM '02. IEEE.