The mixture cure model has been widely applied to survival data in which a fraction of the observations never experience the event of interest, despite long-term follow-up. In this paper, we study the Cox proportional hazards mixture cure model where the covariate effects on the distribution of uncured subjects’ failure time may jump when a covariate exceeds a change point. The nonparametric maximum likelihood estimation is used to obtain the semiparametric estimates. We employ a two-step computational procedure involving the Expectation-Maximization algorithm to implement the estimation. The consistency, convergence rate and asymptotic distributions of the estimators are carefully established under technical conditions and we show that the change point estimator is n consistency. The m out of n bootstrap and the Louis algorithm are used to obtain the standard errors of the estimated change point and other regression parameter estimates, respectively. We also contribute a test procedure to check the existence of the change point. The finite sample performance of the proposed method is demonstrated via simulation studies and real data examples.
GandyAJensenULtkebohmertC.A cox model with a change-point applied to an actuarial problem.Brazil J Probabil Statistics2005;
19: 93–109.
22.
JensenULütkebohmertC.A cox-type regression model with change-points in the covariates.Lifetime Data Analysis2008;
14: 267–285.
23.
LuoXBoyettJM.Estimations of a threshold parameter in cox regression. Commun Stats-Theory Meth1997;
26: 2329–2346.
24.
PonsO.Estimation in a cox regression model with a change-point according to a threshold in a covariate. Ann Stat2003;
31: 442–463.
25.
DupuyJF.Detecting change in a hazard regression model with right-censoring. J Stat Plan Inference2009;
139: 1578–1586.
26.
ChenBEJiangWTuD.A hierarchical Bayes model for biomarker subset effects in clinical trials. Computat Stat Data Analys2014;
71: 324–334.
27.
BouazizONuelG.A change-point model for detecting heterogeneity in ordered survival responses.Stat Meth Med Res2018;
27: 3595–3611.
28.
FarewellVT.The use of mixture models for the analysis of survival data with long-term survivors. Biometrics1982;
38: 1041–1046.
29.
TaylorJMG.Semi-parametric estimation in failure time mixture models. Biometrics1995;
51: 899–907.
30.
PengYDearKBGDenhamJW.A generalized f mixture model for cure rate estimation. Stat Med1998;
17: 813–830.
31.
PengYDearKBG.A nonparametric mixture model for cure rate estimation.Biometrics2000;
56: 237–243.
32.
KlebanovLBYakovlevAY.A new approach to testing for sufficient follow-up in cure-rate analysis. J Stat Plann Inference2007;
137: 3557–3569.
33.
LuW.Maximum likelihood estimation in the proportional hazards cure model. Ann Inst Statist Math2008;
60: 545C574.
34.
MaS.Cure model with current status data. Statist Sinica2009;
19: 233–249.
35.
LuW.Efficient estimation for an accelerated failure time model with a cure fractio. Statist Sinica2010;
20: 661C674.
36.
WangLDuPLiangH.Two-component mixture cure rate model with spline estimated nonparametric component. Biometrics2012;
68: 726–735.
37.
JiangWSunHPengY.Prediction accuracy for the cure probabilities in mixture cure models.Stat Meth Med Res2017;
26: 2029–2041.
38.
ShiHYinG.Landmark cure rate models with time-dependent covariates.Stat Meth Med Res2017;
26: 2042–2054.
39.
BalakrishnanNBaruiSMilienosF.Proportional hazards under Conway–Maxwell-Poisson cure rate model and associated inference. Stat Meth Med Res2017;
26: 2055–2077.
40.
ShiXMaSHuangY. Promoting sign consistency in the cure model estimation and selection. Stat Meth Med Res 2020; 29: 15–28.
41.
OthusMLiYTiwariR.Change-point cure models with application to estimating the change-point effect of age of diagnosis among prostate cancer patients. J Appl Stat2012;
39: 901–911.
42.
KosorokMRSongR.Inference under right censoring for transformation models with a change-point based on a covariate threshold. Ann Stat2007;
35: 957–989.
43.
AndersenPKGillRD.Cox’s regression model for counting processes: a large sample study. Ann Stat1982;
10: 1100–1120.
44.
LouisTA.Finding the observed information matrix when using the em algorithm. J Royal Stat Soc: Ser B (Methodological)1982;
44: 226–233.
45.
DümbgenL.On nondifferentiable functions and the bootstrap.Probabil Theory Related Fields1993;
95: 125–140.
46.
ShaoJ.Bootstrap sample size in nonregular cases. Proc Am Math Soc1994;
122: 1251–1262.
47.
BickelPJSakovA.On the choice of m in the m out of n bootstrap and confidence bounds for extrema. Stat Sinica2008;
18: 967–985.
48.
XuGSenBYingZ.Bootstrapping a change-point cox model for survival data.Electronic J Stat2014;
8: 1345–1379.
49.
WangSZhangJAndrewBL.A Bayesian normal mixture accelerated failure time spatial model and its application to prostate cancer.Stat Meth Med Res2016;
25: 793–806.
50.
HenleySJWardEMScottS, et al.
Annual report to the nation on the status of cancer, Part I: National cancer statistics.Cancer2020;
126: 2225–2249.
51.
HamiltonASPotoskyALGillilandFD, et al.
Racial and ethnic differences in advanced-stage prostate cancer: the prostate cancer outcomes study.J Natl Cancer Inst2001;
93: 388–395.
52.
ColtmanCAHershmanDLCrowleyJJ, et al.
Racial disparities in cancer survival among randomized clinical trials patients of the Southwest Oncology Group. J Natl Cancer Inst2009;
101: 984–992.
53.
ZhangSSunDHeCZ, et al.
A bayesian semi-parametric model for colorectal cancer incidences.Stat Med2006;
25: 285–309.
54.
VicenteGCJorgeLBDipakKD.A new class of regression model for a bounded response with application in the study of the incidence rate of colorectal cancer. Stat Meth Med Res2020;
29: 2015–2033.
55.
HowladerNNooneAMKrapchoM, et al. SEER cancer statistics review, 1975–2017. Bethesda, MD: National Cancer Institute, https://seer.cancer.gov/csr/1975_2017/, based on November 2019 SEER data submission, posted to the SEER web site, April 2020.
56.
FryzlewiczP.Wild binary segmentation for multiple change-point detection. Ann Stat2014;
42: 2243–2281.
57.
ScottAJKnottM.A cluster analysis method for grouping means in the analysis of variance. Biometrics1974; 30: 507–512.
58.
HarchaouiZLévy-LeducC.Multiple change-point estimation with a total variation penalty. J Am Stat Assoc2010;
105: 1480–1493.
59.
LiJJinB.Multi-threshold accelerated failure time model. Ann Stat2018;
46: 2657–2682.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
