EM Algorithm-Based Likelihood Estimation for Some Cure Rate Models

In the recent work of Rodrigues et al. (2009), a flexible cure rate survival model was developed by assuming the number of competing causes of the event of interest to follow the Conway-Maxwell Poisson distribution. This model includes as special cases some of the well-known cure rate models discussed in the literature. As the data obtained from cancer clinical trials are often subject to right censoring, the expectation maximization (EM) algorithm can be used as a powerful and efficient tool for the estimation of the model parameters based on right censored data. In this paper, the cure rate model developed by Rodrigues et al. (2009) is considered and assuming the time-to-event to follow the exponential distribution, exact likelihood inference is developed based on the EM algorithm. The inverse of the observed information matrix is used to compute the standard errors of the maximum likelihood estimates (MLEs). An extensive Monte Carlo simulation study is performed to illustrate the method of inference developed here. Finally, the proposed methodology is illustrated with real data on cutaneous melanoma.

[1]  T. Louis Finding the Observed Information Matrix When Using the EM Algorithm , 1982 .

[2]  Narayanaswamy Balakrishnan,et al.  COM–Poisson cure rate survival models and an application to a cutaneous melanoma data , 2009 .

[3]  D. Cox,et al.  Analysis of Survival Data. , 1985 .

[4]  A Yu Yakovlev,et al.  Stochastic Models of Tumor Latency and Their Biostatistical Applications , 1996 .

[5]  N. Balakrishnan,et al.  Generalized gamma frailty model , 2006, Statistics in medicine.

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

[7]  J. P. Sy,et al.  Estimation in a Cox Proportional Hazards Cure Model , 2000, Biometrics.

[8]  K. Liang,et al.  Asymptotic Properties of Maximum Likelihood Estimators and Likelihood Ratio Tests under Nonstandard Conditions , 1987 .

[9]  Joseph Berkson,et al.  Survival Curve for Cancer Patients Following Treatment , 1952 .

[10]  K. Lange A gradient algorithm locally equivalent to the EM algorithm , 1995 .

[11]  Galit Shmueli,et al.  Conjugate Analysis of the Conway-Maxwell-Poisson Distribution , 2006 .

[12]  R. Maller,et al.  Survival Analysis with Long-Term Survivors , 1996 .

[13]  G. McLachlan,et al.  The EM algorithm and extensions , 1996 .

[14]  W. Nelson Statistical Methods for Reliability Data , 1998 .

[15]  Jim E. Griffin,et al.  A Bayesian Partition Model for Customer Attrition , 2001 .

[16]  T. Minka,et al.  A useful distribution for fitting discrete data: revival of the Conway–Maxwell–Poisson distribution , 2005 .

[17]  A D Tsodikov,et al.  A stochastic model of hormesis. , 1993, Mathematical biosciences.

[18]  Joseph G. Ibrahim,et al.  Bayesian Survival Analysis , 2004 .

[19]  G. McLachlan,et al.  The EM Algorithm and Extensions: Second Edition , 2008 .

[20]  V. Farewell,et al.  The use of mixture models for the analysis of survival data with long-term survivors. , 1982, Biometrics.

[21]  Narayanaswamy Balakrishnan,et al.  Connections of the Poisson weight function to overdispersion and underdispersion , 2008 .

[22]  Paul Janssen,et al.  One-sided tests in shared frailty models , 2008 .

[23]  J J Shuster,et al.  Parametric versus non-parametric methods for estimating cure rates based on censored survival data. , 1992, Statistics in medicine.

[24]  Joseph G. Ibrahim,et al.  Cure rate models: A unified approach , 2005 .