Skip to main content

Stratifying non-small cell lung cancer patients using an inverse of the treatment decision rules: validation using electronic health records with application to an administrative database

Abstract

Background

To validate a stratification method using an inverse of treatment decision rules that can classify non-small cell lung cancer (NSCLC) patients in real-world treatment records.

Methods

(1) To validate the index classifier against the TNM 7th edition, we analyzed electronic health records of NSCLC patients diagnosed from 2011 to 2015 in a tertiary referral hospital in Seoul, Korea. Predictive accuracy, stage-specific sensitivity, specificity, positive predictive value, negative predictive value, F1 score, and c-statistic were measured. (2) To apply the index classifier in an administrative database, we analyzed NSCLC patients in Korean National Health Insurance Database, 2002–2013. Differential survival rates among the classes were examined with the log-rank test, and class-specific survival rates were compared with the reference survival rates.

Results

(1) In the validation study (N = 1375), the overall accuracy was 93.8% (95% CI: 92.5–95.0%). Stage-specific c-statistic was the highest for stage I (0.97, 95% CI: 0.96–0.98) and the lowest for stage III (0.82, 95% CI: 0.77–0.87). (2) In the application study (N = 71,593), the index classifier showed a tendency for differentiating survival probabilities among classes. Compared to the reference TNM survival rates, the index classification under-estimated the survival probability for stages IA, IIIB, and IV, and over-estimated it for stages IIA and IIB.

Conclusion

The inverse of the treatment decision rules has a potential to supplement a routinely collected database with information encoded in the treatment decision rules to classify NSCLC patients. It requires further validation and replication in multiple clinical settings.

Peer Review reports

Background

Lung cancer is one of the most common cancers and is the leading cause of cancer death worldwide [1, 2]. Stratification of lung cancer phenotypes is an essential step for treatment decision making, and TNM classification has an important role for therapeutic and prognostic guidance [3, 4]. As TNM staging is based on numeric dimensions (i.e., tumor size and location) and multinomial dimensions (i.e., invasion of lymph nodes and major organs), the TNM classification method projects the original feature space into a few tens of categories [5,6,7].

Increasingly, routinely collected administrative healthcare databases provide new opportunities in conducting large-scale analysis of actual clinical practice data in real-world settings [8,9,10,11]. However, the routinely collected databases tend to include only codified data elements with existing ontology, and they are prone to misclassification and missing information [12]. Therefore, efficient method to identify decision bases are necessary for studying large scale databases.

Observed, or revealed, decision actions can be used to estimate the information that the decision-makers used [13, 14]. The current study validated a stratification method for non-small lung cancer (NSCLC) patients based on an inverse of the treatment decision rules, which can be applied to codified data elements in administrative databases [11]. We validated the index classifier from two aspects. First, we examined the predictive performance when the index classification was compared to the TNM stage in the electronic health records (EHRs). Second, we examined whether the index classification shows differential survival rates when applied to a population-based administrative database.

Methods

The index classification: an inverse of the treatment decision function

The index classification is based on an inverse of the treatment decision rules given in the recommended treatment regimens in the National Comprehensive Cancer Network (NCCN) guidelines [4, 11]. The inverse of the stage-based treatment decision rules maps from the treatment regimens into patients' status. (Additional file 1: Tables S1, S2 and Methods 1). We only considered a subdomain of treatment patterns that are invertible and had sufficient samples. Seven potential categories for the inverse of the treatment decision function are shown in Table 1. Note that the inverse of the stage-based treatment decision function is pre-specified based on clinical knowledge, without any reference to the validation data.

Table 1 Potential categories for the inverse function of the treatment decision

Validation study for the index classification using electronic health records

We retrospectively validated the index classifier against the reference standard staging information in the EHRs of Asan Medical Center (AMC), a tertiary referral hospital in Seoul, Korea. The study protocol using EHR data was approved by the AMC, and all methods were performed in accordance with the relevant guidelines and regulations. We adopted the estimated sample size of 108 for each category for testing the accuracy of a single modality to detect a pre-specified area under the receiver operating characteristic curve value of 0.8 against a null value of 0.7 with a 95% confidence level and 80% power [15]. We identified a consecutive series of patients who were newly diagnosed with NSCLC and staged between 2011 and 2015, covered under the national health insurance, and received cancer treatment in the hospital with their last visit more than 180 days from the end of observation. We included those aged between 20 and 75 without a preexisting non-pulmonary cancer or a cardiopulmonary comorbidity, who may have been treated according to the treatment guidelines (Additional file 1: Methods 2) [16]. We excluded cases for which complete information for the index classification were not available: initially diagnosed at an external hospital, transferred to an external hospital, missing stage, erroneous stage identified during the manual chart review process, or received no treatment. In the sensitivity analysis, those with cardiopulmonary comorbidities were included in addition to the participants in the main analysis.

