Dynamic nomogram for long-term survival in patients with non-small cell lung cancer after pneumonectomy.

BACKGROUND: The study aims to identify prognostic factors of overall survival (OS) in patients who had pneumonectomy, in order to develop a practical dynamic nomogram model.
METHODS: A total of 2,255 patients with non-small cell lung cancer (NSCLC) who underwent pneumonectomy were identified from 2010-2015 in the Surveillance, Epidemiology, and End Results (SEER) database. The cohort was divided into a training (2011-2015) and a validation [2010] cohort. A nomogram and a risk classification system were constructed from the independent survival factors in multivariable analysis. The predictive accuracy of the nomogram was measured through internal and external validation.
RESULTS: Independent prognostic factors associated with OS were gender, age, pathology, tumor size, N stage, chemotherapy, and radiotherapy. The C-index of the nomogram for OS was 0.675 (95% CI: 0.655-0.694). Similarly, the AUC of the model was 0.733, 0.709, and 0.701 for the 1-, 3-, and 5-year OS, respectively. The calibration curves for survival demonstrated good agreement. Significant statistical differences were found in the OS of patients within different risk groups. An online calculation tool was established for clinical use.
CONCLUSIONS: This novel nomogram was able to provide a reliable prognosis for survival in patients with NSCLC undergoing pneumonectomy. 2021 Journal of Thoracic Disease. All rights reserved.

Introduction

Pneumonectomy was first successfully performed by Dr. Evarts Graham in 1933 to treat lung cancer, and it has been used as a radical surgical treatment ever since (1). Advances in medical technology and instruments have permitted wide acceptance of lobectomy, segmentectomy, and other surgical methods. Recently, sleeve lobectomy and pulmonary arterioplasty have also demonstrated similar oncological benefits as pneumonectomy in selected anatomically feasible patients (2). Pneumonectomy is performed in approximately 8.3% to 15% of all lung cancer surgeries (3,4). The procedure results in comparatively lower postoperative quality of life, and higher rates of complication, morbidity, and mortality than other surgical approaches. However, this operation is still the radical treatment strategy in centrally located or hilum invaded non-small cell lung cancer (NSCLC) (5,6).

The indications for pneumonectomy should be carefully considered because the procedure requires an extensive removal of lung tissue. Accordingly, accurate assessment of risk factors and the prediction of postoperative survival are of great importance for patients who are candidates for pneumonectomy. Due to the low proportion of pneumonectomy in thoracic surgery and the high postoperative risk, any prospective randomized controlled study can be challenging (7). However, retrospective analysis of the national database on prognostic factors may shed some light. This aggressive surgical approach makes it rare to study the prognostic factors on a large scale or to construct prospective studies (8,9). In recent years, nomograms are widely used for predicting prognosis. These predictive models generate individual quantified probabilities of clinical events by integrating prognostic clinicopathological variables. They are useful statistical prognostic models that bring us closer to achieving personalized medicine. Multiple studies have reported that nomogram scoring systems could provide an accurate prognosis of the disease (10-12).

In this study, we used the Surveillance, Epidemiology, and End Results (SEER) Program to analyze the prognostic factors of lung cancer patients who underwent pneumonectomy. A predictive dynamic nomogram was constructed to evaluate the survival and prognosis of patients with different clinical and pathological characteristics in the training set. The evaluation and verification of the model were carried out with data of the training and verification set, respectively.

