Skip to main content
  • Research article
  • Open access
  • Published:

Performance evaluation of case definitions of type 1 diabetes for health insurance claims data in Japan

Abstract

Background

No case definition of Type 1 diabetes (T1D) for the claims data has been proposed in Japan yet. This study aimed to evaluate the performance of candidate case definitions for T1D using Electronic health care records (EHR) and claims data in a University Hospital in Japan.

Methods

The EHR and claims data for all the visiting patients in a University Hospital were used. As the candidate case definitions for claims data, we constructed 11 definitions by combinations of International Statistical Classification of Diseases and Related Health Problems, Tenth Revision. (ICD 10) code of T1D, the claims code of insulin needles for T1D patients, basal insulin, and syringe pump for continuous subcutaneous insulin infusion (CSII). We constructed a predictive model for T1D patients using disease names, medical practices, and medications as explanatory variables. The predictive model was applied to patients of test group (validation data), and performances of candidate case definitions were evaluated.

Results

As a result of performance evaluation, the sensitivity of the confirmed disease name of T1D was 32.9 (95% CI: 28.4, 37.2), and positive predictive value (PPV) was 33.3 (95% CI: 38.0, 38.4). By using the case definition of both the confirmed diagnosis of T1D and either of the claims code of the two insulin treatment methods (i.e., syringe pump for CSII and insulin needles), PPV improved to 90.2 (95% CI: 85.2, 94.4).

Conclusions

We have established a case definition with high PPV, and the case definition can be used for precisely detecting T1D patients from claims data in Japan.

Peer Review reports

Background

Type 1 diabetes (T1D) is a chronic disease caused by the destruction of insulin-producing beta cells of the pancreas [1, 2]. T1D patients need to regularly self-monitor their plasma glucose level and self-inject insulin for all their life [2, 3], and have a higher risk of developing cardiovascular diseases than the general population [2]. Although there are some epidemiological studies on T1D patients in Japan [4, 5], studies on T1D using health insurance claims data are scarce. To assess the prevalence and clinical characteristics of patients with a disease, epidemiological studies using health insurance claims data are valuable, and these studies are common in Japan [6,7,8]. By using nationwide claims data, we can obtain the prevalence and prescription pattern of T1D patients. However, it is known that the name of the disease in the medical records is sometimes not detailed enough, as these are put for either inspection or prescription [9, 10], and relying solely on the International Statistical Classification of Diseases and Related Health Problems, Tenth Revision. (ICD 10) code is not an appropriate way of identifying a patient’s disease. Therefore, we need to decide the case definition for extracting the patient of the disease, and validation study of a case definition for extracting patients of a particular phenotype are often conducted for multiple diseases [11,12,13].

To conduct an epidemiological study on T1D patients, it is necessary to develop a case definition for T1D from the claims data. Although case definition methods for diabetes and type 2 diabetes have been proposed in some studies [14,15,16,17,18,19], proposals for T1D are scarce [20]. Although a review of the medical chart is usually conducted for deciding whether or not a given case is a true case of the target disease in the validation study, this activity is time-consuming and expensive. Moreover, when we review patients who match the case definition, only the positive predictive value (PPV) can be calculated (sensitivity of the algorithm cannot be calculated). Then, we used a newly proposed method called Phevaluator for evaluating the case definition algorithm [21]. Phevaluator is a machine-learning-based method of assessing phenotyping methods. It constructs a predictive model for the disease and calculates the predictive value of being the disease for each individual using the model. Using Phevaluator, we can calculate the performance indexes of the algorithm without reviewing the medical chart.

In this study, we aimed to construct some case definition methods of TID for claims data and evaluate the performance of them using the EHR and claims data of a University Hospital.

Methods

Study population

