A Metaheuristic Approach to Parameter Estimation of a Flexible Parametric Mixture Cure Rate Model with Interval-Censored Data

P. Borges, M. Campos


A flexible parametric mixture cure model, called bi-lognormal cure rate model or simply BLN model, is defined and studied. The BLN model can be effectively used to analyze survival dataset in the presence of long-term survivors, especially when the dataset presents the underlying phenomenon of latent competing risks or when there is evidence that a bimodal hazard function is appropriated to described it, which are advantages over other cure rate models found in the literature. We discuss the maximum likelihood estimation for the model parameters considering interval-censored data through the differential evolution algorithm that is a nature-inspired computing metaheuristic used for global optimization of functions defined in multidimensional spaces. This approach is also used because the likelihood function of the model is multimodal and the direct application of gradient methods in this case is not ideal, since such methods are local search methods with a high chance of getting stuck at a local maximum when the starting point is chosen outside the basin of attraction of a global maximum. In addition, a simulation study was implemented to compare the performance of differential evolution algorithm with the performance of the Newton-Raphson algorithm in terms of bias, root mean square error, and the coverage probability of the asymptotic confidence intervals for the parameters. Finally, an application of the BLN model to real data is presented to illustrate that it can provide a better fit than other mixture cure rate models.


Polyhazard model; Lognormal distribution; differential evolution algorithm; Newton-Raphson algorithm; generalized Turnbull's nonparametric estimator

Full Text:



J. Boag, Maximum likelihood estimates of the proportion of patients cured by cancer therapy, Journal of the Royal Statistical Society B, vol. 11, pp. 15-53, 1949.

J. Berkson and R. Gage, Survival curve for cancer patients following treatment, Journal of the American Statistical Association, vol. 42, pp. 501-515,1952.

V. Farewell, A model for a binary variable with time censored observations, Biometrics, vol. 64, pp. 43-46, 1977.

V. Farewell, The use mixture models for the analysis of survival data with long term survivors, Biometrics, vol. 38, pp. 1041-1046, 1982.

R. Maller and X. Zhou, Testing for the presence of immune or cured individuals in censored survival data, Biometrics, vol. 51, pp. 1197-1205, 1995.

M. Roman, F. Louzada-Neto, V. Cancho, and J. Leite, A new long-term survival distribution for cancer data, Journal of Data Science, vol. 10, pp. 241-258, 2012.

J. Mazucheli, F. Louzada-Neto, and J. Achcar, The polysurvival model with long-term survivors, Brazilian Journal of Probability and Statistics, vol. 26, pp. 313-324, 2012.

D. Cox and D. Oakes, Analysis of Survival Data. London, Chapman and Hall,1984.

A. Yakovlev and A. Tsodikov, Stochastic Models of Tumor Latency and Their Biostatistical Applications. Singapore: World Scientific Publishers, 1996.

M. Chen, J. Ibrahim, and D. Sinha, A new bayesian model for survival data with a surviving fraction, Journal of the American Statistical Association, vol. 94, pp. 909-914, 1999.

F. Cooner, S. Banerjee, B. Carlin, and D. Sinha, Flexible cure rate modelling under latent activation schemes, Journal of the American Statistical Association, vol. 102, pp. 560-572, 2007.

J. Rodrigues, V. Cancho, M. Castro, and F. Louzada-Neto, On the unification of the long-term survival models, Statistics and Probability Letters, vol. 79, pp. 753-759, 2009.

P. Borges, J. Rodrigues, and N. Balakrishnan, Correlated destructive generalized power series cure rate models and associated inference with an application to a cutaneous melanoma data, Computational Statistics and Data Analysis, vol. 56, pp. 1703-1713, 2012.

P. Borges, J. Rodrigues, F. Louzada-Neto, and N. Balakrishnan, A cure rate survival model under a hybrid latent activation scheme, Statistical Methods in Medical Research, vol. 25, pp. 838-856, 2016.

P. Borges, Em algorithm-based likelihood estimation for a generalized gompertz regression model in presence of survival data with long-term survivors: an application to uterine cervical cancer data, Journal of Statistical Computation and Simulation, vol. 87, pp. 1-11, 2017.