The reference stage, which was based on the TNM 7th edition [5] during the study period, was assigned by manual chart review for all of the study participants by a clinician who was blinded to the modeling to minimize potential bias. The stage was then inferred by a researcher who was blinded to the reference standard staging information, using the diagnosis codes, procedure codes, and medication codes.

Overall accuracy was calculated as the proportion of cases with correctly classified stages among the entire study population (Additional file 1: Methods 3). Stage-specific sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), F1 score, and c-statistic were calculated by contrasting each stage with alternative stages [17,18,19]. The 95% confidence intervals for the performance metrics were estimated with 1000 bootstrap resamples.

Stage-specific survival analysis applying the index classification in a population-based administrative database

We applied the index classifier to Korean National Health Insurance database, 2002–2013, retrospectively [20, 21]. We had access to treatment and survival information for the entire cohort of NSCLC patients and to the full claims records for a 2% random sample. We identified a consecutive series of all patients with the diagnosis code for lung cancer with at least one claim code for the treatment of NSCLC between 2004 and 2013. Patients with claim codes for lung cancer between 2002 and 2003 were excluded for a washout period, as they may have been diagnosed with lung cancer before 2002. We considered overall survival using the linked records in the national death registry with any cause. Class-specific survival rates were compared to the reference survival rates in the TNM 7th edition [22].

The proportional hazard assumption was rejected from the independency test for Schoenfeld residuals against time (Additional file 1: Methods 4) [23]. The survival function was estimated with the non-parametric Kaplan–Meier method, and the log-rank test was used to test the null hypothesis of no difference in the survival between groups. Pairwise log-rank tests were performed to explore significantly different survival curves, and the significance level of 0.05 was adjusted with Bonferroni methods. The 95% confidence intervals for the survival estimates were computed with the normal approximation. Analyses were performed using R statistical software version 3.5.3 (R Foundation for Statistical Computing) and SAS version 9.4 (SAS Institute Inc., NC, USA).

We considered the validation study as the main study, and this report follows the STARD reporting guideline for diagnostic studies [24, 25].

Results

Validation study for the index classification using the electronic health records

The study population for the validation study consisted of 1375 NSCLC patients with a mean age of 60.3 (sd = 9.3) and a male-to-female ratio of 1.86. The population selection diagram for the validation study is shown in Fig. 1. The stage distribution was 68.6%, 8.9%, 7.9%, and 14.5% for stage I, II, III, and IV, respectively.

Fig. 1
figure 1

Population selection diagram of the validation study. To validate the index classifier against the TNM 7th edition, a retrospective study was conducted with electronic health records (EHRs) of NSCLC patients in a tertiary referral hospital in Seoul, Korea. We identified a consecutive series of all patients who were newly diagnosed with NSCLC and staged between 2011 and 2015, covered under the national health insurance, and received cancer treatment in the hospital with their last visit more than 180 days from the end of observation. We excluded cases initially diagnosed at an external hospital, transferred to an external hospital, missing stage information, erroneous stage identified during the manual chart review process, received no treatment, missing age information, age of less than 20 years or more than 75 years, preexisting non-pulmonary cancer, or cardiopulmonary comorbidity (eMethods 2)

In the validation study, the overall accuracy of the index classification was 93.8% (95% CI: 92.5–95.0%). Stage I showed the highest stage-specific sensitivity (0.99, 95% CI: 0.98–0.99), PPV (0.97, 95% CI: 0.97–0.98), F1 score (0.98, 95% CI: 0.98–0.99), and c-statistic (0.97, 95% CI: 0.96–0.98), while specificity was highest for stage III (0.99, 95% CI: 0.99–0.99), and NPV was highest for stage IV (0.99, 95% CI: 0.98–0.99) (Table 2). On the other hand, stage-specific sensitivity (0.64, 95% CI: 0.55–0.74), F1 score (0.77, 95% CI: 0.70–0.84), and c-statistic (0.82, 95% CI: 0.77–0.87) were lowest for stage III, and PPV was the lowest for stage II (0.74, 95% CI: 0.67–0.81). A confusion matrix comparing the true stage and the inferred stage in the validation study is shown in Additional file 1: Table S3.

Table 2 Stage-specific predictive performance metric of the index classifier against the reference standard (TNM stage) (N = 1375)

