Impact of grouping complications on mortality in traumatic brain injury: A nationwide population-based study

Abstract

Traumatic brain injury (TBI) is an important health issue with high mortality. Various complications of physiological and cognitive impairment may result in disability or death after TBI. Grouping of these complications could be treated as integrated post-TBI syndromes. To improve risk estimation, grouping TBI complications should be investigated, to better predict TBI mortality. This study aimed to estimate mortality risk based on grouping of complications among TBI patients. Taiwan's National Health Insurance Research Database was used in this study. TBI was defined according to the International Classification of Diseases, Ninth Revision, Clinical Modification codes: 801–804 and 850–854. The association rule data mining method was used to analyze coexisting complications after TBI. The mortality risk of post-TBI complication sets with the potential risk factors was estimated using Cox regression. A total 139,254 TBI patients were enrolled in this study. Intracerebral hemorrhage was the most common complication among TBI patients. After frequent item set mining, the most common post-TBI grouping of complications comprised pneumonia caused by acute respiratory failure (ARF) and urinary tract infection, with mortality risk 1.55 (95% C.I.: 1.51–1.60), compared with those without the selected combinations. TBI patients with the combined combinations have high mortality risk, especially those aged <20 years with septicemia, pneumonia, and ARF (HR: 4.95, 95% C.I.: 3.55–6.88). We used post-TBI complication sets to estimate mortality risk among TBI patients. According to the combinations determined by mining, especially the combination of septicemia with pneumonia and ARF, TBI patients have a 1.73-fold increased mortality risk, after controlling for potential demographic and clinical confounders. TBI patients aged<20 years with each combination of complications also have increased mortality risk. These results could provide physicians and caregivers with important information to increase their awareness about sequences of clinical syndromes among TBI patients, to prevent possible deaths among these patients.

About the Authors:

Chung-Han Ho

Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Writing – original draft

Affiliations Department of Medical Research, Chi-Mei Medical Center, Tainan, Taiwan, Department of Hospital and Health Care Administration, Chia Nan University of Pharmacy and Science, Tainan, Taiwan

ORCID logo http://orcid.org/0000-0001-5925-8477

Fu-Wen Liang

Roles Data curation, Formal analysis, Writing – review & editing

Affiliation: National Cheng Kung University Research Center for Health Data and Department of Public Health, College of Medicine, National Cheng Kung University, Tainan, Taiwan

Jhi-Joung Wang

Roles Project administration, Resources

Affiliation: Department of Medical Research, Chi-Mei Medical Center, Tainan, Taiwan

Chung-Ching Chio

Roles Resources, Software

Affiliation: Department of Neurosurgery, Chi-Mei Medical Center, Tainan, Taiwan

Jinn-Rung Kuo

Roles Conceptualization, Project administration, Writing – original draft

* E-mail: [email protected]

Affiliations Department of Medical Research, Chi-Mei Medical Center, Tainan, Taiwan, Department of Neurosurgery, Chi-Mei Medical Center, Tainan, Taiwan, Department of Biotechnology, Southern Taiwan University of Science and Technology, Tainan, Taiwan

Introduction

Traumatic brain injury (TBI), a major cause of morbidity and mortality in many countries [1,2], is an important global public health problem. TBI is known as a silent disease because its complications are often associated with changes in thinking, feeling, language, or mood [3,4].

Patients with TBI commonly have many sequelae, such as sleep disorders [5], anxiety or depression [6], or post-traumatic stress disorder [7]. Studies have shown that these complications may lead to illness or death [8,9]; however, patients with TBI often pay little attention to these complications, as many of them are not obvious to patients or caregivers. In addition, many studies have suggested that TBI may cause or accelerate the progression of different diseases [9], such as stroke [10,11] and pneumonia [12].

The mortality risk owing to a single sequela after TBI has been analyzed in many studies, as mentioned above, and many researchers have developed different risk prediction models to estimate the probability of death in post-TBI patients [13–15]. However, a major challenge in evaluating the effect of multiple complications among patients with TBI is how to adequately and appropriately present complex postoperative complications among TBI patients. In particular, some complications after TBI are directly caused by brain injury and some are indirectly caused by other complications [9]. For example, stroke is a risk factor of death in all populations [16], and some studies have indicated that TBI is associated with a higher risk of stroke [10,11]. In addition, it has been shown that sleep problems are a common symptom during and after TBI [9,17,18] and that patients with sleep disorders have an increased risk of stroke. Yet it remains unknown how the combined effect of stroke and sleep disorders affects the risk of death among TBI patients. Therefore, diseases and symptoms should be grouped together as a more meaningful risk factor to examine the combined effect of disease and disease, symptoms and symptoms, or disease and symptoms on mortality risk in patients with TBI. The relationships between diseases and symptoms in large databases could be grouped using an association rule approach, to discover potential patterns [19].

