A penalized likelihood approach for mixture cure models

Cure models have been developed to analyze failure time data with a cured fraction. For such data, standard survival models are usually not appropriate because they do not account for the possibility of cure. Mixture cure models assume that the studied population is a mixture of susceptible individuals, who may experience the event of interest, and non-susceptible individuals that will never experience it. Important issues in mixture cure models are estimation of the baseline survival function for susceptibles and estimation of the variance of the regression parameters. The aim of this paper is to propose a penalized likelihood approach, which allows for flexible modeling of the hazard function for susceptible individuals using M-splines. This approach also permits direct computation of the variance of parameters using the inverse of the Hessian matrix. Properties and limitations of the proposed method are discussed and an illustration from a cancer study is presented.

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

[2]  Paul H. C. Eilers,et al.  Flexible smoothing with B-splines and penalties , 1996 .

[3]  David R. Cox,et al.  Regression models and life tables (with discussion , 1972 .

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

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

[6]  Chin-Shang Li,et al.  A semi‐parametric accelerated failure time cure model , 2002, Statistics in medicine.

[7]  Yingwei Peng,et al.  Estimating baseline distribution in proportional hazards cure models , 2003, Comput. Stat. Data Anal..

[8]  A I Goldman,et al.  Survivorship analysis when cure is a possibility: a Monte Carlo study. , 1984, Statistics in medicine.

[9]  Odd O. Aalen,et al.  Modelling Heterogeneity in Survival Analysis by the Compound Poisson Distribution , 1992 .

[10]  Yingwei Peng,et al.  A new estimation method for the semiparametric accelerated failure time mixture cure model , 2007, Statistics in medicine.

[11]  D. Y. Fong,et al.  Estimating the proportion of cured patients in a censored sample , 2005, Statistics in medicine.

[12]  K. Dear,et al.  A Nonparametric Mixture Model for Cure Rate Estimation , 2000, Biometrics.

[13]  J. Boag,et al.  Maximum Likelihood Estimates of the Proportion of Patients Cured by Cancer Therapy , 1949 .

[14]  D. Cox.  Regression Models and Life-Tables , 1972 .

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

[16]  O. Aalen,et al.  Heterogeneity in survival analysis. , 1988, Statistics in medicine.

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

[18]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[19]  K K Yau,et al.  Long-term survivor mixture model with random effects: application to a multi-centre clinical trial of carcinoma. , 2001, Statistics in medicine.

[20]  D Commenges,et al.  A penalized likelihood approach for arbitrarily censored and truncated data: application to age-specific incidence of dementia. , 1998, Biometrics.

[21]  Christopher U. Jones,et al.  Local control of carcinoma of the tonsil by radiation therapy: an analysis of patterns of fractionation in nine institutions. , 1995, International journal of radiation oncology, biology, physics.

[22]  N. Gordon,et al.  Application of the theory of finite mixtures for the estimation of 'cure' rates of treated cancer patients. , 1990, Statistics in medicine.

[23]  K Borch-Johnsen,et al.  Heterogeneity models of disease susceptibility, with application to diabetic nephropathy. , 1994, Biometrics.

[24]  A K Manatunga,et al.  Modelling survival data with a cured fraction using frailty models , 2001, Statistics in medicine.

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

[26]  Daniel Commenges,et al.  A penalized likelihood approach for an illness-death model with interval-censored data: application to age-specific incidence of dementia. , 2002, Biostatistics.

[27]  Benoit Liquet,et al.  Choice between Semi‐parametric Estimators of Markov and Non‐Markov Multi‐state Models from Coarsened Observations , 2007 .

[28]  J. Gamel,et al.  Comparison of parametric and non-parametric survival methods using simulated clinical data. , 1997, Statistics in medicine.

[29]  Standard errors for the cox proportional hazards cure model , 2001 .

[30]  M. E. Ghitany,et al.  Exponential mixture models with long-term survivors and covariates , 1994 .

[31]  Xian Zhou,et al.  A new parametric model for survival data with long-term survivors. , 2004, Statistics in medicine.