- Research article
- Open Access
- Open Peer Review
A simulation model of colorectal cancer surveillance and recurrence
BMC Medical Informatics and Decision Makingvolume 14, Article number: 29 (2014)
Approximately one-third of those treated curatively for colorectal cancer (CRC) will experience recurrence. No evidence-based consensus exists on how best to follow patients after initial treatment to detect asymptomatic recurrence. Here, a new approach for simulating surveillance and recurrence among CRC survivors is outlined, and development and calibration of a simple model applying this approach is described. The model’s ability to predict outcomes for a group of patients under a specified surveillance strategy is validated.
We developed an individual-based simulation model consisting of two interacting submodels: a continuous-time disease-progression submodel overlain by a discrete-time Markov submodel of surveillance and re-treatment. In the former, some patients develops recurrent disease which probabilistically progresses from detectability to unresectability, and which may produce early symptoms leading to detection independent of surveillance testing. In the latter submodel, patients undergo user-specified surveillance testing regimens. Parameters describing disease progression were preliminarily estimated through calibration to match five-year disease-free survival, overall survival at years 1–5, and proportion of recurring patients undergoing curative salvage surgery from one arm of a published randomized trial. The calibrated model was validated by examining its ability to predict these same outcomes for patients in a different arm of the same trial undergoing less aggressive surveillance.
Calibrated parameter values were consistent with generally observed recurrence patterns. Sensitivity analysis suggested probability of curative salvage surgery was most influenced by sensitivity of carcinoembryonic antigen assay and of clinical interview/examination (i.e. scheduled provider visits). In validation, the model accurately predicted overall survival (59% predicted, 58% observed) and five-year disease-free survival (55% predicted, 53% observed), but was less accurate in predicting curative salvage surgery (10% predicted; 6% observed).
Initial validation suggests the feasibility of this approach to modeling alternative surveillance regimens among CRC survivors. Further calibration to individual-level patient data could yield a model useful for predicting outcomes of specific surveillance strategies for risk-based subgroups or for individuals. This approach could be applied toward developing novel, tailored strategies for further clinical study. It has the potential to produce insights which will promote more effective surveillance—leading to higher cure rates for recurrent CRC.
Roughly two-thirds of the more than 140,000 Americans diagnosed with colorectal cancer (CRC) each year will be treated surgically with curative intent . Approximately one-third of these will experience recurrence of the original disease or will develop a new primary (metachronous) CRC [2, 3]. Median survival for those experiencing recurrence of original disease is around two years . Ultimately, nearly 50,000 patients in the U.S. die each year from CRC .
Most patients treated curatively are placed on some type of surveillance program involving periodic follow-up testing to detect preclinical recurrence. For patients who will experience recurrence, prognosis, though generally poor, may be improved if detection occurs prior to symptom onset, particularly if surgical resection of metastatic disease is possible . Reviews of the relatively small number of trials comparing two or more specific follow-up strategies have suggested that more intensive strategies tend to increase survival at five years, detecting recurrence about six months earlier than less intensive strategies [6–8] at a point where salvage surgery with curative intent is more likely to occur (10.7% vs. 5.7%; p = .0002) . However, the strategies tested, the populations studied, and the study periods varied significantly between trials [6–8]. As such, drawing prescriptive conclusions regarding best practices on the basis of these data is difficult.
Meanwhile, in economic terms, CRC surveillance consumes significant resources. U.S. investigators found that five-year follow up can cost as much as US$16,942 per patient . Another group in Europe reported surveillance costs of US$9,011  to US$59,841  per detection of a recurrence resulting in attempted curative salvage surgery.
Not surprisingly, consensus is lacking among expert panels on how best to follow these patients [12–18]. New surveillance trials are in progress, but results may be several years away [19–21]. When complete, these trials will provide valuable information but will have been able to examine only a small fraction of possible surveillance strategies.
Given the potential impact on quantity and quality of life and on health care costs, better tools are needed for informing decisions around postsurgical surveillance of colorectal cancer patients. By simulating the dynamics of recurrence in a population of patients, a realistic computer simulation model could function as a virtual laboratory within which an unlimited number of experiments comparing hypothetical surveillance strategies could be run in silico within a compressed timeframe.
Several models have incorporated dynamics of the adenoma-carcinoma sequence in order to compare the effectiveness of hypothetical screening strategies in patients without history of CRC [22–32]; these include three models used in the National Cancer Institute’s comparative modeling effort CISNET (Cancer Intervention and Surveillance Modeling Network) . Fewer models though have simulated the events following diagnosis and treatment of CRC in order to compare postsurgical surveillance strategies [34–37]. None has captured the dynamics of recurrence in a way that accounts for disease progression during diagnostic delay and that considers the full range of possible metastatic sites.
Capturing the dynamics of CRC recurrence is a major methodological challenge mainly because of the difficulty of estimating parameters describing progression of recurring disease amid the censoring caused by medical and surgical interventions. In order to create a realistic model which allows assessment of any hypothetical surveillance strategy, one must be able to account for disease progression amid diagnostic delay. Here, we describe a new approach to modeling the interaction between natural history of CRC recurrence and early detection of recurrence through surveillance testing—an approach designed to allow the simulation of any potential surveillance strategy. We introduce a basic model we have developed which applies this approach, preliminarily estimate disease progression parameters by calibration based on published outcomes from a classic surveillance trial, and offer a quantitative validation of the model.
Overview of approach
The model itself is comprised of two interacting submodels: a continuous-time disease progression submodel and a discrete time Markov submodel of surveillance testing and re-treatment. In the disease progression submodel, the exact time to earliest recurrence detectability is pre-determined for each simulated patient who will recur based on random draws from an exponential probability distribution. A pair of formulas—both functions of the time to earliest detectability--determines the timing of the transition to unresectability and to the point of symptom onset. Once these pre-determinations are made, the discrete-time Markov surveillance and re-treatment submodel simulates scheduled visits for surveillance testing of asymptomatic patient as frequently as every three months. This submodel references the pre-determined timeline of disease progression to determine whether asymptomatic recurrences are detectable by testing during surveillance visits, and whether recurrences are considered potentially resectable versus unresectable at the time they are discovered. Recurrences may alternatively be detected in the interval between surveillance visits as a result of symptoms which prompt individuals to seek earlier care.
To simulate the impact of any combination and schedule of surveillance tests, the disease progression submodel must be capable—in the extreme—of simulating disease progression in the absence of surveillance. Since most data describing recurrence tends to be “contaminated” by the effect of testing to detect asymptomatic recurrence and by subsequent intervention, it is difficult to directly estimate certain key parameters underlying a disease progression submodel. Therefore, we use a calibration approach to estimate these parameters while using the surveillance and re-treatment submodel described above to control for the effect of a known surveillance regimen on observed disease progression. Once these disease progression parameters are estimated through calibration, the schedule of follow-up tests embodied in the surveillance and re-treatment submodel can be changed to simulate any number of hypothetical follow-up regimens. In the case of the current work, we use the surveillance regimen followed in the intervention arm of a published CRC follow-up trial to calibrate disease progression parameters. We then change the surveillance regimen to be consistent with the control arm of the same study and compare the predicted to the observed outcomes.
Modeling disease progression
Each individual i is randomly assigned a future recurrence status (i.e. will recur at some point or not) and site based on random number draws from a uniform probability distribution with range 0 to 1. Anatomic sites of recurrence are described in terms of local tumor status (none, locoregional, anastomotic only, or intraluminal/metachronous) and metastasis status (none, liver, lung, or multiple/other organs). All recurrences are considered to “begin” at the time of initial treatment since, by definition, there would have been at least some small amount of cancerous tissue which remained despite treatment. Each individual who will recur is assigned a future time Di which represents the earliest point in time at which the recurrence of patient i will be detectable through surveillance testing. Di is assigned randomly to these individuals from a continuous exponential probability distribution with hazard rate rd, leading to a declining incidence of recurrence with increasing time from initial treatment at the population level . The time Ui at which the case is no longer amenable to curative salvage surgery (point of unresectability) is defined as
Where xdu is the baseline (i.e. where Di = 0) interval between Di and Ui and ru is the rate of change in this window. A positive value of Xdu would represent the minimum possible width of the window between earliest detectability and unresectability given ru > 0. Beyond the point of unresectability, the patient will only be eligible for palliation and will be assigned a life expectancy accordingly.
A non-trivial number of recurrences are discovered through work-ups precipitated by patient symptomatic presentations in the intervals between scheduled surveillance visits . This proportion can be expected to grow in inverse proportion to the vigilance of follow-up regimens, and should thus be modeled. The point of recurrence symptom onset Si for patient i is defined as
Where xds is the baseline (Di = 0) interval between Di and Si and rs is the rate of change in this window. A positive value of Xds would represent the minimum width of the window between earliest detectability and symptom onset given rs > 0. The normally distributed error term ϵds (with mean = 0 and standard deviation = σds) is included to ensure sufficient variation in symptom onset so that some patients will develop symptoms before the point of unresectability and some will develop symptoms after that point. Figure 1 illustrates the two possible scenarios for progression of recurring individuals through clinically important points in disease natural history.
Modeling surveillance and re-treatment
We used TreeAge Pro 2012 (Williamstown, MA) to create an individual-based Markov submodel of surveillance testing and re-treatment. Individual-based Markov models allow simulation of phenomena which involve transitions of heterogeneous individuals between states over time. During each discrete three-month time increment—or cycle—each of a series of hypothetical individuals may be in only one state at a time. Within a state, individuals may experience a series of probabilistic events which determine whether they will remain in their current state during the next cycle or transition to a new state. In our model, a series of simulated patients starts each cycle in one of five states: no known recurrence, recurrence curatively treated, recurrence palliatively treated, deceased due to cancer, or deceased due to other causes (See state transition diagram in Figure 2). Three months is used as a cycle length since this represents the shortest recommended inter-visit interval among published recommendations for post-surgical surveillance in CRC patients [12, 15, 16, 40]. Note that only surveillance testing of asymptomatic patients occurs based on discrete time steps. Natural history events are measured in continuous time so that simulated patients will reach the point of detectability (Di), the point of unresectability (Ui), and the point of symptom onset (Si) during the intervals between surveillance visits.
The time horizon of five years for surveillance was chosen because this is the timeframe over which virtually all recurrences manifest , and the majority of surveillance recommendations focus on this five-year horizon [12, 15, 16, 40]. All individuals begin the first cycle—representing the period immediately following curative resection—in the no known recurrence state.
Each of the living states in the surveillance and re-treatment submodel is represented by a decision subtree (see Additional file 1: Figure S1). Whether or not a test is performed during a particular three-month cycle is dictated by a surveillance testing schedule matrix with a row for each three month cycle and a column for each possible testing modality. Using 1’s and 0’s, this matrix can specify any combination and timing of laboratory tests (e.g. carcinoembryonic antigen, or CEA), imaging studies (e.g. abdominal/pelvic CT), and provider visits. Positivity or negativity of each test given during a cycle depends on the current stage of the model compared to Di, the location of any detectable tumors, and the sensitivity and specificity of the test being employed. If any of the tests performed in a given cycle is positive, it is assumed that sufficient further workup using other available tests is undertaken to correctly identify detectable recurrence or otherwise rule out recurrence within that cycle with a combined sensitivity and specificity of 100%. Once recurrence is diagnosed (whether based on detection through scheduled surveillance or stemming from symptoms developing between scheduled visits), patients are assigned a life expectancy. The exact value depends on whether their recurrence has been diagnosed before or after reaching the point of unresectability (Table 1).
Calibration to estimate disease progression parameters
Values for the six disease progression parameters described above (rd, xdu, ru, xds, rs, and σds) are impossible or impractical to observe, and their values can be expected to differ based on disease and patient characteristics of the groups of interest. Therefore, we used a calibration process to estimate values for these parameters. We systematically ran the model with many thousands of different combinations of the unknown parameters in order to find the combination which yielded model outcomes most consistent with a real world data source . As a data source, we used published data from a classic trial of an intensive versus a minimal follow-up strategy for patients who had undergone curative resection for CRC from 1987–1990 at an Italian Center . Among all available trials of surveillance strategies, this study by Pietra et al. offered the most complete data needed for the exercise. Patients in the trial had been treated surgically for Astler-Coller Stage B or C CRC without adjuvant chemotherapy (additional details describing study participants in Table 2). The focus of the trial was on local recurrences; however, data for subjects with metastatic recurrences was also reported. We used data from the intensive follow-up arm to calibrate our model; accordingly, we refer to this group as the calibration group. To account for risk of non-cancer death among individuals who have not been diagnosed with recurrence, we calibrated an additional parameter (m) representing group-specific background (i.e. from non-CRC causes) mortality probability per cycle. We considered this a more reliable method to estimate background mortality in a small group compared to using life table-based estimates.
We selected a set of seven observable outcomes reported by Pietra et al. to serve as calibration targets (note that these are distinct from the seven parameters estimated through calibration). The discrepancies between observed and model-predicted values of these targets were compared for each combination of values of the disease progression and background mortality parameters. The targets used were disease-free survival at five years (DFS5), overall survival at years one through five (OS1 – OS5), and the proportion of patients undergoing eventual curative salvage surgery following recurrence. Group characteristics and recurrence rates, surveillance testing schedule, and observed values for each of the calibration targets are shown for the intensive (calibration) group in the middle column of Table 2.
Searching the parameter space was achieved using a series of increasingly narrow grid searches. Specifically, each of the seven parameters to be calibrated was given an initially broad range (shown in the third column of Table 3) which was divided into three to eight evenly-spaced intervals. Values from across the initial range were systematically sampled so that all possible combinations of parameters were tried (i.e. a grid search of the parameter space ). Each run of the model simulated a cohort of 10,000 patients. Model--generated values of each calibration target were compared to values reported by Pietra et al. Candidate disease progression parameter sets which yielded target values > +/− one percentage point (or > +/− two percentage points if no values varying less than one percentage point were found) from the observed target value for DFS5 were eliminated from consideration. Next, model-generated and observed values for OS5 were compared in this same manner to further reduce the set of candidate parameter values, followed by OS4, OS3, proportion of patients undergoing curative salvage surgery, OS2, and finally OS1. OS2 and OS1 were considered last since these were most subject to influence from small numbers of events. The result was a small group of remaining parameter sets. Based on the range of candidate parameter values suggested by this process, a second round of calibration was performed, this time with a narrower set of candidate parameter values. This process was repeated for subsequent rounds until no additional gains in model fit—as judged by changes in the value of a goodness-of-fit statistic—were achieved. The goodness-of-fit statistic used for this purpose was a simple unweighted sum of squared differences between observed and predicted values for the calibration targets. The lowest goodness-of-fit value from among all sets of parameters from one round of calibration was compared to the lowest value from the previous round. Once no further improvement was seen between rounds, no further rounds of calibration were undertaken. After the final calibration round, the remaining parameter sets collectively contained a small range of candidate values for each parameter. A single best set was chosen which contained the values closest to the midpoint of the final range for each parameter (in the same order of priority which was applied throughout calibration: DFS5, OS5, OS4, OS3, probability of curative salvage operation, OS2, then OS1). This approach of using multiple, increasingly narrow grid searches was preferred on efficiency grounds to an approach of performing a single grid search with a large number of closely-spaced intervals for each parameter. With the latter approach, the number of parameters combinations would have been astronomical. Ultimately, we performed nine successive rounds of calibration, examining a total of 276,960 distinct parameter combinations (10,000 patients simulated for each parameter combination) in order to arrive at the final parameter estimates.
Test sensitivity and specificity, as well as life expectancies, were estimated based on values from published studies (Table 1). Since evidence-based estimates were obtainable, calibration of these parameters was not necessary. Sensitivity and specificity values shown are for individual tests in isolation. It was assumed that metastases to multiple locations or to single sites besides the liver or lungs would be inoperable. It was also assumed that lung metastases were 50% as likely as liver metastases to be operable [55–57].
One-way sensitivity analyses were performed on all non-calibrated model parameters according to the ranges shown in Table 1. The primary outcomes of interest in sensitivity analyses were OS5 and proportion of subjects undergoing curative salvage surgery for recurrence of CRC. The same randomization algorithm seed value was used throughout sensitivity analysis to minimize the impact of stochastic variability on comparisons with baseline model outputs.
Once values for the six parameters describing the progression of CRC recurrence and the parameter representing background mortality for those without known recurrence had been estimated by calibration, the model was used to predict the outcomes that would be expected in the minimal follow-up arm of the Pietra trial; we will refer to this arm as the validation group. The primary purpose of this validation was to assess the model’s ability to predict cancer-related outcomes for a group with known disease characteristics and non-cancer mortality risk under a specified surveillance regimen. Keeping all calibrated disease progression parameters constant, the follow-up regimen modeled was modified to be consistent with that undergone by subjects in the validation group (see Table 2). Because the Pietra trial was not sufficiently large that subject characteristics were evenly distributed between groups through randomization, the proportion of patients relapsing and the distribution of metastatic sites in the model (see Table 2) were adjusted in validation to match what was reported for the minimal follow-up group. Specifically, probability of recurrence was adjusted from 39.4% to 40.4%, and distribution of metastatic disease between liver/lung/other was adjusted from 26.7%/0.0%/73.3% to 14.2%/4.8%/81.0%. Pietra et al. reported a substantial difference in DFS5 between the intensive and minimal surveillance groups (68% and 53% respectively). Only 1.0 percentage point of this difference was explained by a difference in recurrence rate. Since eventual cancer-related deaths do not affect disease-free survival calculation (rather, it is the diagnosis of recurrence which precedes death that counts), it was assumed that the remaining difference in DFS5 was attributable to differences in background mortality. We accordingly adjusted the five-year background cumulative probability of death from other causes in the validation group by 14.0%, equating to a 3.4 percentage point increment in annual mortality from other causes in the validation group compared to the calibration group. No adjustments were made to the structure of the model itself.
Final calibrated parameter estimates are shown in the fourth column of Table 3. Figure 3 compares observed values from the intensive follow-up (calibration) group in the Pietra trial to the model-generated outputs for the same group using the final best-fitting parameter set. As seen from the figure, predicted overall survival is slightly overestimated early in the five-year time horizon, while five-year disease-free survival and the rate of curative salvage surgery are slightly underestimated by the model. Additional file 2: Figure S2 depicts the disease progression milestones for 20 simulated patients with recurrence based on the calibrated parameters.
In order to examine the sensitivity of key model results to changes in non-calibrated parameters, we conducted a deterministic sensitivity analysis, varying all non-calibrated parameters in the fitted model across the ranges shown in the third column of Table 1. The outcomes of interest for sensitivity analysis were the proportion of patients undergoing eventual salvage surgery (denominator = all patients initially treated for CRC, whether or not they eventually recurred) and OS5.
Figure 4a depicts sensitivity analysis results for the outcome representing the proportion of patients undergoing curative salvage surgery. The most influential parameter was test sensitivity of the CEA assay, with a 15 percentage point increase in sensitivity to 79% resulting in a 7.9% increase in proportion undergoing curative salvage surgery, and a 15 percentage point decrease to 49% resulting in a 7.4% decrease in proportion undergoing curative salvage surgery. Varying all other parameters over the specified ranges yielded variation in curative salvage surgeries of within 5% of the base case.
Figure 4b shows sensitivity analysis results for the outcome of overall survival at 5 years. Life expectancy values had the greatest impact on this outcome, though no parameter proved highly influential.
In order to assess the predictive value of the model for a similar population under a different follow-up regimen, we ran the model for the minimal follow-up arm of the Pietra trial and compared model-predicted with observed outcomes. This comparison is shown in Figure 5. The largest difference between observed and predicted outcomes for the validation group was in the proportion of subjects who would eventually undergo curative salvage surgery for recurrence. The model overestimated this proportion at 10%, while only 6% of patients in reality underwent curative salvage surgery in the validation group.
We have introduced a new approach for simulating colorectal cancer surveillance and recurrence and have described a model which applies this approach. We estimated a best-fitting disease progression parameter set for this model by comparing iteratively generated model outputs to observed outcomes in one arm of a classic randomized controlled trial comparing two follow-up regimens (calibration). We then validated the model by demonstrating that it closely approximated several key outcomes for a separate group of patients followed less intensively in the same trial (validation). Application of this approach in subsequent based on larger, individual-level data sets could facilitate more personalized post-surgical surveillance of CRC survivors, ultimately translating into fewer low-yield tests and more lives saved through early detection of recurrence.
The approach of developing a generic disease progression submodel and calibrating its parameters based on outcomes data was taken because of the lack of biological data needed to directly estimate these parameters. The natural history behavior of recurring colorectal cancer is difficult to observe directly and is extremely heterogeneous. This heterogeneity likely stems, at least partly, from the greater genotypic diversity of genomically unstable tumor tissue which has been exposed to the stresses of treatment.
Though the specific values for the parameter estimates pertain to the patients in the Pietra trial and should not be construed to broadly represent disease behavior in the context of more recent treatment norms, they are consistent with generally-recognized patterns of CRC recurrence. Specifically, the estimates imply that those who recur later after initial treatment (higher values of Di) tend to exhibit a more indolent course of disease (as represented by a positive value of rds) and are more likely to be considered curable when diagnosed (as represented by a positive value of rdu) [53, 58–60].
In sensitivity analysis, the proportion of patients undergoing curative salvage surgery was not highly sensitive to any single parameter. This is partly explained by the fact that approximately six of ten individuals did not recur; including these individuals in the denominator somewhat attenuates the variation in this outcome. However, even if only patients who recurred are included in the denominator, no parameter influenced the outcome by more than 20%. An interesting pattern did emerge. Probability of salvage surgery was most directly influenced by changes in sensitivity of tests which were not specific to an anatomic site (i.e. CEA and clinical exam/interview). In the simplest terms, these “tests” cast a wide net, with the potential to generate a true positive for any recurrence, regardless of site. This observation reinforces the rationale behind the current pursuit of newer, highly sensitive circulating markers for recurrence.
Paradoxically, proportion of curative salvage surgeries moved in inverse proportion to the specificity of those tests which targeted an anatomic site (e.g. hepatic ultrasound). This finding suggests that false positives for such tests actually provide an indirect benefit by precipitating further workups which diagnose true recurrences in other locations. This paradox underscores the importance of considering quality-of-life impact of testing as well as cost considerations. With these issues in mind, our model has been designed so that future versions may be used to compare the costs per quality-adjusted life year of alternative surveillance regimens. Specifically, the surveillance and re-treatment submodel has the ability to track costs due to surveillance testing, work-up of positive surveillance tests, and subsequent curative or palliative treatment for each simulated individual at each point in time. Quality-adjustment of life expectancy will be achieved by applying health state utility parameters to assigned life expectancies.
The fact that overall survival of CRC survivors in general is only modestly affected by changes in many of the parameters associated with diagnostic test characteristics and differential survival suggests an important implication. Because most survivors will not recur, and many of those who do recur have a narrow window within which early detection may change prognosis, developing tools to better identify those at greatest risk and to individualize their surveillance will likely yield greater benefit than applying generalized surveillance recommendations to all survivors. With calibration to larger, individual-based data sets covering patients with a range of risk factor profiles, it is hoped that future iterations of the model will provide a means to compare the effectiveness of proposed surveillance regimens for risk-based subgroups of patients and to even suggest optimized regimens for individuals.
Multiple models have been developed to evaluate CRC screening strategies in healthy populations. These models simulate the sequence by which benign polyps transform to adenocarcinomas, and by which these adenocarcinomas grow and invade healthy tissue [22–32]. Some of these screening models have spawned research questions involving disease natural history [28, 61], have informed development of U.S. Preventive Services Task Force (USPSTF) guidelines on colorectal cancer screening , and have been applied by the Centers for Medicare and Medicaid Services (CMS) to compare the effectiveness of CRC screening strategies [63, 64].
Only a few models have addressed the issue of surveillance of curatively treated CRC patients for early detection of recurrence. Two of these have not accounted for disease progression during diagnostic delay in a manner that would allow realistic assessment of novel surveillance strategies [34, 35]. A separate group of investigators  used a sophisticated mechanistic, organ-level simulation to compare different treatment and surveillance strategies in patients with hepatic metastases. While the model does account for disease progression amid diagnostic delay, it does not provide information useful in determining optimal follow-up for non-hepatic recurrences. Also, the approach taken may be less feasible than a more population-based approach when considering metastases to other sites about which less detailed data are available.
A final, recent example of a model examining CRC surveillance is provided by Erenay et al. . These investigators used longitudinal patient data to estimate natural history parameters related to formation and progression of metachronous colorectal tumors in patients who were status post curative resection. Though well-done, this model does not address recurrence of original cancers (which represent over 90% of relapses ). Furthermore, it is not validated against any other data source.
The approach we have described uses calibration of a disease progression submodel interacting with a surveillance and re-treatment submodel to estimate disease progression parameters while controlling for the impact of surveillance on disease natural history. The effect is to estimate group-specific parameters for a model describing disease progression in the absence of surveillance. Once calibrated, novel combinations and schedules of surveillance tests can be coupled with the disease progression model to predict outcomes for a group with similar underlying characteristics. Further calibration of the model to patient-level outcomes data containing individual risk factors will enable prediction of outcomes for populations or individuals with specific risk profiles, paving the way for more personalized and effective surveillance and the possibility for earlier detection of recurrence at a stage where cure is more likely.
The chief limitation of this work stems from the small size of the patient groups used in calibration and validation, and from the lack of individual-level data. With larger groups, the probability that underlying disease behavior and background mortality would be similar between populations used for calibration and populations for whom outcome predictions were being sought would be greater. We chose the Pietra trial because of the wealth of outcomes data it provided and the between-study-arm similarity in subject characteristics compared to other available published studies. In validation, we attempted to alleviate limitations related to group size by assigning the appropriate background mortality rates and distribution of metastatic sites to match those reported for the validation group. In a scenario of making actual population-level predictions, one might not have the luxury of knowing these characteristics of the group of interest. However, for our purposes, these adjustments can be seen as strengthening the internal validation of the model by controlling for known differences between groups. In the future, individual-level data describing more diverse groups will be used to provide more generalizable parameter estimates and a higher-resolution understanding of the effects of different risk factors. With larger sample sizes, far less variability can be assumed in metastasis pattern and background mortality between calibration populations and populations for which predictions are to be made.
A minor limitation relates to the fact that we did not explicitly model the occurrence of second primary (metachronous) CRC’s. Because no subjects in the calibration group were diagnosed with metachronous CRC, calibrating the model to simultaneously simulate incidence of metachronous CRC and recurrence of previously-treated disease was not a feasible option. Fortunately, the impact of this limitation on the integrity of the model validation should be minimal since only one individual in the validation group was reported to have experienced metachronous CRC. Because metachronous CRC (normally representing between 1.6%  and 7.4% [7, 37, 65] of CRC recurrences) behaves as a primary cancer and is associated with higher probability of cure compared to recurrent disease, it will be important to explicitly model separate natural history processes for these two types of events in future iterations of the model.
We do not consider the fact that the study is based on data gathered in the 1980’s and 1990’s to be a limitation. The purpose of this work was not to make prescriptive assertions about survivorship care of CRC patients today, but rather to demonstrate a method and preliminarily validate it using publicly-available data. That said, some useful generalizations are possible regarding the impact of diagnostic test characteristics and the importance of developing the ability to better customize follow-up of CRC survivors as elaborated above.
Scheduled post-operative surveillance of patients curatively treated for colorectal cancer allows the detection of asymptomatic recurrence which can create the possibility of cure for some patients. We have described and preliminarily validated a simulation model, based on a novel approach, which allows comparison of important patient outcomes under hypothetical combinations and schedules of common surveillance tests.
Next steps include adapting the model to a discrete event simulation framework which will allow calibration to individual-level patient data in which subject risk factors as well as precise dates of testing events and clinical milestones are known. Calibration and validation with larger, individual-level data sets representing patients with a more diverse range of demographic and disease characteristics is planned. Use of richer data sets (e.g. from combined clinical trials, or based on cancer registry data) will enable the use of future versions of the model to compare the effectiveness and cost-effectiveness of published surveillance recommendations for subgroups based on cancer stage, location (e.g. rectum versus colon), age, comorbidity status, presence of molecular markers, adjuvant treatment received, and other considerations. In addition, the model could be adapted to utilize an optimization algorithm to generate candidate follow-up strategies for such subgroups—allowing the more efficient design of meaningful clinical trials. Future adaptations could also power interactive clinical decision aids for planning the management of individual patients, an approach that itself could be validated in a randomized trial.
JR is a Preventive Medicine physician and health services researcher in the Department of Family Medicine and Community Health at Case Western Reserve University School of Medicine. His main research focus involves simulation models of preventive interventions in cancer.
KMA is a practicing colorectal surgeon working at the University Hospital of North Norway. He serves as a research manager at the Norwegian National Centre for Integrative Care and Telemedicine, and was the principal investigator in a recently published randomized trail assessing colorectal cancer surveillance.
CYK is a researcher at the Massachusetts General Hospital Institute for Technology Assessment. He is an expert in methods for parameter estimation using calibration in complex models involving cancer.
NJM is the Dr. Lester E. Coleman, Jr. Chair Professor of Cancer Research and Therapeutics and Chief of the Division of Hematology and Oncology at University Hospitals Case Medical Center/Case Western Reserve University in Cleveland, OH. He serves as Associate Director of for Clinical Research at the Case Comprehensive Cancer Center, is Co-Chair of the National Cancer Institute Gastrointestinal Cancer Steering Committee, and serves on the Board of Directors of the American Society of Clinical Oncology (ASCO).
MWK is Chair of the Department of Quantitative Health Sciences at Cleveland Clinic in Cleveland, OH. His research is primarily focused on the development, validation, and use of prediction models in cancer.
QH is a graduate student in Biostatistics at Case Western Reserve University.
XA is a graduate student in Biostatistics at Case Western Reserve University.
GSC is a Professor of Medicine-Gastroenterology at University Hospitals Case Medical Center/Case Western Reserve University in Cleveland, OH. He directs the Office of Comparative Effectiveness Research at the Case Clinical and Translational Science Collaborative and co-leads the Cancer Prevention, Control and Population Research Program at Case Comprehensive Cancer Center.
Jemal A, Siegel R, Ward E, Hao Y, Xu J, Thun MJ: Cancer statistics, 2009. CA Cancer J Clin. 2009, 59 (4): 225-249. 10.3322/caac.20006.
Bohm B, Schwenk W, Hucke HP, Stock W: Does methodic long-term follow-up affect survival after curative resection of colorectal carcinoma?. Dis Colon Rectum. 1993, 36 (3): 280-286. 10.1007/BF02053511.
Scholefield JH, Steele RJ: Guidelines for follow up after resection of colorectal cancer. Gut. 2002, 51 (Suppl 5): V3-V5.
Cunningham D, Atkin W, Lenz HJ, Lynch HT, Minsky B, Nordlinger B, Starling N: Colorectal cancer. Lancet. 2010, 375 (9719): 1030-1047. 10.1016/S0140-6736(10)60353-4.
Kjeldsen BJ, Kronborg O, Fenger C, Jorgensen OD: A prospective randomized study of follow-up after radical surgery for colorectal cancer. Br J Surg. 1997, 84 (5): 666-669. 10.1002/bjs.1800840523.
Jeffery M, Hickey BE, Hider PN: Follow-up strategies for patients treated for non-metastatic colorectal cancer. Cochrane Database Syst Rev. 2007, 1: CD002200-
Tjandra JJ, Chan MK: Follow-up after curative resection of colorectal cancer: a meta-analysis. Dis Colon Rectum. 2007, 50 (11): 1783-1799. 10.1007/s10350-007-9030-5.
Figueredo A, Rumble RB, Maroun J, Earle CC, Cummings B, McLeod R, Zuraw L, Zwaal C: Follow-up of patients with curatively resected colorectal cancer: a practice guideline. BMC Cancer. 2003, 3: 26-
Virgo KS, Vernava AM, Longo WE, McKirgan LW, Johnson FE: Cost of patient follow-up after potentially curative colorectal cancer treatment. JAMA. 1995, 273 (23): 1837-1841. 10.1001/jama.1995.03520470045030.
Bleeker WA, Mulder NH, Hermans J, Otter R, Plukker JT: Value and cost of follow-up after adjuvant treatment of patients with Dukes’ C colonic cancer. Br J Surg. 2001, 88 (1): 101-106. 10.1046/j.1365-2168.2001.01638.x.
Audisio RA, Setti-Carraro P, Segala M, Capko D, Andreoni B, Tiberio G: Follow-up in colorectal cancer patients: a cost-benefit analysis. Ann Surg Oncol. 1996, 3 (4): 349-357. 10.1007/BF02305664.
Desch CE, Benson AB 3rd, Somerfield MR, Flynn PJ, Krause C, Loprinzi CL, Minsky BD, Pfister DG, Virgo KS, Petrelli NJ: Colorectal cancer surveillance: 2005 update of an American Society of Clinical Oncology practice guideline. J Clin Oncol. 2005, 23 (33): 8512-8519. 10.1200/JCO.2005.04.0063.
Engstrom PF, Arnoletti JP, 3rd Benson AB, Chen YJ, Choti MA, Cooper HS, Covey A, Dilawari RA, Early DS, Enzinger PC, Fakih MG, Fleshman J, Fuchs C, Grem JL, Kiel K, Knol JA, Leong LA, Lin E, Mulcahy MF, Rao S, Ryan DP, Saltz L, Shibata D, Skibber JM, Sofocleous C, Thomas J, Venook AP, Willett C: NCCN Clinical Practice Guidelines in Oncology: rectal cancer. J Natl ComprCanc Netw. 2009, 7 (8): 838-881.
Engstrom PF, Arnoletti JP, 3rd Benson AB, Chen YJ, Choti MA, Cooper HS, Covey A, Dilawari RA, Early DS, Enzinger PC, Fakih MG, Fleshman J, Fuchs C, Grem JL, Kiel K, Knol JA, Leong LA, Lin E, Mulcahy MF, Rao S, Ryan DP, Saltz L, Shibata D, Skibber JM, Sofocleous C, Thomas J, Venook AP, Willett C: NCCN Clinical Practice Guidelines in Oncology: colon cancer. J Natl ComprCanc Netw. 2009, 7 (8): 778-831.
Winawer S, Fletcher R, Rex D, Bond J, Burt R, Ferrucci J, Ganiats T, Levin T, Woolf S, Johnson D, Kirk L, Litin S, Simmang C: Colorectal cancer screening and surveillance: clinical guidelines and rationale-Update based on new evidence. Gastroenterology. 2003, 124 (2): 544-560. 10.1053/gast.2003.50044.
Anthony T, Simmang C, Hyman N, Buie D, Kim D, Cataldo P, Orsay C, Church J, Otchy D, Cohen J, Perry WB, Dunn G, Rafferty J, Ellis CN, Rakinic J, Fleshner P, Stahl T, Gregorcyk S, Ternent C, 3rd Kilkenny JW, Whiteford M: Practice parameters for the surveillance and follow-up of patients with colon and rectal cancer. Dis Colon Rectum. 2004, 47 (6): 807-817. 10.1007/s10350-004-0519-x.
Rex DK, Kahi CJ, Levin B, Smith RA, Bond JH, Brooks D, Burt RW, Byers T, Fletcher RH, Hyman N, Johnson D, Kirk L, Lieberman DA, Levin TR, O'Brien MJ, Simmang C, Thorson AG, Winawer SJ: Guidelines for colonoscopy surveillance after cancer resection: a consensus update by the American Cancer Society and the US Multi-Society Task Force on Colorectal Cancer. Gastroenterology. 2006, 130 (6): 1865-1871. 10.1053/j.gastro.2006.03.013.
Davila RE, Rajan E, Baron TH, Adler DG, Egan JV, Faigel DO, Gan SI, Hirota WK, Leighton JA, Lichtenstein D, Qureshi WA, Shen B, Zuckerman MJ, VanGuilder T, Fanelli RD: ASGE guideline: colorectal cancer screening and surveillance. Gastrointest Endosc. 2006, 63 (4): 546-557. 10.1016/j.gie.2006.02.002.
Wille-Jorgensen P, Laurberg S, Pahlman L, Carriquiry L, Lundqvist N, Smedh K, Svanfeldt M, Bengtson J: An interim analysis of recruitment to the COLOFOL trial. Colorectal Dis. 2009, 11 (7): 756-758. 10.1111/j.1463-1318.2008.01668.x.
University hospitals Southampton NHS Foundation Trust: Follow-Up Study of Patients Who Have Undergone Surgery for Stage I, Stage II, or Stage III Colorectal Cancer. The FACS study. 2011,http://www.clinicaltrials.gov/ct2/show/NCT00560365?term=FACS&rank=3,
Grossmann EM, Johnson FE, Virgo KS, Longo WE, Fossati R: Follow-up of colorectal cancer patients after resection with curative intent-the GILDA trial. Surg Oncol. 2004, 13 (2–3): 119-124.
Clemen RT, Lacke CJ: Analysis of colorectal cancer screening regimens. Health Care Manag Sci. 2001, 4 (4): 257-267. 10.1023/A:1011886010184.
Frazier AL, Colditz GA, Fuchs CS, Kuntz KM: Cost-effectiveness of screening for colorectal cancer in the general population. JAMA. 2000, 284 (15): 1954-1961. 10.1001/jama.284.15.1954.
Gyrd-Hansen D, Sogaard J, Kronborg O: Colorectal cancer screening: efficiency and effectiveness. Health Econ. 1998, 7 (1): 9-20. 10.1002/(SICI)1099-1050(199802)7:1<9::AID-HEC304>3.0.CO;2-H.
Haug U, Brenner H: A simulation model for colorectal cancer screening: potential of stool tests with various performance characteristics compared with screening colonoscopy. Cancer Epidemiol Biomarkers Prev. 2005, 14 (2): 422-428. 10.1158/1055-9965.EPI-04-0411.
Lejeune C, Arveux P, Dancourt V, Fagnani F, Bonithon-Kopp C, Faivre J: A simulation model for evaluating the medical and economic outcomes of screening strategies for colorectal cancer. Eur J Cancer Prev. 2003, 12 (1): 77-84. 10.1097/00008469-200302000-00012.
Lejeune C, Arveux P, Dancourt V, Bejean S, Bonithon-Kopp C, Faivre J: Cost-effectiveness analysis of fecal occult blood screening for colorectal cancer. Int J Technol Assess Health Care. 2004, 20 (4): 434-439.
Loeve F, Boer R, van Oortmarssen GJ, van Ballegooijen M, Habbema JD: The MISCAN-COLON simulation model for the evaluation of colorectal cancer screening. Comput Biomed Res. 1999, 32 (1): 13-33. 10.1006/cbmr.1998.1498.
Neilson AR, Whynes DK: Cost-effectiveness of screening for colorectal cancer: a simulation model. IMA J Math Appl Med Biol. 1995, 12 (3–4): 355-367.
Song K, Fendrick AM, Ladabaum U: Fecal DNA testing compared with conventional colorectal cancer screening methods: a decision analysis. Gastroenterology. 2004, 126 (5): 1270-1279. 10.1053/j.gastro.2004.02.016.
Vijan S, Hwang EW, Hofer TP, Hayward RA: Which colon cancer screening test? A comparison of costs, effectiveness, and compliance. Am J Med. 2001, 111 (8): 593-601. 10.1016/S0002-9343(01)00977-9.
Roberts S, Wang L, Klein R, Neiss R, Dittus R: Development of a simulation model of colorectal cancer. ACM Trans Modeling Comp Simulation. 2007, 18 (1): 4.1-4.30.
Kuntz KM, Lansdorp-Vogelaar I, Rutter CM, Knudsen AB, van Ballegooijen M, Savarino JE, Feuer EJ, Zauber AG: A systematic comparison of microsimulation models of colorectal cancer: the role of assumptions about adenoma progression. Med Decis Making. 2011, 31 (4): 530-539. 10.1177/0272989X11408730.
Renehan AG, O’Dwyer ST, Whynes DK: Cost effectiveness analysis of intensive versus conventional follow up after curative resection for colorectal cancer. BMJ. 2004, 328 (7431): 81-10.1136/bmj.328.7431.81.
Borie F, Combescure C, Daures JP, Tretarre B, Millat B: Cost-effectiveness of two follow-up strategies for curative resection of colorectal cancer: comparative study using a Markov model. World J Surg. 2004, 28 (6): 563-569.
Gazelle GS, Hunink MG, Kuntz KM, McMahon PM, Halpern EF, Beinfeld M, Lester JS, Tanabe KK, Weinstein MC: Cost-effectiveness of hepatic metastasectomy in patients with metastatic colorectal carcinoma: a state-transition Monte Carlo decision analysis. Ann Surg. 2003, 237 (4): 544-555.
Erenay FS, Alagoz O, Banerjee R, Cima RR: Estimating the unknown parameters of the natural history of metachronous colorectal cancer using discrete-event simulation. Med Decis Making. 2011, 31 (4): 611-624. 10.1177/0272989X10391809.
Sargent DJ, Patiyil S, Yothers G, Haller DG, Gray R, Benedetti J, Buyse M, Labianca R, Seitz JF, O’Callaghan CJ, Francini G, Grothey A, O'Connell M, Catalano PJ, Kerr D, Green E, Wieand HS, Goldberg RM, de Gramont A: End points for colon cancer adjuvant trials: observations and recommendations based on individual patient data from 20,898 patients enrolled onto 18 randomized trials from the ACCENT Group. J Clin Oncol. 2007, 25 (29): 4569-4574. 10.1200/JCO.2006.10.4323.
Graham RA, Wang S, Catalano PJ, Haller DG: Postsurgical surveillance of colon cancer: preliminary cost analysis of physician examination, carcinoembryonic antigen testing, chest x-ray, and colonoscopy. Ann Surg. 1998, 228 (1): 59-63. 10.1097/00000658-199807000-00009.
National Comprehensive Cancer Network: NCCN Clinical Practice Guidelines in Oncology: Colon Cancer. 2014,http://www.nccn.org/professionals/physician_gls/pdf/colon.pdf,
Tan E, Gouvas N, Nicholls RJ, Ziprin P, Xynos E, Tekkis PP: Diagnostic precision of carcinoembryonic antigen in the detection of recurrence of colorectal cancer. Surg Oncol. 2009, 18 (1): 15-24. 10.1016/j.suronc.2008.05.008.
McCormack PM, Ginsberg KB, Bains MS, Burt ME, Martini N, Rusch VW, Ginsberg RJ: Accuracy of lung imaging in metastases with implications for the role of thoracoscopy. Ann Thorac Surg. 1993, 56 (4): 863-865. 10.1016/0003-4975(93)90344-H. discussion 865–866
McIntosh J, Sylvester PA, Virjee J, Callaway M, Thomas MG: Pulmonary staging in colorectal cancer–is computerised tomography the answer?. Ann R Coll Surg Engl. 2005, 87 (5): 331-333. 10.1308/003588405X60579.
Valk PE, Abella-Columna E, Haseman MK, Pounds TR, Tesar RD, Myers RW, Greiss HB, Hofer GA: Whole-body PET imaging with [18 F]fluorodeoxyglucose in management of recurrent colorectal cancer. Arch Surg. 1999, 134 (5): 503-511. 10.1001/archsurg.134.5.503. discussion 511–503
Valls C, Lopez E, Guma A, Gil M, Sanchez A, Andia E, Serra J, Moreno V, Figueras J: Helical CT versus CT arterial portography in the detection of hepatic metastasis of colorectal carcinoma. AJR Am J Roentgenol. 1998, 170 (5): 1341-1347. 10.2214/ajr.170.5.9574613.
Valls C, Andia E, Sanchez A, Guma A, Figueras J, Torras J, Serrano T: Hepatic metastases from colorectal cancer: preoperative detection and assessment of resectability with helical CT. Radiology. 2001, 218 (1): 55-60. 10.1148/radiology.218.1.r01dc1155.
Glover C, Douse P, Kane P, Karani J, Meire H, Mohammadtaghi S, Allen-Mersh TG: Accuracy of investigations for asymptomatic colorectal liver metastases. Dis Colon Rectum. 2002, 45 (4): 476-484. 10.1007/s10350-004-6224-y.
Konopke R, Bunk A, Kersting S: The role of contrast-enhanced ultrasound for focal liver lesion detection: an overview. Ultrasound Med Biol. 2007, 33 (10): 1515-1526. 10.1016/j.ultrasmedbio.2007.04.009.
Staib L, Schirrmeister H, Reske SN, Beger HG: Is (18)F-fluorodeoxyglucose positron emission tomography in recurrent colorectal cancer a contribution to surgical decision making?. Am J Surg. 2000, 180 (1): 1-5. 10.1016/S0002-9610(00)00406-2.
Winawer SJ, Stewart ET, Zauber AG, Bond JH, Ansel H, Waye JD, Hall D, Hamlin JA, Schapiro M, O’Brien MJ, Sternberg SS, Gottlieb LS: A comparison of colonoscopy and double-contrast barium enema for surveillance after polypectomy. National Polyp Study Work Group. N Engl J Med. 2000, 342 (24): 1766-1772. 10.1056/NEJM200006153422401.
World Health Organization: Life tables for WHO member states. 1990, [http://www.who.int/gho/mortality_burden_disease/life_tables/life_tables/en/]
Pietra N, Sarli L, Costi R, Ouchemi C, Grattarola M, Peracchia A: Role of follow-up in management of local recurrences of colorectal cancer: a prospective, randomized study. Dis Colon Rectum. 1998, 41 (9): 1127-1133. 10.1007/BF02239434.
O’Connell MJ, Campbell ME, Goldberg RM, Grothey A, Seitz JF, Benedetti JK, Andre T, Haller DG, Sargent DJ: Survival following recurrence in stage II and III colon cancer: findings from the ACCENT data set. J Clin Oncol. 2008, 26 (14): 2336-2341. 10.1200/JCO.2007.15.8261.
Kong CY, McMahon PM, Gazelle GS: Calibration of disease simulation model using an engineering approach. Value Health. 2009, 12 (4): 521-529. 10.1111/j.1524-4733.2008.00484.x.
Simmonds PC, Primrose JN, Colquitt JL, Garden OJ, Poston GJ, Rees M: Surgical resection of hepatic metastases from colorectal cancer: a systematic review of published studies. Br J Cancer. 2006, 94 (7): 982-999. 10.1038/sj.bjc.6603033.
Pfannschmidt J, Dienemann H, Hoffmann H: Surgical resection of pulmonary metastases from colorectal cancer: a systematic review of published series. Ann Thorac Surg. 2007, 84 (1): 324-338. 10.1016/j.athoracsur.2007.02.093.
Sorbye H, Pfeiffer P, Cavalli-Bjorkman N, Qvortrup C, Holsen MH, Wentzel-Larsen T, Glimelius B: Clinical trial enrollment, patient characteristics, and survival differences in prospectively registered metastatic colorectal cancer patients. Cancer. 2009, 115 (20): 4679-4687. 10.1002/cncr.24527.
Hughes KS, Rosenstein RB, Songhorabodi S, Adson MA, Ilstrup DM, Fortner JG, Maclean BJ, Foster JH, Daly JM, Fitzherbert D, Sugarbaker P, Iwatsuki S, Starzl T, Ramming KP, Longmire WP, O’Toole K, Petrelli NJ, Herrera L, Cady B, McDermott W, Nims T, Enker WE, Coppa GF, Blumgart LH, Bradpiece H, Urist M, Aldrete JS, Schlag P, Hohenberger P, Steele G, et al: Resection of the liver for colorectal carcinoma metastases. A multi-institutional study of long-term survivors. Dis Colon Rectum. 1988, 31 (1): 1-4. 10.1007/BF02552560.
Olak J, Wexler MJ, Rodriguez J, McLean AP: Hepatic resection for metastatic disease. Can J Surg. 1986, 29 (6): 435-439.
Petrelli NJ, Nambisan RN, Herrera L, Mittelman A: Hepatic resection for isolated metastasis from colorectal carcinoma. Am J Surg. 1985, 149 (2): 205-209. 10.1016/S0002-9610(85)80065-9.
Telford JJ, Levy AR, Sambrook JC, Zou D, Enns RA: The cost-effectiveness of screening for colorectal cancer. CMAJ. 2010, 182 (12): 1307-1313. 10.1503/cmaj.090845.
Pignone M, Saha S, Hoerger T, Mandelblatt J: Cost-effectiveness analyses of colorectal cancer screening: a systematic review for the U.S. Preventive Services Task Force. Ann Intern Med. 2002, 137 (2): 96-104. 10.7326/0003-4819-137-2-200207160-00007.
Knudsen AB, Lansdorp-Vogelaar I, Rutter CM, Savarino JE, van Ballegooijen M, Kuntz KM, Zauber AG: Cost-effectiveness of computed tomographic colonography screening for colorectal cancer in the medicare population. J Natl Cancer Inst. 2010, 102 (16): 1238-1252. 10.1093/jnci/djq242.
Lansdorp-Vogelaar I, Kuntz KM, Knudsen AB, Wilschut JA, Zauber AG, van Ballegooijen M: Stool DNA testing to screen for colorectal cancer in the Medicare population: a cost-effectiveness analysis. Ann Intern Med. 2010, 153 (6): 368-377. 10.7326/0003-4819-153-6-201009210-00004.
Kjeldsen BJ, Kronborg O, Fenger C, Jorgensen OD: The pattern of recurrent colorectal cancer in a prospective randomised study and the characteristics of diagnostic tests. Int J Colorectal Dis. 1997, 12 (6): 329-334. 10.1007/s003840050118.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6947/14/29/prepub
JR, NJM, MWK, and GSC were supported for this work in part by a National Cancer Institute Program Grant (5P30 CA043703-21). JR was supported by a Mentored Research Scholar Grant in Applied and Clinical Research, MRSG-13-315-01-CPHPS from the American Cancer Society. The funders played no role in the study design, in the collection, analysis, and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.
The authors confirm that neither the submitted manuscript nor any similar manuscript, in whole or in part, is under consideration, in press, or published elsewhere.
The authors declare that they have no competing interests.
JR and KMA conceived the original idea for the work. JR developed the model with input from KMA, CYK, NJM, GSC, and MWK. JR, KMA, XA, and QH developed parameters estimates based on review of the literature. JR and CYK planned and performed model calibration and validation. JR wrote the manuscript with editorial input from all other authors. All authors read and approved the final manuscript.