To the best of our knowledge, few studies have attempted to explore the effect of grouping complications among post-TBI patients. Therefore, this study aimed not only to explore the associations among post-TBI complications but also to estimate mortality risk by grouping TBI complications in post-TBI patients, using a nationwide administrative database.

Methods

Data source

Taiwan launched a government-run single-payer National Health Insurance program in 1995. The medical claims database, the National Health Insurance Research Database (NHIRD), was established for research purposes. In this study, data were obtained from the NHIRD between January 1997 and December 2013. The NHIRD covers 99% of inpatient and outpatient medical benefit claims for Taiwan’s 23 million residents. The database comprises detailed information regarding clinical visits for each insured person, including date of diagnosis; diagnostic codes according to the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM); payments for consultations; and prescription details [20,21].

Ethics statement

Encrypted personal identifiers are used in the NHIRD to protect patients’ privacy and prevent the possibility of an ethical violation according to regulations of the Bureau of National Health Insurance in Taiwan. The anonymous identification numbers have linked claims information, such as patient sex, date of birth, medical services received, and prescription information, which can be used by researchers. Therefore, informed consent was not required for the present study, which was approved for exemption by the Institutional Review Board (IRB) of Chi Mei Medical Center (IRB: 10410-E01).

Patient selection and definition

Diagnostic codes were accessed via the inpatient and outpatient claims databases of the NHIRD. A total of 139,254 patients with TBI were enrolled in this study. TBI patients were defined as inpatients with a diagnosis code of TBI (ICD-9-CM codes 801–804 and 850–854), and who had intensive care unit (ICU) records between 2001 and 2012. We excluded patients who died within one month after TBI diagnosis owing to severe TBI. In addition, to ensure that all selected patients had new-onset TBI, we also excluded those with a TBI diagnosis before December 31, 2000.

Demographic variables (age and sex) and clinical characteristics, including length of stay (LOS), length of ICU stay, length of noninvasive ventilation (NV), and Charlson comorbidity index (CCI) score, were used as potential confounders to estimate the mortality risk among post-TBI patients. The CCI score is commonly used as a surrogate of patients’ disease severity [22,23]. Comorbidities measured by the CCI score are based on medical records for the one year prior to the date of TBI diagnosis; the CCI score has been used in many studies of TBI [24,25]. The length of ICU stay in TBI patients is often associated with a severe or poor clinical condition [26,27]. Therefore, the length of hospital stay, length of ICU stay, and the length of NV days among patients with TBI were further classified into four groups, according to their statistical distribution: a short-term group (minimum to first quartile), mild-term group (first quartile to median), long-term group (median to third quartile), and prolonged group (third quartile to maximum). This approach has been used in previous studies [28].

Measurements

The main event in this study was mortality in post-TBI patients with more than 30 days’ survival. Follow-up time was calculated from the date of TBI admission to the date of death, withdrawal from the insurance system, or the end of the study (31 December, 2013), whichever came first. The last two conditions were defined as censors in this study. Mortality was defined as having a death record in the inpatient claim dataset; a ICD-9-CM diagnosis code of 798 (sudden death, cause unknown) in the emergency department, according to the outpatient claim dataset; or withdrawal from the insurance program and without reenrollment within 180 days of the last health care visit.

The main aim of this study was to examine the grouping of complications in post-TBI patients using the association rule method. Complications were based on ICD-9-CM diagnosis codes after the date of TBI admission, obtained from inpatient and outpatient claims datasets. However, only complications that occurred within the first year after the date of TBI admission were selected to estimate the mortality risk among post-TBI patients. The primary reason for selecting only complications within the first year was to avoid inclusion of diseases that may be unrelated to TBI, such as the common cold. In addition, we excluded patients with certain chronic diseases before TBI admission, such as hypertension (ICD-9-CM: 410), diabetes mellitus (DM) (ICD-9-CM: 250), and hyperlipidemia (ICD-9-CM: 272).

The association rule, a data mining approach, was used to evaluate multiple complications, to improve the accuracy of mortality risk prediction. The association rule method was first introduced by Piatetsky-Shaprio [19], and it is often used in large databases for pattern mining in different fields, such as computer science [29] and for medical applications [30,31]. In this study, frequent item set mining, an association rule method, was used to investigate the grouping of post-TBI complications for estimating mortality risk. The support (%) of complications set X was defined as the proportion of diagnoses among all health care visits after the date of TBI admission from inpatient and outpatient datasets. Support is a measure of how frequently the rule occurs in the database. The resultant integrated post-TBI risk profile will provide useful information for TBI patients and their families.