We used the data of a University Hospital in Japan from 2009 to 2019. Of those, the data from 2009 to 2014 were used as training data for constructing a predictive model for the disease, and the remaining were used as test data for evaluating performances of candidate case definitions. However, only patients who did not visit the hospital from 2009 to 2014 were used for test group (validation data) to separate patients in the training data and the test data. Electronic healthcare records (EHR) data and health insurance claims data were used in this study. EHR data were used for determining the “true” T1D patients in the medical chart review, as described below. From the claims data, data on age, sex, diseases, medications, and medical practices for all the visiting patients were used.

Extraction of cases of T1D from EHR data

We extracted “true” T1D patients from EHR data. First, we extracted all possible TID patients from visiting patients from 2009 to 2014. All the “possible” patients are those who met one or more of the following five criteria.

  1. Patients who were diagnosed with TID or insulin-dependent diabetes.

  2. Patients who met all of the criteria (a), (b), and (c).

    1. (a)

      Those who were prescribed insulin treatment.

    2. (b)

      Those who had serum C-peptide immunoreactivity (CPR) less than 0.6 ng/ml at least once.

    3. (c)

      Those who had earlier been diagnosed with ketoacidosis.

  3. Patients whose insulin autoantibody (anti-glutamic acid decarboxylase antibody; GAD or anti-insulinoma‐associated protein-2 antibody; IA2) was positive.

  4. Patients who were introduced as definitely the patients with TID by diabetologists.

  5. Patients whose serum CPR values were less than 0.2 ng/ml at least once.

Medical chart review was conducted against all the possible TID patients by three diabetologists, and each patient was either classified as “true” TID patients or not. Then, those who fell within the criteria were classified as TID patients.

Phevaluator

Phevaluator is a machine-learning-based method of evaluating case definitions [21]. In this method, a predictive model is constructed in the training data for classifying the target disease patients and the other population. Then, we apply the predictive model to the patients of test group (validation data), and calculate the predictive probability of being a patient with the disease for each patient. We explain Phevaluator using the example where the number of patients in the validation data are 4 in Table 1. Let patients A and C be tested positive according to a candidate case definition, and patients B and D be tested negative according to the same case definition. From the predictive values of being a true patient or not for A and C, we calculate the cumulative probabilities of true positivity (TPs) and cumulative probabilities of false positivity (FPs). Similarly, from the predictive values of being a true patient or not for B and D, we calculate the cumulative probabilities of true negativity (TNs) and the cumulative probabilities of false negativity (FNs). Then, sensitivity, specificity, PPV, and negative predictive value (NPV) can be calculated as follows: sensitivity: \(TPs/(TPs+FNs)\); specificity: \(TNs/(FPs+TNs)\); PPV: \(TPs/(TPs+FPs)\); NPV: \(TNs/(FNs+TNs)\). We also calculated the F-score by \(2\times sensitivity\times PPV/(sensitivity+PPV)\).

Table 1 Schematic table of the calculation method of performance indexes by Phevaluator

Case definitions

As possible case definitions, we evaluated the performances of multiple definitions. We constructed definitions using the following four codes: Confirmed disease name of TID, ICD 10 code: of E10; Claims code of insulin needles for T1D patients: claims code of 114,010,970; Claims code of basal insulin (long-acting insulin analog, intermediate-acting insulin, biphasic insulin); and Claims code of syringe pump for continuous subcutaneous insulin infusion (CSII): claims codes of 114004810 and 114022010. The codes are closely associated with T1D. As concerns the method of glycemic control, there are two types of methods: multiple insulin injections and CSII [22]; we used the claims codes of both types of methods. Also, in Japan, medical remuneration points of prescription of injector needles for TID patients are set as same as those for hemophilia or other patients and are higher than those for other diseases. TID patients who self-inject insulin must have this claims code, and we used it. Basal insulin consists of a long-acting insulin analog, intermediate-acting insulin, and biphasic insulin. The list of claims codes for basal inulin is shown in the Additional file 1. As the codes associated with CSII, we used the claims code for the syringe pump for intermittent infusion. Through combinations of these four items, we tested the performance of 11 types of case definitions.

Statistical analysis

