Augmenting electronic health record data with social and environmental determinant of health measures to understand regional factors associated with asthma exacerbations ========================================================================================================================================================================= * Alana Schreibman * Kimberly Lactaoen * Jaehyun Joo * Patrick K. Gleeson * Gary E. Weissman * Andrea J. Apter * Rebecca A. Hubbard * Blanca E. Himes ## Abstract Electronic health records (EHRs) provide rich data for diverse populations but often lack information on social and environmental determinants of health (SEDH) that are important for the study of complex conditions such as asthma, a chronic inflammatory lung disease. We integrated EHR data with seven SEDH datasets to conduct a retrospective cohort study of 6,656 adults with asthma. Using Penn Medicine encounter data from January 1, 2017 to December 31, 2020, we identified individual-level and spatially-varying factors associated with asthma exacerbations. Black race and prescription of an inhaled corticosteroid were strong risk factors for asthma exacerbations according to a logistic regression model of individual-level risk. A spatial generalized additive model (GAM) identified a hotspot of increased exacerbation risk (mean OR = 1.41, SD 0.14, p < 0.001), and inclusion of EHR-derived variables in the model attenuated the spatial variance in exacerbation odds by 34.0%, while additionally adjusting for the SEDH variables attenuated the spatial variance in exacerbation odds by 66.9%. Additional spatial GAMs adjusted one variable at a time revealed that neighborhood deprivation (OR = 1.05, 95% CI: 1.03, 1.07), Black race (OR = 1.66, 95% CI: 1.44, 1.91), and Medicaid health insurance (OR = 1.30, 95% CI: 1.15, 1.46) contributed most to the spatial variation in exacerbation odds. In spatial GAMs stratified by race, adjusting for neighborhood deprivation and health insurance type did not change the spatial distribution of exacerbation odds. Thus, while some EHR-derived and SEDH variables explained a large proportion of the spatial variance in asthma exacerbations across Philadelphia, a more detailed understanding of SEDH variables that vary by race is necessary to address asthma disparities. More broadly, our findings demonstrate how integration of information on SEDH with EHR data can improve understanding of the combination of risk factors that contribute to complex diseases. **Author summary** Electronic health records constitute an important source of data for understanding the health of large and diverse real-world populations, however, they do not routinely capture socioeconomic and environmental factors known to affect health outcomes. We show how electronic health record data can be augmented to include individual measures of air pollution exposures, neighborhood socioeconomic status, and the natural and built environment using patients’ residential addresses to study asthma exacerbations, episodes of worsening disease that remain a major public health challenge in the United States. We found that on an individual patient-level, Black race and prescription of an inhaled corticosteroid were the factors most strongly associated with asthma exacerbations. In contrast, neighborhood deprivation, race, and health insurance type accounted for the most spatial variation in exacerbation risk across Philadelphia. Our findings provide insight into factors that contribute to asthma disparities in our region and present a framework for future efforts to expand the scope of electronic health record data. ## Introduction Electronic health records (EHRs) are a source of rich patient-level data for large and diverse populations that can be used for research due to their widespread availability [1]. However, EHR data often contain incomplete or low-quality measures of social and environmental determinants of health (SEDH), limiting their utility for the study of complex diseases [2]. Recent efforts to address this limitation have included developing methodologies to integrate external data on the physical, built, and social environment via linkage with patient addresses [3–5], with high-resolution geospatial datasets providing the closest estimate of individual exposures [6]. Integrated EHR and SEDH datasets can be used to understand both individual-level outcomes and patterns of risk across a spatial region, thereby providing insights for both precision medicine and precision public health efforts. Because many environmental exposures pose a greater risk to select groups of people, integrating external information on SEDH with EHR data is also helpful to address health disparities according to race, ethnicity, and socioeconomic status. Asthma, a chronic disease that is characterized by inflammation and reversible narrowing of the airways, affects over 25 million people or approximately 8% of the United States population [7]. Racial and ethnic disparities in its morbidity and mortality are well known, and those living in poverty are also more likely to have asthma [7–11]. The clinical goals of asthma management are to control patients’ symptoms and minimize long-term risk of lung function decline [12]. This includes preventing asthma exacerbations, episodes of worsening disease which require treatment with systemic steroids [13]. However, despite guideline-directed clinical management, asthma exacerbations remain common, contributing to asthma-related morbidity and mortality as well as higher health care costs and utilization [14,15]. Risk factors for asthma exacerbations in adults include female sex [16], obesity [17], current or past smoking [18], comorbid allergic rhinitis or chronic obstructive pulmonary disease (COPD) [17,19], and a history of previous exacerbations [20]. Observational studies have also found independent associations between asthma exacerbations and exposure to particulate matter (PM), gaseous air pollutants, and to mixtures of pollutants such as traffic-related air pollution (TRAP) [21]. Similarly, associations between asthma exacerbations and living in poverty [22], substandard housing conditions, such as the presence of mold and pests [23,24], and neighborhood “greenness” [25,26] have been documented. Because asthma exacerbations result from complex interactions among various biological, social, and environmental factors that vary across individuals and geographically, creating generalizable models of exacerbations remains an unachieved goal. Further, because in the United States there is a high correlation between minoritized race or ethnicity, poverty, and harmful environmental exposures, disentangling relationships among key variables is difficult [9,27]. Approaches that model many asthma-related variables in specific regions may lead to the identification of actionable strategies to reduce exacerbations locally, and using EHR data to create these models has the advantage of providing health information for the specific catchment region served by a given healthcare system. Few studies have linked EHRs with a diverse set of SEDH variables to study asthma exacerbations using patient-level data [28–30], and to our knowledge, none have used geospatial analysis techniques to understand the contribution of these factors to the spatial distribution of exacerbation risk. Here, we show how EHR data can be extended to include individualized measures of air pollution exposures, socioeconomic status indices, and measures of the natural and built environment to identify local factors associated with asthma exacerbation risk. ## Methods A retrospective cohort analysis was performed using de-identified EHR data from Penn Medicine, a large health system that serves the greater Philadelphia area, from encounters dated 1/1/2017-12/31/2020. An overview of our study design is shown in Fig 1. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/11/02/2024.10.24.24316063/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2024/11/02/2024.10.24.24316063/F1) Figure 1. Overview of study design. Graphical overview of study design, including processing and linkage of electronic health record (EHR) and social and environmental determinants of health (SEDH) data, cohort selection including spatial filtering by assessing the representation of the EHR cohort compared to the underlying population, and patient-level and geospatial analyses on the expanded EHR dataset. ### Ethics statement Our study was approved by the University of Pennsylvania Institutional Review Board (IRB) under protocol number 824789. Formal consent was not obtained, as a waiver of HIPAA Authorization was granted for the conduct of this research. ### Study population Patient-level encounter data was obtained for adults (i.e. age ≥ 18 years) who had at least one encounter with an International Classification of Diseases (ICD)-10 code for asthma (J45*) and who were prescribed a short-acting β2-agonist (SABA) (S1 Table). The most recent residential address for each patient was obtained and geocoded using previously described methods [3]. Demographic, comorbidity, and medication data for encounters during the study period was extracted and used to compute several variables, hereafter referred to as “EHR-derived variables.” These included: age at first encounter, sex, race, ethnicity, body mass index (BMI), health insurance type, smoking status, chronic obstructive pulmonary disease (COPD), allergic rhinitis, a modified Elixhauser score [31], inhaled corticosteroid (ICS) prescription, and years followed (defined as the number of years between first and last encounter). Additional details, including inclusion and exclusion criteria and definitions of the EHR-derived variables, are provided in S1 Text. ### Outcome Asthma exacerbations were defined as encounters with an oral corticosteroid (OCS) prescription (S1 Table) and either 1) a primary asthma diagnosis code (ICD-10, J45*) for encounters with primary diagnosis codes listed or 2) a nonprimary asthma ICD-10 code for encounters without a primary diagnosis listed but only one or two ICD-10 codes listed. A count of exacerbations during the study period was computed for each patient. ### SEDH data Seven external datasets were integrated with our EHR dataset via linkage with patient geocodes, creating variables that are hereafter referred to as “SEDH variables.” Additional details on data processing are reported in S1 Text, and dataset sources and spatiotemporal dimensions are summarized in S2 Table. All processed SEDH data is available in Sensor-based Analysis of Pollution in the Philadelphia Region with Information on Neighborhoods and the Environment (SAPPHIRINE), a web application that integrates spatially distributed high-resolution social and environmental data in the greater Philadelphia region to facilitate the conduct of local health studies [32]. #### Air pollution exposures Average pollutant exposures were assigned to each patient using high-resolution (∼1×1 km2) geophysical model estimates. 2017-2019 NO2 estimates (in parts per billion by volume, or ppbv) and 2017-2020 PM2.5 estimates (in μg/m3) were downloaded from resources reported in Cooper et al. and van Donkelaar et al., respectively, temporally averaged, and linked to the study cohort using bilinear interpolation [33,34]. Exposure to other toxic airborne chemicals from point sources was estimated using the EPA Toxics Release Inventory (TRI) [35]. Total toxic air release exposure by patient was computed as the sum of toxic releases (in kilograms) over the study period within a 1-km circular buffer of each patient’s residential address. Exposure to air pollution from mobile sources was estimated as the sum of the daily vehicle distance traveled (DVDT), a metric computed using traffic data published by the Pennsylvania Department of Transportation, within a 300-m circular buffer [36]. #### Neighborhood socioeconomic environment Socioeconomic disadvantage for each person was summarized as the Area Deprivation Index (ADI), a validated index computed for each Census block group based on several ACS variables [37]. A higher ADI score indicates greater disadvantage. 2018 ADI data was extracted from the Neighborhood Atlas and was assigned to each patient by block group [38]. #### Built and natural environment Asthma-related housing code violation data (i.e. pests, water damage, and indoor air contamination) was obtained from a Philadelphia Department of Licenses and Inspections dataset reported by OpenDataPhilly (S3 Table) [39]. For each block group, the number of violations during the study period per 100 people (based on the 2019 ACS population estimates) was computed and assigned to each patient. Vegetation density was summarized as the normalized difference vegetation index (NDVI), an index with values ranging from −1.0 to 1.0 where higher positive values represent higher vegetation density. NDVI was computed in Google Earth Engine using surface reflectance images from the Landsat 8 satellite during the study period and assigned to each patient as the mean value within a 300-m circular buffer [40]. ### Statistical analysis Analyses were conducted in R 4.2 [41]. #### Study area To minimize bias in our geospatial analyses introduced by uneven spatial density across the Penn Medicine catchment area, we restricted our study region to spatial areas in which the geospatial representativeness of our EHR cohort was adequate compared to the underlying population. Following methods described in Xie et al. [42], we computed a spatial representation ratio (SRR), defined as the ratio between the proportion of our EHR cohort living in each census block group and the proportion of the Philadelphia population living in that block group, as reported by the 2019 American Community Survey (ACS). We defined our study region as contiguous census tracts (i.e. adjacent or separated by a non-residential area such as a park) with a mean SRR of 0.5 or greater. Only patients who resided in this study region were included in analyses. #### Modeling individual-level risk factors Chi-squared and Kruskal-Wallis rank sum tests were used in bivariate analyses to assess associations between patient characteristics and asthma exacerbation level (i.e. 0, 1-2, 3-4, 5+) during the study period, and to compare the characteristics of complete cases to individuals with missing data. To identify patient-level factors associated with asthma exacerbations, we fit logistic regression models with asthma exacerbations as a binary case-control (0 vs >0) outcome. This approach was chosen to match the dichotomous outcome used in spatial analyses. Logistic regression models were initially adjusted for EHR-derived variables only, and then adjusted for both EHR-derived and SEDH variables (EHR & SEDH-adjusted). Years followed was included as a covariate in both models to account for variation in the length of available follow-up across patients. Model fit was assessed using the Akaike information criterion (AIC). We checked for multicollinearity by computing Pearson’s correlation coefficients between all EHR and SEDH variables, while selecting White race and Private health insurance type as reference levels for the two nominal categorical variables, and we computed variance inflation factors (VIF) for all independent variables. #### Sensitivity analysis for individual-level risk factors We conducted sensitivity analyses of individual-level risk factors by fitting negative binomial regression models with asthma exacerbations represented as a count outcome and comparing results to those of logistic regression models. The negative binomial regression models were adjusted first for EHR-derived variables only and then for EHR & SEDH variables, while including years followed as an offset in each model. #### Modeling spatial risk factors To estimate local odds of exacerbation as a function of location, spatial generalized additive models (GAMs) were fit with a binary case-control outcome (0 vs >0 exacerbations) on a grid of points across the study region, using the R *MapGAM* package and previously described methods (S1 Text) [3,43]. Maps of spatial effect predictions were created for the smoothed spatial term of each model, where the spatial odds ratio (OR) at each point represented the ratio between the odds of exacerbation at that point and the median odds across all points. First, a univariable model adjusted only for years followed, hereafter referred to as “unadjusted model”, was fit to identify hotspots and coldspots across the study region. Next, multivariable models adjusted for 1) only EHR-derived variables and 2) EHR & SEDH variables were compared to the unadjusted model by computing the mean and standard deviation (SD) of the spatial ORs at points within any overlapping hotspots, and by computing the percent difference between the variance in ORs across the full study region in the adjusted models and the unadjusted model. To understand the contribution of individual variables to the observed spatial effects, we fit additional models adjusted for each EHR-derived or SEDH variable one variable at a time (all models were also adjusted for years followed) and computed the percent difference between the variance in ORs in these models and the unadjusted model. Model fit was assessed using the AIC. #### Stratified analysis We conducted a stratified analysis to further evaluate the association between race and asthma exacerbations. Chi-square and Wilcoxon rank sum tests were used to assess bivariate relationships between race and other variables. To test whether any EHR or SEDH variables that were correlated with race influenced spatial patterns of asthma exacerbation risk independently of race, spatial GAMs adjusted one variable at a time were fit on each race stratum using the same approach described above. Before fitting the models, SRR selection was repeated for each stratum to account for the uneven geographic distribution of race across the initial study region. ## Results Following selection of patients based on inclusion/exclusion criteria and spatial filtering (S1 Fig), the retrospective study cohort consisted of 6,656 asthma patients, 2,329 of whom had one or more exacerbations (Table 1), with residence in 249 census tracts in Philadelphia (S2 Fig). The spatial distribution of all processed SEDH datasets within the study region is shown in S3 Fig. The EHR-derived variables years followed, age, race, BMI, health insurance type, smoking status, COPD, allergic rhinitis, Elixhauser comorbidity score, and ICS, as well as the SEDH variable ADI were the most significantly associated with exacerbations according to bivariate analyses (p < 0.001) (Table 1). Patients with more exacerbations during the study period were followed for more years, more likely to be aged 35-54, of Black race, or class 2 or class 3 obese, and more likely to have Medicaid health insurance, a history of smoking, COPD, allergic rhinitis, higher Elixhauser comorbidity scores, an ICS prescription, or live in neighborhoods with higher ADI. The proportion of patients in the study cohort who were prescribed controller medications including ICS, leukotriene modifiers, long-acting β2-agonists (LABA), and biologic therapies was positively associated with exacerbation count (S4 Table). Bivariate analyses comparing the distribution of characteristics between the study cohort and patients excluded due to missing EHR data found statistically significant differences in years followed, age, sex, health insurance type, COPD, allergic rhinitis, Elixhauser comorbidity score, and ICS (p < 0.001) (S5 Table). View this table: [Table 1.](http://medrxiv.org/content/early/2024/11/02/2024.10.24.24316063/T1) Table 1. Patient characteristics by exacerbation count levels. Shown are the characteristics of patients according to their number of exacerbations during the study period. For each exacerbation level, the number and percentage of patients are shown for categorical variables, and the Median and Interquartile Range (IQR) are shown for continuous variables. Based on Pearson’s correlation coefficients, there were no strong correlations between any EHR-derived or SEDH variables, except for NO2 and PM2.5 (ρ = 0.77) (S4 Fig). Moderate correlations (|ρ| > 0.50) were observed between age and Medicare health insurance, Black race and ADI, and NO2 and NDVI. Furthermore, adjusted generalized VIFs for all variables did not exceed 2.00, suggesting that all variables could be included in multivariable models (S6 Table). ### Individual risk factors associated with asthma exacerbations ORs for the EHR-derived and SEDH variables included in the multivariable logistic regression model are summarized in Table 2. The EHR-derived variables years followed, age 35-54, Black race, and ICS had the strongest positive associations with having at least one exacerbation during the study period (p < 0.001) (Table 2). These effects persisted after additionally adjusting for SEDH variables, and of the SEDH variables in the logistic regression model, only NO2 exposure was positively associated with exacerbations (p = 0.0059). Inclusion of the SEDH variables did not improve model fit as determined by AIC (8,320 for both the EHR-adjusted and EHR & SEDH-adjusted models). Sensitivity analyses showed that the risk factors identified by logistic and negative binomial regression models were mostly consistent, with Black race and ICS prescription having the strongest effects between both models (p < 10−4), and variables such as Medicaid health insurance and Elixhauser comorbidity score of 1-9 having statistically significant p-values in both models though smaller in the negative binomial regression (p < 0.001) than the logistic regression (p < 0.05) (S7 Table). However, the negative binomial regression model did not identify statistically significant associations between asthma exacerbations and age 35-54 or NO2 exposure as did the logistic regression, although the directions of effect were consistent in both models. In addition, the negative binomial regression model identified an association with the SEDH variable housing code violations (p = 0.045) that was not present in the logistic regression model. View this table: [Table 2.](http://medrxiv.org/content/early/2024/11/02/2024.10.24.24316063/T2) Table 2. Individual-level asthma exacerbation risk factors in multivariable logistic regression models. Shown are the adjusted odds ratios (ORs), 95% confidence intervals (CIs), and p-values for logistic regression models of asthma exacerbations as a dichotomous outcome adjusted for EHR-derived variables only and for both EHR and SEDH variables. ### Spatial risk factors associated with exacerbations Maps of ORs across the study region for the unadjusted, EHR-adjusted, and EHR & SEDH-adjusted spatial GAMs are shown in Fig 2. In the unadjusted model, the global test of the null hypothesis that asthma exacerbations were not associated with geographic location was significant (p < 0.001) (Fig 2A). Local tests identified a statistically significant hotspot of exacerbations (p < 0.01) in West and South Philadelphia with a mean spatial OR of 1.41 (SD 0.14). In the model adjusted for EHR-derived variables (Fig 2B) the global test statistic remained significant (p < 0.001). Local tests again identified a hotspot (p < 0.01) in West and South Philadelphia which had a decreased mean spatial OR of 1.27 (SD 0.055). In this EHR-adjusted model, the variance in spatial ORs across the study region was 34.0% lower than the variance of the unadjusted model (S5 Fig). In the model adjusted for both EHR-derived and SEDH variables, the global test statistic remained significant (p < 0.001, Fig 2C) and local tests identified a hotspot (p < 0.01) in West Philadelphia that overlapped geographically with the one in the other models, although of a smaller area and with a smaller mean spatial OR of 1.24 (SD 0.042) compared to the EHR-adjusted model. The variance in spatial ORs in the EHR & SEDH-adjusted model was strongly attenuated (66.9% lower than that of the unadjusted model), suggesting that these variables partially explained the spatial correlation of exacerbations (S5 Fig). The ORs for the other terms included in the spatial GAMs (i.e., the EHR-derived and SEDH variables) are summarized in Table 3. All variables (i.e., years followed, age 35-54, Black race, class 3 obesity, Medicaid health insurance, Elixhauser score 1-9, ICS) that were significant in the multivariable logistic regression models (Table 2) were also significant in the spatial GAMs (Table 3), except for NO2 exposure, which was significant in the logistic regression model but not in the spatial GAM. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/11/02/2024.10.24.24316063/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2024/11/02/2024.10.24.24316063/F2) Figure 2. Spatial odds ratios (ORs) of exacerbations before and after adjusting for EHR-derived and SEDH variables. (A) Unadjusted spatial GAM (adjusted only for years followed). (B) Spatial GAM adjusted for EHR-derived variables only. (C) Spatial GAM adjusted for both EHR-derived and SEDH variables. Base maps were created using the Stamen Design from Stadia Maps. View this table: [Table 3.](http://medrxiv.org/content/early/2024/11/02/2024.10.24.24316063/T3) Table 3. Spatial asthma exacerbation risk factors in multivariable spatial GAMs. Shown are the adjusted odds ratios (ORs), 95% confidence intervals (CIs), and p-values for spatial GAMs of asthma exacerbations as a dichotomous outcome adjusted for EHR-derived variables only and for both EHR-derived and SEDH variables. In the spatial GAM models adjusted one variable at a time, ADI, race, and health insurance type most strongly attenuated the variation in spatial ORs (Fig 3), individually accounting for 55.2%, 38.5%, and 26.5%, respectively, of the variation in the unadjusted model (S5 Fig). In these models, each variable was positively associated with exacerbations (p < 0.001): ADI with OR = 1.05 (95% CI: 1.03, 1.07); Black race with OR = 1.66 (95% CI: 1.44, 1.91); and Medicaid health insurance with OR = 1.30 (95% CI: 1.15, 1.46) (S8 Table). The spatial distribution of these variables followed similar patterns: high ADI and high density of Black patients and patients with Medicaid insurance cooccurred in West, North, and South Philadelphia (Fig 3). No other EHR (S6 Fig) or SEDH (S7 Fig) variable attenuated the hotspot area or its effect size, however, adjustment for NO2 levels resulted in an expansion of a coldspot area (S7A Fig). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/11/02/2024.10.24.24316063/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2024/11/02/2024.10.24.24316063/F3) Figure 3. Spatial distribution of individual variables that most strongly attenuated the spatial odds ratios (ORs) of exacerbations along with corresponding spatial GAM results adjusted for these individual variables. Spatial distribution in the study region of (A) the area deprivation index (ADI), (B) race, and (C) health insurance type of patients. Corresponding spatial GAMs adjusted only for years followed and (D) ADI, (E) race, or (F) health insurance type. ### Stratified analysis by race Bivariate analysis comparing the distribution of all EHR and SEDH characteristics between patients of Black and White race found statistically significant differences for all variables except ethnicity and ICS (Table 4). Given that ADI, race, and health insurance type had the strongest relationship with asthma exacerbations in spatial analyses, and their spatial distributions were similar, we performed stratified spatial analyses by race to determine whether ADI and health insurance status remained significantly associated with asthma exacerbations within race-stratified groups. After applying SRR inclusion criteria in strata according to race, an additional 265 Black patients were excluded from the stratified spatial analysis, resulting in a sample size of 4,363 patients (Fig 4A). The unadjusted spatial GAM for patients of Black race had a significant global test statistic (p < 0.001) and local tests identified hotspots in West and South Philadelphia consistent with results of the full cohort (Fig 4B). In contrast to spatial analyses in the full cohort, adjusting for ADI and health insurance separately for Black patients did not attenuate the spatial variance in ORs, instead increasing the variance by 4.70% and 2.79%, respectively, compared to the unadjusted model (Fig 4C and 4D). In patients of White race, 73 additional patients were removed after applying SRR inclusion criteria, resulting in 1,383 patients included in the spatial model for which the global test statistic was no longer significant (p = 0.098, S8 Fig). View this table: [Table 4.](http://medrxiv.org/content/early/2024/11/02/2024.10.24.24316063/T4) Table 4. Patient characteristics by race. Shown are the number and percentage of patients in each level for categorical variables, and the Median and Interquartile Range (IQR) for continuous variables in patients of White race versus Black race. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/11/02/2024.10.24.24316063/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2024/11/02/2024.10.24.24316063/F4) Figure 4. Spatial odds ratios (ORs) of exacerbations among Black patients along with the effects of ADI and health insurance type on this distribution. (A) SRR values for the updated study region used in spatial GAMs for patients of Black race only (SRR = 1 indicates no representativeness bias). (B) Unadjusted spatial GAM (adjusted only for years followed) for patients of Black race (N = 4,363). Spatial GAMs adjusted additionally for (C) area deprivation index (ADI) and (D) health insurance type. Base maps were created using the Stamen Design from Stadia Maps. ## Discussion Our analysis of individual-level risk factors found that Black race and ICS prescription had the strongest positive associations with asthma exacerbations, as determined by individual-level logistic regression and spatial GAM models. In our cohort of Penn Medicine asthma patients, 66% with no exacerbations were Black compared to 83% with 5+ exacerbations, consistent with known racial disparities in asthma and observations in past Penn Medicine cohorts [44,45]. With regard to the observed association between exacerbations and ICS prescription, international asthma management guidelines underwent a major shift in 2019, recommending ICS as part of the first-line treatment for all asthma patients [46]. This shift was not reflected in our cohort, with only 76.0% of patients prescribed ICS. We observed exacerbations in patients in our cohort regardless of whether they had an ICS prescription, however, the strong association between ICS and exacerbations suggests that patients with more severe asthma were prescribed ICS more frequently than those with milder asthma. Our sensitivity analysis found that logistic regression identified positive associations between age 35-54 and NO2 exposure that were not observed in a negative binomial regression model, suggesting that these factors were associated with risk of having at least one exacerbation compared to none, but not with having a higher count of exacerbations. Our spatial analysis of patient-level data revealed several important insights. First, we observed that asthma exacerbation risk across Philadelphia was spatially correlated and identified a hotspot with 41% higher odds of exacerbation compared to the median across the study region. This finding is consistent with the results of community-based pediatric asthma screening in Philadelphia, which have found that local asthma prevalence can vary significantly from regional or national estimates [47]. Our findings are also consistent with studies that have assessed spatial heterogeneity of pediatric and adult asthma data in other United States metropolitan areas, albeit with less granular spatial resolution. Zárate et al. observed statistically significant spatial patterning of asthma-related emergency department visits across census tracts in Central Texas [48], Grunwell et al. identified a group of contiguous census tracts in the State of Georgia with high rates of admission to the pediatric intensive care unit for asthma [49], and Harris et al. and Corburn et al. identified a statistically significant cluster of zip codes in St. Louis, Missouri and of census tracts in New York City, respectively, with elevated pediatric asthma hospitalization rates [50,51]. Our observation of a West and South Philadelphia hotspot of exacerbation risk is also consistent with analyses of past Penn Medicine cohorts [3,4] that focused on validating methods for augmenting EHR datasets rather than identifying factors associated with the hotspots. The factors we found to be associated with asthma exacerbations according to EHR- and EHR & SEDH-adjusted spatial GAM models were largely consistent with individual-level logistic regression findings, but the spatial analysis provided an improved context to understand the individual-level results. Adjusting for all EHR and SEDH variables decreased the variance in spatial ORs by 66.9%, indicating that these variables together accounted for a large proportion of the spatial variance in exacerbation odds. By adjusting the spatial model one variable at a time, we found that ADI, race, and health insurance type most attenuated the hotspot area and effect size by reducing the variance of spatial ORs, suggesting that these variables were the most influential in determining the spatial distribution of exacerbation risk. Our findings are consistent with known asthma disparities by race/ethnicity and socioeconomic status [7–11], as well as past neighborhood-level analyses of pediatric asthma: Harris et al. and Corburn et al. found that asthma hospitalization hotspots in St. Louis and New York City had greater proportions of non-White residents and greater rates of poverty, unemployment, high-density housing, and lack of access to a household vehicle, although they did not test for statistical significance [50,51]. Grunwell et al. found a statistically significant difference between hotspot and non-hotspot census tracts in Georgia for poverty, unemployment, and high-density housing, but not for race [49]. Nearly all patient characteristics in our cohort, including health insurance type and ADI, differed significantly between Black and White patients (Table 4). Most notably, 8.5% of White patients had Medicaid health insurance compared to 39% of Black patients, and the median ADI for White patients was 2.60 (IQR 1.50-4.30) compared to 8.40 for Black patients (IQR 6.50-9.40), making it difficult to assess confounding in our spatial models. Although the VIF indicated that all variables could be included in a multivariable model without substantially inflating variance, collinearity between race and ADI as well as race and health insurance type (S4 Fig) may help explain why ADI and health insurance type were statistically significant in bivariable spatial models but not in the EHR & SEDH-adjusted spatial GAM nor in the EHR & SEDH-adjusted logistic regression (Tables 2 and 3). We attempted to overcome some of these limitations by stratifying our analysis by race. We found that, unlike in the full cohort, adjusting for ADI and health insurance type for Black patients did not attenuate the variance in spatial ORs, suggesting that in our cohort the association between these variables and asthma exacerbations was confounded by race. Our results are consistent with previous observations that racial disparities in asthma control persist even after accounting socioeconomic status [52], but that it is difficult to separate out the effects of socioeconomic status from the effects of race [53]. Confounding of the asthma-socioeconomic status relationship by race has also been observed in past neighborhood-level analyses of asthma morbidity. Zárate et al. found that spatial patterning of asthma-related emergency department visits in Central Texas was partially explained by socioeconomic characteristics in White patients, but not in Black or Hispanic patients [48]. More broadly, understanding the relative contributions of many social determinants of health to asthma is made difficult by their unequal distribution across racial/ethnic groups in the United States [10]. Our health insurance type variable serves as a proxy for individual-level socioeconomic status, and the relationship we measured between it and asthma exacerbation risk is potentially mediated by several pathways including increased rates of smoking [54], psychosocial stress [55], and obesity [56], or an unmeasured variable that varies with socioeconomic status and is also associated with minoritized racial and ethnicity groups in the United States [57]. On the other hand, low neighborhood-level socioeconomic status, which we measured with ADI, is associated with differential exposure to indoor and outdoor air pollution, psychosocial stress from neighborhood violence, and community norms surrounding health behaviors, all of which have been linked to asthma exacerbations [27]. Due to geographic segregation by race that resulted from structural racism and is common in many US cities, including in Philadelphia as we observed in our study, race and neighborhood-level SES are also highly correlated [58]. Thus, our inability to identify an association between asthma exacerbations and ADI when restricting our analysis to Black patients may be due to a restriction of the range of people across the ADI spectrum relative to the range observed in all patients [48]. Future work is needed to understand what specific SEDH variables that covary with race, ADI, and health insurance type are the primary drivers of local disparities in asthma exacerbation risk across Philadelphia. Our results demonstrate that integration of diverse SEDH datasets with the EHR and the use of both spatial and non-spatial modeling approaches are helpful to understand factors contributing to complex health conditions in real world populations. In both our logistic regression models and spatial GAMs, model fit was not improved by adjusting for SEDH variables in addition to EHR-derived variables. However, in our spatial models, adjusting for SEDH variables resulted in twice as much reduction of the initial variance in spatial ORs compared to adjusting for EHR variables only. Our non-spatial and spatial models also identified different sets of factors associated with exacerbations. For example, ICS prescription was found to have a strong positive association with exacerbations in logistic regression models but did not contribute to the spatial distribution of ORs; conversely, ADI attenuated the variance in spatial ORs more strongly than all other variables tested but was not significant in either logistic regression or negative binomial models. These findings present a framework for future efforts to expand the scope of EHR data, which, especially as the spatial resolution of SEDH datasets continues to increase, will allow for improved individualized exposure estimates. In the future, integration of SEDH data into the EHR may be helpful to tailor asthma management strategies and for health systems to create population-level interventions to improve health of their patients. This study is strengthened by the high spatial resolution of both our SEDH data and our analysis. Integrating the highest resolution SEDH data available during the time of the study period allowed us to most closely approximate individual-level exposures, and analyzing EHR data at a fine resolution allowed us to understand local health patterns that may not be visible at the census tract or zip-code level. Additional strengths included accounting for many variables and increasing the likelihood that Penn Medicine was patients’ primary care provider by applying SRR restriction. Our study is also subject to limitations, including some related to use of EHR data, such as missingness, entry error, and phenotype misclassification. Additionally, the geocoded addresses used in our study reflect residence at the time of the data pull, but they do not account for residential mobility during the study period, nor does residence information provide a full assessment of environmental exposures. ## Conclusion By integrating seven datasets containing information on SEDH with an EHR dataset to create individualized exposure assessments, we identified non-spatial and spatial factors associated with asthma exacerbations. Race and prescription of an ICS were most strongly associated with exacerbations in individual-level models. Race also accounted for the most spatial variation in exacerbation odds, along with ADI and health insurance type. Because these three variables had similar spatial distributions, understanding which contributes most to disparities in asthma exacerbations requires additional study of people living in the region identified as a hotspot. Our findings demonstrate how integrating diverse data types and geospatial modeling approaches with EHR data are helpful to understand complex diseases locally. ## Data Availability The social and environmental determinants of health datasets that support the findings of this study are publicly available in Sensor-based Analysis of Pollution in the Philadelphia Region with Information on Neighborhoods and the Environment (SAPPHIRINE), offered as a web application ([http://sapphirine.org](http://sapphirine.org)) and R package ([https://github.com/HimesGroup/sapphirine](https://github.com/HimesGroup/sapphirine)). Based on ethical and legal considerations, such as that the data was not collected with informed consent, the electronic health record (EHR) data used in this study cannot be shared widely. To request access to EHR data at the Penn Medicine hospital system, reach out to the PennDnA office ([https://www.med.upenn.edu/penndna/](https://www.med.upenn.edu/penndna/)). The EHR data cleaning methodology is provided as supplementary files to support reproducibility despite being unable to share the data itself. [http://sapphirine.org](http://sapphirine.org) ## Supporting information **S1 Text. Supplementary Methods.** **S1 Figure. Flowchart of patient cohort selection.** Overview of steps followed to select final patient cohort (N = 6,656) from EHR data on all Penn Medicine patients with at least one asthma ICD code (N = 86,787). **S2 Figure. Selection of study region using the spatial representation ratio (SRR).** (A) SRR values, defined as the cohort population residing in a block group divided by the underlying population as reported by the 2019 American Community Survey, for all of Philadelphia and for the selected study region (inset box). SRR = 1 indicates no representativeness bias. Density plots of the study cohort (B) before and (C) after filtering for the study region. Base maps were created using the Stamen Design from Stadia Maps. **S3 Figure. Spatial distribution of SEDH datasets that were integrated with EHR data.** The following maps are shown for the spatial area which comprised our study region: (A) raster of NO2 pollution levels, (B) raster of PM2.5 pollution levels, (C) point sites of toxic releases and the total summed emissions at each site, (D) line segments of roadways and the daily vehicle miles traveled (DVDT) on each, (E) housing violations per block group, normalized by the underlying 2019 American Community Survey population, (F) raster of the normalized difference vegetation index (NDVI). A map of area deprivation index (ADI), which most strongly reduced the spatial variance of odds of exacerbation risk in our spatial GAMs, is shown in Figure 3A. Base maps were created using the Stamen Design from Stadia Maps. **S4 Figure. Pairwise correlation between all EHR and SEDH variables.** Measures in each box correspond to Pearson’s correlation coefficients. For nominal categorical variables, reference levels are as follows: White (race), Private (health insurance type). **S5 Figure. Influence of individual EHR and SEDH variables on the odds ratio (OR) of the unadjusted spatial GAM.** Percent reduction in the variance of ORs across the study region for one-variable-at-a-time adjusted spatial GAMs compared to the unadjusted model. OR changes for models adjusted one variable at a time (in addition to years followed) are shown in blue. Multivariable models (both EHR-adjusted and EHR & SEDH-adjusted) are shown in red for comparison. **S6 Figure. Spatial GAMs adjusted one-at-a-time for the EHR-derived variables that did not greatly reduce variance.** Spatial odds ratios (ORs) of exacerbation are shown after adjusting for years followed and one-at-a-time for the following variables whose percent reduction in variance of ORs was less than 25: (A) age, (B) sex, (C) ethnicity, (D) BMI, (E) smoking status, (F) COPD, (G) allergic rhinitis, (H) Elixhauser comorbidity score, (I) ICS. Base maps were created using the Stamen Design from Stadia Maps. **S7 Figure. Spatial GAMs adjusted one-at-a-time for the SEDH variables that did not greatly reduce variance.** Spatial odds ratios (ORs) of exacerbation are shown after adjusting for years followed and one-at-a-time for the following variables whose reduction in variance of ORs was less than 25: (A) NO2, (B) PM2.5, (C) toxic releases exposure, (D) vehicular traffic, (E) housing violations, (F) normalized difference vegetation index (NDVI). Base maps were created using the Stamen Design from Stadia Maps. **S8 Figure. Spatial odds ratios (ORs) of exacerbations among White patients along with the effects of ADI and health insurance type on this distribution.** (A) SRR values for the updated study region used in spatial GAMs for patients of White race only (SRR = 1 indicates no representativeness bias). (B) Unadjusted spatial GAM adjusted only for years followed for patients of White race (N = 1,383). Spatial GAMs adjusted additionally for (C) area deprivation index (ADI) and (D) health insurance type. Base maps were created using the Stamen Design from Stadia Maps. **S1 Table**. **Generic medication names included in medication classes.** The following generic drug names recorded in the EHR during the study period were used for asthma and exacerbation phenotyping as well as used as independent variables in select models (i.e., ICS). Instances in which these drugs were listed as investigational or nasal formulations were not included. **S2 Table. Sources and spatiotemporal dimensions of geospatial datasets merged with EHR data.** **S3 Table**. **Asthma-related housing code violations extracted from the Philadelphia Department of Licenses and Inspections.** **S4 Table**. **Patient medications by exacerbation count levels.** Shown are the number and percentage of patients receiving each of the medication types listed according to their number of exacerbations during the study period. **S5 Table. Characteristics of complete cases and patients excluded due to missingness.** Shown are the number and percentage of patients in each level for categorical variables and the Median and Interquartile Range (IQR) for the Years followed variable in complete cases versus those excluded due to missingness in the sex, ethnicity, health insurance type, BMI, and smoking status variables. **S6 Table. Adjusted Generalized Variance Inflation Factors (GVIFs) for each EHR and SEDH variable included in the EHR & SEDH-adjusted negative binomial and logistic regression models.** **S7 Table**. **Individual-level asthma exacerbation risk factors in multivariable negative binomial regression models.** Shown are the adjusted incidence rate ratios (IRRs), 95% confidence intervals (CIs), and p-values for negative binomial models of asthma exacerbations as a count outcome adjusted for EHR-derived variables only and for both EHR-derived and SEDH variables. **S8 Table. Spatial GAMs of asthma exacerbations adjusted for individual risk factors that most changed risk**. Shown are the adjusted odds ratios (ORs), 95% confidence intervals (CIs), and p-values for spatial GAMs of asthma exacerbations as a dichotomous outcome adjusted for years followed and one-at-a-time for ADI, race, and health insurance type, the three variables whose percent reduction in variance of ORs was greater than 25. ## Acknowledgements We would like to thank Sunil Thomas from the University of Pennsylvania Penn Data Store for extracting the EHR data used for this project. ## Footnotes * Updated formatting of some sections. * Received October 24, 2024. * Revision received November 1, 2024. * Accepted November 2, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Mooney SJ, Westreich DJ, El-Sayed AM. Commentary: Epidemiology in the era of big data. Epidemiology. 2015;c: 390–394. doi:10.1097/EDE.0000000000000274 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/EDE.0000000000000274&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25756221&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 2. 2.Cook LA, Sachs J, Weiskopf NG. The quality of social determinants data in the electronic health record: a systematic review. J Am Med Inform Assoc JAMIA. 2021;29: 187–196. doi:10.1093/jamia/ocab199 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jamia/ocab199&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34664641&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 3. 3.Xie S, Greenblatt R, Levy MZ, Himes BE. Enhancing Electronic Health Record Data with Geospatial Information. AMIA Summits Transl Sci Proc. 2017;2017: 123–132. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28815121&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 4. 4.Xie S, Himes BE. Approaches to link geospatially varying social, economic, and environmental factors with electronic health record data to better understand asthma exacerbations. AMIA Annu Symp Proc. 2018;2018: 1561–1570. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30815202&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 5. 5.Cui Y, Eccles KM, Kwok RK, Joubert BR, Messier KP, Balshaw DM. Integrating multiscale geospatial environmental data into large population health studies: challenges and opportunities. Toxics. 2022;10: 403. doi:10.3390/toxics10070403 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/toxics10070403&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35878308&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 6. 6.Schreibman A, Xie S, Hubbard RA, Himes BE. Linking ambient NO2 pollution measures with electronic health record data to study asthma exacerbations. AMIA Summits Transl Sci Proc. 2023;2023: 467–476. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37350870&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 7. 7.CDC. Asthma: most recent national asthma data. Atlanta, GA: US Department of Health and Human Services, CDC; 2020. Available: [https://www.cdc.gov/asthma/most\_recent\_national\_asthma\_data.htm](https://www.cdc.gov/asthma/most\_recent_national_asthma_data.htm) 8. 8.Leong AB, Ramsey CD, Celedón JC. The challenge of asthma in minority populations. Clin Rev Allergy Immunol. 2012;43: 156–183. doi:10.1007/s12016-011-8263-1 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12016-011-8263-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21538075&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 9. 9.Forno E, Celedón JC. Health disparities in asthma. Am J Respir Crit Care Med. 2012;185: 1033–1035. doi:10.1164/rccm.201202-0350ED [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.201202-0350ED&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22589306&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 10. 10.Grant T, Croce E, Matsui EC. Asthma and the social determinants of health. Ann Allergy Asthma Immunol Off Publ Am Coll Allergy Asthma Immunol. 2022;128: 5–11. doi:10.1016/j.anai.2021.10.002 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.anai.2021.10.002&link_type=DOI) 11. 11.Lovinsky-Desir S, Riley IL, Bryant-Stephens T, De Keyser H, Forno E, Kozik AJ, et al. Research priorities in pediatric asthma morbidity: addressing the impacts of systemic racism on children with asthma in the United States: an official American Thoracic Society workshop report. Ann Am Thorac Soc. 2024;21: 1349–1364. doi:10.1513/AnnalsATS.202407-767ST [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1513/AnnalsATS.202407-767ST&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=39352175&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 12. 12.Global Initiative for Asthma (GINA). Global strategy for asthma management and prevention (2024 update). Bethesda, MD: Global Initiative for Asthma (GINA); 2024. 13. 13.Fuhlbrigge A, Peden D, Apter AJ, Boushey HA, Camargo CA, Gern J, et al. Asthma outcomes: exacerbations. J Allergy Clin Immunol. 2012;129: S34–S48. doi:10.1016/j.jaci.2011.12.983 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2011.12.983&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22386508&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 14. 14.Krishnan V, Diette GB, Rand CS, Bilderback AL, Merriman B, Hansel NN, et al. Mortality in patients hospitalized for asthma exacerbations in the United States. Am J Respir Crit Care Med. 2006;174: 633–638. doi:10.1164/rccm.200601-007OC [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.200601-007OC&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16778163&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000240664100004&link_type=ISI) 15. 15.Ivanova JI, Bergman R, Birnbaum HG, Colice GL, Silverman RA, McLaurin K. Effect of asthma exacerbations on health care costs among asthmatic patients with moderate and severe persistent asthma. J Allergy Clin Immunol. 2012;129: 1229–1235. doi:10.1016/j.jaci.2012.01.039 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2012.01.039&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22326484&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 16. 16.Schatz M, Clark S, Camargo CA. Sex differences in the presentation and course of asthma hospitalizations. Chest. 2006;129: 50–55. doi:10.1378/chest.129.1.50 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1378/chest.129.1.50&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16424412&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 17. 17.Schatz M, Mosen DM, Kosinski M, Vollmer WM, Magid DJ, O’Connor E, et al. Predictors of asthma control in a random sample of asthmatic patients. J Asthma. 2007;44: 341–345. doi:10.1080/02770900701344421 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/02770900701344421&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17530535&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000247078700017&link_type=ISI) 18. 18.Pedersen SE, Bateman ED, Bousquet J, Busse WW, Yoxall S, Clark TJ. Determinants of response to fluticasone propionate and salmeterol/fluticasone propionate combination in the Gaining Optimal Asthma controL study. J Allergy Clin Immunol. 2007;120: 1036–1042. doi:10.1016/j.jaci.2007.07.016 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2007.07.016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17935765&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250973400009&link_type=ISI) 19. 19.Price D, Zhang Q, Kocevar VS, Yin DD, Thomas M. Effect of a concomitant diagnosis of allergic rhinitis on asthma-related health care use by adults. Clin Exp Allergy. 2005;35: 282–287. doi:10.1111/j.1365-2222.2005.02182.x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2222.2005.02182.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15784104&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000227837300007&link_type=ISI) 20. 20.Miller MK, Lee JH, Miller DP, Wenzel SE. Recent asthma exacerbations: A key predictor of future exacerbations. Respir Med. 2007;101: 481–489. doi:10.1016/j.rmed.2006.07.005 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.rmed.2006.07.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16914299&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000244829500014&link_type=ISI) 21. 21.Guarnieri M, Balmes JR. Outdoor air pollution and asthma. The Lancet. 2014;383: 1581–1592. doi:10.1016/S0140-6736(14)60617-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(14)60617-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24792855&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 22. 22.Weitzman M, Gortmaker S, Sobol A. Racial, social, and environmental risks for childhood asthma. Am J Dis Child. 1990;144: 1189–1194. doi:10.1001/archpedi.1990.02150350021016 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archpedi.1990.02150350021016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2239856&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 23. 23.Black PN, Udy AA, Brodie SM. Sensitivity to fungal allergens is a risk factor for life-threatening asthma. Allergy. 2000;55: 501–504. doi:10.1034/j.1398-9995.2000.00293.x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1034/j.1398-9995.2000.00293.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10843433&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000086710800013&link_type=ISI) 24. 24.Gelber LE, Seltzer LH, Bouzoukis JK, Pollart SM, Chapman MD, Platts-Mills TAE. Sensitization and exposure to indoor allergens as risk factors for asthma among patients presenting to hospital. Am Rev Respir Dis. 1993;147: 573–578. doi:10.1164/ajrccm/147.3.573 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/ajrccm/147.3.573&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8442589&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 25. 25.Lovasi GS, O’Neil-Dunne JPM, Lu JWT, Sheehan D, Perzanowski MS, MacFaden SW, et al. Urban tree canopy and asthma, wheeze, rhinitis, and allergic sensitization to tree pollen in a New York City birth cohort. Environ Health Perspect. 2013;121: 494–500. doi:10.1289/ehp.1205513 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1289/ehp.1205513&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23322788&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000323706100035&link_type=ISI) 26. 26.Lovasi GS, Quinn JW, Neckerman KM, Perzanowski MS, Rundle A. Children living in areas with more street trees have lower prevalence of asthma. J Epidemiol Community Health. 2008;62: 647–649. doi:10.1136/jech.2007.071894 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiamVjaCI7czo1OiJyZXNpZCI7czo4OiI2Mi83LzY0NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzExLzAyLzIwMjQuMTAuMjQuMjQzMTYwNjMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 27. 27.Gold DR, Wright R. Population disparities in asthma. Annu Rev Public Health. 2005;26: 89–113. doi:10.1146/annurev.publhealth.26.021304.144528 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.publhealth.26.021304.144528&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15760282&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228981500005&link_type=ISI) 28. 28.Fecho K, Ahalt SC, Appold S, Arunachalam S, Pfaff E, Stillwell L, et al. Development and application of an open tool for sharing and analyzing integrated clinical and environmental exposures data: asthma use case. JMIR Form Res. 2022;6: e32357. doi:10.2196/32357 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2196/32357&link_type=DOI) 29. 29.Beck AF, Huang B, Ryan PH, Sandel MT, Chen C, Kahn RS. Areas with high rates of police-reported violent crime have higher rates of childhood asthma morbidity. J Pediatr. 2016;173: 175–182.e1. doi:10.1016/j.jpeds.2016.02.018 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpeds.2016.02.018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26960918&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 30. 30.Rasmussen SG, Ogburn EL, McCormack M, Casey JA, Bandeen-Roche K, Mercer DG, et al. Association between unconventional natural gas development in the Marcellus Shale and asthma exacerbations. JAMA Intern Med. 2016;176: 1334–1343. doi:10.1001/jamainternmed.2016.2436 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamainternmed.2016.2436&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27428612&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 31. 31.van Walraven C, Austin PC, Jennings A, Quan H, Forster AJ. A modification of the Elixhauser comorbidity measures into a point system for hospital death using administrative data. Med Care. 2009;47: 626. doi:10.1097/MLR.0b013e31819432e5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/MLR.0b013e31819432e5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19433995&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 32. 32.Christie C, Xie S, Diwadkar AR, Greenblatt RE, Rizaldi A, Himes BE. Consolidated Environmental and Social Data Facilitates Neighborhood-Level Health Studies in Philadelphia. AMIA Annu Symp Proc. 2022;2021: 305–313. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35308932&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 33. 33.Cooper MJ, Martin RV, McLinden CA, Brook JR. Inferring ground-level nitrogen dioxide concentrations at fine spatial resolution applied to the TROPOMI satellite instrument. Environ Res Lett. 2020;15: 104013. doi:10.1088/1748-9326/aba3a5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1088/1748-9326/aba3a5&link_type=DOI) 34. 34.van Donkelaar A, Hammer MS, Bindle L, Brauer M, Brook JR, Garay MJ, et al. Monthly global estimates of fine particulate matter and their uncertainty. Environ Sci Technol. 2021;55: 15287–15300. doi:10.1021/acs.est.1c05309 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/acs.est.1c05309&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34724610&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 35. 35.US EPA. TRI basic data files: calendar years 1987-present. 3 Mar 2013 [cited 25 Jul 2023]. Available: [https://www.epa.gov/toxics-release-inventory-tri-program/tri-basic-data-files-calendar-years-1987-present](https://www.epa.gov/toxics-release-inventory-tri-program/tri-basic-data-files-calendar-years-1987-present) 36. 36.RMSTRAFFIC (Traffic Volumes). 2023 [cited 25 Jul 2023]. Available: [https://data-pennshare.opendata.arcgis.com/datasets/rmstraffic-traffic-volumes/explore](https://data-pennshare.opendata.arcgis.com/datasets/rmstraffic-traffic-volumes/explore) 37. 37.Kind AJH, Buckingham WR. Making neighborhood-disadvantage metrics accessible — the neighborhood atlas. N Engl J Med. 2018;378: 2456–2458. doi:10.1056/NEJMp1802313 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMp1802313&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29949490&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 38. 38.University of Wisconsin School of Medicine Public Health. 2018 Area Deprivation Index v3.1. Available: [https://www.neighborhoodatlas.medicine.wisc.edu/](https://www.neighborhoodatlas.medicine.wisc.edu/). Accessed July 2022. 39. 39.Licenses and Inspections code violations. In: OpenDataPhilly [Internet]. [cited 25 Jul 2023]. Available: [https://opendataphilly.org/datasets/licenses-and-inspections-code-violations/](https://opendataphilly.org/datasets/licenses-and-inspections-code-violations/) 40. 40.Ermida SL, Soares P, Mantas V, Göttsche F-M, Trigo IF. Google Earth Engine open-source code for land surface temperature estimation from the Landsat series. Remote Sens. 2020;12: 1471. doi:10.3390/rs12091471 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/rs12091471&link_type=DOI) 41. 41.R Core Team (2020). R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; Available: [https://www.R-project.org/](https://www.R-project.org/) 42. 42.Xie SJ, Kapos FP, Mooney SJ, Mooney S, Stephens KA, Chen C, et al. Geospatial divide in real-world EHR data: analytical workflow to assess regional biases and potential impact on health equity. AMIA Summits Transl Sci Proc. 2023;2023: 572–581. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37350875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 43. 43.Bai L, Bartell S, Bliss R, Vieira V. Mapping smoothed effect estimates from individual-level spatial data. 2022. Available: [https://search.r-project.org/CRAN/refmans/MapGAM/html/MapGAM-package.html](https://search.r-project.org/CRAN/refmans/MapGAM/html/MapGAM-package.html) 44. 44.Moorman JE, Akinbami LJ, Bailey CM, Zahran HS, King ME, Johnson CA, et al. National surveillance of asthma: United States, 2001-2010. Vital Health Stat 3. 2012; 1–58. 45. 45.Greenblatt RE, Zhao EJ, Henrickson SE, Apter AJ, Hubbard RA, Himes BE. Factors associated with exacerbations among adults with asthma according to electronic health record data. Asthma Res Pract. 2019;5: 1. doi:10.1186/s40733-019-0048-y [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40733-019-0048-y&link_type=DOI) 46. 46.Global Initiative for Asthma (GINA). Global strategy for asthma management and prevention (2019 update). Bethesda, MD: Global Initiative for Asthma (GINA); 2019. 47. 47.Bryant-Stephens T, West C, Dirl C, Banks T, Briggs V, Rosenthal M. Asthma prevalence in philadelphia: description of two community-based methodologies to assess asthma prevalence in an inner-city population. J Asthma. 2012;49: 581–585. doi:10.3109/02770903.2012.690476 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3109/02770903.2012.690476&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22747076&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 48. 48.Zárate Rebecca A, Bhavnani D, Chambliss S, Hall EM, Zigler C, Cubbin C, et al. Neighborhood-level variability in asthma-related emergency department visits in Central Texas. J Allergy Clin Immunol. 2024;154: 933–939. doi:10.1016/j.jaci.2024.05.024 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2024.05.024&link_type=DOI) 49. 49.Grunwell JR, Opolka C, Mason C, Fitzpatrick AM. Geospatial analysis of social determinants of health identifies neighborhood hot spots associated with pediatric intensive care use for life-threatening asthma. J Allergy Clin Immunol Pract. 2022;10: 981–991.e1. doi:10.1016/j.jaip.2021.10.065 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaip.2021.10.065&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34775118&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 50. 50.Harris KM. Mapping inequality: Childhood asthma and environmental injustice, a case study of St. Louis, Missouri. Soc Sci Med. 2019;230: 91–110. doi:10.1016/j.socscimed.2019.03.040 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.socscimed.2019.03.040&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30999144&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) 51. 51.Corburn J, Osleeb J, Porter M. Urban asthma and the neighbourhood environment in New York City. Health Place. 2006;12: 167–179. doi:10.1016/j.healthplace.2004.11.002 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.healthplace.2004.11.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16338632&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234734900005&link_type=ISI) 52. 52.Haselkorn T, Lee JH, Mink DR, Weiss ST. Racial disparities in asthma-related health outcomes in severe or difficult-to-treat asthma. Ann Allergy Asthma Immunol. 2008;101: 256–263. doi:10.1016/S1081-1206(10)60490-5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1081-1206(10)60490-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18814448&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259056300007&link_type=ISI) 53. 53.Boudreaux ED, Emond SD, Clark S, Camargo CA. Acute asthma among adults presenting to the emergency department: the role of race/ethnicity and socioeconomic status. Chest. 2003;124: 803–812. doi:10.1378/chest.124.3.803 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1378/chest.124.3.803&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12970001&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185344100009&link_type=ISI) 54. 54.Hiscock R, Bauld L, Amos A, Fidler JA, Munafò M. Socioeconomic status and smoking: a review. Ann N Y Acad Sci. 2012;1248: 107–123. doi:10.1111/j.1749-6632.2011.06202.x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1749-6632.2011.06202.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22092035&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305677500006&link_type=ISI) 55. 55.Matthews KA, Gallo LC, Taylor SE. Are psychosocial factors mediators of socioeconomic status and health connections? Ann N Y Acad Sci. 2010;1186: 146–173. doi:10.1111/j.1749-6632.2009.05332.x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1749-6632.2009.05332.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20201872&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000277908000009&link_type=ISI) 56. 56.Wang Y, Beydoun MA. The obesity epidemic in the United States—gender, age, socioeconomic, racial/ethnic, and geographic characteristics: a systematic review and meta-regression analysis. Epidemiol Rev. 2007;29: 6–28. doi:10.1093/epirev/mxm007 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/epirev/mxm007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17510091&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248364800002&link_type=ISI) 57. 57.Williams DR. Race, socioeconomic status, and health: the added effects of racism and discrimination. Ann N Y Acad Sci. 1999;896: 173–188. doi:10.1111/j.1749-6632.1999.tb08114.x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1749-6632.1999.tb08114.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10681897&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F02%2F2024.10.24.24316063.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000085238100015&link_type=ISI) 58. 58.Charles CZ. The dynamics of racial residential segregation. Annu Rev Sociol. 2003;29: 167–207. doi:10.1146/annurev.soc.29.010202.100002 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.soc.29.010202.100002&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185094800008&link_type=ISI)