Statistical analysis

The frequency with percentage and median with interquartile range (IQR) were presented for discrete variables and continuous variables, respectively. The Pearson’s chi-squared test was used to examine distribution differences by age group, sex, comorbidities, categorized CCI score, LOS, length of ICU stay, and length of NV days between patients with TBI who were alive and those who died. The Wilcoxon rank-sum test was used to compare LOS, length of ICU stay, and length of NV days at first TBI diagnosis between deceased and alive TBI patients.

In addition, the Kaplan-Meier failure plot was applied to describe the cumulative incidence rate of mortality, and the log-rank test was used to compare the risk difference between patients with grouping of complications and those without grouping. A Cox proportional regression model was used to estimate the mortality risk for different complication groupings, adjusted for potential confounding factors including age, sex, CCI score, LOS, length of ICU stay, and length of NV days. The proportional hazards assumption was assessed using Schoenfeld residuals for each variable of interest.

To test the hypothesis that the effect of grouping post-TBI complications on mortality varied by age and sex in patients with TBI, the model was also stratified by age and sex groups. SAS version 9.4 software (SAS Institute, Inc, Cary, NC, USA) was used to perform all statistical analyses. Kaplan-Meier curves were generated using Stata version 12 (StataCorp, College Station, TX, USA). All significance levels were set at a p-value <0.05.

Results

A total of 139,254 patients were enrolled in this study. The distribution of the baseline characteristics between TBI patients who died before the end of the study and those who were still alive is presented in Table 1. The median follow-up time for the enrolled participants was 4.39 years (IQR = 1.76–7.63 years). The distribution of all characteristics, including age, sex, and comorbid conditions, was significantly different between the deceased group and the survival group. Compared with TBI patients who were still alive, those patients who died were significantly older and had higher CCI scores, longer lengths of stay, longer lengths of ICU stay, and a longer length of NV days. The median LOS for deceased TBI patients was 16 days (IQR = 8.0–30.0 days) and that for the survival group was 10 days (IQR = 6.0–17.0 days). The median length of ICU stay for TBI patients who died and those who survived was 5 days (IQR = 2.0–15.0 days) and 3 days (IQR = 1.0–5.0 days), respectively. Furthermore, as predicted, TBI patients who died before the end of study had more comorbid conditions. Among deceased TBI patients, 17.7% had cerebrovascular disease, 13.2% had chronic pulmonary disease, 12.3% had ulcer disease, and 23.8% had DM. Among TBI patients who survived, 5.0% of them had cerebrovascular disease, 4.2% had chronic pulmonary disease, 5.1% had ulcer disease, and 9.4% had DM. In addition, TBI patients aged over 65 years had the highest risk of death (HR = 27.64, 95% CI = 24.64–29.71). Male TBI patients had 1.14 times (95% CI = 1.11–1.17) greater risk for death than their female counterparts. Higher CCI scores were associated with higher risk of death; patients with a baseline CCI score of 2 had 4.59 times (95% CI = 4.42–4.78) greater risk of death and those with a baseline score of 3 had 6.69 times (95% CI = 6.49–6.70) greater risk, compared with patients with a baseline score of 0. For each one-day increase in length of ICU stay, there was a 5% (IQR = 2.00–15.00) increase in the mortality risk. Patients with one additional day of ventilator use had a 4% (95% CI = 1.043–1.045) increased risk of death.

[Figure omitted. See PDF.]

Table 1. Baseline information of 139,254 patients with TBI from 2001–2012.

https://doi.org/10.1371/journal.pone.0190683.t001

The five most common complications occurring within one year of TBI diagnosis are shown in Table 2. A total 3,853,713 diagnostic records of post-TBI patients within one year were included in this study. Intracerebral hemorrhage (2.66%) was the most common complication within the first year after TBI, followed by acute upper respiratory infections of multiple or unspecified sites (1.75%), dizziness and giddiness (1.21%), constipation (1.03%), and urinary tract infection (0.92%).

[Figure omitted. See PDF.]

Table 2. Top five complications within one year among post-TBI patients.

https://doi.org/10.1371/journal.pone.0190683.t002

