Lognormal lifetimes and likelihood-based inference for flexible cure rate models based on COM-Poisson family

Recently, a new cure rate survival model has been proposed by considering the Conway-Maxwell Poisson distribution as the distribution of the competing cause variable. This model includes some of the well-known cure rate models discussed in the literature as special cases. Cancer clinical trials often lead to right censored data and so the EM algorithm can be used as an efficient tool for the estimation of the model parameters based on right censored data. By considering this Conway-Maxwell Poisson-based cure rate model and by assuming the lognormal distribution for the time-to-event variable, the steps of the EM algorithm are developed here for the estimation of the parameters of different cure rate survival models. The standard errors of the MLEs are obtained by inverting the observed information matrix. An extensive Monte Carlo simulation study is performed to illustrate the method of inference developed. Model discrimination between different cure rate models is addressed by the likelihood ratio test as well as by Akaike and Bayesian information criteria. Finally, the proposed methodology is illustrated with a real data on cutaneous melanoma.

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

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

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

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

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

[6]  R. Rigby,et al.  Generalized additive models for location, scale and shape , 2005 .

[7]  K. Yau,et al.  Mixture cure model with random effects for clustered interval‐censored survival data , 2011, Statistics in medicine.

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

[9]  Anthony Y. C. Kuk,et al.  A mixture model combining logistic regression with proportional hazards regression , 1992 .

[10]  Narayanaswamy Balakrishnan,et al.  Destructive weighted Poisson cure rate models , 2011, Lifetime data analysis.

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

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

[13]  C. C. Kokonendji,et al.  Semiparametric estimation for count data through weighted distributions , 2009 .

[14]  Debanjan Mitra,et al.  Left truncated and right censored Weibull data and likelihood inference with an illustration , 2012, Comput. Stat. Data Anal..

[15]  T. Mathew,et al.  Inferences on the means of lognormal distributions using generalized p-values and generalized confidence intervals , 2003 .

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

[17]  Heleno Bolfarine,et al.  Cure rate survival models with missing covariates: a simulation study , 2013 .

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

[19]  Narayanaswamy Balakrishnan,et al.  Correlated destructive generalized power series cure rate models and associated inference with an application to a cutaneous melanoma data , 2012, Comput. Stat. Data Anal..

[20]  J. Duistermaat,et al.  Distributions: Theory and Applications , 2010 .

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

[22]  Debanjan Mitra,et al.  Likelihood inference for lognormal data with left truncation and right censoring with an illustratio , 2011 .

[23]  Nicholas T. Longford.  Inference with the Lognormal Distribution , 2008 .

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

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

[26]  E. Crow,et al.  Lognormal Distributions: Theory and Applications , 1987 .

[27]  Joseph G. Ibrahim,et al.  A New Bayesian Model For Survival Data With a Surviving Fraction , 1999 .

[28]  Yu Shen,et al.  A Semiparametric Regression Cure Model for Interval-Censored Data , 2009, Journal of the American Statistical Association.

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

[30]  J M Taylor,et al.  Semi-parametric estimation in failure time mixture models. , 1995, Biometrics.

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

[32]  Narayanaswamy Balakrishnan,et al.  EM Algorithm-Based Likelihood Estimation for Some Cure Rate Models , 2012 .

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

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

[35]  Josemar Rodrigues,et al.  On the unification of long-term survival models , 2009 .

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

[37]  R. Rigby,et al.  Generalized Additive Models for Location Scale and Shape (GAMLSS) in R , 2007 .

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

[39]  Xue Li,et al.  Statistical inference for the common mean of two log-normal distributions and some applications in reliability , 2006, Comput. Stat. Data Anal..

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

[41]  J W Denham,et al.  A generalized F mixture model for cure rate estimation. , 1998, Statistics in medicine.