A LASSO-based survival prediction model for patients with synchronous colorectal carcinomas based on SEER
Introduction
Colorectal cancer (CRC) is the third common cancer and ranks third as a cause of cancer-related death for males in America (1,2). The definition of multiple primary colorectal carcinomas (MPCC) is the presence of 2 or more primary invasive adenocarcinomas diagnosed in patients. Synchronous colorectal carcinomas (SCC) is identified as the second invasive adenocarcinomas diagnosis within 6 months after the first invasive adenocarcinomas diagnosis (3). Metachronous colorectal carcinomas (MCC) is identified as the second invasive adenocarcinomas diagnosis after more than 6 months after the first invasive adenocarcinomas diagnosis (4). Among patients suffering from CRC, SCC contribute 1% to 8% (5). Patients with SCC must be fully studied since it’s not a rare incident of SCC among CRC (6).
Currently, the literature of SCC is mostly small series (<80 patients) which described epidemiology and clinicopathology (3-5,7-9). In these studies, a part of them indicated there was no appreciable differences between patients with synchronous tumors and single neoplasm in survival when compared to individuals who had single neoplasms. In contrast, the other part of them indicated individuals who had metachronous carcinomas have been observed to show a poor clinical outcome after the development of the second carcinoma (7). The prognosis of SCC still controversial. What’s more, there are few reports on the impact factors of synchronous colorectal carcinoma’s overall survival (OS) and formulation of prognostic models (8). So, we evaluated the impact factors of SCC on the OS and made a prognostic model with a large cohort of patients.
In this study, our aim was to develop and validate a nomogram based on treatment variables, surgical variables, clinical characteristics and tumor characteristics to predict the survival of SCC patients. The data was obtained from the population-based Surveillance, Epidemiology, and End Results (SEER) database which contains a large sample size and has a long follow-up time. Since the prediction model is associated with the first and second-time treatment variables, surgical variables and tumor characteristics, it may exist a multicollinearity problem between first and second-time variables. Furthermore, because of incorporating too many variables, there may be over-fitting in predicting model. For these reasons, we selected the least absolute shrinkage and selection operator (LASSO) method to deal with the above concerns. In order to determine whether the prediction model fitted with LASSO-based Cox regression was better, we compared the LASSO model (fit by Cox regression after variables selection by LASSO Cox regression) to the COX model (fit by Cox regression), TNM model (established in first time T, N, M grade) and TTNNMM model (established in first and second times T, N, M grade). We present the following article in accordance with the STROBE reporting checklist (available at https://tcr.amegroups.com/article/view/10.21037/tcr-20-1860/rc).
Methods
Data source
We identified the survivors from the National Cancer Institute’s Surveillance, Epidemiology, and End Results (SEER) 13-registry database by analyzing patients diagnosed from 2004 to 2013. The SEER has developed and maintained high-quality, validated data on causes of death among cancer survivors, providing insight into relative and cause-specific deaths in this population (10,11). Data was retrieved using SEER*Stat 8.3.5 (Surveillance Research Program, National Cancer Institute, Bethesda, MD). Our study was determined and it exempted the data from Colorectal Surgery Union Hospital in Fujian. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).
Patients
Patients over 18 years old who were diagnosed with colorectal carcinoma between 2004 and 2013 with surgery were initially analyzed. Patients diagnosed with single colorectal carcinoma were excluded to explore SCC. Patients who survival time was unknown were excluded to explore the epidemiology, pathogenesis and factors that influenced the survival of SCC. Patients with an unknown grade of the tumor, unknown T stage, unknown N stage, and unknown M stage were excluded for further comparison of the feasibility of the TNM model and TTNNMM model. Patients with unknown prognostic characteristics (including race, tumor size and location) were also excluded. The clinicopathologic variables were then collected from the SEER 13 database, including gender, race, sex, delta t, months survived and first and second times age of diagnosis, marital, location of tumor, TNM staging (12), histologic grade, number of lymph nodes examined, number of positive lymph nodes, tumor size, radiation sequence, chemotherapy, and surgical related variables. Then, patients who had CRC more than 3 times or multiple metachronous primary carcinomas were excluded. Lastly, we excluded patients who survived less than 1 month or other variable unknown (Figure 1).
In the construction of the survival predicting model, the internal cohort included patients from SEER database, while the external validation cohort consisted of patients from Colorectal Surgery Union Hospital in Fujian database.
Statistical analysis
Statistical analysis was carried out with R software (version 3.4.2; http://www.Rproject.org) and SPSS (Statistical Product and Service Solutions, version 22.0). The packages in R used in this study are as follows. Statistical significance levels were all two-sided, with statistical significance set at 0.05.
Variables selection and models constructing
The LASSO method, which is suitable for the regression of high-dimensional data (13,14), was used to select the most useful predictive variables from the primary data set. The “glmnet” package was used to perform the LASSO Cox regression model analysis (15). To compare the differences in the Cox method and LASSO-based Cox method, we separately used the Cox regression or LASSO-based Cox regression to construct models. The COX model was selected and constructed using the internal cohort by backward Cox analysis using Akaike information criterion (AIC) selection criteria and the best model was selected with the least AIC (16-18). The LASSO model was selected and constructed using the internal cohort by LASSO-based Cox regression. The TNM model was established using the internal cohort by the first time T, N and M staging and that of TTNNMM model was established using the internal cohort by the first and second times T, N and M staging (19).
Compare models
The LASSO model using LASSO-based Cox regression and COX model using backward Cox analysis was first internally validated in the internal cohort using a bootstrap method (1,000 bootstraps resamples) and then externally validated in the external cohorts. The 3- and 5-year OS calibration of the LASSO model and COX model were performed by comparing the observed survival with the predicted survival in the internal and external cohorts. Then for survival testing with the LASSO model and COX model of the specificity, time-dependent receiver operating characteristic (timeROC) curves were estimated for two cohorts by inverse probability of censoring weighting estimators (KM-weight) at 3-, and 5-year (20,21). Sequential AUCs were compared among the four models using identically and independently distributed representations of the AUC estimators (22). Also, the overall prognostic performance of the four models was assessed using the Bayesian information criterion (BIC) via bootstrap-resampling analysis. Lastly, four models were evaluated with AIC, C-index (23), the net reclassification improvement (NRI) (24,25) and integrated discrimination improvement (IDI) (26).
Clinical use
The net benefit and clinical usefulness of the four models above were estimated with decision curve analysis (DCA) throughout the whole cohort (27).
Nomogram for a visualization model
For the purpose of illustration and clinical applicability, we created a nomogram based on the LASSO model. In the nomogram, model-based score points for each predictor variable category were displayed, which has to be summarized for any individual patient. From the resulting total number of points, the corresponding predicted survival probabilities from the nomogram could be easily read.
Results
Clinical characteristics
From the data obtained from 2004 to 2013, 4,616 patients with SCC in the SEER database were found. Patient characteristics are shown in Table 1. There are significant correlations in age and slight correlation in pN, pM, examined lymph nodes, Surg Prim Site, and chemotherapy between the twice synchronous colorectal. Patients with SCC were mostly older (>65 years), more often men (54.1%), and likely to have a depth of invasion by T3 (56.8% and 41.5% by the first and second times). Tumors were mostly situated in the cecum, ascending colon and sigmoid colon.
Table 1
Characteristics | First time, total N=4,616 (n, %) | Second time, total N=4,616 (n, %) | rho |
---|---|---|---|
Age (years) | 0.999 | ||
0–49 | 369 (8.0) | 368 (8.0) | |
50–59 | 619 (13.4) | 618 (13.4) | |
60–64 | 454 (9.8) | 450 (9.7) | |
65–69 | 573 (12.4) | 577 (12.5) | |
70–74 | 652 (14.1) | 646 (14.0) | |
75–79 | 721 (15.6) | 727 (15.7) | |
80–84 | 656 (14.2) | 654 (14.2) | |
85+ | 573 (12.4) | 576 (12.5) | |
Site | 0.087 | ||
Large intestine, NOS | 34 (0.7) | 38 (0.8) | |
Rectum | 468 (10.1) | 599 (13.0) | |
Rectosigmoid junction | 228 (4.9) | 286 (6.2) | |
Sigmoid colon | 816 (17.7) | 786 (17.0) | |
Descending colon | 287 (6.2) | 372 (8.1) | |
Splenic flexure | 163 (3.5) | 165 (3.6) | |
Transverse colon | 531 (11.5) | 583 (12.6) | |
Hepatic flexure | 268 (5.8) | 221 (4.8) | |
Ascending colon | 807 (17.5) | 766 (16.6) | |
Cecum | 1,014 (22.0) | 800 (17.3) | |
pT | 0.205 | ||
Tis/T0 | 95 (2.1) | 253 (5.5) | |
T1 | 624 (13.5) | 1,229 (26.6) | |
T2 | 702 (15.2) | 909 (19.7) | |
T3 | 2,620 (56.8) | 1,917 (41.5) | |
T4 | 575 (12.5) | 308 (6.7) | |
pN | 0.638 | ||
N0 | 2,746 (59.5) | 3,032 (65.7) | |
N1 | 1,163 (25.2) | 983 (21.3) | |
N2 | 707 (15.3) | 601 (13.0) | |
pM | 0.460 | ||
M0 | 4,083 (88.5) | 4,067 (88.1) | |
M1 | 533 (11.5) | 549 (11.9) | |
Examined lymph nodes | 0.770 | ||
1–14 | 1,652 (35.8) | 1,813 (39.3) | |
15–39 | 2,558 (55.4) | 2,417 (52.4) | |
40+ | 406 (8.8) | 386 (8.4) | |
Grade | 0.306 | ||
Unknown | 233 (5.0) | 499 (10.8) | |
Grade I | 360 (7.8) | 483 (10.5) | |
Grade II | 3,006 (65.1) | 2,924 (63.3) | |
Grade III | 880 (19.1) | 629 (13.6) | |
Grade IV | 137 (3.0) | 81 (1.8) | |
Tumor size | 0.172 | ||
1–29 mm | 860 (18.6) | 1,445 (31.3) | |
30–99 mm | 3,035 (65.7) | 2,238 (48.5) | |
100–900 mm | 192 (4.2) | 96 (2.1) | |
Unknown | 529 (11.5) | 837 (18.1) | |
Surg prim site | 0.559 | ||
Radical resection | 4,159 (90.1) | 4,158 (90.1) | |
Combined organ resection | 403 (8.7) | 361 (7.8) | |
Partial resection | 29 (0.6) | 63 (1.4) | |
Surgical resection | 13 (0.3) | 22 (0.5) | |
Radical resection + ostomy | 12 (0.3) | 12 (0.3) | |
Chemotherapy | 0.796 | ||
No/unknown | 3,063 (66.4) | 3,237 (70.1) | |
Yes | 1,553 (33.6) | 1,379 (29.9) | |
Sex | |||
Male | 2,498 (54.1) | ||
Female | 2,118 (45.9) | ||
Race | |||
White | 3,773 (81.7) | ||
Black | 489 (10.6) | ||
Other | 354 (7.7) | ||
Marital | |||
Single | 581 (12.6) | ||
Unmarried | 4 (0.1) | ||
Marry | 2,425 (52.5) | ||
Separated | 43 (9.0) | ||
Divorce | 386 (8.4) | ||
Widowed | 1,000 (21.7) | ||
Unknown | 177 (3.8) |
NOS, nitric oxide synthase
Predictive variable selection
It was showed that 31 variables were reduced to 11 potential predictors on the basis of 4,616 patients by LASSO-based Cox regression in the internal cohort (Figure 2A-2C) or were reduced to 16 potential predictors by Cox regression base minute AIC (Figure 2D).
Results of the selected variables with Cox regression and LASSO-based Cox regression are listed in Table 2. Table 2 indicates that the age of the first time SCC diagnosis, sex, first time size, first time surgery, second time marital, second time grade, second time chemotherapy, first and second times pT, pN, pM, regional nodes examined and site of disease was significantly associated with OS by Cox regression. Table 2 also indicates that the age of first time SCC diagnosis, sex, second-time chemotherapy, first and second times pT, pN, pM, regional nodes examined were significantly associated with OS by LASSO-based Cox regression.
Table 2
Variable | COX | LASSO | |||||||
---|---|---|---|---|---|---|---|---|---|
adj.p | HR | 95% CI, lower | 95% CI, upper | adj.p | HR | 95% CI, lower | 95% CI, upper | ||
First time characteristics | |||||||||
Age (years) | <0.001 | 2.293 | 2.091 | 2.514 | <0.001 | 2.327 | 2.129 | 2.544 | |
Sex | <0.001 | 0.794 | 0.727 | 0.868 | <0.001 | 0.816 | 0.749 | 0.889 | |
Site | 0.023 | 1.092 | 1.012 | 1.178 | |||||
pT | <0.001 | 1.253 | 1.188 | 1.323 | <0.001 | 1.247 | 1.182 | 1.316 | |
pN | <0.001 | 1.733 | 1.490 | 2.015 | <0.001 | 1.688 | 1.450 | 1.964 | |
pM | <0.001 | 2.444 | 2.114 | 2.824 | <0.001 | 2.422 | 2.096 | 2.798 | |
Examined lymph nodes | <0.001 | 0.691 | 0.597 | 0.799 | 0.035 | 0.778 | 0.617 | 0.982 | |
Tumor size | 0.134 | 1.132 | 0.962 | 1.332 | |||||
Surg prim site | 0.013 | 1.681 | 1.118 | 2.527 | |||||
Second time characteristics | |||||||||
Marital | 0.046 | 1.055 | 1.001 | 1.111 | |||||
Site | 0.043 | 0.926 | 0.860 | 0.998 | |||||
Grade | 0.091 | 1.204 | 0.971 | 1.493 | |||||
pT | <0.001 | 1.312 | 1.197 | 1.438 | <0.001 | 1.349 | 1.236 | 1.473 | |
pN | 0.07 | 1.159 | 0.988 | 1.359 | 0.034 | 1.190 | 1.013 | 1.397 | |
pM | 0.003 | 1.253 | 1.079 | 1.457 | 0.003 | 1.255 | 1.081 | 1.458 | |
Examined lymph nodes | 0.312 | 0.887 | 0.704 | 1.119 | |||||
Chemotherapy | <0.001 | 0.715 | 0.642 | 0.795 | <0.001 | 0.721 | 0.649 | 0.802 |
LASSO, least absolute shrinkage and selection operator; SCC, synchronous colorectal carcinomas; OS, overall survival.
Results from the relation between first and second times pT, pN, pM, grade and regional nodes are listed in Table 1.
Development of COX model and LASSO model
The multivariable regression model for age, sex, marital, race, site, pT, pN, pM, radiation, chemotherapy, surgery, nodes examined, etc. were included in the regression after variables were selected by the LASSO-based Cox regression or Cox regression. We showed hazard ratios with 95% CIs for covariates which are included in Table 2.
Apparent performance of the LASSO model or COX model in the internal cohort
The calibration curves of the LASSO model and COX model for the probability of OS in 3- and 5-year between prediction and observation in the internal cohort (Figure 3A-3D) were plotted to assess the calibration of the COX model and LASSO model, which were accompanied with the Hosmer-Lemeshow test (A significant test statistic implies that the model calibrates perfectly).
Validation of the LASSO model and COX model
The external validation was tested in the external cohort. The LASSO model was formed in the internal cohort and was applied to all the patients of the external cohort. The calibration curves in 3- and 5-year (Figure 4A,4B) were derived on the basis of the regression analysis.
C-index and AIC
To quantify the discrimination performance of the LASSO model, COX model, TNM model, and TTNNMM model, Harrell’s C-index and AIC were applied (Table 3). The C-index for the LASSO model, COX model, TNM model and TTNNMM model were 0.710 (95% CI: 0.703 to 0.717), 0.712 (95% CI: 0.705 to 0.719), 0.637 (95% CI: 0.631 to 0.644) and 0.651 (95% CI: 0.644 to 0.657), which were confirmed to be 0.710, 0.712, 0.637 and 0.651 via bootstrapping validation. The AIC for the LASSO model, COX model, TNM model, and TTNNMM model were 33,431, 33,420, 34,043, 33,994. According to C-index and AIC, there are no significant difference between COX model and LASSO model.
Table 3
Model | AIC | C-index | Concordance |
---|---|---|---|
LASSO | 33,431 | 0.710 | 0.710 (0.703–0.717) |
COX | 33,420 | 0.712 | 0.712 (0.705–0.719) |
TNM | 34,043 | 0.637 | 0.637 (0.631–0.644) |
TTNNMM | 33,994 | 0.651 | 0.651 (0.644–0.657) |
LASSO model: the model fit by Cox regression after variables selection by LASSO Cox regression; COX model: the model fit by Cox regression; TNM model: the model established in first time T, N, M grade; TTNNMM model: the model established in first and second times T, N, M grade. AIC, akaike information criterion; LASSO, least absolute shrinkage and selection operator.
Predictive accuracy of COX model and LASSO model
According to the survROC curves for 1-, 3-, 5-year OS for LASSO model, the COX model, TNM model, and TTNNMM model (Figure 5A-5D), the ROC curve (a general measure of predictiveness) was found to be greater in 3- and 5-year.
Whether apparent different performance of the LASSO and COX model
TimeAUC
Time-dependent ROC curves were generated to compare the sequential trends of the LASSO, COX, TNM and TTNNMM model for OS. The time-dependent ROC curve of the LASSO model was continuously superior to that of the COX model, TNM model and TTNNMM model (Figure 6).
BIC
The prognostic performances of the LASSO, COX, TNM, and TTNNMM model were compared using BIC, which is not only a measure of the goodness of fit of an estimated statistical model but also accurately considers the number of parameters included in the model. As shown in Figure 7, there was no significant difference between the COX and LASSO model after the bootstrap analysis (BIC 4.49, 95% CI: –2.92 to 11.91) but there was a significant difference between the TNM and LASSO model (BIC 1,178.76, 95% CI: 1,171.15–1,186.37), also TTNNMM and LASSO model (BIC 1,098.57, 95% CI: 1,092.05–1,105.09).
NRI and IDI
The discriminant ability for LASSO model, COX model, TNM model, and TTNNMM model was calculated using NRI and IDI (Table 4). Compared to the TNM model and TTNNMM model, LASSO model was found to be a higher discriminant and possess reclassification indices (IDI 0.072 and 0.064; P<0.001; NRI 0.525 and 0.466) (Table 4). In addition, compared to the COX model, the LASSO model doesn’t significantly decrease the discriminant and reclassification indices (IDI −0.002, P=0.058; NRI −0.009) (Table 4).
Table 4
Model | NRI | IDI | ||||||
---|---|---|---|---|---|---|---|---|
Value | 95% CI, lower | 95% CI, upper | Value | 95% CI, lower | 95% CI, upper | P value | ||
LASSO & COX | −0.009 | −0.054 | 0.012 | −0.002 | −0.006 | 0.000 | 0.058 | |
LASSO & TNM | 0.525 | 0.464 | 0.589 | 0.072 | 0.061 | 0.087 | <0.001 | |
LASSO & TTNNMM | 0.466 | 0.415 | 0.538 | 0.064 | 0.052 | 0.079 | <0.001 |
LASSO model: the model fit by Cox regression after variables selection by LASSO Cox regression; COX model: the model fit by Cox regression; TNM model: the model established in first time T, N, M grade; TTNNMM model: the model established in first and second times T, N, M grade. NRI, net reclassification improvement; IDI, integrated discrimination improvement; LASSO, least absolute shrinkage and selection operator.
Clinical use
DCA was conducted to determine the clinical usefulness of the LASSO model by quantifying the net benefits at different threshold probabilities. We also plotted the decision curve for the four models in 3- and 5-year (Figure 8A,8B).
Visualization of SCC survival prediction model
Survival prediction model of the nomogram was established based on factors selected by LASSO-based Cox regression (Figure 9). The nomogram showed that first time age had the most contribution to prognosis, followed by first- and second-times T stage, N stage, metastases and examined lymph nodes. Sex had a modest effect on survival. Each subtype of the variables was assigned a score. A straight line can be drawn down at each time point on the total point scale to determine the estimated probability of survival, according to the total number of points. For each predictor, the points assigned on the 0–10 scale at the top are read and then these points are added. The number on the “Total Points” scale were found and then the corresponding predictions of 3-, and 5-year risk are recorded.
Discussion
In this study, we developed and validated a prognostic model about SCC which based on large data from SEER by LASSO-based Cox regression. Our results show that the OS is associated with age, sex, second-time chemotherapy and first and second times pT, pN, pM.
Notably, the second time examined lymph nodes didn’t show enough predictive strength on the basis of Cox regression, which makes a common strategy to exclude this variable for model development. However, it may be a result of nuances in the data set or confounding by other predictors that reject important predictors (17,28), for which no significant statistical association with OS does not definitively imply that examined lymph nodes are unimportant. In addition, more lymph nodes examined may mean a better quality of operation. Therefore, we kept the second time examined lymph nodes as a candidate factor in the process of model development. For the same reason, we kept the first-time tumor size and the second-time grade and pN in the COX model.
Grade, size, surgery and marital which may be multi-collinearity bias with pT, pN, pM, and age were not included in the LASSO model. Grade, size, surgery may be associated with TNM grading (12). Besides, the old aged were more likely to be widowed. Also, because of overfitting site was not included in the LASSO model.
From Table 1, we found that patients with SCC were generally older (>65 years), more often male, and likely with the depth of invasion of T3. There may be less estrogen to protect in male (29) and a high probability of microsatellite instability (MSI) in older patients (30). Besides, it may be as a result of tumor biologic characteristics with the depth of invasion by T3. Therefore, patients who are older (>65 years), men, and depth of invasion by T3 should closely monitor the postoperative enteroscopy for early detection of SCC.
Tables 3,4 and Figures 3-8 show that there was no significant difference between the LASSO model and COX model (NRI, IDI, c index, ROC, AIC, and BIC) but the LASSO model was obviously better than the TNM model and TTNNMM model (NRI, IDI, c index, AIC, ROC, and BIC). Compared to the COX model, the LASSO model significantly reduced the variables included which minimized overfitting and collinearity. Moreover, Figure 6 shows that the LASSO model performs better than the Cox model in the timeAUC. Although the LASSO model included fewer variables, the LASSO model performed better in the timeAUC compared to the COX model.
The most important and final argument for the use of the nomogram is based on the need to interpret the individual need for additional treatment or care. However, the clinical consequences of a particular level of discrimination or degree of miscalibration cannot be adequately assessed by the risk-prediction discrimination, performance, and calibration (17,31). Therefore, in order to justify the clinical usefulness, it is crucial to ascertain whether the LASSO model-assisted decisions can improve patient outcomes. With this aim, in this study, the application of the DCA instead of the multi-institutional prospective for the validation of the model was performed. This novel method offers an insight into the clinical consequences on the basis of threshold probability, from which the net benefit could be derived (17,32). Through the decision curve plot, we can ascertain whether the probability of threshold of a patient or doctor is 5% using the LASSO model is more beneficial than either the TNM model or TTNNMM model and not inferior compared to the COX model (Figure 8A,8B).
There are some limitations in the present work that should be discussed. The collection of the SEER database is retrospective. There is a lack of molecular data and data for biological prognostic factors that might also influence the prognosis of SCC patients. In recent years, increased research with gene markers, such as MSI, SSA, BRAF V600E associated with SCC has been proposed (3,8). Moreover, there might be some increase in the bias that we excluded all patients who had missing data from the collected variables. The study didn’t incorporate detailed chemotherapy and radiation methods due to the lack of adequate information and large bias of the information. Finally, although this nomogram performed well in both internal and external cohorts, due to the influence of deaths related to the operation, the data should be used with caution when predicting 1-year risk. Even so, although it didn’t include the genomic characteristics, excluded patients who had missing data and was retrospect to analysis, it was the first model to perform a prognosis OS of SCC.
Based on the database content, the main influencing factors were screened for the LASSO model. Due to the limitations of the database, some important factors weren’t covered. In the future, we hope to have relevant data to incorporate it into our research.
In conclusion, this study presents a prognosis nomogram that incorporates both the first-time and the second-time variables and can be conveniently used to facilitate the prediction of OS in patients with SCC.
Acknowledgments
We thank all the staff in the Department of Colorectal Surgery of FMUUH and Dr. Jerry for polishing the manuscript.
Funding: This study was supported by National Clinical Key Specialty Construction Project (General Surgery) of China and Wu Jieping Medical Foundation (No. 320.2710.1833).
Footnote
Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at https://tcr.amegroups.com/article/view/10.21037/tcr-20-1860/rc
Peer Review File: Available at https://tcr.amegroups.com/article/view/10.21037/tcr-20-1860/prf
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tcr.amegroups.com/article/view/10.21037/tcr-20-1860/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin 2018;68:7-30. [Crossref] [PubMed]
- Huang YQ, Liang CH, He L, et al. Development and Validation of a Radiomics Nomogram for Preoperative Prediction of Lymph Node Metastasis in Colorectal Cancer. J Clin Oncol 2016;34:2157-64. [Crossref] [PubMed]
- Hu H, Chang DT, Nikiforova MN, et al. Clinicopathologic features of synchronous colorectal carcinoma: A distinct subset arising from multiple sessile serrated adenomas and associated with high levels of microsatellite instability and favorable prognosis. Am J Surg Pathol 2013;37:1660-70. [Crossref] [PubMed]
- Brandariz L, Alegre C, Rueda D, et al. New Perspectives in Multiple Primary Colorectal Cancer: A Surgical Approach. Digestion 2016;94:57-65. [Crossref] [PubMed]
- van Leersum NJ, Aalbers AG, Snijders HS, et al. Synchronous colorectal carcinoma: a risk factor in colorectal cancer surgery. Dis Colon Rectum 2014;57:460-6. [Crossref] [PubMed]
- Papadopoulos V, Michalopoulos A, Basdanis G, et al. Synchronous and metachronous colorectal carcinoma. Tech Coloproctol 2004;8:s97-s100. [Crossref] [PubMed]
- Fante R, Roncucci L, Di Gregorio C, et al. Frequency and clinical features of multiple tumors of the large bowel in the general population and in patients with hereditary colorectal carcinoma. Cancer 1996;77:2013-21. [Crossref] [PubMed]
- Lam AK, Chan SS, Leung M. Synchronous colorectal cancer: clinical, pathological and molecular implications. World J Gastroenterol 2014;20:6815-20. [Crossref] [PubMed]
- Liu YL, Xu HT, Jiang SX, et al. Prognostic significance of lymph node status in patients with metastatic colorectal carcinoma treated with lymphadenectomy. J Surg Oncol 2014;109:234-8. [Crossref] [PubMed]
- Howlader N, Ries LA, Mariotto AB, et al. Improved estimates of cancer-specific survival rates from population-based data. J Natl Cancer Inst 2010;102:1584-98. [Crossref] [PubMed]
- Mariotto AB, Noone AM, Howlader N, et al. Cancer survival: an overview of measures, uses, and interpretation. J Natl Cancer Inst Monogr 2014;2014:145-86. [Crossref] [PubMed]
- Amin MB, Greene FL, Edge SB, et al. The Eighth Edition AJCC Cancer Staging Manual: Continuing to build a bridge from a population-based to a more "personalized" approach to cancer staging. CA Cancer J Clin 2017;67:93-9.
- Sauerbrei W, Royston P, Binder H. Selection of important variables and determination of functional form for continuous predictors in multivariable model building. Stat Med 2007;26:5512-28. [Crossref] [PubMed]
- Tibshirani R. Regression shrinkage and selection via the lasso: A retrospective. J R Statist Soc B 2011;73:273-82. [Crossref]
- Tibshirani R. The lasso method for variable selection in the Cox model. Stat Med 1997;16:385-95. [Crossref] [PubMed]
- Venables WN, Ripley BD. Modern Applied Statistics with S. Statistics & Computing 2002;52:704-5.
- Collins GS, Reitsma JB, Altman DG, et al. Transparent Reporting of a Multivariable Prediction Model for Individual Prognosis or Diagnosis (TRIPOD): the TRIPOD Statement, 2015.
- Sauerbrei W, Boulesteix AL, Binder H. Stability investigations of multivariable regression models derived from low- and high-dimensional data. J Biopharm Stat 2011;21:1206-31. [Crossref] [PubMed]
- Rao SJ. Regression Modeling Strategies: With Applications to Linear Models, Logistic Regression, and Survival Analysis by Frank E. Harrell. Publications of the American Statistical Association 2005;98:257-8. [Crossref]
- Blanche P, Dartigues JF, Jacqmin-Gadda H. Estimating and comparing time-dependent areas under receiver operating characteristic curves for censored event times with competing risks. Stat Med 2013;32:5381-97. [Crossref] [PubMed]
- Heagerty PJ, Lumley T, Pepe MS. Time-dependent ROC curves for censored survival data and a diagnostic marker. Biometrics 2000;56:337-44. [Crossref] [PubMed]
- Rodríguez-Álvarez MX, Meira-Machado L, Abu-Assi E, et al. Nonparametric estimation of time-dependent ROC curves conditional on a continuous covariate. Stat Med 2016;35:1090-102. [Crossref] [PubMed]
- Hanley JA, McNeil BJ. A method of comparing the areas under receiver operating characteristic curves derived from the same cases. Radiology 1983;148:839-43. [Crossref] [PubMed]
- Pencina MJ, D'Agostino RB Sr, Steyerberg EW. Extensions of net reclassification improvement calculations to measure usefulness of new biomarkers. Stat Med 2011;30:11-21. [Crossref] [PubMed]
- Tangri N, Stevens LA, Griffith J, et al. A predictive model for progression of chronic kidney disease to kidney failure. JAMA 2011;305:1553-9. [Crossref] [PubMed]
- Chambless LE, Cummiskey CP, Cui G. Several methods to assess improvement in risk prediction models: extension to survival analysis. Stat Med 2011;30:22-38. [Crossref] [PubMed]
- Vickers AJ, Cronin AM, Elkin EB, et al. Extensions to decision curve analysis, a novel method for evaluating diagnostic tests, prediction models and molecular markers. BMC Med Inform Decis Mak 2008;8:53. [Crossref] [PubMed]
- Collins GS, Reitsma JB, Altman DG, et al. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD Statement. Eur J Clin Invest 2015;45:204-14. [Crossref] [PubMed]
- Fernandez E, La Vecchia C, Braga C, et al. Hormone replacement therapy and risk of colon and rectal cancer. Cancer Epidemiol Biomarkers Prev 1998;7:329-33. [PubMed]
- Seo HM, Chang YS, Joo SH, et al. Clinicopathologic characteristics and outcomes of gastric cancers with the MSI-H phenotype. J Surg Oncol 2009;99:143-7. [Crossref] [PubMed]
- Van Calster B, Vickers AJ. Calibration of risk prediction models: impact on decision-analytic performance. Med Decis Making 2015;35:162-9. [Crossref] [PubMed]
- Balachandran VP, Gonen M, Smith JJ, et al. Nomograms in oncology: more than meets the eye. Lancet Oncol 2015;16:e173-80. [Crossref] [PubMed]