Most patients included in this study had more than one complication. Using the frequent item set mining approach, we found five sets of the most common post-TBI grouping combinations that were highly correlated (Table 3). The most common set was pneumonia with acute respiratory failure (ARF) and urinary tract infection (UTI) (N = 16,859, 0.055%), followed by the set of unspecified septicemia with pneumonia and ARF (N = 12,988, 0.042%), the set of unspecified septicemia with ARF and UTI (N = 9,778, 0.032%), the set of unspecified septicemia with pneumonia and UTI (N = 8,423, 0.027%), and the combination of pneumonia with chronic respiratory failure and UTI (N = 6,535, 0.021%).

[Figure omitted. See PDF.]

Table 3. Top five grouped complications within the first year after TBI.

https://doi.org/10.1371/journal.pone.0190683.t003

Kaplan-Meier failure estimates of mortality between TBI patients with and without grouping of post-TBI complications are shown in Fig 1. Consistently, the survival probability for TBI patients with any set of grouped complications was significantly lower than those without grouped complications (log-rank test, p < .0001) over the follow-up period.

[Figure omitted. See PDF.]

Fig 1. Kaplan-Meier failure estimates of mortality between TBI patients with and without grouped post-TBI complications.

https://doi.org/10.1371/journal.pone.0190683.g001

In the Cox regression model (Table 4), each combination of post-TBI complications was a significant predictor of mortality. After adjusting for age, sex, CCI score, LOS, length of ICU stay, and length of ventilation days, most sets of grouped complications were significantly associated with increased risk of death, except for the set comprising pneumonia, chronic respiratory failure, and UTI.

[Figure omitted. See PDF.]

Table 4. Hazard ratios of mortality among TBI patients with grouping of post-TBI complications, compared with those without grouping.

https://doi.org/10.1371/journal.pone.0190683.t004

The results of stratified analyses by age and sex are presented in Table 5. For each age group, after controlling for potential confounders, the grouping of complications was associated with increased risk of death, and the effect was stronger in younger patients with Combinations 2, 3, and 4. However, elderly TBI patients (aged 65 years or more) with Combination 5 had a significantly lower risk of death compared with those of the same age who did not have that set of complications (HR = 0.92, 95% CI = 0.87–0.98). In addition, grouping of complications increased the risk of death in both men and women, and the adverse effect was stronger in female TBI patients.

[Figure omitted. See PDF.]

Table 5. Hazard ratios of mortality among TBI patients, with and without combinations of three complications, stratified by age and sex.

https://doi.org/10.1371/journal.pone.0190683.t005

Discussion

To the best of our knowledge, this is the first study to report the risk of mortality based on grouping of post-TBI complications among TBI patients using population-based administrative data. The main findings of this study were as follows. (1) The five most common post-TBI complications were related to a clinical diagnosis of TBI; these were intracerebral hemorrhage, acute upper respiratory infections, dizziness and giddiness, constipation, and urinary tract infection. (2) The main combinations of post-TBI complications were commonly found in respiratory care centers or respiratory care wards; these included pneumonia, acute respiratory failure, urinary tract infection, unspecified septicemia, and chronic respiratory failure. (3) The adverse effects were greater when the grouping of TBI complications included unspecified septicemia. (4) The effect of increased mortality risk owing to grouping of post-TBI complications was greater in younger and in female TBI patients.

Consistent with our previous studies, older age, male sex, higher CCI scores, longer length of hospital stay, longer ICU stay, longer ventilator use, and comorbidities of DM or cerebrovascular, cardiac, renal, or liver diseases were associated with significantly higher mortality (Table 1) [32,33]. Our results showed that intracerebral hemorrhage (2.66%), a type of cerebrovascular accident, was the most common diagnosis among post-TBI patients (Table 2). Intracerebral hemorrhage is one of the diagnosis codes used in the definition of stroke in previous studies [34]. Stroke is the second leading cause of death worldwide, and it is also the main cause of disability in adults [35]. In addition, most studies have indicated that TBI patients have a higher risk of stroke [11,36]. Hence, the results of the present study provide powerful evidence for the high incidence of intracerebral hemorrhage among post-TBI patients.

Recently, co-existence of two or more chronic diseases in an individual, also known as multiple morbidity, complex morbidity, or complex chronic disease, has been proved to have a much greater impact on health care systems than previously thought [37], such as by reducing life quality and life expectancy [38,39] and increasing drug–drug interaction side effects [40]. However, when discussing the effects of complications on mortality, potential prognostic factors are usually evaluated individually. The mortality risk based on grouping of post-TBI complications has not been investigated. In particular, most studies of TBI have used the neural network approach, a machine learning method, to predict potential risk factors [41–43]; however, grouping complications using association rule mining may offer a novel approach to evaluating the combined effects of risks in the real world. Therefore, this issue is worth further investigation, as our results have shown.