We constructed the predictive model for classifying TID patients with non-TID patients using the claims data from 2009 to 2014, and the data on age, sex, and diseases, medications, and medical practices were used as explanatory variables. Regarding diseases, we distinguished between suspected diseases and confirmed disease and classified each of the ICD10 codes of diseases based on the first three digits of the ICD10 codes. For medications, we classified claims codes for medications by their generic names.

Regarding the outcome variable, the patients who were definitely patients of TID should be used as the cases in the evaluation by Phevaluator [21]. Therefore, from the patients who were classified as TID, we excluded patients whose pancreas were transplanted because their insulin secretion ability was probably boosted by the transplantation. We also excluded the patients who had either confirmed or suspected type 2 diabetes (T2D). Also, as concerns the controls, we needed to extract those who were definitely not patients with TID [21]. Then, the controls were selected from the patients whose medical charts were not reviewed. Also, those who had either suspected or confirmed TID were excluded. Furthermore, we randomly sampled the control patients to adjust the ratio of the cases and the controls [21].

Explanatory variables, except for age and sex, were transformed into dummy variables based on whether a patient had a code during the periods or not. However, the variables used for extracting the true patients need to be excluded from the explanatory variables [21]. Therefore, we excluded variables of confirmed disease names of T1D, T2D, and ketoacidosis from the explanatory variables. Also, we excluded variables of suspected disease names related to T1D and T2D. Moreover, claims code for test for the insulin receptor autoantibody and CPR, and insulin medications were excluded from the explanatory variables.