In the sensitivity analysis, 198 NSCLC patients with cardiopulmonary comorbidity were included in the analysis (N = 1573). The overall accuracy of the index classification was 93.5% (95% CI: 92.3–94.7%). Sensitivity, specificity, PPV, NPV, and c-statistic were all within the confidence limits of the main analysis. Results for the sensitivity analysis are shown in Additional file 1: Tables S4, S5.

Stage-specific survival analysis applying the index classification in a population-based administrative database

The study population for the application study consisted of 71,593 NSCLC patients with a mean age of 64.3 (sd = 10.2) and a male-to-female ratio of 72:28. Among a total of 166,203 patients who had a diagnosis code and at least one procedure or medication code for the treatment of NSCLC in the national health insurance database, 71,593 patients were eligible for analysis. Of these patients, 32.9% (23,571/71,593) received chemotherapy alone, 29.9% (21,378/71,593) received chemotherapy and radiotherapy, and 24.5% (17,543/71,593) underwent surgical resection alone.

When applying the index classifier for NSCLC patients, the independency test for Schoenfeld residuals against time rejected the proportional hazard assumption (p = 0.043). The global null hypothesis of no difference in the survival rate was rejected using the log-rank test (p < 0.001). Pairwise log-rank tests corrected with Bonferroni methods rejected most null hypotheses, except those comparing class 2 vs. class 5, and class 3 vs. class 4 (Additional file 1: Table S6). Class-specific survival curves (A) and the inferred stage-specific (B) survival curves are shown in Fig. 2. The class 1 corresponds to multiple stages (i.e., IA, IB, IIA) and stage IIIA corresponds to multiple classes (i.e., 3, 4, 5), so the comparability is limited. Those who received surgery only (class 1) showed lower survival for the initial months, but higher long-term survival compared with those who received combined modality (class 2, 3, 4, 5). Within stage IIIA, those with neoadjuvant chemotherapy and surgery (class 5) showed long-term survival similar to that of stage IIB (class 2). When comparing the reference survival rates by the TNM classification, the class-specific survival rates based on the index algorithm under-estimated the survival probability for stage IIIB and IV, and over-estimated it for stage IIB (Table 3).

Fig. 2
figure 2

Class-specific (A) and stage-specific (B) survival curves with 95% confidence intervals and risk tables. A The index classification using an inverse of treatment decision function was applied to NSCLC cancer patients in Korean National Health Insurance database (N = 71,593), and class-specific survival function was estimated with the non-parametric Kaplan–Meier method. B The cancer stages (TNM 7th edition) were inferred from the index classification, and stage-specific survival function was estimated with the non-parametric Kaplan–Meier method. Stage I and IIA are aggregated in the survival curve because the index classification cannot distinguish them

Table 3 Two-year and five-year survival by TNM classification (reference) vs index classification (N = 71,593)

Discussion

We proposed a stratification method using an inverse of the treatment decision rules that can classify lung cancer patients with real-world treatment records that are commonly available in routinely collected administrative databases. To validate the index classifier, we first evaluated the predictive performance of the index classifier against the reference standard staging information in an EHR database of a tertiary referral hospital. In addition, we examined whether the index classification showed differential survival rates when applied to a nationally representative administrative database that covers the hospital in the validation study. To our knowledge, this is the first study to apply and validate a stratification method based on an inverse of the treatment decision rules using both EHRs and a nationwide administrative database.

The index classification method considered treatment pattern-based categories as the baseline stratifying variable, rather than as a response or outcome variable. Stratifying phenotypes based on potential treatment options can be a practical way of defining patient categories. Modern classification of diseases started from Carl Linnaeus's classification (1707–1778), which influenced the international classification of disease and cause of death (ICD) [26, 27], which is an ontology that defines formal semantic relationships between concepts [27]. The TNM stage classification of cancer is a nomenclature describing the anatomic extent of cancer with limited predictive ability [28]. On the other hand, treatment decisions are made based on causal inference, or counterfactual prediction of an outcome given a treatment [29]. The classification of diseases has a major role in treatment decisions, but when its predictive ability is limited, additional qualitative clinical judgement is required [4, 28].

The index classification method used in the current study supplements the disease ontology information using both treatment information and treatment decision rules. The treatment decision rules provide the structural linkages between the disease information and the treatment information by the totality of the evidence in the literature. Utilization of treatment patterns enables the classification to be more interpretable and transparent to clinicians, and therefore, satisfies an essential need to deliver clinical impact with a data-driven system [30]. Treatment decision-making could also be seen as the ultimate goal of diagnostic development, which is denoted by the term "theragnostics" in certain fields [31, 32].