Our findings indicated that the most common set of post-TBI complications was pneumonia with acute respiratory failure (ARF) and urinary tract infection (UTI) (0.055%); other combinations were also highly related to pneumonia. Patients with TBI have been shown to have a greater risk of pneumonia [44,45]. In addition, TBI patients might have airway obstruction, aspiration, or hypoxia [46], which could lead to pneumonia, as could using intubation or mechanical ventilation [47]. Our findings when grouping complications were consistent with those other studies, as pneumonia, ARF, and UTI are well-known post-TBI complications [48].

Septicemia, a common post-TBI complication, might be associated with other neurological outcomes, according to different levels of disease severity [49]. Studies with rat models have indicated that sepsis would increase the post-TBI mortality rate [50,51]. In our study, Combinations 2 and 3, which include septicemia, presented higher mortality risk (Combination 2: HR = 1.73, 95% CI = 1.68–1.78; Combination 3: HR = 1.53, 95% CI = 1.48–1.58) than other groupings. Therefore, it is important that intensive chest care, a timely fever workup, and adequate antibiotic treatment are provided for patients with TBI, to prevent septicemia.

Moreover, the risk of death among TBI patients according to age and sex varied, as the prognosis differed according to these factors [52,53]. Age is directly associated with post-TBI mortality [54,55]. Thus, a stratified analysis by different ages and sex was necessary. In our study, the effect of increased mortality risk was stronger among younger TBI patients. A possible reason for this effect is that younger adults might receive lower levels of critical care consultation, which may lead to higher mortality in this population [56]. This implies that TBI patients who are younger or female and who have the identified grouping combinations are more vulnerable (Table 5). Therefore, we emphasize that all health caregivers in daily practice should be alert when dealing with high-risk TBI patients who have lung or urinary system illnesses, to prevent the occurrence of pneumonia, respiratory failure, urinary tract infection, and septicemia after TBI.

Some limitations should be noted. First, all diagnoses were based on claims data and ICD-9-CM diagnosis codes, which might result in disease misclassifications, missing data, or reporting discrepancies. However, all diagnostic records were coded by medical professionals, and many studies have validated the accuracy of the main diagnosis codes in the NHIRD [57–60]. Second, the level of disease severity, imaging information, and behavioral information such as smoking or alcohol consumption, were unavailable in the database. Although the CCI score, the length of ICU stay, and the length of NV days were used to estimate the severity level among TBI patients, validation using a clinical severity scale such as the Glasgow Coma scale or Barthel scale, as well as the severity variables mentioned above, should be considered in future research. In addition, the causes of TBI such as vehicular accidents or falls, which might be related to disease severity, were unknown in our study. Finally, although we excluded patients who died within one month after TBI diagnosis, multiple injuries may also affect the mortality risk of patients as they might have higher levels of clinical severity.

In conclusion, we used population-based health insurance claims data to demonstrate that grouping of post-TBI complications was associated with increased mortality risk. The findings of this study could improve the limitation of using a single disease to estimate mortality risk in future TBI research. These results can provide physicians and caregivers with important information to increase awareness about certain combinations of clinical syndromes among TBI patients.

Acknowledgments

This study is based in part on data obtained from the National Health Insurance Research Database provided by the Bureau of National Health Insurance of the Taiwan Department of Health and managed by the National Health Research Institutes. The interpretation and conclusions contained herein do not represent those of the Bureau of National Health Insurance, the Taiwan Department of Health, or the National Health Research Institutes. The authors would also like to thank Yi-Chen Chen, Yi-Ting Chen, Chia-Li Chang, and Hsin-Min Wu for their assistance in helping this project to move forward.

Citation: Ho C-H, Liang F-W, Wang J-J, Chio C-C, Kuo J-R (2018) Impact of grouping complications on mortality in traumatic brain injury: A nationwide population-based study. PLoS ONE 13(1): e0190683. https://doi.org/10.1371/journal.pone.0190683

References

1. Tagliaferri F, Compagnone C, Korsic M, Servadei F, Kraus J. A systematic review of brain injury epidemiology in Europe. Acta neurochirurgica. 2006;148(3):255–68. pmid:16311842

2. Alaranta H, Koskinen S, Leppänen L, Palomäki H. Nationwide epidemiology of hospitalized patients with first-time traumatic brain injury with special reference to prevention. Wiener medizinische Wochenschrift (1946). 1999;150(22):444–8.

3. Goldstein M. Traumatic brain injury: a silent epidemic. Annals of neurology. 1990;27(3):327-. pmid:2327740