The gradient-boosting decision tree was used in constructing the predictive model. However, if we use all the explanatory variables, the size of the data would become very large. Therefore, we calculated in advance the relative risk of each explanatory variable to the outcome variable and used the top 500 variables for the predictive model. The area under the curve (AUC) of the predictive model was calculated by tenfold cross-validation within the training data. A predictive model of the test data was applied to each patient, and the predictive value for T1D was calculated. Also, confidence intervals of the performance indexes were calculated by bootstrap sampling. All statistical analyses were conducted using R version 3.6.3 (https://cran.r-project.org/bin/windows/base/old/3.6.3/).

Results

Figure 1 shows the flowchart of the training data in this study. Finally, 296 patients were used as cases, and 69023 patients were used as controls.

Fig. 1
figure 1

Flowchart of subjects used for the training data

Table 2 shows the baseline characteristics of the analyzed data. More than 7000 variables in total were used in the construction of the predictive model.

Table 2 Baseline characteristics of the analyzed data

Table 3 shows the variable importance of gradient boosting decision trees, and the result of the top 10 variables are shown. The variable importance of the claims code of insulin needles for T1D patients was significantly higher than other variables, and it was suggested that the claims code is crucial for the classification of T1D patients. Also, the mean AUC of the predictive model by 10-cross validation was 0.935. Therefore, a predictive model with high classification ability was constructed using the analyzed data.

Table 3 Variable importance of gradient boosting decision trees (top 10 variables)

Table 4 shows the result of the evaluation of candidate case definitions for T1D. The results of NPV and specificity are not shown in Table 4 because the values were almost 100% for all the case definitions. The sensitivity and PPV of confirmed disease T1D were relatively low. On the other hand, sensitivity and PPV for the code of insulin needles for T1D patients outperformed those of the confirmed disease, and the F-score was the highest among the candidate definitions. By combining the confirmed T1D and the claims codes of the syringe pump for CSII, PPV became the highest, but the sensitivity dropped. The 9th case definition had the highest F-score among the case definitions whose PPV was approximately 90%.

Table 4 The result of the evaluation of candidate case definitions for type 1 diabetes

Discussion

We used Phevaluator for evaluating the performance of the case definitions in order to extract true TID patients from the medical records. The predictive model with high classification ability was constructed, and the performance indexes were considered to be accurately estimated. As a result of variable importance showed, the claims code of insulin needles for T1D patients was found to have the highest predictive ability. T1D patients need to inject insulin regularly [2, 3], and injection needles are considered to be vital for classifying T1D patients and non-T1D patients.

As a result of performance evaluations, PPV and sensitivity varied depending on the case definitions. Both sensitivity and PPV for the case definition of only using the ICD 10 code of T1D were not high among the case definitions for extracting the true T1D patients. It was suggested that more than half of the number of patients with T1D patients were not attached to the ICD 10 code of T1D, and there is a possibility that many T1D patients treated with other disease names except the ICD 10 code of E10. It was also demonstrated in a previous study that the sensitivity of 1 diabetes was relatively low for extracting “true” T1D patients [20]. On the other hand, the PPV increased by using the claims code of insulin needles for T1D patients, and the sensitivity was almost unchanged. It is considered that a large part of patients with confirmed diagnosis of T1D use insulin needles for T1D patients because decrease in the sensitivity was small. Furthermore, by using both the confirmed diagnosis of T1D and the claims code of insulin needles as the case definition, PPV increased further. As one possible reason, the claims code of insulin needles is used not only for T1D patients, but the same claims code is used for injector needles for hemophilia patients. Therefore, by restricting the cases to those who have ICD10 codes of T1D, PPV is considered to have improved. On the other hand, PPV did not improve by adding the code of basal insulin for the case definition. Basal insulins are required for T1D patients who cannot secrete insulin [23], and the sensitivity was highest when using only basal insulin. As the result shows, the number of extracted cases by the case definition remained almost unchanged by adding the codes of basal insulin compared with the definition of only using the insulin needles, and it is considered that patients who have the code of insulin needles tend to have a claims code of basal insulin too. Although PPV increased when using the claims codes of the syringe pump for CSII, sensitivity tended to decrease because the rate of T1D patients using CSII is lower compared with injection treatment [24]. However, the F-score of the confirmed diagnosis of TID with the claims code of either of the treatment methods has the highest F-score among the case definitions whose PPV were sufficiently large. Taking into account that a certain percentage of T1D patients use CSII in Japan [24], this case definition is considered to be useful in identifying true T1D patients.

This study has some limitations. First, we used the data of only one site for evaluating the case definitions. It is considered that some patients receive insulin medications in other clinics or hospital, and this might have affected the result. Similar validation studies need to be conducted using nationwide data for confirming the results of this study. As another limitation, although whether a candidate T1D patient is a “true” case or control was judged by diabetologists in the medical chart review for the classification of TID patients, there is still a possibility that some misdiagnosis occurred in the classification. Finally, we could not obtain a case definition whose PPV and sensitivity were both sufficiently high values, and we still need to seek claims codes for improving the sensitivity of the case definition. However, the PPV was high enough, and it was suggested that we could precisely identify true T1D patients from claims data. The actual conditions of prescription patterns, comorbidities, or medical expenditures for T1D patients is uncertain in Japan at this moment, and epidemiological studies on T1D using claims data need to be conducted by using the proposed case definition. This is the first study that derived a case definition of T1D from claims data in Japan, and further studies for case definition and epidemiological studies of T1D are needed.

Conclusions

As a result of the performance evaluation of the case definitions for T1D, it was suggested that the ICD10 code of T1D should not be used for assessing the true patients with T1D. The F-score was highest when using both the confirmed diagnosis of T1D and either of the claims codes of two insulin treatment methods (i.e., syringe pump for CSII and insulin needles) among the case definitions whose PPV were sufficiently large. Therefore, the proposed case definition can be used for precisely detecting T1D patients from claims data in Japan.

Availability of data and materials

The datasets analyzed during the current study are not publicly available due hospital data were used but are available from the corresponding author on reasonable request.

Abbreviations

ICD10:

International Statistical Classification of Diseases and Related Health Problems, Tenth Revision

T1D:

Type 1 diabetes

EHR:

Electronic health care records

PPV:

Positive predictive value

CSII:

Continuous subcutaneous insulin infusion

CPR:

C-peptide immunoreactivity

AUC:

Area under the curve

T2D:

Type 2 diabetes

References

  1. Butalia S, Kaplan GG, Khokhar B, Rabi DM. Environmental risk factors and type 1 diabetes: past, present, and future. Can J Diabetes. 2016;40(6):586–93.

    Article  Google Scholar 

  2. Atkinson MA, Eisenbarth GS, Michels AW. Type 1 diabetes. Lancet. 2014;383(9911):69–82.

    Article  Google Scholar 

  3. Janež A, Guja C, Mitrakou A, et al. Insulin therapy in adults with type 1 diabetes mellitus: a narrative review. Diabetes Ther. 2020;11(2):387–409.

    Article  Google Scholar 

  4. Onda Y, Sugihara S, Ogata T, et al. Incidence and prevalence of childhood-onset Type 1 diabetes in Japan: the T1D study. Diabet Med. 2017;34(7):909–15.

    Article  CAS  Google Scholar 

  5. Kawasaki E, Matsuura N, Eguchi K. Type 1 diabetes in Japan. Diabetologia. 2006;49(5):828–36.

    Article  CAS  Google Scholar 

  6. Mahlich J, Tsukazawa S, Wiegand F. Estimating prevalence and healthcare utilization for treatment-resistant depression in japan: a retrospective claims database study. Drugs Real World Outcomes. 2018;5(1):35–43.

    Article  Google Scholar 

  7. Nishimura R, Kato H, Kisanuki K, et al. Treatment patterns, persistence and adherence rates in patients with type 2 diabetes mellitus in Japan: a claims-based cohort study. BMJ Open. 2019;9(3):e025806.

    Article  Google Scholar 

  8. Sato K, Ohno T, Ishii T, Ito C, Kaise T. The prevalence, characteristics, and patient burden of severe asthma determined by using a Japan health care claims database. Clin Ther. 2019;41(11):2239–51.

    Article  Google Scholar 

  9. Khan A, Ramsey K, Ballard C, et al. Limited accuracy of administrative data for the identification and classification of adult congenital heart disease. J Am Heart Assoc. 2018;7(2):e007378.

    Article  Google Scholar 

  10. Oake J, Aref-Eshghi E, Godwin M, et al. Using electronic medical record to identify patients with dyslipidemia in primary care settings: international classification of disease code matters from one region to a national database. Biomed Inform Insights. 2017;9:1178222616685880.

    Article  Google Scholar 

  11. van Mourik MS, van Duijn PJ, Moons KG, Bonten MJ, Lee GM. Accuracy of administrative data for surveillance of healthcare-associated infections: a systematic review. BMJ Open. 2015;5(8):e008424.

    Article  Google Scholar 

  12. Lee CK, Ha HJ, Oh SJ, et al. Nationwide validation study of diagnostic algorithms for inflammatory bowel disease in Korean National Health Insurance Service database. J Gastroenterol Hepatol. 2020;35(5):760–8.

    Article  Google Scholar 

  13. Quan H, Khan N, Hemmelgarn BR, et al. Validation of a case definition to define hypertension using administrative data. Hypertension. 2009;54(6):1423–8.

    Article  CAS  Google Scholar 

  14. Richesson RL, Rusincovitch SA, Wixted D, et al. A comparison of phenotype definitions for diabetes mellitus. J Am Med Inform Assoc. 2013;20(e2):e319–26.

    Article  Google Scholar 

  15. Esteban S, Rodríguez Tablado M, Peper FE, et al. Development and validation of various phenotyping algorithms for Diabetes Mellitus using data from electronic health records. Comput Methods Programs Biomed. 2017;152:53–70.

    Article  Google Scholar 

  16. Chen G, Khan N, Walker R, Quan H. Validating ICD coding algorithms for diabetes mellitus from administrative data. Diabetes Res Clin Pract. 2010;89(2):189–95.

    Article  Google Scholar 

  17. Khokhar B, Jette N, Metcalfe A, et al. Systematic review of validated case definitions for diabetes in ICD-9-coded and ICD-10-coded data in adult populations. BMJ Open. 2016;6(8):e009952.

    Article  Google Scholar 

  18. Zheng T, Xie W, Xu L, et al. A machine learning-based framework to identify type 2 diabetes through electronic health records. Int J Med Inform. 2017;97:120–7.

    Article  Google Scholar 

  19. Kagawa R, Kawazoe Y, Ida Y, et al. Development of type 2 diabetes mellitus phenotyping framework using expert knowledge and machine learning approach. J Diabetes Sci Technol. 2017;11(4):791–9.

    Article  Google Scholar 

  20. Klompas M, Eggleston E, McVetta J, Lazarus R, Li L, Platt R. Automated detection and classification of type 1 versus type 2 diabetes using electronic health record data. Diabetes Care. 2013;36(4):914–21.

    Article  Google Scholar 

  21. Swerdel JN, Hripcsak G, Ryan PB. PheValuator: development and evaluation of a phenotype algorithm evaluator. J Biomed Inform. 2019;97:103258.

    Article  Google Scholar 

  22. Ross LJ, Neville KA. Continuous subcutaneous insulin infusion versus multiple daily injections for type 1 diabetes. J Paediatr Child Health. 2019;55(6):718–22.

    Article  Google Scholar 

  23. Matejko B, Kukułka A, Kieć-Wilk B, Stąpór A, Klupa T, Malecki MT. Basal insulin dose in adults with type 1 diabetes mellitus on insulin pumps in real-life clinical practice: a single-center experience. Adv Med. 2018;2018:1473160.

    Article  Google Scholar 

  24. Murata T, Aoki Y, Kato Y, et al. The percentage of continuous subcutaneous insulin infusion usage among adult type 1 diabetes mellitus patients in Japan: a cross-sectional study at national hospital organization hospitals. J Diabetes Sci Technol. 2017;11(5):1055–6.

    Article  Google Scholar 

Download references

Acknowledgements

Enago has proofread the manuscript.

Funding

This study was supported by grants from the Ministry of Health, Labour and Welfare (Nos. H29-Juinkankitou-Ippan-004 and H28-Jyunkankitou-Ippan-006).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: NN, NT. Data curation: CN, SK, KA, SM, MM, YM, NN. Formal analysis: CN, TO. Methodology: CN, NN, TO. Funding acquisition: NN, NT. Writing—original draft: TO. Writing—review and editing: NN, NT, SK, KA, SM, MM, YM, CN, TK, TO. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Tasuku Okui.

Ethics declarations

Ethics approval and consent to participate

This study was performed in accordance with the Declaration of Helsinki and was approved by the ethical committees of Faculty of Medicine of the Kyushu University (Number 2019–587). Also, we used existing and de-identified clinical data of a hospital, and obtaining informed consents from all participants is difficult. In this case, informed consent from all the participants are not essential according to Ethical Guidelines for Medical and Health Research Involving Human Subjects in Japan (https://www.mhlw.go.jp/file/06-Seisakujouhou-12600000-Seisakutoukatsukan/0000168764.pdf). In addition, need of informed consent from participants was waived by the ethics committee of Faculty of Medicine of the Kyushu University. Instead, we made an opportunity for declining participation in this study by disclosing information about this study in our homepage (https://www.mic.hosp.kyushu-u.ac.jp/files/2020/05/8adc802e9d2b750023a9dda50f8e4845.pdf).

Competing interests

The authors declare that they have no competing interest.

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.

List of basal insulin used in the analysis.

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

Okui, T., Nojiri, C., Kimura, S. et al. Performance evaluation of case definitions of type 1 diabetes for health insurance claims data in Japan. BMC Med Inform Decis Mak 21, 52 (2021). https://doi.org/10.1186/s12911-021-01422-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12911-021-01422-z

Keywords