key: cord-0028066-01qrprml authors: Pelletier, Jonathan H.; Rakkar, Jaskaran; Simon, Dennis; Au, Alicia K.; Fuhrman, Dana Y.; Clark, Robert S.B.; Kochanek, Patrick M.; Horvat, Christopher M. title: Association between pediatric TBI mortality and median family income in the United States: A retrospective cohort study date: 2022-01-05 journal: Lancet Reg Health Am DOI: 10.1016/j.lana.2021.100164 sha: 1a48f4748eb00c7e468d2405d3c86fa0625d9d52 doc_id: 28066 cord_uid: 01qrprml BACKGROUND: There are regional disparities in pediatric traumatic brain injury (TBI) mortality across the United States, but the factors underlying these differences are unclear. METHODS: We performed a retrospective cross-sectional analysis of the Pediatric Health Information System database including inpatient hospital encounters for children less than 18 years old with a primary diagnosis of TBI between 2010–2019. FINDINGS: Lower median family income was associated with pediatric TBI mortality. Encounters from zip-codes with a median family income of <$20,000 had a 3.1% (29/950) mortality, as opposed to 1.3% (29/2,267) mortality for zip-codes with a median family income of >$80,000 (p = 0.00096). In multivariable logistic regression, every $10,000 of income was associated with an odds ratio of mortality of 0.94 (95% confidence interval 0.90 – 0.98). 82.5% (397/481) of ballistic TBI injuries were caused by a firearm. Lower income was associated with a higher proportion of ballistic TBI injuries (2.5% [24/950] for <$20,000 versus 0.3% [7/2,267] for >$80,000, p < 0.0001). In multivariable logistic regression, ballistic TBI injuries were associated with an odds ratio of mortality of 5.19 (95% confidence interval 4.00 – 6.73). United States regional variation in pediatric TBI mortality was linearly associated with the percentage of ballistic TBI (adjusted r-squared 0.59, p = 0.0097). INTERPRETATION: Children from lower income zip-codes are more likely to sustain a ballistic TBI, and more likely to die. Further work is necessary to determine causal factors underlying these associations and to design interventions that prevent these injuries and/or improve outcomes. Pediatric traumatic brain injury (TBI) has a worldwide incidence of 47-280 per 100,000 children, with male predominance after age three years. 1 The estimated annual incidence of pediatric TBI in the United States (US) is ~70/100,000 children, with a bimodal age distribution, peaking at ages less than five years old and in the late teenage years. 1, 2 While US pediatric TBI mortality decreased between 2000-2010, associated with decreased automobile-related TBI, mortality rates between 2010-2017 increased slightly, predominantly driven by increases in suicides and firearm violence. 3 US pediatric and adult TBI mortality varies by region, with higher mortality in more rural areas and the southeastern US. [3] [4] [5] However, it is unclear whether these outcomes are related to differences in TBI mechanisms, 3, 6 access to pediatric trauma centres, 7, 8 or other factors. The relationship between socioeconomic factors and pediatric TBI mortality is also unclear. Single-center and county-level studies have come to varying conclusions regarding the influence of socioeconomic status on pediatric trauma outcomes. 9, 10 Lastly, the effects of the rise in the incidence of TBI caused by firearm violence 3 on pediatric TBI mortality are also not well described. Given this uncertainty, we sought to analyze sociodemographic and regional factors associated with pediatric TBI mortality using a large multi-centre database of US children's hospitals. 11 American, American Indian or Alaska Native, Asian, Native Hawaiian or Pacific Islander, and other, as required by the United States Office of Management and Budget) 13 , ethnicity (Hispanic or Not Hispanic), age, sex, Rural-Urban Commuter Area (RUCA) codes, 14 zipcode median income, admission diagnosis, intensive care unit (ICU) admission, use of invasive mechanical ventilation (IMV), and complex chronic conditions (CCC). 15, 16 To ensure that changes in admission numbers reflected seasonal trends, rather than database expansion (as the PHIS database has grown over the past decade) 11 we restricted the analysis to hospitals providing data since 2010. Patients age zero to eighteen years were eligible for inclusion if they were discharged between January 1st, 2010 and December 31st, 2019 from one of the PHIS centres with a primary encounter diagnosis of TBI, according to International Classification of Diseases (ICD) version 9 or version 10 codes, as defined by the US Military Health System. 17 Diagnoses were filtered to exclude sequela and subsequent encounters. Hospitals were grouped according to their US Census Division. 18 Outcomes included the number of monthly and annual admissions over time, and in-hospital mortality, defined in PHIS. 11 This study was granted an exemption by the University of Pittsburgh Institutional Review Board, as it was a secondary analysis of an anonymized database. This manuscript follows the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guidelines for cross-sectional studies. 19 Cohort demographics were described using summary statistics and t-tests, chi-square tests, or Fisher exact tests as appropriate. Hospital charges were adjusted by the Centers for Medicare & Medicaid Services (CMS) wage/price index according to hospital zip code in PHIS. 11 Charge-over-time analyses were adjusted for the quarterly Gross Domestic Product (GDP) provided by the Bureau of Economic Analysis and expressed in 2010 dollars. 20 Patients were classified as having ballistic TBI if their encounter included a secondary diagnostic code for ballistic injuries. For analyses incorporating family income, each encounter was mapped to patient-level zip-code 2010 median family income using PHIS. Admission numbers were transformed into time-series data based on the date of admission and mapped to US Census Region. Chi-square testing was used to compare mortality proportions across hospitals, regions, and time. Seasonality testing was performed using Webel and Ollech's method. 21 Admission season was defined as winter (December-February), spring (March-May), summer (June-August), or autumn (September-November) per United States weather patterns. Quantification of seasonal trends in admission rates between 2010-2019 were displayed using LOcally Estimated Scatterplot Smoothing (LOESS). 22 For analyses considering RUCA codes, 14 zip-code level RUCA codes were extracted from PHIS and categorized according to population density and geographic isolation to match the categorizations of the US Census Bureau. 23, 24 To understand whether income and ballistic injuries were associated with TBI mortality after adjustment for confounders, we trained a stepwise multivariable logistic regression model using data available at the time of admission, including demographics, geography, seasonality, and admission diagnosis based on imaging findings using k-fold crossvalidation. 25 Age was included as a spline with knots at 1, 5, 10, and 15 years, as these groups were previously known to display differing mortality trends. 3 The linearity assumption for income and mortality was verified by visual inspection before model development. To understand whether regional disparities in income and ballistic injuries were associated with previously-known regional differences in pediatric TBI mortality, 3, 4 we conducted linear regression between these variables and pediatric TBI mortality across the US Census Bureau regions. We conducted two sensitivity analyses. In the first sensitivity analysis, we included TBI cases with codes that were definite or possible for TBI as defined by the US Military Health System (probable TBI). 17 In the second sensitivity analysis, as a proxy for severe TBI, we included only patients with coding for TBI and the use of invasive mechanical ventilation (IMV) and without a pre-existing complex chronic condition (TBI + IMV), as supplied by PHIS (as PHIS does not include Glasgow Coma Scale score data). 11 All statistical analyses were performed using RStudio version 1.4.1103 (RStudio, Boston, MA) and R versions 4.0.5 and 4.1.1 (R Foundation for Statistical Computing, Vienna, Austria) with the following packages: caret, cowplot, ggrepel, gifski, kableExtra, knitr, lattice, lubridate, precrec, seastests, splines tidyverse, usmap. The funders had no role in study design, data collection, analysis, interpretation, or writing of the report. TBI admissions met criteria for seasonality using Webel and Ollech's method. As shown in Supplemental Figure 1 , there was an overall summer-predominance of admissions, with a median (IQR) 504 (480 -526) in July compared to 344 (326 -358) in January (p = 0.00018). TBI admissions typically peaked in the afternoon and evening hours, as shown in Supplemental Figure 2 . Overall survival to discharge was 97.3% (49,491/50,872). There was a higher percentage of deaths in the winter compared with the summer, though the absolute number of deaths per season was similar (311/10,032 [3.1%] Figure 1 , median family income was inversely related to TBI mortality. Zip-codes with median family income < $20,000 were associated with an overall 3.1% (29/950) mortality, compared with 1.3% (29/2,267) mortality among zip-codes with a median family income greater than $80,000 (p = 0.00096). When stratified according to United States Department of Agriculture RUCA codes, this inverse association between median family income and mortality was present in urban areas and isolated small rural towns, but not in large rural cities or small rural towns, as shown in Supplemental Figure 3 . When stratified according to United States Office of Management and Budget race categories, there were significant differences in mortality according to income for patients of white race (3.0% [6/194] for income <$20,000 versus 1.1% [17/1,608] for income >$80,000, p = 0.041), but not for patients of other races, as shown in Supplemental Figure 4 . As shown in Figure 2 , there was more than 3-fold difference in TBI mortality across the United States Census Divisions (4.1% in the East South Central Division vs. 1.2% in the Pacific Division, p < 0.0001). Also shown in Figure 2 , there was 10-fold difference in TBI mortality at the hospital level (range 0.6%−6.1%, p < 0.0001). The 2010 median income ranged from $33,249 in the East South Central region to $51,964 in the New England region. As shown in Figure 3 , regional variation in overall TBI mortality was not significantly associated with 2010 median income (adjusted r-squared 0.20, p = 0.13). The Middle Atlantic region had substantially lower TBI mortality than other regions with similar income. After removing the Middle Atlantic region, median income was inversely associated with TBI mortality (adjusted r-squared 0.66, p = 0.0085, Supplemental Figure 5 ). There were 481/50,872 (0.9%) of TBI admissions with coding for gunshot wound or explosive device. Of these, 397/481 (82.5%) were injured by a firearm, 82/481 (17.0%) were injured by an air gun, and 2/481 (0.4%) were injured by an explosive. As shown in Figure 1 , family income was inversely associated with ballistic TBI. Zip-codes with median family income < $20,000 were associated with an overall 2.5% (24/950) ballistic TBI, compared with 0.3% (7/2,267) ballistic TBI among zip-codes with a median family income greater than $80,000 (p < 0.0001). Compared to the entire TBI cohort, patients with ballistic injuries were older, median (IQR) age 11.1 (4.8-14.7) vs. 5.7 (1.1-12.2), p < 0.0001 (Supplemental Table 1 As shown in Figure 4 , the regional percentage of ballistic TBI was inversely correlated with the 2010 median income for the region (adjusted r-squared 0.64, p = 0.0061). The Middle Atlantic region had notably lower ballistic TBI than other regions with comparable income. Regional variation in overall TBI mortality was directly associated with the percentage of TBI patients secondary to ballistic injuries (adjusted r-squared 0.59, p = 0.0097, Figure 5 ). The Middle Atlantic Region had mortality consistent with other regions with low proportions of ballistic TBI. Ballistic TBI admissions did not meet the criteria for seasonality using Webel and Ollech's method. The number of ballistic TBI admissions did not significantly change over time (Supplemental Figure 6 , 44 in 2010 versus 56 in 2019, p = 0.26). The percentage of TBI patients secondary to ballistic injuries ranged from 0.3% in the New England region to 2% in the East South Central region. A multivariable logistic regression model to elucidate factors associated with pediatric TBI mortality is shown in Table 2 . The linearity assumption for family income is shown in Supplemental Figure 7 . The spline relationship between age and mortality is shown in Supplemental Figure 8 . After adjustment, lower median family income was significantly associated with pediatric TBI mortality (odds ratio 0.94, 95% CI 0.90 -0.98 for every $10,000 of family income). Ballistic injury was also significantly associated with mortality (odds ratio 5.19, 95% CI 4.00 -6.73). Regional variation in TBI mortality remained significant, with patients in the East South Central United States having an odds ratio (95% CI) of mortality of 2.09 (1.61 -2.71) compared to patients in the Pacific region. The demographics for the sensitivity analyses including patients with probable TBI and restricting the cohort to TBI patients undergoing invasive mechanical ventilation (TBI + IMV) are shown in Supplemental Table 1 and stratified by mortality in Supplemental Table 2 and Supplemental Table 3 , respectively. Regional-and hospital-level variation for these sensitivity analyses is shown in Supplemental Figure 9 for the probable TBI group and Supplemental Figure 10 for the TBI + IMV group. Multivariable logistic regression for factors associated with mortality in the probable TBI group is shown in Supplemental Table 4 and in the TBI + IMV group is shown in Supplemental Table 5 . Briefly, median family income was inversely associated with mortality in the Probable TBI cohort (OR 0.94, 95% CI 0.90-0.98 per $10,000), the point estimate was similar in the smaller TBI + IMV cohort, but the association did not reach statistical significance (OR 0.95, 95% CI 0.87-1.05 per $10,000). Ballistic injury was significantly associated with mortality in both cohorts (Probable TBI OR 5.38, 95% CI 4.16-6.96, TBI + IMV OR 10.68, 95% CI 6. 45-17.70 ). This study reports the novel finding that zip-code level income is associated with pediatric TBI mortality in the United States. While previous studies have come to conflicting conclusions concerning the association between socioeconomic status and pediatric trauma outcomes, 9,10 the present study is in agreement with adult data showing lower income zip-codes have higher in-hospital trauma mortality. 26 Figure 1 shows that ballistic injuries increased significantly in zip-codes with the lowest 2010 median family incomes. Zip-codes with a 2010 median household income less than $40,000 had a 36% higher pediatric TBI mortality rate than those with a median income greater than $40,000. In multivariable logistic regression analysis, every $10,000 of family income was associated with an odds ratio of mortality of 0.94 (Table 2 ). In an additional novel analysis, this mortality difference was especially notable in areas classified as "urban" by RUCA codes, as shown in Supplemental Figure 3 , arguing that associations between income and mortality do not simply mirror distance to trauma centres. 7, 8, 26 While the point-estimate was similar, the association between income and mortality was not significant in the sensitivity analysis of patients receiving invasive mechanical ventilation (Supplemental Table 5 ). It is unclear whether this represents true differences in care-seeking patterns or is due to reduced power or possible model overfitting, as this cohort was only 10.3% of the size of the main analysis. Additionally, the association between income and mortality was present among patients of white race, but not other United States Office of Management and Budget race categories (Supplemental Figure 4) . Further work is needed to understand whether this represents differences in care-seeking patterns or TBI mechanisms among these populations, or reflects structural racism. [27] [28] [29] Ballistic TBI (82% of which was caused by firearms) represents a small, but particularly severe phenotype. This finding is in keeping with recently published trauma literature. 30 While representing less than 1% of pediatric TBI, these patients were nearly nine times more likely to die than the overall cohort. A public health intervention aimed at reducing firearm-related TBI in children 31 could also reduce pediatric TBI mortality. Given the disproportionate incidence and associated increased mortality of ballistic TBI in communities with lower household-level income (Figure 1 ), such an intervention would also align with the American Academy of Pediatrics' goals of reducing health disparities. 32 Interestingly, while overall rates of pediatric TBI decreased in the database between 2010-2019, ballistic TBIs did not (Supplemental Figure 6 ). In some states, physicians have been restricted from inquiring about gun ownership and offering counselling on safe gun storage for children. 33, 34 Our data suggest that this may undermine physician efforts to reduce mortality in pediatric TBI. The present study also offers insight into regional differences in pediatric TBI mortality. As in previous studies, we found a greater than three-fold difference in TBI mortality across the US Census Divisions (Figure 2 ). [3] [4] [5] Previous work has shown these differences to persist after adjustment for hospital characteristics. 5 We demonstrate that this difference is robust to adjustment for sociodemographic factors (e.g. age, sex, RUCA codes). Notably, we also show correlations between income, ballistic TBI, and TBI mortality (Figures 3-5) . Taken together, these trends suggest that actions to reduce poverty and prevent ballistic injuries may reduce disparities in pediatric TBI mortality. Caution is needed when interpreting the meaning of these findings. While the present study identified ten-fold interhospital variation in TBI mortality, 35 it does not necessarily imply discrepancies in inpatient care. Previous work has noted disparities in pediatric TBI mortality persist after adjustment for hospital characteristics. 5 While the practice patterns of high-performing centres should be analyzed as part of a positive-deviance approach to quality improvement, 36 survival in TBI is dependent upon a complex chain of healthcare delivery. Many previous studies have demonstrated that rapid access to pediatric trauma centres reduces mortality. 7, 8, [37] [38] [39] [40] While we did perform multivariable adjustment for RUCA codes, we were unable to precisely determine important factors such as transport time, which may be more important in improving pediatric TBI mortality than changes to inpatient practice. Studies from developing nations have shown that improving firstresponder care in underserved areas can reduce trauma mortality. 41 Additional analyses, including more granular data and additional, non-children's hospital trauma centres, are needed before designing public health interventions. Our study also shows the seasonality of pediatric TBI, with greater incidence in the summer, but stable absolute mortality, as shown in Supplemental Figure 1 . Despite the differing mechanisms, the overall summer predominance of TBI is similar to that for trauma overall. 15, 42 This summer-predominant seasonality of TBI is likely a reflection of seasonal changes in recreational behaviour, as sports represent a significant aetiology of pediatric TBI. 3, 6 The summertime increase in overall incidence coupled with no change in absolute mortality seen in Supplemental Figure 1 may suggest that there is an increase in milder sports-related TBI during the warmer months. Another likely cause of increased incidence of TBI in the summer and spring is increased automobile traffic during these months, 43 as automobile accidents are a predominant cause of pediatric TBI. 3, 6 This study has several limitations. While including over 50,000 encounters from 38 hospitals across the US, these data may not be fully representative of the US pediatric population. The median income of the included sample was $40,643, compared favourably to a nationwide average of $49,445 in 2010. 44 However, PHIS predominantly includes large children's hospitals, many of which are dedicated pediatric trauma centres. Because access to pediatric trauma care has previously been shown to reduce mortality, 7,8 the present study may underestimate mortality rates among children hospitalized with TBI. Pediatric trauma care is known to be inconsistently structured across the United States, and the distance to trauma care is related to mortality. 26, 45 Because income data, but not geocoded zip-code data, are available in PHIS, it was not possible to conduct an analysis regarding the distance between a patient's home and care center. We attempted to account for this relationship through the use of RUCA codes, but recognize that a geocoding approach would have been more precise. Though PHIS is quality controlled, 11 all database analysis is subject to misclassification bias from errors in diagnostic coding. Additionally, we were unable to completely discern the factors contributing to mortality for each TBI. We attempted to adjust for the intrinsic heterogeneity of TBI by classifying diagnoses according to type (e.g. concussion versus subarachnoid haemorrhage), but this approach is limited in its granularity. While we excluded subsequent encounters for the same injury based on ICD codes, approximately 1% of patients had more than one TBI encounter. The clinical information in PHIS is limited, thus we could not fully control for severity of illness when assessing interhospital variability in mortality. Lastly, though using IMV is commonly used as a surrogate for severe TBI, it is imperfect. 46 In conclusion, lower family income and ballistic TBI are associated with pediatric TBI mortality in the United States. These associations are significant after multivariable adjustment for a variety of confounders. Pediatric TBI mortality varies three-fold regionally, and ten-fold across hospitals in the PHIS database. This variability is correlated with regional income and rates of ballistic TBI. Further work is needed to develop targeted interventions to address these disparities. Refer to Web version on PubMed Central for supplementary material. National Institutes of Health. Evidence before this study Mortality in pediatric traumatic brain injury (TBI) varies across the United States (US), with more rural areas and the Southeastern US experiencing higher mortality. However, it is unclear whether these differences are related to differences in TBI mechanisms, pre-hospital care, trauma care, or sociodemographic factors. For the first time, this study demonstrates an association between median family income (assessed at the zip-code level) and pediatric TBI mortality in the United States. This association between income and TBI mortality was present when stratified according to rural / urban status, and in multivariable logistic regression after adjustment for confounders. A second association between the proportion of ballistic (predominantly firearm) TBI and mortality was also demonstrated. On a regional basis, higher levels of ballistic TBI were linearly associated with increases in TBI mortality. Median family income and ballistic injuries are associated with pediatric TBI mortality. Further research is necessary to determine whether interventions to reduce poverty and firearm violence may reduce pediatric TBI mortality. Variation in TBI Mortality. Panel A shows hospitals represented according to their US Census region. Alaska and Hawaii were excluded as there are no PHIS member hospitals in those states. For each region, the mortality is normalized to the total number of admissions in that region and expressed as a percentage according to the scale on the right. Panel B is a barplot of TBI Mortality by Hospital. Individual hospitals are displayed on the x-axis, with % of TBI mortality on the y-axis. The numbers within the bars represent the average mortality at that hospital over the 10-year study period. Correlation between regional median household income and % TBI mortality. The 2010 median family income for the patient's zip code (grouped by US Census Division) is plotted on the x-axis and the % TBI mortality is plotted on the y-axis. The solid line represents the line of best fit using linear regression. The gray shaded region represents the 95% confidence interval of the model. Correlation between regional median household income and % ballistic TBI. The 2010 median family income for the patient's zip code (grouped by US Census Division) is plotted on the x-axis and the % ballistic TBI is plotted on the y-axis. The solid line represents the line of best fit using linear regression. The gray shaded region represents the 95% confidence interval of the model. Correlation between % ballistic TBI and % TBI mortality. The % ballistic TBI (grouped by US Census Division) is plotted on the x-axis and the % TBI mortality is plotted on the y-axis. The solid line represents the line of best fit using linear regression. The gray shaded region represents the 95% confidence interval of the model. Epidemiology of Global Pediatric Traumatic Brain Injury: Qualitative Review Incidence of Pediatric Traumatic Brain Injury and Associated Hospital Resource Utilization in the United States Traumatic brain injury mortality among U.S. children and adolescents ages 0-19 years Geographical Disparity and Traumatic Brain Injury in America: Rural Areas Suffer Poorer Outcomes Facility characteristics and inhospital pediatric mortality after severe traumatic brain injury Pediatric Traumatic Brain Injury: Characteristic Features, Diagnosis, and Management Pediatric trauma centers: Correlation of ACS-verified trauma centers with CDC statewide pediatric mortality rates The association of trauma center closures with increased inpatient mortality for injured patients Family Income is Not Associated with Outcomes in Pediatric Patients with Critical Traumatic Injury A Population-Based Analysis of Socioeconomic Status and Insurance Status and Their Relationship With Pediatric Trauma Hospitalization and Mortality Rates PHIS Corticosteroids and mortality in children with bacterial meningitis Census.gov. Census.gov USDA ERS -Rural-Urban Commuting Area Codes Trends in US Pediatric Hospital Admissions in 2020 Compared With the Decade Before the COVID-19 Pandemic Pediatric complex chronic conditions classification system version 2: updated for ICD-10 and complex medical technology dependence and transplantation Surveillance Case Definitions | Health.mil Census Divisions | Monitoring References | National Centers for Environmental Information (NCEI) The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement: guidelines for reporting observational studies Adjusting Health Expenditures for Inflation: A Review of Measures for Health Services Research in the United States Condensing information from multiple seasonality tests with random forests SAS/STAT ® 14.1 User's Guide Rural Urban Commuting Area Codes Data Metropolitan and Micropolitan Akaike H Information Theory and an Extension of the Maximum Likelihood Principle Associations of Distance to Trauma Care, Community Income, and Neighborhood Median Age With Rates of Injury Mortality Racism as a Determinant of Health: A Systematic Review and Meta-Analysis Intersection of Bias, Structural Racism, and Social Determinants With Health Care Inequities STRUCTURAL RACISM AND HEALTH INEQUITIES: Old Issues, New Directions Evaluation of Injury Severity and Resource Utilization in Pediatric Firearm and Sharp Force Injuries Firearm Violence Among Youth: Public Health Strategies for Prevention The Impact of Racism on Child and Adolescent Health Patients and Guns: Florida Physicians Are Not Asking Gag Laws" and Gun Safety Variation in Pediatric Traumatic Brain Injury Outcomes in the United States Positive deviance: a different approach to achieving patient safety Effect of predicted travel time to trauma care on mortality in major trauma patients in Nova Scotia Pediatric trauma mortality: an ecological analysis evaluating correlation between injury-related mortality and geographic access to trauma care in the United States in 2010 A National Analysis of Pediatric Trauma Care Utilization and Outcomes in the United States Access to Pediatric Trauma Care in the United States Rural pre-hospital trauma systems improve trauma outcome in low-income countries: a prospective study from North Iraq and Cambodia Zucker-braun NS. Seasonal, Weather, and Temporal Factors in the Prediction of Admission to a Pediatric Trauma Center Figure 2 -Travel on U.S. Highways By Month -December 2019 -Policy | Federal Highway Administration Income, Poverty and Health Insurance Coverage in the United States: 2010 -Income & Wealth Pediatric Traumatic Brain Injury Is Inconsistently Regionalized in the United States Early Use of Antiseizure Medication in Mechanically Ventilated Traumatic Brain Injury Cases: A Retrospective Pediatric Health Information System Database Study. Pediatr. Crit. Care Med Publish Ahead of Print Zip-Code Median Income (S) Author manuscript; available in PMC Author manuscript; available in PMC Admit Season (Reference: Summer) Admit Season: Spring Admit Season: Autumn Admit Season: Winter Admission Dx (Reference: Concussion) Admission Dx: Skull Fracture without ICH Admission Dx: Skull Fracture with ICH Admission Dx: Cranial Nerve or Carotid Artery Admission Dx: Epidural Hemorrhage Admission Dx: Subdural Hemorrhage 60.66 (19 Admission Dx: Subarachnoid Hemorrhage Admission Dx: Intracerebral / Intracerebellar Hemorrhage Admission Dx: Unspecified Intracranial Hemorrhage Admission Dx: Cerebral / Cerebellar Contusion or Laceration 60 Admission Dx: Diffuse Axonal Injury Admission Dx: Traumatic Cerebral Edema Author manuscript; available in PMC