4. Faul M, Xu L, Wald MM, Coronado VG. Traumatic brain injury in the United States: emergency department visits, hospitalizations, and deaths. Atlanta (GA): Centers for Disease Control and Prevention. National Center for Injury Prevention and Control. 2010;2:1–9.

5. Castriotta RJ, Wilde MC, Lai JM, Atanasov S, Masel BE, Kuna ST. Prevalence and consequences of sleep disorders in traumatic brain injury. Journal of clinical sleep medicine: JCSM: official publication of the American Academy of Sleep Medicine. 2007;3(4):349.

6. Hiott DW, Labbate L. Anxiety disorders associated with traumatic brain injuries. NeuroRehabilitation. 2002;17(4):345–55. pmid:12547982

7. Hoffman JM, Dikmen S, Temkin N, Bell KR. Development of posttraumatic stress disorder after mild traumatic brain injury. Archives of physical medicine and rehabilitation. 2012;93(2):287–92. pmid:22289239

8. Vitaz TW, Jenks J, Raque GH, Shields CB. Outcome following moderate traumatic brain injury. Surgical neurology. 2003;60(4):285–91. pmid:14505837

9. Masel BE, DeWitt DS. Traumatic brain injury: a disease process, not an event. Journal of neurotrauma. 2010;27(8):1529–40. pmid:20504161

10. Feigin VL, Barker-Collo S, Krishnamurthi R, Theadom A, Starkey N. Epidemiology of ischaemic stroke and traumatic brain injury. Best Practice & Research Clinical Anaesthesiology. 2010;24(4):485–94.

11. Chen YH, Kang JH, Lin HC. Patients with traumatic brain injury: population-based study suggests increased risk of stroke. Stroke. 2011;42(10):2733–9. pmid:21799162

12. Howle AA, Nott MT, Baguley IJ. Aspiration pneumonia following severe traumatic brain injury: prevalence and risk factors for long-term mortality. Brain Impairment. 2011;12(03):179–86.

13. Collaborators MCT, Perel P, Arango M, Clayton T, Edwards P, Komolafe E, et al. Predicting outcome after traumatic brain injury: practical prognostic models based on large cohort of international patients. Bmj. 2008;336(7641):425–9. pmid:18270239

14. Hukkelhoven CW, Steyerberg EW, Habbema JDF, Farace E, Marmarou A, Murray GD, et al. Predicting outcome after traumatic brain injury: development and validation of a prognostic score based on admission characteristics. Journal of neurotrauma. 2005;22(10):1025–39. pmid:16238481

15. Perel P, Edwards P, Wentz R, Roberts I. Systematic review of prognostic models in traumatic brain injury. BMC medical informatics and decision making. 2006;6(1):38.

16. Johnston SC, Mendis S, Mathers CD. Global variation in stroke burden and mortality: estimates from monitoring, surveillance, and modelling. The Lancet Neurology. 2009;8(4):345–54. pmid:19233730

17. Castriotta RJ, Lai JM. Sleep disorders associated with traumatic brain injury. Archives of physical medicine and rehabilitation. 2001;82(10):1403–6. pmid:11588744

18. Webster JB, Bell KR, Hussey JD, Natale TK, Lakshminarayan S. Sleep apnea in adults with traumatic brain injury: a preliminary investigation. Archives of physical medicine and rehabilitation. 2001;82(3):316–21. pmid:11245752

19. Piatetsky-Shapiro G. Discovery, analysis, and presentation of strong rules. Knowledge discovery in databases. 1991:229–38.

20. Cheng T- M. Taiwan’s new national health insurance program: genesis and experience so far. Health Affairs. 2003;22(3):61–76. pmid:12757273

21. Cheng T-M. Taiwan’s national health insurance system: high value for the dollar. Six Countries, Six Reform Models: Their Healthcare Reform: Experience of Israel, the Netherlands, New Zealand, Singapore, Switzerland and Taiwan Hackensack, NJ: World Scientific. 2009:171–204.

22. Deyo RA, Cherkin DC, Ciol MA. Adapting a clinical comorbidity index for use with ICD-9-CM administrative databases. Journal of clinical epidemiology. 1992;45(6):613–9. pmid:1607900

23. Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. Journal of chronic diseases. 1987;40(5):373–83. pmid:3558716

24. Shi H- Y, Hwang S- L, Lee K- T, Lin C- L. In-hospital mortality after traumatic brain injury surgery: a nationwide population-based comparison of mortality predictors used in artificial neural network and logistic regression models: Clinical article. Journal of neurosurgery. 2013;118(4):746–52. pmid:23373802