A practical advantage of the classification using an inverse of the treatment decision rules is that it can be applied to most routinely collected healthcare databases that have codified information of diagnoses, procedures, and treatments. Therefore, the classification can be applied to data generated while delivering care for building a “learning system” [33,34,35]. When the treatment decision rules are used for the baseline stratifying variable to be examined against future health outcomes, the results may feedback to update the treatment decision rules, enabling a continuous learning system. The index classification mitigates some known challenges in analyzing observational routinely-collected health data [12, 36, 37] by utilizing data for major procedures, such as surgery, chemotherapy, and radiotherapy [38]. Modeling the inverse of the decision function is a way of estimating the decision bases (i.e., patients' status, including stage and performance status) when the direct measurement is not available. This is in line with the estimation method of unmeasured decision bases using observed, or revealed, decision actions in economics [13] and statistics [14]. Treatment decisions can essentially encode such information, and therefore, it can be used as a proxy for such information [39], given a quantitative evaluation regarding the degree of bias and error. Therefore, it is an important first step to evaluate the predictive performance of the index classifier against the reference standard staging information in the EHRs.

Another advantage of the classification using an inverse of the treatment decision rules is that it may be considered as a phenotypical classification differentiating survival rates. In this sense, the treatment pattern for a patient can be considered as a decision function that maps from both the patient's performance status and the cancer stage into the treatment. Then, the survival prognosis from the patient's performance status and the cancer stage may be approximated by the classifier based on an inverse of the treatment decision rules. This approach may be useful when information on a patient's performance status is not readily available.

In the validation study using the EHRs, the index classifier showed high specificity for all stages of NSCLC. However, the sensitivity was particularly low for stage III (Table 2). There is especially high uncertainty when the treatment involves radiotherapy, which can be used for different purposes: it can be used as adjuvant or neoadjuvant therapy for surgical resection in early-stage disease, or it can be used as palliative treatment in combination with chemotherapy in advanced stage disease.

Differential class-specific survival rates have been used as a justification for the TNM classification of lung cancer [28, 40]. The methodological principles for development and validation of newer edition of the TNM classification are to make an ordinal classification system such that within-group variance is minimized while between-group variance is maximized, while not compromising the former edition. For lung cancer, the prognosis for survival was chosen as the measure for the clustering of the stage groups, while having the discriminatory power maintained across different factors, such as time, geographical location, diagnostic method, histology, and patient characteristics. Therefore, we examined whether the index classification results in differential class-specific survival rates in the underlying population. When the index classifier was applied to the nationwide administrative database, the class-specific survival rates based on the index algorithm under-estimated the survival probability for TNM stage IIIB and IV and over-estimated it for stage IIB. The discrepancies may rise from misclassifications of patients who received non-standard treatments due to unique clinical situations. The initially low survival of those who underwent surgery only (class 1) can be a result of immediate postoperative complications [41]. Observation of stage IIIA patients with neoadjuvant chemotherapy and surgery (class 5) showing long-term survival similar to stage IIB patients (class 2) is hypothesis-generating: this may be related to non-significant differences between adjuvant and neoadjuvant chemotherapy [42,43,44]. Additional studies need to confirm the long-term outcomes beyond five years among cancer survivors [45].

This study has several limitations. First, this study considered only a subdomain of treatment patterns that are invertible. When considering the decision functions with many-to-many relationships, inversion can be erroneous. The index classification showed some discrepancies compared to the TNM stage classification, and therefore, potential biases and errors should be addressed appropriately when it is used as a proxy variable. The present study was not meant to develop a precise prediction model, which needs to incorporate all important predictive factors, such as histology [46], genotype [47], environmental factors [48], and behavioral factors, such as smoking status [49]. Variations in healthcare practice need to be accounted for appropriately, and additional validation may be required before application to a new setting. Regional variations, such as the proportion of smokers [50, 51], genotype distribution [52, 53], and screening policy [54, 55], should also be considered when developing a precise prediction model. The utility of the index classification can be decreased in settings where treatment decision rules are less likely to be strictly followed, such as in budget-constrained settings.

Despite these limitations, the stratification of NSCLC patients using an inverse of the treatment decision rules can supplement a routinely collected database with information encoded in treatment decision rules to classify NSCLC patients by providing phenotypical classification that differentiates survival prognosis. Clinical prognostic value of the index classification needs to be validated and replicated in multiple clinical settings.

Conclusion

We validated a stratification method of NSCLC patients using an inverse of the treatment decision rules in comparison to the TNM stage classification and demonstrated that the index classification has a potential to differentiate survival probabilities when applied to a population-based administrative database. The inverse of the treatment decision rules has a potential to supplement a routinely collected database with information encoded in treatment decision rules to classify NSCLC patients. It requires further validation and replication in multiple clinical settings.

Availability of data and materials

The data used in the validation study were provided by the Asan Medical Center (AMC). The data used in the application study were provided by the National Health Information Database (NHIS-2017-4-003) managed by the National Health Insurance Service (NHIS) of Korea. The data from AMC that are analysed during the current study are not publicly available due the institution's privacy regulations but are available from the corresponding author on reasonable request. The data from NHIS that are analysed during the current study are available from NHIS, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of NHIS.

References

  1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  Google Scholar 

  2. Jung K-W, Won Y-J, Kong H-J, Lee ES. Community of population-based regional cancer registries. Cancer statistics in Korea: incidence, mortality, survival, and prevalence in 2015. Cancer Res Treat Off J Korean Cancer Assoc. 2018;50(2):303–16.

    Google Scholar 

  3. Brierley JD, Gospodarowicz MK, Wittekind C. TNM classification of malignant tumours. John Wiley & Sons; 2016.

  4. Ettinger DS, Bepler G, Bueno R, Chang A, Chang JY, Chirieac LR, et al. Non-small cell lung cancer: clinical practice guidelines in oncology™. JNCCN J Natl Compr Cancer Netw. 2006;4(6):548–82.

    Article  CAS  Google Scholar 

  5. Goldstraw P, Crowley J, Chansky K, Giroux DJ, Groome PA, Rami-Porta R, et al. The IASLC lung cancer staging project: proposals for the revision of the TNM stage groupings in the forthcoming (seventh) edition of the TNM Classification of malignant tumours. J Thorac Oncol Off Publ Int Assoc Study Lung Cancer. 2007;2(8):706–14.

    Google Scholar 

  6. Goldstraw P, Chansky K, Crowley J, Rami-Porta R, Asamura H, Eberhardt WEE, et al. The IASLC lung cancer staging project: proposals for revision of the TNM stage groupings in the forthcoming (Eighth) edition of the TNM classification for lung cancer. J Thorac Oncol Off Publ Int Assoc Study Lung Cancer. 2016;11(1):39–51.

    Google Scholar 

  7. Detterbeck FC, Boffa DJ, Kim AW, Tanoue LT. The eighth edition lung cancer stage classification. Chest. 2017;151(1):193–203.

    Article  Google Scholar 

  8. Fleurence RL, Curtis LH, Califf RM, Platt R, Selby JV, Brown JS. Launching PCORnet, a national patient-centered clinical research network. J Am Med Inform Assoc JAMIA. 2014;21(4):578–82.

    Article  Google Scholar 

  9. Hripcsak G, Ryan PB, Duke JD, Shah NH, Park RW, Huser V, et al. Characterizing treatment pathways at scale using the OHDSI network. Proc Natl Acad Sci U S A. 2016;113(27):7329–36.

    Article  CAS  Google Scholar 

  10. Sherman RE, Anderson SA, Dal Pan GJ, Gray GW, Gross T, Hunter NL, et al. Real-world evidence—what is it and what can it tell us? N Engl J Med. 2016;375(23):2293.

    Article  Google Scholar 

  11. Kim M-H, Park S, Park YR, Ji W, Kim S, Choo M, et al. Developing computational phenotyping algorithm to classify lung cancer stages using common data elements [Internet]. In: American Medical Informatics Association 2019 Annual Symposium; 2019 Nov 18; Washington, D.C. Available from: https://www.amia.org/amia2019

  12. Haut ER, Pronovost PJ, Schneider EB. Limitations of administrative databases. JAMA. 2012;307(24):2589.

    Article  CAS  Google Scholar 

  13. Samuelson PA. Consumption theory in terms of revealed preference. Economica. 1948;15(60):243–53.

    Article  Google Scholar 

  14. Swartz RJ, Cox DD, Cantor SB, Davies K, Follen M. Inverse decision theory: characterizing losses for a decision rule with applications in cervical cancer screening. J Am Stat Assoc. 2006;101(473):1–8.

    Article  CAS  Google Scholar 

  15. Hajian-Tilaki K. Sample size estimation in diagnostic test studies of biomedical informatics. J Biomed Inform. 2014;48:193–204.

    Article  Google Scholar 

  16. Driessen EJM, Schulkes KJG, Dingemans A-MC, Van Loon JGM, Hamaker ME, Aarts MJ, et al. Patterns of treatment and survival among older patients with stage III non-small cell lung cancer. Lung Cancer Amst Neth. 2018;116:55–61.

    Article  Google Scholar 

  17. Hand DJ, Till RJ. A simple generalisation of the area under the ROC curve for multiple class classification problems. Mach Learn. 2001;45(2):171–86.

    Article  Google Scholar 

  18. Davis J, Goadrich M. The relationship between precision-recall and ROC curves. In: Proceedings of the 23rd international conference on Machine learning. ACM; 2006. pp. 233–40.

  19. Pencina MJ, D’Agostino RB. Evaluating discrimination of risk prediction models: the C statistic. JAMA. 2015;314(10):1063–4.

    Article  CAS  Google Scholar 

  20. Kim L, Kim J-A, Kim S. A guide for the utilization of health insurance review and assessment service national patient samples. Epidemiol Health. 2014;30:e2014008.

    Article  Google Scholar 

  21. Lee J, Lee JS, Park SH, Shin SA, Kim K. Cohort profile: the national health insurance service-national sample cohort (NHIS-NSC), South Korea. Int J Epidemiol. 2017;46(2):e15. https://doi.org/10.1093/ije/dyv319

    Article  Google Scholar 

  22. Choi CM, Kim HC, Jung CY, Cho DG, Jeon JH, Lee JE, Ahn JS, Kim SJ, Kim Y, Choi YD, Suh YG, Kim JE, Lee B, Won YJ, Kim YC. Report of the Korean Association of Lung Cancer Registry (KALC-R), 2014. Cancer Res Treat. 2019;51(4):1400–10. https://doi.org/10.4143/crt.2018.704.

    Article  Google Scholar 

  23. Stensrud MJ, Hernán MA. Why test for proportional hazards? JAMA. 2020;323(14):1401–2. https://doi.org/10.1001/jama.2020.1267.

    Article  Google Scholar 

  24. Bossuyt PM, Reitsma JB, Bruns DE, Gatsonis CA, Glasziou PP, Irwig L, et al. STARD 2015: an updated list of essential items for reporting diagnostic accuracy studies. BMJ. 2015;28(351): h5527.

    Article  Google Scholar 

  25. Cohen JF, Korevaar DA, Altman DG, Bruns DE, Gatsonis CA, Hooft L, et al. STARD 2015 guidelines for reporting diagnostic accuracy studies: explanation and elaboration. BMJ Open. 2016;6(11):e012799.

    Article  Google Scholar 

  26. Knibbs GH. The international classification of disease and causes of death and its revision. Med J. 1929;1:2–12.

    Google Scholar 

  27. Haendel MA, Chute CG, Robinson PN. Classification, ontology, and precision medicine. N Engl J Med. 2018;379(15):1452–62.

    Article  CAS  Google Scholar 

  28. Detterbeck FC, Chansky K, Groome P, Bolejack V, Crowley J, Shemanski L, et al. The IASLC Lung cancer staging project: methodology and validation used in the development of proposals for revision of the stage classification of NSCLC in the forthcoming (Eighth) edition of the TNM classification of lung cancer. J Thorac Oncol Off Publ Int Assoc Study Lung Cancer. 2016;11(9):1433–46.

    Google Scholar 

  29. Hernán MA, Hsu J, Healy B. A second chance to get causal inference right: a classification of data science tasks. Chance. 2019;32(1):42–9.

    Article  Google Scholar 

  30. Kelly CJ, Karthikesalingam A, Suleyman M, Corrado G, King D. Key challenges for delivering clinical impact with artificial intelligence. BMC Med. 2019;17(1):195.

    Article  Google Scholar 

  31. Ozdemir V, Williams-Jones B, Glatt SJ, Tsuang MT, Lohr JB, Reist C. Shifting emphasis from pharmacogenomics to theragnostics. Nat Biotechnol. 2006;24(8):942–6.

    Article  CAS  Google Scholar 

  32. Chen XS. Introducing theranostics journal—from the editor-in-chief. Theranostics. 2011;8(1):1–2.

    Article  Google Scholar 

  33. Smith M, Halvorson G, Kaplan G. What’s needed is a health care system that learns: recommendations from an IOM report. JAMA. 2012;308(16):1637–8.

    Article  CAS  Google Scholar 

  34. Miriovsky BJ, Shulman LN, Abernethy AP. Importance of health information technology, electronic health records, and continuously aggregating data to comparative effectiveness research and learning health care. J Clin Oncol Off J Am Soc Clin Oncol. 2012;30(34):4243–8.

    Article  Google Scholar 

  35. Bindman AB. The agency for healthcare research and quality and the development of a learning health care system. JAMA Intern Med. 2017;177(7):909–10.

    Article  Google Scholar 

  36. Sarrazin MSV, Rosenthal GE. Finding pure and simple truths with administrative data. JAMA. 2012;307(13):1433–5.

    Article  Google Scholar 

  37. Hripcsak G, Albers DJ. Next-generation phenotyping of electronic health records. J Am Med Inform Assoc JAMIA. 2013;20(1):117–21.

    Article  Google Scholar 

  38. Quan H, Parsons GA, Ghali WA. Validity of procedure codes in international classification of diseases, 9th revision, clinical modification administrative data. Med Care. 2004;42(8):801–9.

    Article  Google Scholar 

  39. Wickens MR. A note on the use of proxy variables. Econometrica. 1972;40(4):759–61.

    Article  Google Scholar 

  40. Groome PA, Bolejack V, Crowley JJ, Kennedy C, Krasnik M, Sobin LH, et al. The IASLC lung cancer staging project: validation of the proposals for revision of the T, N, and M descriptors and consequent stage groupings in the forthcoming (seventh) edition of the TNM classification of malignant tumours. J Thorac Oncol Off Publ Int Assoc Study Lung Cancer. 2007;2(8):694–705.

    Google Scholar 

  41. Cao C, Wang D, Chung C, Tian D, Rimner A, Huang J, et al. A systematic review and meta-analysis of stereotactic body radiation therapy versus surgery for patients with non-small cell lung cancer. J Thorac Cardiovasc Surg. 2019;157(1):362–73.

    Article  Google Scholar 

  42. Watanabe S-I, Nakagawa K, Suzuki K, Takamochi K, Ito H, Okami J, et al. Neoadjuvant and adjuvant therapy for stage III non-small cell lung cancer. Jpn J Clin Oncol. 2017;47(12):1112–8.

    Article  Google Scholar 

  43. MacLean M, Luo X, Wang S, Kernstine K, Gerber DE, Xie Y. Outcomes of neoadjuvant and adjuvant chemotherapy in stage 2 and 3 non-small cell lung cancer: an analysis of the national cancer database. Oncotarget. 2018;9(36):24470–9.

    Article  Google Scholar 

  44. Tao X, Yuan C, Zheng D, Ye T, Pan Y, Zhang Y, et al. Outcomes comparison between neoadjuvant chemotherapy and adjuvant chemotherapy in stage IIIA non-small cell lung cancer patients. J Thorac Dis. 2019;11(4):1443–55.

    Article  Google Scholar 

  45. Sugimura H, Yang P. Long-term survivorship in lung cancer: a review. Chest. 2006;129(4):1088–97.

    Article  Google Scholar 

  46. Noguchi M, Morikawa A, Kawasaki M, Matsuno Y, Yamada T, Hirohashi S, et al. Small adenocarcinoma of the lung. Histol Charact Progn Cancer. 1995;75(12):2844–52.

    CAS  Google Scholar 

  47. Li B, Cui Y, Diehn M, Li R. Development and validation of an individualized immune prognostic signature in early-stage nonsquamous non-small cell lung cancer. JAMA Oncol. 2017;3(11):1529–37.

    Article  Google Scholar 

  48. Samet JM, Avila-Tang E, Boffetta P, Hannan LM, Olivo-Marston S, Thun MJ, et al. Lung cancer in never smokers: clinical epidemiology and environmental risk factors. Clin Cancer Res Off J Am Assoc Cancer Res. 2009;15(18):5626–45.

    Article  Google Scholar 

  49. Alexandrov LB, Ju YS, Haase K, Van Loo P, Martincorena I, Nik-Zainal S, et al. Mutational signatures associated with tobacco smoking in human cancer. Science. 2016;354(6312):618–22.

    Article  CAS  Google Scholar 

  50. Wakelee HA, Chang ET, Gomez SL, Keegan TH, Feskanich D, Clarke CA, et al. Lung cancer incidence in never smokers. J Clin Oncol Off J Am Soc Clin Oncol. 2007;25(5):472–8.

    Article  Google Scholar 

  51. Cho J, Choi SM, Lee J, Lee C-H, Lee S-M, Kim D-W, et al. Proportion and clinical features of never-smokers with non-small cell lung cancer. Chin J Cancer. 2017;36(1):20.

    Article  Google Scholar 

  52. Takano T, Fukui T, Ohe Y, Tsuta K, Yamamoto S, Nokihara H, et al. EGFR mutations predict survival benefit from gefitinib in patients with advanced lung adenocarcinoma: a historical comparison of patients treated before and after gefitinib approval in Japan. J Clin Oncol Off J Am Soc Clin Oncol. 2008;26(34):5589–95.

    Article  CAS  Google Scholar 

  53. Yatabe Y, Kerr KM, Utomo A, Rajadurai P, Tran VK, Du X, et al. EGFR mutation testing practices within the Asia Pacific region: results of a multicenter diagnostic survey. J Thorac Oncol Off Publ Int Assoc Study Lung Cancer. 2015;10(3):438–45.

    CAS  Google Scholar 

  54. Sone S, Takashima S, Li F, Yang Z, Honda T, Maruyama Y, et al. Mass screening for lung cancer with mobile spiral computed tomography scanner. Lancet Lond Engl. 1998;351(9111):1242–5.

    Article  CAS  Google Scholar 

  55. National Lung Screening Trial Research Team, Aberle DR, Adams AM, Berg CD, Black WC, Clapp JD, et al. Reduced lung-cancer mortality with low-dose computed tomographic screening. N Engl J Med. 2011;365(5):395–409.

Download references

Acknowledgements

Part of the validation study for the index classification using electronic health records was presented as a poster at the American Medical Informatics Association 2019 Annual Symposium. The authors appreciate the reviewers for their helpful comments.

Funding

This research was supported by a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (grant number: HI18C1735). No sponsor took any role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Author information

Authors and Affiliations

Authors

Contributions

Conception and design of the work: M-HK; SP; YRP; WJ; S-GK; MC; S-SH; JCL; HRK; C-MC. Data collection: SP; YRP; WJ; S-GK; MC; C-MC. Data analysis and interpretation: M-HK; S-GK. Drafting the article: M-HK; SP. Critical revision of the article: M-HK; SP; YRP; WJ; S-GK; MC; S-SH; JCL; HRK; C-MC. Final approval of the version to be published: M-HK; SP; YRP; WJ; S-GK; MC; S-SH; JCL; HRK; C-MC. All authors read and approved the final manuscript.

Authors’ information

Chang-Min Choi (Asan Medical Center) had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Min-Hyung Kim (Harvard T.H. Chan School of Public Health), Seul-Gi Kim (Asan Medical Center), Wonjun Ji (Asan Medical Center), and Minji Choo (Asan Medical Center) were responsible for the conduct and analysis of the validation study using the EHR data. Sojung Park (Ewha Womans University Medical Center) and Seul-Gi Kim (Asan Medical Center) were responsible for the conduct and analysis of the application study using NHIS data. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Chang-Min Choi.

Ethics declarations

Ethics approval and consent to participate

The study protocol using EHR data was approved by the AMC (IRB no. AMC 2018–1515). The data were deidentified and consent requirements were waived by the Internal Review Board (IRB). The study protocol was approved by the NHI, Korea (IRB no. 2015–0236). The data were deidentified and the consent requirements were waived by the IRB.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflicts of interest related to this study.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

eTable 1. A simplified summary of the treatment guideline used in this study. eTable 2. TNM staging of non-small cell lung cancer (7th edition). eTable 3. Confusion matrix comparing the true stage and the inferred stage in the validation study (N = 1375). eTable 4. Sensitivity analysis for the confusion matrix comparing the true stage and the inferred stage in the validation study, including the lung cancer patients with cardiopulmonary comorbidities (N = 1573). eTable 5. Sensitivity analysis for the stage-specific predictive performance metrics of the index classifier against the reference standard (TNM stage), including the lung cancer patients with cardiopulmonary comorbidity (N = 1573). eTable 6. Pairwise log-rank tests corrected with Bonferroni methods. eMethods 1. Explanation of the treatment guideline used in this study. eMethods 2. Rationale for the population selection. eMethods 3. Flow Diagram for the Analytic Steps in the Validation Study for the Index Classification using the Electronic Health Records. eMethods 4. Flow Diagram for the Analytic Steps for the Stage-specific Survival Analysis Applying the Index Classification in a Population-based Administrative Database.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kim, MH., Park, S., Park, Y.R. et al. Stratifying non-small cell lung cancer patients using an inverse of the treatment decision rules: validation using electronic health records with application to an administrative database. BMC Med Inform Decis Mak 23, 3 (2023). https://doi.org/10.1186/s12911-022-02088-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12911-022-02088-x

Keywords