Unified multivariate survival model with a surviving fraction: an application to a Brazilian customer churn data

In this paper we propose a new lifetime model for multivariate survival data in presence of surviving fractions and examine some of its properties. Its genesis is based on situations in which there are m types of unobservable competing causes, where each cause is related to a time of occurrence of an event of interest. Our model is a multivariate extension of the univariate survival cure rate model proposed by Rodrigues et al. [37]. The inferential approach exploits the maximum likelihood tools. We perform a simulation study in order to verify the asymptotic properties of the maximum likelihood estimators. The simulation study also focus on size and power of the likelihood ratio test. The methodology is illustrated on a real data set on customer churn data.

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

[2]  L G Hanin Iterated birth and death process as a model of radiation cell survival. , 2001, Mathematical biosciences.

[3]  Francisco Louzada-Neto,et al.  Influence diagnostics for polyhazard models in the presence of covariates , 2008, Stat. Methods Appl..

[4]  M. Kendall Theoretical Statistics , 1956, Nature.

[5]  Jorge Alberto Achcar,et al.  Bayesian inference for polyhazard models in the presence of covariates , 2001 .

[6]  Christophe Mues,et al.  Mixture cure models in credit scoring: If and when borrowers default , 2012, Eur. J. Oper. Res..

[7]  Jorge Alberto Achcar,et al.  Mixture hazard models for lifetime data , 2002 .

[8]  Heleno Bolfarine,et al.  Modeling the presence of immunes by using the exponentiated-Weibull model , 2001 .

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

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

[11]  R. Kay The Analysis of Survival Data , 2012 .

[12]  F. Louzada,et al.  On the Distribution of the Minimum or Maximum of a Random Number of i.i.d. Lifetime Random Variables , 2012 .

[13]  D. Clayton A model for association in bivariate life tables and its application in epidemiological studies of familial tendency in chronic disease incidence , 1978 .

[14]  Sudipto Banerjee,et al.  Analysis of cure rate survival data under proportional odds model , 2011, Lifetime data analysis.

[15]  Daniel R. Jeske,et al.  Statistical Inference: An Integrated Approach , 2003 .

[16]  F. Louzada,et al.  The FGM Long-Term Bivariate Survival Copula Model: Modeling, Bayesian Estimation, and Case Influence Diagnostics , 2013 .

[17]  Francisco Louzada,et al.  The exponential‐Poisson model for recurrent event data: An application to a set of data on malaria in Brazil , 2015, Biometrical journal. Biometrische Zeitschrift.

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

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

[20]  Francisco Louzada,et al.  The complementary Weibull geometric distribution , 2014 .

[21]  Heleno Bolfarine,et al.  Comparing several accelerated life models , 1993 .

[22]  Francisco Louzada-Neto,et al.  Extended Hazard Regression Model for Reliability and Survival Analysis , 1997, Lifetime data analysis.

[23]  K. Yamaguchi Accelerated Failure-Time Regression Models with a Regression Model of Surviving Fraction: An Application to the Analysis of “Permanent Employment” in Japan , 1992 .

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

[25]  Francisco Louzada,et al.  A New Long-Term Survival Distribution for Cancer Data , 2012, Journal of Data Science.

[26]  Francisco Louzada,et al.  A multiple time scale survival model with a cure fraction , 2012 .

[27]  D. Oakes A Model for Association in Bivariate Survival Data , 1982 .

[28]  Francisco Louzada-Neto,et al.  The complementary exponential power lifetime model , 2011, Comput. Stat. Data Anal..

[29]  Dirk Van den Poel,et al.  Customer base analysis: partial defection of behaviourally loyal clients in a non-contractual FMCG retail setting , 2005, Eur. J. Oper. Res..

[30]  Francisco Louzada,et al.  On estimation and influence diagnostics for a bivariate promotion lifetime model based on the FGM copula: a fully bayesian computation , 2013 .

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

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

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

[34]  Vitor A. A. Marchi,et al.  The Complementary Exponentiated Exponential Geometric Lifetime Distribution , 2013 .

[35]  F Louzada-Neto,et al.  Polyhazard Models for Lifetime Data , 1999, Biometrics.

[36]  Francisco Louzada-Neto,et al.  A general hazard model for lifetime data in the presence of cure rate , 2011 .

[37]  N. Chatterjee,et al.  A Bivariate Cure‐Mixture Approach for Modeling Familial Association in Diseases , 2001, Biometrics.

[38]  F. Louzada,et al.  The complementary exponential power series distribution , 2013 .

[39]  D. Karlis An EM algorithm for multivariate Poisson distribution and related models , 2003 .

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

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