25. Thompson HJ, Dikmen S, Temkin N. Prevalence of comorbidity and its association with traumatic brain injury and outcomes in older adults. Research in gerontological nursing. 2012;5(1):17. pmid:22165997

26. Rothen HU, Stricker K, Einfalt J, Bauer P, Metnitz PG, Moreno RP, et al. Variability in outcome and resource use in intensive care units. Intensive care medicine. 2007;33(8):1329–36. pmid:17541552

27. Arabi Y, Venkatesh S, Haddad S, AL SHIMEMERI A, Al Malik S. A prospective study of prolonged stay in the intensive care unit: predictors and impact on resource utilization. International Journal for Quality in Health Care. 2002;14(5):403–10. pmid:12389806

28. Weissman C. Analyzing the impact of long-term patients on ICU bed utilization. Intensive care medicine. 2000;26(9):1319–25. pmid:11089759

29. Pasquier N, Bastide Y, Taouil R, Lakhal L. Discovering frequent closed itemsets for association rules. Database Theory—ICDT’99: Springer; 1999. p. 398–416.

30. Chen T- J, Chou L- F, Hwang S- J. Application of a data-mining technique to analyze coprescription patterns for antacids in Taiwan. Clinical Therapeutics. 2003;25(9):2453–63. pmid:14604744

31. Tu C- Y, Chen T- J, Chou L- F. Application of Frequent Itemsets Mining to Analyze Patterns of One-Stop Visits in Taiwan. PloS one. 2011;6(7):e14824. pmid:21747926

32. Liao JC, Ho CH, Liang FW, Wang JJ, Lin KC, Chio CC, et al. One-year mortality associations in hemodialysis patients after traumatic brain injury—an eight-year population-based study. PLoS One. 2014;9(4):e93956. pmid:24714730

33. Wee HY, Ho CH, Liang FW, Hsieh KY, Wang CC, Wang JJ, et al. Increased risk of new-onset depression in patients with traumatic brain injury and hyperlipidemia: the important role of statin medications. J Clin Psychiatry. 2016;77(4):505–11. pmid:27035519

34. Wu M- P, Lin H- J, Weng S- F, Ho C- H, Wang J- J, Hsu Y- W. Insomnia Subtypes and the Subsequent Risks of Stroke Report From a Nationally Representative Cohort. Stroke. 2014;45(5):1349–54. pmid:24699057

35. O'Donnell MJ, Xavier D, Liu L, Zhang H, Chin SL, Rao-Melacini P, et al. Risk factors for ischaemic and intracerebral haemorrhagic stroke in 22 countries (the INTERSTROKE study): a case-control study. The Lancet. 2010;376(9735):112–23.

36. Yasseen B, Colantonio A, Ratcliff G. Prescription medication use in persons many years following traumatic brain injury. Brain injury. 2008;22(10):752–7. pmid:18787984

37. van Oostrom SH, Picavet HS, de Bruin SR, Stirbu I, Korevaar JC, Schellevis FG, et al. Multimorbidity of chronic diseases and health care utilization in general practice. BMC Fam Pract. 2014;15:61. pmid:24708798

38. DuGoff EH, Canudas-Romo V, Buttorff C, Leff B, Anderson GF. Multiple chronic conditions and life expectancy: a life table analysis. Med Care. 2014;52(8):688–94. pmid:25023914

39. Wang L, Palmer AJ, Cocker F, Sanderson K. Multimorbidity and health-related quality of life (HRQoL) in a nationally representative population sample: implications of count versus cluster method for defining multimorbidity on HRQoL. Health Qual Life Outcomes. 2017;15(1):7. pmid:28069026

40. Calderon-Larranaga A, Poblador-Plou B, Gonzalez-Rubio F, Gimeno-Feliu LA, Abad-Diez JM, Prados-Torres A. Multimorbidity, polypharmacy, referrals, and adverse drug events: are we doing things well? Br J Gen Pract. 2012;62(605):e821–6. pmid:23211262

41. Rughani AI, Dumont TM, Lu Z, Bongard J, Horgan MA, Penar PL, et al. Use of an artificial neural network to predict head injury outcome. Journal of neurosurgery. 2010;113(3):585–90. pmid:20020844

42. Shi H- Y, Hwang S- L, Lee K- T, Lin C- L. In-hospital mortality after traumatic brain injury surgery: a nationwide population-based comparison of mortality predictors used in artificial neural network and logistic regression models. Journal of neurosurgery. 2013;118(4):746–52. pmid:23373802

43. Eftekhar B, Ketabchi E, Ardebili HE, Mohammad K, Ghodsi M. Comparison of artificial neural network and logistic regression models for prediction of mortality in head trauma based on initial clinical data. BMC Medical Informatics and Decision Making. 2005;5(1):3.