We present the following article in accordance with the TRIPOD reporting checklist (available at http://dx.doi.org/10.21037/jtd-20-3203).

Methods

Study population

Data in the study were acquired from the SEER database of the National Cancer Institute in the United States. The registry data program collects and publishes cancer incidence and survival data from population-based cancer registries covering approximately 34.6% of the United States population (https://seer.cancer.gov/data/). The SEER*Stat software (version 8.3.8) was used to access the Incidence - SEER 18 Regs Custom Data (with additional treatment fields), Nov 2018 Sub (1975–2016 varying). Patients in the database were included if they had: (I) histology codes (International Classification of Diseases for Oncology, third edition, ICD-O-3) for adenocarcinoma: 8140, 8144, 8230, 8244, 8250-8255, 8260, 8310, 8323, 8333, 8470, 8480, 8481, 8490, 8550, 8551, 8574; squamous cell carcinoma: 8050, 8052, 8070-8074, 8083, 8084, 8123; neuroendocrine carcinoma: 8013, 8240, 8246, 8249; or other NSCLC: 8012; 8022; 8030-8033; 8035; 8046; 8051; 8082; 8200; 8430; 8560 (13); (II) “Primary Site-labeled” of any ICD-O-3 value ranging from C34.0 to C34.9 for a primary tumor site in the lung; (III) “Histology” with positive confirmation by diagnostic pathology; and (IV) “Therapy. Rx Sum--Surg Prim Site (1998+)” with values “55, 56, 65, 66, 70” for pneumonectomy. Furthermore, we selected patients who were diagnosed from 2010 to 2015 using the 7th AJCC Staging System. Patients with a diagnosis only confirmed by autopsy, 30-day operative mortality, or distant metastasis (M1) were excluded. Complete-case analysis is used to handle the missing data. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

Variables

Baseline demographic variables included age, gender, race, year of diagnosis, marriage, vital status, cause of death, and survival months. Tumor descriptors included laterality, histopathologic subtypes, grading, tumor size, T stage, number of lymph nodes dissection (LND), number of lymph nodes positive (LNP), N stage, TNM stage, surgery type, and additional therapy (e.g., chemotherapy and radiotherapy). According to the SEER variables dictionary, radical pneumonectomy was defined as pneumonectomy with mediastinal lymph node dissection, and extended pneumonectomy was radical pneumonectomy with the dissection of surrounding structures such as the diaphragm, pleura, and chest wall. The primary endpoint was overall survival (OS). The variable “Survival months” was used to identify survival time. OS was calculated by the difference in months from the diagnosis to death due to any reason listed under the variable “Vital status recode”. This study was a registry population-based research, so it did not involve blinding of research analysis and outcome or adjusted association between each candidate predictor and outcome.

Predictive model

According to the statement of transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD), the authors set eligible patients diagnosed in 2011–2015 as the internal training cohort, and the patients diagnosed in 2010 as the external validation cohort (14). The training cohort was used to establish the predictive model and to develop the nomogram. The validation cohort was used to validate this model. The Cox proportional hazards model was used for the univariate analysis of the OS in the training cohort. Any univariate result with a P value less than 0.05 was entered into a multivariate analysis. A hazard ratio (HR) and a corresponding 95% confidence interval (CI) were calculated. Based on results of the multivariate analysis, a nomogram for OS was constructed from its independent prognostic factors. Harrell’s concordance index (C-index) and the area under the time-dependent receiver operating characteristic curve (AUC) were used to measure the performance of the nomogram. Internal calibration plots with 400 bootstrap resamples from the training cohort was used to evaluate the nomogram performance, which compared the predicted and observed probabilities of 1-, 3-, and 5-year OS. The “rms” package was used to perform the external validation in the test cohort with the same predictors for the nomogram model from the training set. We divided the training group into high-risk group and low-risk group according to the median value of the total score for each patient, and the OS differences between low- and high-risk groups in the validation and total cohort were also evaluated. The authors analyzed the risk classification system based on the nomogram model to distinguish high-risk and low-risk groups in the enrolled patients with stage I–III NSCLC after pneumonectomy. An easy-to-use web nomogram model calculation was developed at https://shinyapps.io/ and the most updated version will also be available online.

Statistical analyses

Differences in continuous data were analyzed by the Mann-Whitney U test or Student’s t-test. Categorical variables were analyzed by Chi-squared or Fisher’s exact tests, as appropriate. Kaplan-Meier graphs were created for survival analysis. Statistical analyses were performed by SPSS (Version 24.0; IBM Corporation, Armonk, NY, USA). R Project (Version 4.0.3) was used to build and evaluate the performance of the nomogram. “Rms” package was used to build and evaluate the nomogram model. “Shiny” package was used to develop the dynamic version. All statistical tests were two-tailed, and a P value less than 0.05 was considered as statistically significant.

Results

Patient characteristics

Pneumonectomy accounts for about 3.9% (2,583/64,874) of all lung malignant tumor resections performed in the SEER database within the same study period. The proportion of pneumonectomy was reduced from 4.7% to 3.1% during 2011–2015. A total of 2,255 patients who met the inclusion criteria were included in the study. Among these patients, 62.5% were male, and the median age was 63 years old (IQR, 56–70 years old). Squamous cell carcinoma was the most common pathological type, accounting for 48.0%. The median tumor size was 5.0 cm (IQR, 3.2–7.0 cm). Lymph node metastasis was found in 59.6% of patients. Pneumonectomy was performed on the left side in 60.1% of the patients. Among all the patients, 70.5% underwent radical pneumonectomy including mediastinal lymph node dissection. Chemotherapy was given in 53.5% of the patients and 20.6% received radiotherapy. The training group contained 1,846 (81.9%) patients and the verification group contained 409 (18.1%) patients. Clinicopathological characteristics of the patients between the training group and the verification group are shown in .

Table 1

Patient characteristics in the study (N=2,255)

Variables	Total cohort (N=2,255)	Training cohort (N=1,846)	Validation cohort (N=409)	P
Gender, n (%)				0.773
Female	846 (37.5)	690 (37.4)	156 (38.1)
Male	1,409 (62.5)	1,156 (62.6)	253 (61.9)
Age				0.699
Median (IQR)	63 [56–70]	63 [56–70]	63 [56–70]
Mean (± SD)	62.73 (±10.46)	62.69 (±10.52)	62.91 (±10.23)
Race, n (%)				0.796
Caucasoid	1,926 (85.4)	1,575 (85.3)	351 (85.8)
Other	329 (14.6)	271 (14.7)	58 (14.2)
Marriage, n (%)				0.007
Married	1,341 (59.5)	1,122 (60.8)	219 (53.5)
Other	914 (40.5)	724 (39.2)	190 (46.5)
Pathology, n (%)				0.308
NEU	143 (6.3)	122 (6.6)	21 (5.1)
ADC	800 (35.5)	666 (36.1)	134 (32.8)
SQC	1,083 (48.0)	875 (47.4)	208 (50.9)
Other	229 (10.2)	183 (9.9)	46 (11.2)
Grading, n (%)				0.642
I	169 (7.5)	141 (7.6)	28 (6.8)
II	804 (35.7)	650 (35.2)	154 (37.7)
III/IV	1,066 (47.3)	873 (47.3)	193 (47.2)
Unknown	216 (9.6)	182 (9.9)	34 (8.3)
Laterality, n (%)				0.718
Left	1,355 (60.1)	1,106 (59.9)	249 (60.9)
Right	900 (39.9)	740 (40.1)	160 (39.1)
Tumor size (cm), n (%)				0.584
≤3	496 (22.0)	398 (21.6)	98 (24.0)
≤5	702 (31.1)	572 (31.0)	130 (31.8)
≤7	486 (21.6)	401 (21.7)	85 (20.8)
>7	512 (22.7)	423 (22.9)	89 (21.8)
Unknown	59 (2.6)	52 (2.8)	7 (1.7)
T stage (7th), n (%)				0.669
T1	230 (10.2)	181 (9.8)	49 (12.0)
T2	922 (40.9)	763 (41.3)	159 (38.9)
T3	702 (31.1)	575 (31.1)	127 (31.1)
T4	379 (16.8)	310 (16.8)	69 (16.9)
Unknown	22 (1.0)	17 (0.9)	5 (1.2)
N stage (7th), n (%)				0.547
N0	894 (39.6)	735 (39.8)	159 (38.9)
N1	853 (37.8)	705 (38.2)	148 (36.2)
N2	476 (21.1)	378 (20.5)	98 (24.0)
N3	15 (0.7)	13 (0.7)	2 (0.5)
Unknown	17 (0.8)	15 (0.8)	2 (0.5)
TNM stage (7th), n (%)				0.515
I	381 (16.9)	307 (16.6)	74 (18.1)
II	799 (35.4)	667 (36.1)	132 (32.3)
III	1,046 (46.4)	848 (45.9)	198 (48.4)
Unknown	29 (1.3)	24 (1.3)	5 (1.2)
Number of LND				0.060
Median (IQR)	14 [9–21]	14 [9–21]	13 [7–19]
Mean (± SD)	15.49 (±10.28)	15.68 (±10.34)	14.61 (±9.98)
Number of LNP				0.012
Median (IQR)	1 (0–3)	1 (0–3)	1 (0–3)
Mean (± SD)	1.93 (±3.14)	1.84 (±3.05)	2.32 (±3.53)
Surgery type, n (%)				0.166
Pneumonectomy	590 (26.2)	468 (25.3)	122 (29.8)
Radical Pneumo	1,589 (70.5)	1,314 (71.2)	275 (67.2)
Extended Pneumo	76 (3.4)	64 (3.5)	12 (2.9)
Chemotherapy				0.604
Yes	1,206 (53.5)	992 (53.7)	214 (52.3)
No/unknown	1,049 (46.5)	854 (46.3)	195 (47.7)
Radiotherapy				0.621
Yes	465 (20.6)	377 (20.4)	88 (21.5)
No	1,790 (79.4)	1,469 (79.6)	321 (78.5)

ADC, adenocarcinoma; SQC, squamous cell carcinoma; NEU, neuroendocrine carcinoma; Pneumo, pneumonectomy; LND, lymph nodes dissection; LNP, lymph nodes positive.

Survival analyses and prediction model

The 1-, 3-, and 5-year OS of pneumonectomy were 75.9%, 54.9%, and 44.4%, respectively. The estimated median OS time was 46.0 months for the whole cohort. The median follow-up time was 44.5 months. As shown in , the univariate Cox regression model revealed that gender, age, pathology type, grading, laterality, tumor size, N stage, number of LNP, chemotherapy, and radiotherapy were significantly associated with OS. Multivariate analysis further confirmed that gender, age, pathology type, laterality, tumor size, N stage, chemotherapy, and radiotherapy were independent factors of long-term survival. The prognostic significant factors were included in the construction of the predictive model for the nomogram (). An online version of the nomogram to assist researchers and clinicians could be accessed at https://thoracic.shinyapps.io/nomogrampneumonectomy/. Predicted survival probability across time could be easily determined by inputting clinical features and reading output figures and tables generated by the webserver (). The C-index of OS, which indicates discrimination ability, was 0.675 (95% CI: 0.655–0.694). Similarly, the AUC of the prediction model was 0.733, 0.709, and 0.701 for the 1-, 3-, and 5-year OS, respectively (). These findings indicated that the nomogram had moderate predictive ability for OS. As shown in and 3, the internal () and external () calibration plots of the predictive model for the 1-, 3-, and 5-year OS were in good agreement.

Table 2

OS univariate and multivariate analysis of prognostic factors in the training cohort.

Variables	Univariate analysis of OS			Multivariate analysis of OS
	HR (95% CI)	P		HR (95% CI)	P
Gender
Male vs. female	1.392 (1.204–1.609)	<0.001		1.201 (1.021–1.414)	0.027
Age	1.026 (1.019–1.033)	<0.001		1.025 (1.017–1.034)	<0.001
Race
Other vs. Caucasoid	1.021 (0.844–1.236)	0.828
Marriage
Other vs. married	1.120 (0.977–1.286)	0.105
Histology		<0.001
ADC vs. NEC	2.437 (1.632–3.640)	<0.001		1.995 (1.226–3.247)	0.005
SQC vs. NEC	2.691 (1.811–3.999)	<0.001		2.131 (1.305–3.480)	0.002
Other vs. NEC	2.811 (1.819–4.345)	<0.001		2.217 (1.303–3.774)	0.003
Grading		<0.001
II vs. I	1.606 (1.153–2.238)	0.005		0.977 (0.681–1.402)	0.899
III & IV vs. I	2.193 (1.587–3.030)	<0.001		1.281 (0.896–1.833)	0.175
Laterality
Right vs. left	1.379 (1.204–1.580)	<0.001		1.472 (1.267–1.709)	<0.001
Tumor size		<0.001
≤5 vs. ≤3 cm	1.232 (1.003–1.513)	0.046		1.165 (0.931–1.456)	0.182
≤7 vs. ≤3 cm	1.322 (1.062–1.644)	0.012		1.254 (0.988–1.591)	0.063
>7 vs. ≤3 cm	1.944 (1.581–2.390)	<0.001		1.792 (1.427–2.250)	<0.001
N stage		<0.001
N1 vs. N0	1.468 (1.252–1.722)	<0.001		1.559 (1.291–1.883)	<0.001
N2 vs. N0	1.719 (1.435–2.060)	<0.001		1.906 (1.519–2.391)	<0.001
N3 vs. N0	2.809 (1.445–5.459)	0.002		4.482 (2.128–9.438)	<0.001
Number of LND	0.998 (0.991–1.004)	0.487
Number of LNP	1.033 (1.019–1.048)	<0.001		1.016 (0.995–1.037)	0.127
Surgery type		0.168
Radical Pneumo vs. Pneumo	0.911 (0.781–1.063)	0.237
Extended Pneumo vs. Pneumo	1.215 (0.843–1.751)	0.297
Chemotherapy
Yes vs. no/unknown	0.769 (0.672–0.880)	<0.001		0.558 (0.471–0.660)	<0.001
Radiotherapy
Yes vs. no	1.244 (1.061–1.458)	0.007		1.428 (1.168–1.747)	0.001

OS, overall survival; HR, hazard ratio; CI, confidence interval; ADC, adenocarcinoma; SQC, squamous cell carcinoma; NEC, neuroendocrine carcinoma; Pneumo, pneumonectomy; LND, lymph nodes dissection; LNP, lymph nodes positive.

Figure 1

The predictive model developed in this study. Nomogram (A) and online calculation tool (B) predicting the OS for patients undergoing pneumonectomy diagnosed NSCLC. OS, overall survival; NSCLC, non-small cell lung cancer.

Figure 2

The ROC curve and calibration plots of the predictive model in the training cohort. The area under the ROC curve of 1-, 3-, and 5-year OS in the training cohort were 0.733, 0.709 and 0.701, respectively (A,B,C); Calibration plots comparing actual and predicted overall survival probabilities at 1-, 3- and 5-year follow-up in the training set (D). Perfect prediction would correspond to a slope of 1 (diagonal 45-degree gray line). ROC, receiver operating characteristic.

Figure 3

Calibration plots of the predictive model in the validation cohort. Calibration plots comparing actual and predicted overall survival probabilities at 1-, 3- and 5-year (A,B,C) follow-up in the validation set.

Risk classification system

The novel risk classification system placed patients into the low-risk (934/1,846; score ≤203) or high-risk group (912/1,846; score >203) by the median value of the total points for each patient in the training cohort (). In the total cohort, the medians OS of patients in the low- and high-risk groups were 75.0 and 26.0 months, respectively (Log-rank P<0.001, ). The Kaplan-Meier curves showed that the OS of the two groups can be significantly differentiated by the risk classification system both in the training and validation cohorts (Log-rank P<0.001, ). There were also significant survival differences between the low- and high-risk groups in stage I, II, or III lung cancer (Log-rank P<0.05, ). The high-risk patients in stage I and II still had worse survival outcome than those in stage II and III with low-risk (3-year mortality: stage I/high-risk vs. stage II/low-risk: 46.0% vs. 29.4%, Log-rank P=0.009; stage II/high-risk vs. stage III/low-risk: 57.0% vs. 39.6%, Log-rank P<0.001, ). However, the prognosis difference between patients in stage I/high-risk and stage III/low-risk groups did not show significant statistical differences (5-year mortality: stage I/high-risk vs. stage III/low-risk: 52.9% vs. 51.8%, Log-rank P=0.456, ).

Figure 4

The OS Kaplan-Meier curves of patients in low- and high-risk groups. In the total (A), training (B), and validation (C) cohort, patients in low-risk group had significantly superior OS than the high-risk group (Log-rank P<0.001). OS, overall survival.

Figure 5

The cumulative overall mortality of patients in different stages. Patients with low-risk had significantly better survival than patients with high-risk in TNM stage I (A), II (B), and III (C) (Log-rank P<0.05). The patients in stage I and II with high-risk had significantly worse survival than stage II and III with low-risk (Stage I/high-risk vs. Stage II/low-risk, Log-rank P=0.009; Stage II/high-risk vs. Stage III/low-risk, Log-rank P<0.001), but the patients in stage I/high-risk and stage III/low-risk groups had similar cumulative overall mortality without significant statistical difference (Log-rank P=0.456) (D).

Discussion

This study focused on the construction of a survival prediction model for patients who had undergone pneumonectomy, which systematically analyzed the clinicopathological factors of patients identified in the SEER database. We obtained multiple independent prognostic variables to build a reliable predictive nomogram, a convenient online tool, and an effective risk classification system.

In recent years, there have been improvements in anesthesia, surgical techniques, and perioperative nursing, each contributing to a reduction in mortality associated with pneumonectomy (15,16). Pneumonectomy is a high-risk surgical procedure accounting for about 15% of lung cancer surgeries (4). During the period of this study, the proportion of pneumonectomy was only 3.9% among all lung cancer resections and with a reduction tendency during 2011–2015, which was lower than the data from the European Society of Thoracic Surgeons (ESTS) (7.4%) between 2010 and 2013 (17). Similar studies from English and Danish registries data showed a significant decrease in pneumonectomy proportion over time (18,19). Pneumonectomy will still be used as an inevitable and effective method of treatment for some patients with centrally located or locally advanced lung cancer (20).

Previously, researchers have created the Thoracic Surgery Scoring System (Thoracoscore) using the French Epithor database, to evaluate the prognosis of patients undergoing thoracic surgery (21,22). This system has later been included in the guidelines for the British Thoracic Society (23). However, one study found that Thoracoscore was not reliable in predicting mortality after pneumonectomy because of the under-representation of this high-risk procedure in the database (6%) (24). Several nomograms for lung cancer have been developed, but few systemically predicted outcomes after pneumonectomy. For instance, Cheng and colleagues constructed a nomogram prediction model for the prognosis of patients undergoing pneumonectomy from a cohort of 100 patients (25). This model incorporated circulating blood biomarkers and tumor characteristics, and demonstrated prognostic superiority over the pTNM staging system. But, the calibration curve of the study did not show good agreement in the discovery and validation cohorts, which is likely due to its small sample size. Compared with the current only nomogram study on patients post pneumonectomy, our research has a bigger sample size, and the nomogram produced good results in both training and validation cohorts. For patients undergoing pneumonectomy in the same cancer stage, the risk classification system can significantly distinguish between the high- and low-risk groups. Considering the small proportion of this highly risk procedure, it was difficult to conduct randomized studies. Therefore, more researches based on the prospective designed database including multicenter was expected to in-depth evaluate the effects of different pathological indicators on short- and long-term prognosis.

In this study, an effective risk classification system was constructed by nomogram prediction model. The prognosis of patients in each stage of lung cancer in high-risk group was significantly worse than that in low-risk group. There was no significant difference in long-term overall mortality between patients in stage I/high-risk group and stage III/low-risk group. Based on our risk classification system, we could more accurately identify high-risk patients in different stages of lung cancer and provide more rigorous follow-up and care. The authors also provided an online version of the dynamic nomogram which was easier-to-use for physicians to predict the long-term survival of patients with different clinicopathological features.

Pneumonectomy laterality has an impact on the long-term survival of patients. In this study, we found that patients with a right-side pneumonectomy presented worse prognosis than those with a left-side pneumonectomy. This observation may be related to several factors. First, 55% of ventilation is performed by the right lung; right-side pneumonectomy will thus proportionately affect the overall lung function (1). Secondly, mediastinal displacement with the asymmetric position of the heart to the left also results in relief of the left side after pneumonectomy. The right side lacks a corresponding resistance to mediastinal displacement, which may have a negative impact on hemodynamic stability and pose a possible life-threatening risk (26). In addition, some researchers also reported that the proportion of serious postoperative complications is higher on the right side (acute respiratory distress syndrome and bronchopleural fistula), which leads to higher morbidity and mortality. This is due to a lack of retraction on the right stump as opposed to the left, which makes the right stump more prone to inflammation and breakdown (27).

In this study, we found that the use of chemotherapy improved OS of patients while the radiotherapy showed a negative effect. Possible explanation may be related to the fact that chemotherapy serves as a more systemic treatment, while radiotherapy is a local treatment. In addition, radiotherapy can increase the burden on cardiopulmonary function, leading to serious complications and reduced long-term survival. Researchers also found that chemotherapy improves the prognosis of patients with pneumonectomy (28). Kim et al. demonstrated that forced expiratory volume in 1 second (FEV1) at the beginning of postoperative radiotherapy—that is closely related to the extent of resection—was the only significant prognostic factor for OS (29). Given that patients undergoing pneumonectomy usually have reduced lung function, Karnofsky performance status scores and cardiopulmonary tolerance should be carefully evaluated and optimized when considering radiotherapy or chemotherapy.

However, the current study has several limitations. First, this study was a retrospective study with common defects. Second, only the variables available in the SEER database could be analyzed with no access to more detailed information such as the Charlson comorbidity index, pulmonary function evaluation, postoperative complications, the chemotherapy regimen, and the surgery approaches (e.g., open or video-assisted thorascopic surgery). The missing information also brought effects on the accuracy of the nomogram. Third, due to the lack of tumor specific site and invading information, we are unable to fully utilize the 8th staging pT criteria to evaluate the primary tumor in the prediction model.

Pneumonectomy accounted for about 3.9% of lung cancer operations in the same period. The nomogram was developed from a number of independent prognostic clinicopathological variables and demonstrated its good utility in C-index and AUC results. The predictive model also had a satisfactory fit to the calibration curves with both internal and external verification. The nomogram and online calculator have essential practical significance in clinical settings, which can help physicians to evaluate the long-term survival of patients and optimize the personalized treatments for lung cancer.

The article’s supplementary files as