J. Balka, A. F. Desmond, and P. D. McNicholas, Review and implementation of cure models based on first hitting times for wiener processes, Lifetime Data Analysis, vol. 15, pp. 147-176, 2009.

F. Louzada-Neto, Polyhazard regression models for lifetime data, Biometrics, vol. 55, pp. 1281-1285, 1999.

J. Kalbeisch and R. Prentice, The Statistical Analysis of Failure Time Data. New York, Wiley, 1980.

J. Berger and D. Sun, Bayesian analysis for the poly-weibull distribution, Journal of the American Statistical Association, vol. 88, pp. 1412-1418, 1993.

S. Basu, A. Basu, and C. Mukhopadhyay, Bayesian analysis for masked system failure data using non-identical weibull models, Statistics and Probability Letters, vol. 78, pp. 255-275, 1999.

L. Kuo and T. Yang, Bayesian reliability modeling for masked system lifetime, Statistics and Probability Letters, vol. 47, pp. 229-241, 2000.

J. Mazucheli, F. Louzada-Neto, and J. Achcar, Bayesian inference for polyhazard models in the presence of covariates, Computational Statistics and Data Analysis, vol. 38, pp. 1-14, 2001.

F. Louzada-Neto, C. Andrade, and F. Almeida, On the non-identifiability

problem arising on the poly-weibull model, Communications in Statistics -

Simulation and Computation, vol. 33, pp. 541-552, 2004.

J. Fachini, E. Ortega, and F. Louzada-Neto, lnfluence diagnostics for polyhazard models in the presence of covariates, Statistical Methods and Applications, vol. 17, pp. 413-433, 2008.

R. Tsai and L. Hotta, Polyhazard models with dependent causes, Brazilian Journal of Probability and Statistics, vol. 27, pp. 357-376, 2013.

R. Tsai and L. Hotta, Polyhazard models with dependent causes, fitting distributions with the polyhazard model with dependence, Communications in Statistics - Theory and Methods, vol. 44, pp. 1886-1895, 2015.

N. Demiris, D. Lunn, and L. Sharples, Survival extrapolation using the polyweibull model, Statistical Methods in Medical Research, vol. 24, pp. 287-301, 2015.

J. Sun, The statistical analysis of interval-censored failure time data. New York, Chapman and Hall, 2006.

E. Hashimoto, E. Ortega, and G. Cordeiro, A new long-term survival model with interval-censored data, Sankhya, vol. 77, pp. 207-239, 2015.

S. Pal and N. Balakrishnan, Likelihood inference for com-poisson cure rate model with interval-censored data and weibull lifetimes, Statistical Methods in Medical Research, vol. 26, pp. 2093-2113, 2017.

S. Scolas, A. EL Ghouch, C. Lengrand, and A. Oulhaj, Variable selection in a flexible parametric mixture cure model with interval-censored data, Statistics in Medicine, vol. 35, pp. 1210-1225, 2016.

V. Calsavara, A. Rodrigues, R. Rocha, V. Tomazella, and F. Louzada-Neto, Defective regression models for cure rate modeling with interval-censored data, Biometrical Journal, vol. 61, pp. 841-859, 2019.

V. Chan and W. Meeker, A failure-time model for infant mortality and wearout failure models, IEEE Transactions on Reliability, vol. 48, pp. 377-387, 1999.

P. Garthwaite, I. Jolliffe, and B. Jones, Statistical Inference, 2nd ed. New York, Oxford University Press, 2002.

A. Khuri, Advanced Calculus with Applications, 2nd ed. New York, John Wiley and Sons, 2003.

R. Storn and K. Price, Differential evolution - a simple and efficient heuristic for global optimization over continuous spaces, Journal of Global Optimization, vol. 11, pp. 341-359, 1997.

DOI: https://doi.org/10.5540/tcam.2023.024.03.00535

Article Metrics

Metrics Loading ...

Metrics powered by PLOS ALM


  • There are currently no refbacks.

Trends in Computational and Applied Mathematics

A publication of the Brazilian Society of Applied and Computational Mathematics (SBMAC)


Indexed in:




Desenvolvido por:

Logomarca da Lepidus Tecnologia