44. Hamele M, Stockmann C, Cirulis M, Riva-Cambrin J, Metzger R, Bennett TD, et al. Ventilator-associated pneumonia in pediatric traumatic brain injury. Journal of neurotrauma. 2016;33(9):832–9. pmid:26203702

45. Jovanovic B, Milan Z, Markovic-Denic L, Djuric O, Radinovic K, Doklestic K, et al. Risk factors for ventilator-associated pneumonia in patients with severe traumatic brain injury in a Serbian trauma centre. International Journal of Infectious Diseases. 2015;38:46–51. pmid:26166697

46. Lockey D, Coats T, Parr M. Aspiration in severe trauma: a prospective study. Anaesthesia. 1999;54(11):1097–8. pmid:10540100

47. Drakulovic MB, Torres A, Bauer TT, Nicolas JM, Nogué S, Ferrer M. Supine body position as a risk factor for nosocomial pneumonia in mechanically ventilated patients: a randomised trial. The Lancet. 1999;354(9193):1851–8.

48. Launey Y, Nesseler N, Le Cousin A, Feuillet F, Garlantezec R, Mallédant Y, et al. Effect of a fever control protocol-based strategy on ventilator-associated pneumonia in severely brain-injured patients. Critical Care. 2014;18(6):1.

49. Jawad I, Lukšić I, Rafnsson SB. Assessing available information on the burden of sepsis: global estimates of incidence, prevalence and mortality. population. 2012;2(3):4.

50. Venturi L, Miranda M, Selmi V, Vitali L, Tani A, Margheri M, et al. Systemic sepsis exacerbates mild post-traumatic brain injury in the rat. Journal of neurotrauma. 2009;26(9):1547–56. pmid:19257801

51. Shultz SR, McDonald SJ, Haar CV, Meconi A, Vink R, van Donkelaar P, et al. The potential for animal models to provide insight into mild traumatic brain injury: translational challenges and strategies. Neuroscience & Biobehavioral Reviews. 2016.

52. Roozenbeek B, Maas AI, Menon DK. Changing patterns in the epidemiology of traumatic brain injury. Nature Reviews Neurology. 2013;9(4):231–6. pmid:23443846

53. Liesemer K, Riva-Cambrin J, Bennett KS, Bratton SL, Tran H, Metzger RR, et al. Use of Rotterdam CT Scores for Mortality Risk Stratification in Children With Traumatic Brain Injury. Pediatric critical care medicine: a journal of the Society of Critical Care Medicine and the World Federation of Pediatric Intensive and Critical Care Societies. 2014.

54. Sarani B, Temple-Lykens B, Kim P, Sonnad S, Bergey M, Pascual JL, et al. Factors associated with mortality and brain injury after falls from the standing position. Journal of Trauma and Acute Care Surgery. 2009;67(5):954–8.

55. Sanchez-Olmedo J, Flores-Cordero J, Rincon-Ferrari M, Perez-Ale M, Munoz-Sanchez M, Dominguez-Roldan J, et al., editors. Brain death after severe traumatic brain injury: the role of systemic secondary brain insults. Transplantation proceedings; 2005: Elsevier.

56. Thompson HJ, Rivara FP, Jurkovich GJ, Wang J, Nathens AB, MacKenzie EJ. Evaluation of the effect of intensity of care on mortality after traumatic brain injury. Critical care medicine. 2008;36(1):282. pmid:18007264

57. Cheng C- L, Lee C- H, Chen P- S, Li Y- H, Lin S- J, Yang Y- HK. Validation of acute myocardial infarction cases in the national health insurance research database in taiwan. Journal of Epidemiology. 2014;24(6):500–7. pmid:25174915

58. Hsieh C- Y, Chen C- H, Li C- Y, Lai M- L. Validating the diagnosis of acute ischemic stroke in a National Health Insurance claims database. Journal of the Formosan Medical Association. 2015;114(3):254–9. pmid:24140108

59. Cheng C- L, Chien H- C, Lee C- H, Lin S- J, Yang Y- HK. Validity of in-hospital mortality data among patients with acute myocardial infarction or stroke in National Health Insurance Research Database in Taiwan. International journal of cardiology. 2015;201:96–101. pmid:26292275

60. Hung L- C, Sung S- F, Hsieh C- Y, Hu Y- H, Lin H- J, Chen Y- W, et al. Validation of a novel claims-based stroke severity index in patients with intracerebral hemorrhage. Journal of Epidemiology. 2016.