{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,11,10]],"date-time":"2026-11-10T18:44:25Z","timestamp":1794336265232,"version":"4.9.0"},"reference-count":30,"publisher":"Oxford University Press (OUP)","issue":"1","license":[{"start":{"date-parts":[[2026,1,1]],"date-time":"2026-01-01T00:00:00Z","timestamp":1767225600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/pages\/standard-publication-reuse-rights"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,1,20]]},"abstract":"<jats:title>Summary<\/jats:title>\n                  <jats:p>In non-small cell lung cancer (NSCLC) clinical trials, tumor burden (TB) is a key longitudinal biomarker for assessing treatment effects. Typically, standard-of-care (SOC) therapies and some novel interventions initially decrease TB; however, many patients subsequently experience an increase\u2014indicating disease progression\u2014while others show a continuous decline. In patients with an eventual TB increase, the change point marks the onset of progression and must occur before the time of the event. To capture these dynamics, we propose a cure-rate joint model that integrates time-to-event and longitudinal TB data and distinguishes change-point and non-change-point trajectory groups. For the change-point group, our approach flexibly estimates an individualized change point by leveraging time-to-event information. We use a Monte Carlo Expectation-Maximization (MCEM) algorithm for efficient parameter estimation. Simulation studies show that our model can recover diverse disease progression patterns and yield robust marginal TB outcome estimates while accommodating censoring complexity. When applied to a Phase 3 NSCLC trial comparing cemiplimab monotherapy to SOC, the treatment group demonstrates delayed tumor regrowth and consistently lower TB over time, highlighting the clinical utility of our approach. The implementation code is publicly available on https:\/\/github.com\/quyixiang\/JoCuR.<\/jats:p>","DOI":"10.1093\/biostatistics\/kxag025","type":"journal-article","created":{"date-parts":[[2026,6,25]],"date-time":"2026-06-25T11:48:04Z","timestamp":1782388084000},"source":"Crossref","is-referenced-by-count":1,"title":["Cure rate joint model for time-to-event data and longitudinal tumor burden with potential change points"],"prefix":"10.1093","volume":"27","author":[{"ORCID":"https:\/\/orcid.org\/0009-0009-7485-7895","authenticated-orcid":false,"given":"Yixiang","family":"Qu","sequence":"first","affiliation":[{"name":"Department of Biostatistics, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Campus Box 7420, 135 Dauer Drive, Chapel Hill, NC 27599-7420,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-6112-9030","authenticated-orcid":false,"given":"Ethan M","family":"Alt","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Campus Box 7420, 135 Dauer Drive, Chapel Hill, NC 27599-7420,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Weibin","family":"Zhong","sequence":"additional","affiliation":[{"name":"Biostatistics and Data Management, Regeneron Pharmaceuticals, 777 Old Saw Mill River Road, Tarrytown, NY 10591,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jeen","family":"Liu","sequence":"additional","affiliation":[{"name":"Biostatistics and Data Management, Regeneron Pharmaceuticals, 777 Old Saw Mill River Road, Tarrytown, NY 10591,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Chenguang","family":"Wang","sequence":"additional","affiliation":[{"name":"Biostatistics and Data Management, Regeneron Pharmaceuticals, 777 Old Saw Mill River Road, Tarrytown, NY 10591,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Joseph G","family":"Ibrahim","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Campus Box 7420, 135 Dauer Drive, Chapel Hill, NC 27599-7420,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,7,28]]},"reference":[{"key":"2026102110234656000_kxag025-B1","doi-asserted-by":"publisher","first-page":"94","DOI":"10.1186\/s12874-020-00976-2","article-title":"Bayesian joint modelling of longitudinal and time to event data: a methodological review","volume":"20","author":"Alsefri","year":"2020","journal-title":"BMC Med Res Methodol"},{"key":"2026102110234656000_kxag025-B2","doi-asserted-by":"publisher","first-page":"e70021","DOI":"10.1002\/sim.70021","article-title":"Jointly modeling time-to-event and longitudinal data with individual-specific change points: a case study in modeling tumor burden","volume":"44","author":"Alt","year":"2025","journal-title":"Stat Med"},{"key":"2026102110234656000_kxag025-B3","doi-asserted-by":"publisher","first-page":"2113","DOI":"10.1002\/sim.8889","article-title":"Change-point joint model for identification of Plateau of activity in early phase trials","volume":"40","author":"Altzerinakou","year":"2021","journal-title":"Stat Med"},{"key":"2026102110234656000_kxag025-B4","doi-asserted-by":"publisher","first-page":"334","DOI":"10.1093\/ije\/dyu262","article-title":"Joint modelling of repeated measurement and time-to-event data: an introductory tutorial","volume":"44","author":"Asar","year":"2015","journal-title":"Int J Epidemiol"},{"key":"2026102110234656000_kxag025-B5","doi-asserted-by":"publisher","first-page":"1256","DOI":"10.1177\/0962280219853599","article-title":"Joint longitudinal and time-to-event cure models for the assessment of being cured","volume":"29","author":"Barbieri","year":"2020","journal-title":"Stat Methods Med Res"},{"key":"2026102110234656000_kxag025-B6","doi-asserted-by":"publisher","first-page":"827","DOI":"10.1097\/EDE.0000000000000723","article-title":"Bayesian piecewise linear mixed models with a random change point: an application to BMI rebound in childhood","volume":"28","author":"Brilleman","year":"2017","journal-title":"Epidemiology"},{"key":"2026102110234656000_kxag025-B7","doi-asserted-by":"publisher","first-page":"3502","DOI":"10.1177\/0962280218808821","article-title":"Joint longitudinal and time-to-event models for multilevel hierarchical data","volume":"28","author":"Brilleman","year":"2019","journal-title":"Stat Methods Med Res"},{"key":"2026102110234656000_kxag025-B8","doi-asserted-by":"publisher","first-page":"1255","DOI":"10.1016\/j.cmpb.2012.08.013","article-title":"smcure: an r-package for estimating semiparametric mixture cure models","volume":"108","author":"Cai","year":"2012","journal-title":"Comput Methods Programs Biomed"},{"key":"2026102110234656000_kxag025-B9","doi-asserted-by":"publisher","first-page":"165","DOI":"10.1177\/1536867X1301300112","article-title":"Joint modeling of longitudinal and survival data","volume":"13","author":"Crowther","year":"2013","journal-title":"Stata J"},{"key":"2026102110234656000_kxag025-B10","doi-asserted-by":"publisher","first-page":"305","DOI":"10.1111\/biom.12537","article-title":"Using the SAEM algorithm for mechanistic joint models characterizing the relationship between nonlinear PSA kinetics and survival in prostate cancer patients","volume":"73","author":"Desm\u00e9e","year":"2017","journal-title":"Biometrics"},{"key":"2026102110234656000_kxag025-B11","doi-asserted-by":"publisher","first-page":"5786","DOI":"10.1002\/sim.3380","article-title":"A random change point model for assessing variability in repeated measures of cognitive function","volume":"27","author":"Dominicus","year":"2008","journal-title":"Stat Med"},{"key":"2026102110234656000_kxag025-B12","doi-asserted-by":"publisher","first-page":"299","DOI":"10.1111\/j.1467-985X.2007.00514.x","article-title":"A joint latent class changepoint model to improve the prediction of time to graft failure","volume":"171","author":"Garre","year":"2008","journal-title":"J R Stat Soc Ser A Stat Soc"},{"key":"2026102110234656000_kxag025-B13","doi-asserted-by":"publisher","first-page":"2181","DOI":"10.1002\/sim.6141","article-title":"Joint modeling of survival and longitudinal non-survival data: current methods and issues. Report of the DIA Bayesian joint modeling working group","volume":"34","author":"Gould","year":"2015","journal-title":"Stat Med"},{"key":"2026102110234656000_kxag025-B14","doi-asserted-by":"publisher","first-page":"91","DOI":"10.1016\/S0167-9473(02)00148-2","article-title":"Bayesian and profile likelihood change point methods for modeling cognitive function over time","volume":"42","author":"Hall","year":"2003","journal-title":"Comput Stat Data Anal"},{"key":"2026102110234656000_kxag025-B15","doi-asserted-by":"publisher","first-page":"e076316","DOI":"10.1136\/bmj-2023-076316","article-title":"The estimands framework: a primer on the ich e9 (r1) addendum","volume":"384","author":"Kahan","year":"2024","journal-title":"BMJ"},{"key":"2026102110234656000_kxag025-B16","doi-asserted-by":"publisher","first-page":"1452","DOI":"10.1111\/bcp.15200","article-title":"Modelling the association between biomarkers and clinical outcome: an introduction to nonlinear joint models","volume":"88","author":"Kerioui","year":"2022","journal-title":"Br J Clin Pharmacol"},{"key":"2026102110234656000_kxag025-B17","doi-asserted-by":"publisher","first-page":"236","DOI":"10.2307\/2533329","article-title":"Change points in the series of T4 counts prior to AIDS","volume":"51","author":"Kiuchi","year":"1995","journal-title":"Biometrics"},{"key":"2026102110234656000_kxag025-B18","doi-asserted-by":"publisher","first-page":"615","DOI":"10.2307\/2290194","article-title":"Hierarchical Bayes models for the progression of HIV infection using longitudinal CD4 T-cell numbers","volume":"87","author":"Lange","year":"1992","journal-title":"J Am Stat Assoc"},{"key":"2026102110234656000_kxag025-B19","doi-asserted-by":"publisher","first-page":"547","DOI":"10.1093\/biostatistics\/3.4.547","article-title":"The joint modeling of a longitudinal disease progression marker and the failure time process in the presence of cure","volume":"3","author":"Law","year":"2002","journal-title":"Biostatistics"},{"key":"2026102110234656000_kxag025-B20","doi-asserted-by":"publisher","first-page":"637","DOI":"10.1111\/j.1467-9868.2008.00704.x","article-title":"Fully exponential Laplace approximations for the joint modelling of survival and longitudinal data","volume":"71","author":"Rizopoulos","year":"2009","journal-title":"J R Stat Soc Ser B Stat Methodol"},{"key":"2026102110234656000_kxag025-B21","doi-asserted-by":"publisher","first-page":"592","DOI":"10.1016\/S0140-6736(21)00228-2","article-title":"Cemiplimab monotherapy for first-line treatment of advanced non-small-cell lung cancer with pd-l1 of at least 50%: a multicentre, open-label, global, phase 3, randomised, controlled trial","volume":"397","author":"Sezer","year":"2021","journal-title":"Lancet"},{"key":"2026102110234656000_kxag025-B22","doi-asserted-by":"publisher","first-page":"434","DOI":"10.1016\/j.immuni.2018.03.014","article-title":"Regulation and function of the pd-l1 checkpoint","volume":"48","author":"Sun","year":"2018","journal-title":"Immunity"},{"key":"2026102110234656000_kxag025-B23","doi-asserted-by":"publisher","first-page":"232","DOI":"10.1002\/sim.4107","article-title":"Joint modeling of survival time and longitudinal data with subject-specific changepoints in the covariates","volume":"30","author":"Tapsoba","year":"2011","journal-title":"Stat Med"},{"key":"2026102110234656000_kxag025-B24","first-page":"809","article-title":"Joint modeling of longitudinal and time-to-event data: an overview","author":"Tsiatis","year":"2004","journal-title":"Stat Sin"},{"key":"2026102110234656000_kxag025-B25","author":"U.S. Food and Drug Administration","year":"2018"},{"key":"2026102110234656000_kxag025-B26","author":"U.S. Food and Drug Administration","year":"2021"},{"key":"2026102110234656000_kxag025-B27","doi-asserted-by":"publisher","first-page":"684","DOI":"10.1016\/j.csda.2012.07.024","article-title":"Change point models for cognitive tests using semi-parametric maximum likelihood","volume":"57","author":"van den Hout","year":"2013","journal-title":"Comput Stat Data Anal"},{"key":"2026102110234656000_kxag025-B28","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1155\/2012\/640153","article-title":"Analysis of longitudinal and survival data: joint modeling, inference methods, and issues","volume":"2012","author":"Wu","year":"2012","journal-title":"J Probab Stat"},{"key":"2026102110234656000_kxag025-B29","doi-asserted-by":"crossref","first-page":"330","DOI":"10.2307\/2533118","article-title":"A joint model for survival and longitudinal data measured with error","volume":"53","author":"Wulfsohn","year":"1997","journal-title":"Biometrics"},{"key":"2026102110234656000_kxag025-B30","first-page":"835","article-title":"Joint longitudinal-survival-cure models and their application to prostate cancer","author":"Yu","year":"2004","journal-title":"Stat Sin"}],"container-title":["Biostatistics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/biostatistics\/article-pdf\/27\/1\/kxag025\/70230057\/kxag025.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/biostatistics\/article-pdf\/27\/1\/kxag025\/70230057\/kxag025.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,10,21]],"date-time":"2026-10-21T14:23:53Z","timestamp":1792592633000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/biostatistics\/article\/doi\/10.1093\/biostatistics\/kxag025\/8743781"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,1]]},"references-count":30,"journal-issue":{"issue":"1","published-print":{"date-parts":[[2026,1,20]]}},"URL":"https:\/\/doi.org\/10.1093\/biostatistics\/kxag025","relation":{},"ISSN":["1465-4644","1468-4357"],"issn-type":[{"value":"1465-4644","type":"print"},{"value":"1468-4357","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,1]]},"published":{"date-parts":[[2026,1]]},"article-number":"kxag025"}}