Correlated destructive generalized power series cure rate models and associated inference with an application to a cutaneous melanoma data

In this paper, we propose a new cure rate survival model, which extends the model of Rodrigues et al. (2011) by incorporating a structure of dependence between the initiated cells. To create the structure of the correlation between the initiated cells, we use an extension of the generalized power series distribution by including an additional parameter @r (the inflated-parameter generalized power series (IGPS) distribution, studied by Kolev and Minkova (2000)). It has a natural interpretation in terms of both a ''zero-inflated'' proportion and a correlation coefficient. In our approach, the number of initiated cells is assumed to follow the IGPS distribution. The IGPS distribution is a natural choice for modeling correlated count data that exhibit overdispersion. The primary advantage of this distributional assumption is that the correlation structure induced by the additional parameter @r results in a natural characterization of the association between the initiated cells. Moreover, it provides a simple and realistic interpretation for the biological mechanism of the occurrence of the event of interest as it includes a process of destruction of tumor cells after an initial treatment or the capacity of an individual exposed to irradiation to repair initiated cells that result in cancer being induced. This means that what is recorded is only the undamaged portion of the original number of initiated cells not eliminated by the treatment or repaired by the repair system of an individual. Parameter estimation of the proposed model is then discussed through the maximum likelihood estimation procedure. Finally, we illustrate the usefulness of the proposed model by applying it to real cutaneous melanoma data.

[1]  Joseph G Ibrahim,et al.  A General Class of Bayesian Survival Models with Zero and Nonzero Cure Fractions , 2005, Biometrics.

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

[3]  S T Rachev,et al.  A stochastic model of radiation carcinogenesis: latent time distributions and their properties. , 1993, Mathematical biosciences.

[4]  Leda D. Minkova,et al.  Inflated-Parameter Family of Generalized Power Series Distributions and Their Application in Analysis of Overdispersed Insurance Data , 2001 .

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

[6]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[7]  M. V. Aarset The null distribution for a test of constant versus «bathtub» failure rate , 1985 .

[8]  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.

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

[10]  Gleb Haynatzki,et al.  A new statistical model of tumor latency time , 2000 .

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

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

[13]  Ana Maria Barral,et al.  Immunological Studies in Malignant Melanoma : Importance of TNF and the Thioredoxin System , 2001 .

[14]  C W Chen,et al.  A stochastic two-stage carcinogenesis model: a new approach to computing the probability of observing tumor in animal bioassays. , 1991, Mathematical biosciences.

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

[16]  Ram C. Tripathi,et al.  Modified Power Series Distribution , 2014 .

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

[18]  Chin-Shang Li,et al.  Identifiability of cure models , 2001 .

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

[20]  David Collett Modelling Survival Data in Medical Research , 1994 .

[21]  J. Ibrahim,et al.  Bayesian Inference for Multivariate Survival Data with a Cure Fraction , 2002 .

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

[23]  M. May Bayesian Survival Analysis. , 2002 .

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

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

[26]  Bradley P Carlin,et al.  Flexible Cure Rate Modeling Under Latent Activation Schemes , 2007, Journal of the American Statistical Association.

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

[28]  Bradley P Carlin,et al.  Parametric Spatial Cure Rate Models for Interval‐Censored Time‐to‐Relapse Data , 2004, Biometrics.

[29]  N. Draper,et al.  Applied Regression Analysis , 1966 .

[30]  P. C. Consul,et al.  New class of location-parameter discrete probability distributions and their characterizations , 1990 .

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

[32]  D. Collet Modelling Survival Data in Medical Research , 2004 .

[33]  Vernon T. Farewell,et al.  Mixture models in survival analysis: Are they worth the risk? , 1986 .

[34]  Leda D. Minkova A GENERALIZATION OF THE CLASSICAL DISCRETE DISTRIBUTIONS , 2002 .

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