Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Associations between long-term fine particulate matter exposure and hospital procedures in heart failure patients

  • Samantha Catalano,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing

    Current address: Donald and Barbara Zucker School of Medicine at Hofstra/Northwell, Hempstead, New York, United States of America

    Affiliation Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

  • Joshua Moyer,

    Roles Data curation

    Affiliation Center for Public Health and Environmental Assessment, US Environmental Protection Agency, Chapel Hill, North Carolina, United States of America

  • Anne Weaver,

    Roles Data curation, Writing – review & editing

    Affiliation Center for Public Health and Environmental Assessment, US Environmental Protection Agency, Chapel Hill, North Carolina, United States of America

  • Qian Di,

    Roles Data curation

    Affiliation Research Center for Public Health, School of Medicine, Tsinghua University, Beijing, China

  • Joel D. Schwartz,

    Roles Data curation, Resources, Writing – review & editing

    Affiliation Harvard TH Chan School of Public Health, Boston, Massachusetts, United States of America

  • Michael Catalano,

    Roles Writing – review & editing

    Affiliation Division of Cardiovascular Surgery, Department of Surgery, Hospital of the University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

  • Cavin K. Ward-Caviness

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Supervision, Visualization, Writing – original draft, Writing – review & editing

    ward-caviness.cavin@epa.gov

    Affiliation Center for Public Health and Environmental Assessment, US Environmental Protection Agency, Chapel Hill, North Carolina, United States of America

Abstract

Background

Ambient fine particulate matter (PM2.5) contributes to global morbidity and mortality. One way to understand the health effects of PM2.5 is by examining its impact on performed hospital procedures, particularly among those with existing chronic disease. However, such studies are rare. Here, we investigated the associations between annual average PM2.5 and hospital procedures among individuals with heart failure.

Methods

Using electronic health records from the University of North Carolina Healthcare System, we created a retrospective cohort of 15,979 heart failure patients who had at least one of 53 common (frequency > 10%) procedures. We used daily modeled PM2.5 at 1x1 km resolution to estimate the annual average PM2.5 at the time of heart failure diagnosis. We used quasi-Poisson models to estimate associations between PM2.5 and the number of performed hospital procedures over the follow-up period (12/31/2016 or date of death) while adjusting for age at heart failure diagnosis, race, sex, year of visit, and socioeconomic status.

Results

A 1 μg/m3 increase in annual average PM2.5 was associated with increased glycosylated hemoglobin tests (10.8%; 95% confidence interval = 6.56%, 15.1%), prothrombin time tests (15.8%; 95% confidence interval = 9.07%, 22.9%), and stress tests (6.84%; 95% confidence interval = 3.65%, 10.1%). Results were stable under multiple sensitivity analyses.

Conclusions

These results suggest that long-term PM2.5 exposure is associated with an increased need for diagnostic testing on heart failure patients. Overall, these associations give a unique lens into patient morbidity and potential drivers of healthcare costs linked to PM2.5 exposure.

Introduction

Long-term exposure to ambient particulate matter smaller than 2.5 μm in diameter (PM2.5) has been shown to increase mortality and morbidity. PM2.5 can come from a variety of sources such as combustion, blown dust, industry activities, and traffic. PM2.5 is also formed via secondary chemical reactions of primary emitted pollutants [1]. Despite global efforts to reduce particulate matter concentrations, ambient air pollution remains a substantial contributor to the global burden of disease [2]. There are several proposed underlying mechanisms through which PM2.5 causes adverse health effects. The larger sized fractions of PM2.5 are not readily transported into the circulatory system and are instead most often linked to the triggering of systemic inflammation from the lungs and other intermediate pathways, e.g. lipid oxidation, that can cause extra-pulmonary effects. For example, at the molecular level, PM2.5 inhalation can increase the intracellular levels of reactive oxygen species, which causes pulmonary oxidative stress; this can instigate systemic inflammatory responses at the physiological level [3, 4]. The smallest sized fractions of PM2.5, often called ultrafine particulates, have shown an ability to translocate out of the lungs where they may exert direct toxicity on the heart, kidneys, and other organ systems [5].

These contributions of PM2.5 to global morbidity and mortality are exacerbated by aging populations and an increasing prevalence of chronic conditions such as heart failure (HF) [6]. It is estimated that by 2030, 1 in every 33 people in the United States will suffer from HF [7]. Moreover, the total direct medical costs of HF are projected to increase from $21 billion in 2012 to $53 billion in 2030 [7].

Hospital procedures—which we defined as any billable task or service a medical practitioner may provide to a patient—contribute to healthcare costs. Further, an increased number of non-elective hospital procedures may indicate worsening health in the individual and population. Long-term exposure to particulate matter is associated with multiple health outcomes, but it has never been investigated in relation to the performance of hospital procedures on a patient population. Quantifying the relationship between elevated concentrations of ambient particulate matter and hospital procedures would provide further evidence of the impact of poor air quality on individual health. Since hospital procedures are an indicator of patient morbidity and a driver of medical costs, and air pollution has been found to have a close temporal association with hospitalizations, studies investigating the associations between air pollution and hospital procedures may also present a new lens into understanding air pollution-associated health effects and the impacts of environmental exposures on the economics of the healthcare industry [8].

Here we estimate the associations between the frequencies of specific hospital procedures performed and the annual-average PM2.5 exposures among HF patients. Long-term PM2.5 exposure has been associated with increased mortality, hospitalizations, and readmissions among HF patients [9, 10], but this will be the first study to examine the impact of air pollution on performed hospital procedures in this understudied patient population. We hypothesize that long-term exposure to air pollution is associated with increased performance of hospital procedures among HF patients.

Methods

Study cohort

This study used the US Environmental Protection Agency (EPA) Clinical and Archived records Research for Environmental Studies (CARES) resource, which has been used in previous studies of environmental health among HF patients [9, 10]. EPA CARES is a resource that merges electronic health records from the University of North Carolina Healthcare System (UNCHCS) with environmental exposure data. This study was approved by the institutional review board of the University of North Carolina–Chapel Hill (IRB 17–0150), and informed consent was waived for this analysis of existing health records. For this study, we used electronic health records from patients with a diagnosis of HF between July 1, 2004 and December 31, 2016, which has been used in previous studies of mortality and hospital admissions [9, 10]. Patients for this HF cohort had to be aged 20 or older at the time of first HF diagnosis (to remove cases likely due to congenital disease) and had to have a valid address listed in their electronic health records for their date of HF diagnosis (to determine environmental exposures) [9]. We further required that individuals have at least one hospital visit prior to their HF diagnosis to limit the impact of those who were diagnosed at another hospital system and later seen within a UNCHCS hospital or clinic [10]. HF was defined according to the International Classification of Diseases, Ninth Revision (ICD‐9) codes 428.x and the International Classification of Diseases, Tenth Revision (ICD‐10) codes I50.x [9]. Electronic health records for this population included demographics, clinical variables such as address history, dates of procedures, and dates of diagnoses (via ICD-9 and ICD-10 codes). For this study environmental exposures were linked to individuals based on their date of HF diagnosis and the associated address from their electronic health record allowing us to see how environmental exposures at this critical time period impacted long-term health as reflected in hospital procedures. Census block group-level socioeconomic data from the 2010 census—chosen as it was the approximate midpoint of the study observation period—was used to evaluate area-level socioeconomic status based on geocoded primary address, as information on individual-level socioeconomic status is missing from the electronic health record. Individuals were considered successfully geocoded with a minimum of zip code level geocoding.

Individuals who did not receive any of the common (defined as >10% in frequency, as detailed in the Procedures section of the Methods below) procedures over their follow-up time may have represented a distinct set of (potentially healthier) patients or those who received healthcare at multiple hospital systems not in communication with the UNCHCS and, thus, they were removed from primary analyses—though we conducted sensitivity analyses to examine the impact of this choice.

Air pollution data

PM2.5 concentrations were estimated at 1x1 km resolution from 2000–2016 using a hybrid model that incorporated aerosol optical depth, ground-based monitoring, chemical transport models, land use, and meteorology [11]. Study participants were linked to PM2.5 data based on their geocoded primary residence at the date of HF diagnosis. Daily PM2.5 was averaged to create an annual average estimate based on the 365 days of exposure prior to their HF diagnosis.

Procedures

Out of all the observed procedures, we removed any “procedure” that was simply a hospital visit or a follow-up of a previous procedure. The data was further limited to only those procedures that were seen in at least 10% of the population to allow for a sufficient sample size for primary and subgroup analyses.

To eliminate potential duplicate procedure entries, we aggregated procedures to the day of their occurrence—effectively changing them into an indicator of whether the procedure was performed on that patient on a particular day—and merged procedures that changed names over the study period using the Current Procedure Terminology (CPT) code, which was static in our dataset over the study timeframe. We observed that tests classified as stress tests, including cardiovascular stress tests with and without report, technetium tc-99m sestamibi, and regadenoson injections, generally co-occurred. We thus combined these stress tests into one procedure category termed “stress tests” and treated them as a single procedure for the analyses. After all data merging and cleaning, there were 53 distinct procedures (S1 Table in S1 File). These procedures were primarily routine and diagnostic, likely due to the decision to only include procedures seen in at least 10% of the population.

Statistical analysis

We used a quasi-Poisson regression model (which allows for over or under dispersion) to model the association between annual average PM2.5 exposure at the time of HF diagnosis and the procedure count for each of the 53 procedures over the follow-up observation period (12/31/2016 or date of death) [12]. Confounders were selected based on previous studies done in this population and included individual level confounders as well as area-level (census block group) confounders to control for urbanicity and area-level socioeconomic status [9, 10]. Individual-level confounders were age at HF diagnosis; race; sex; and smoking status. Census block group confounders were based on the 2010 census and were median income; median household value; year of HF diagnosis; and the percentage of the census block that is urban, on public assistance, or below the poverty line. The screening mammography was the only sex-specific procedure and thus not adjusted for sex. The primary exposure for all models was the 365-day average (annual average) PM2.5 based on the day of HF diagnosis, and the log of the follow-up time was included as an offset in the models. The outcome was the total number of days on which a given procedure was performed for each patient over the follow-up time.

We further examined the robustness of the associations to confounder inclusion by additionally adjusting for access to healthcare and healthy foods, both of which may be related to socioeconomic status and potentially confound responses to air pollution. To perform this adjustment, we utilized county-level data on the number of primary care physicians, dentists, mental health providers, and non-primary care physicians per 100,000 people, along with average annual healthcare cost and food environment index from the 2015 County Health Rankings [13]. The food environment index is a scaled index ranging from 0 (worst) to 10 (best) that incorporates both limited access to healthy foods and food insecurity using data from the United States Department of Agriculture [14, 15].

We then performed a series of sensitivity analyses to understand the robustness of the significant associations from our primary analyses. We examined a Poisson mixed-effects model with a random intercept for residential zip code. The model was adjusted for the same variables as the original quasi-Poisson regression (including the offset term). To improve convergence, continuous variables were standardized to a mean of zero and a standard deviation of one.

Excluding patients who did not receive any procedure may have introduced bias if the patients who were excluded substantially differed from those who were included. We examined this potential source of bias by including all patients and using a zero-inflated Poisson model to account for the number of patients who did not receive procedures. The confounder adjustment for the zero-inflated Poisson model remained the same as the primary analysis.

To examine if results differed among those exposed to PM2.5 less than the current national ambient air quality standard, we restricted our study cohort to those with an annual average PM2.5 < 12 μg/m3 and reran our primary analysis on this restricted cohort [16].

We also examined population strata based on age, sex, race, urbanicity, and median household income to observe if there were strata-specific associations. All stratified models were run only on those procedures identified as Bonferroni significant in the primary analyses. The subgroups for age included less than 50 years, 50 to 65 years, and greater than or equal to 65 years. Racial subgroups were defined as the Black and White patients, as there were not enough in the “other” racial category for analyses. Urbanicity was determined using the percentage of the census block group based on the 2010 census; urban areas were defined as 100% urban (which captured slightly more than the top third of the distribution), rural areas were defined as 0% urban (which captured slightly less than the bottom third of the distribution), and all other areas were defined as 1%-99% urban, which we call suburban (though many are quite rural or urban). The subgroups for median household income were above and below the median income of $49,355.

To acknowledge the competing risk of death, we performed a sensitivity analysis where observations were weighted using stabilized inverse probability weights for the risk of death, as done in similar studies [10, 17]. The confounder adjustment remained the same as for the primary analysis.

All analyses were run in R version 4.0.2 [18]. Associations were converted to the percent change in outcome by exponentiating the regression coefficient, subtracting one, and multiplying the result by 100. All results are given as the percent change and associated 95% confidence interval (CI) per 1 μg/m3 increase in PM2.5. Statistical significance for the associations was set at p < 0.05/53, reflecting a Bonferroni multiple test adjustment for the number of procedures examined.

Results

The final study cohort consisted of 15,979 patients with an average follow-up time of 2.94 years (Table 1, S1 Fig in S1 File). Procedure frequencies per patient were consistent across the major demographic categories of age, race, and sex (Table 2).

We also examined the 4,941 patients out of the original patient cohort who did not receive any of the 53 procedures (S2 Table in S1 File). This patient population was slightly older, less sick in terms of co-existing conditions, and consisted of fewer smokers and former smokers. While no major differences in their demographics were observed, we cannot rule out unobserved systematic differences between these populations.

In our primary analyses three procedures were significant after multiple test correction: glycosylated hemoglobin tests (10.8% increase; 95% CI = 6.56%, 15.1%); prothrombin time tests (15.8% increase; 95% CI = 9.07%, 22.9%); and stress tests (6.84% increase; 95% CI = 3.65%, 10.1%) (Fig 1, Table 3, S3 Table in S1 File). Associations for all three Bonferroni significant procedures were attenuated after further adjustment for county-level access to healthcare and healthy foods, though associations with prothrombin time tests remained significant (1.46% increase; 95% CI = 0.71%, 2.22%; S4 Table in S1 File).

thumbnail
Fig 1. Associations between annual average PM2.5 and hospital procedures among heart failure patients.

Shown here are associations between annual average PM2.5 and the three Bonferroni significant procedures. Significance was determined according to the primary analysis approach (quasi-Poisson). We also show the percent change and 95% confidence interval associated with the sensitivity analysis for the Poisson Mixed-Effects model (which included a random intercept for zip code) and Zero-Inflated Poisson model (which included all individuals as opposed to just those with at least 1 procedure). See Methods for complete description of all analysis approaches. Percent change and 95% confidence interval presented per 1 μg/m3 increase in PM2.5. Gly Hemoglobin = Glycosylated Hemoglobin Test; Prothrombin = Prothrombin Time Test.

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

As there may be clustering of environmental and sociodemographic characteristics by zip code, we ran a Poisson mixed-effects model with a random intercept for zip code. Associations with the three significant procedures remained unchanged (Fig 1). Associations also remained similar when the patients who did not undergo each of the 53 procedures examined were added into the analysis cohort (Fig 1, S4 Table in S1 File). Given the strong mortality risk among HF patients, the risk of death may present an important competing risk to account for. Using inverse probability weighting, as described in the Methods, we adjusted for the competing risk of death, and observed no differences in the associations as compared to the primary analyses described above (S2 Fig in S1 File).

To understand if associations differed in areas with lower ambient PM2.5 concentrations, we examined associations for all individuals whose annual average PM2.5 exposure was less than the current PM2.5 annual average standard (National Ambient Air Quality Standard) of 12 μg/m3. Among this population, associations between PM2.5 and prothrombin time tests remained the same (15.1% increase; 95% CI = 8.07%, 22.7%), while associations between PM2.5 and glycosylated hemoglobin tests were attenuated (5.98% increase; 95% CI = 1.33%, 10.8%). No association between annual average PM2.5 and stress tests were seen amongst those with exposures less than 12 μg/m3 (S4 Table in S1 File).

We examined potential strata-specific association differences based on age, sex, race, urbanicity, and median household income of the census block group for the three Bonferroni significant procedures. For glycosylated hemoglobin tests, associations were stronger among Black study participants than White participants—and associations were stronger in residents of the suburban and urban census block groups than in those of the rural census block group. However, the wide, overlapping confidence intervals of these associations preclude any strong statements on strata-specific associations (S3-S5 Figs in S1 File).

Discussion

We examined associations between annual average PM2.5 exposure and the occurrence of hospital procedures in HF patients and observed three statistically significant associations, all of which are related to diagnostic testing for cardiorespiratory health.

Exposure to particulate matter is associated with altered glucose metabolism and increased risk of diabetes [19]. We observed statistically significant associations between glycosylated hemoglobin (HbA1C) tests and annual average PM2.5 exposure. Glycosylated hemoglobin tests are reflective of a patient’s blood sugar concentrations throughout the previous 3 months, and they may be used to both screen for diabetes in high-risk individuals and monitor progression of diabetes [20]. Elevated ambient concentrations of PM2.5 may contribute to elevated blood glucose and prevalence of diabetes, and, thus, an increased need for glycosylated hemoglobin tests [21, 22].

The other two significant procedures—stress tests and prothrombin time tests—are related to cardiovascular testing. In this study, “stress tests” included cardiovascular stress tests with and without interpretation and report, technetium sestamibi radiological imaging, and regadenoson injections, thus encompassing both exercise stress tests and pharmacologic stress tests [23]. Exercise and pharmacologic stress tests can help physicians understand the cause and severity of a HF patient’s condition or evaluate and modify their treatment plans. They also are often used in the diagnosis, workup, and management planning of ischemic heart disease and valvular heart disease [24, 25]. Since PM2.5 exposure is a risk factor for cardiovascular morbidity and mortality, these stress tests may be more necessary for patients with increased morbidity due to higher long-term exposure to PM2.5 [26].

Prothrombin time tests measure how long blood takes to clot and are often used to help detect and diagnose a bleeding disorder or excessive clotting disorder. Hematological aberrations are associated with long-term exposure to PM2.5, the mechanism of which may be related to increased inflammatory cytokine levels, oxidative stress, platelet activation, stimulated coagulation pathway, and reduced fibrinolysis [27]. Thus, disrupted hemostasis due to long-term PM2.5 exposure may increase the need for prothrombin time tests.

Associations were robust to a range of modeling scenarios (S4 Table in S1 File). Associations with prothrombin time tests were unchanged when restricted to exposures less than the current national standards, while associations with glycosylated hemoglobin were approximately halved and associations with stress tests attenuated towards the null. As each test may point towards a different underlying biological mechanism (hemostasis, cardiovascular function, and glucose metabolism), these differences may give clues as to the different mechanisms of increased morbidity that may be seen among HF patients as the concentration (and potentially composition) of ambient PM2.5 changes. Importantly, all associations were substantially attenuated when additionally adjusting for county-level indicators of access to healthcare and healthy food options, suggesting that these variables may play a role in moderating associations with annual average PM2.5 (S4 Table in S1 File). Formal examinations of the social and built environment factors that may moderate these associations were outside the scope of this manuscript and would be best performed with more comprehensive and geospatially resolved maps of these factors.

The strengths of our approach include the examination of HF patients, a vulnerable patient population that is more likely to have severe health outcomes from environmental exposures. Despite this vulnerability, HF patients have been understudied with respect to their environmental health risks. Another strength of this study was the large sample size and up to 12 years of follow-up. This study also used high spatial resolution (1x1 km) of the air pollution exposure data and patient address data, including street level geocoding for most patients and no coarser than zip code level geocoding.

Our approach also had a variety of limitations. One limitation is that the electronic health records do not include all possible confounders; in particular, individual-level socioeconomic status was missing. We used data from the US census to adjust for area-level socioeconomic status, which is a widely utilized approach when individual-level data is missing. As we only used records from UNCHCS, procedures that were not performed at an affiliated hospital or clinic were not captured in our data unless those records were transferred to the UNCHCS electronic health record system; further, these associations may not generalize other healthcare systems, such as those with a more rural catchment area. However, previous associations utilizing just the UNCHCS HF patient population have shown concordance with studies done utilizing the entire southeastern United States, so this analysis potentially has at least regional generalizability [10]. Finally, while we used air pollution exposure data with a high spatial resolution, the exposures were ultimately determined using the primary addresses of the patients. While this is standard in the field, reflecting that many people spend most of their time at their primary residence, it may not completely capture the PM2.5 exposure of individuals.

Moving forward, it will be important to investigate short-term exposures alongside long-term exposures. Future studies should also evaluate other patient populations and consider exposures beyond PM2.5 total mass, including the examination of different air pollutants such as ozone, ultrafine particulates/particle number count, and the chemical composition of PM2.5. Finally, this study did not have the scope to connect performed hospital procedures with the broader economics of hospital system operations. However, given that PM2.5 is estimated to cost the US economy billions of dollars annually, quantifying the economic impact of PM2.5 through the lens of performed hospital procedures might provide unique insights that can help quantify the economic impacts of PM2.5 exposure as well as the health, wellness, and morbidity impacts [28].

Conclusions

We observed that elevated long-term exposure to PM2.5 is associated with the increased performance of glycosylated hemoglobin tests, prothrombin time tests, and stress tests among HF patients. HF patients presented a unique patient group to examine this important question, as they have both increased environmental sensitivity and high rates of hospital utilization. This study adds to the field by quantifying PM2.5 impacts on performed hospital procedures, thereby allowing researchers to better understand air pollution health effects and helping to provide key information needed to estimate the total burden of PM2.5 on both patients and hospital systems.

Acknowledgments

This work does not necessarily represent the views or policy of the US EPA. Any mention of trade names does not constitute endorsement.

References

  1. 1. US Environmental Protection Agency. Integrated Science Assessment for Particulate Matter, EPA/600/R-19/188. US Environmental Protection Agency; 2020. Available from: https://cfpub.epa.gov/ncea/isa/recordisplay.cfm?deid=347534
  2. 2. Cohen AJ, Brauer M, Burnett R, Anderson HR, Frostad J, Estep K, et al. Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the Global Burden of Diseases Study 2015. The Lancet. 2017 May;389(10082):1907–191. pmid:28408086
  3. 3. Feng S, Gao D, Liao F, Zhou F, Wang X. The health effects of ambient PM2.5 and potential mechanisms. Ecotoxicology and Environmental Safety. 2016 Jun;128:67–74. pmid:26896893
  4. 4. Brook RD, Rajagopalan S, Pope CA, Brook JR, Bhatnagar A, Diez-Roux AV, et al. Particulate Matter Air Pollution and Cardiovascular Disease. Circulation. 2010 Jun;121(21):2331–78.
  5. 5. Miller MR, Raftis JB, Langrish JP, McLean SG, Samutrtai P, Connell SP, et al. Inhaled Nanoparticles Accumulate at Sites of Vascular Disease. ACS Nano. 2017 Apr 26;11(5):4542–52. pmid:28443337
  6. 6. Virani SS, Alonso A, Aparicio HJ, Benjamin EJ, Bittencourt MS, Callaway CW, et al. Heart Disease and Stroke Statistics—2021 Update: A Report From the American Heart Association. Circulation. 2021 Jan 27;143(8). Available from: https://www.ahajournals.org/doi/10.1161/CIR.0000000000000950 pmid:33501848
  7. 7. Heidenreich PA, Albert NM, Allen LA, Bluemke DA, Butler J, Fonarow GC, et al. Forecasting the Impact of Heart Failure in the United States: A Policy Statement From the American Heart Association. Circulation: Heart Failure. 2013 Apr 24;6(3):606–619. Available from: https://www.ahajournals.org/doi/pdf/10.1161/HHF.0b013e318291329a pmid:23616602
  8. 8. Shah AS, Langrish JP, Nair H, McAllister DA, Hunter AL, Donaldson K, et al. Global association of air pollution and heart failure: a systematic review and meta-analysis. The Lancet. 2013 Sep;382(9897):1039–1048. pmid:23849322
  9. 9. Ward‐Caviness CK, Weaver AM, Buranosky M, Pfaff ER, Neas LM, Devlin RB, et al. Associations Between Long‐Term Fine Particulate Matter Exposure and Mortality in Heart Failure Patients. Journal of the American Heart Association. 2020 Mar 17;9(6). pmid:32172639
  10. 10. Ward‐Caviness CK, Danesh Yazdi M, Moyer J, Weaver AM, Cascio WE, Di Q, et al. Long‐Term Exposure to Particulate Air Pollution Is Associated With 30‐Day Readmissions and Hospital Visits Among Patients With Heart Failure. Journal of the American Heart Association. 2021 May 4. pmid:33942627
  11. 11. Di Q, Amini H, Shi L, Kloog I, Silvern R, Kelly J, et al. An ensemble-based model of PM2.5 concentration across the contiguous United States with high spatiotemporal resolution. Environment International. 2019 Sep;130:104909. pmid:31272018
  12. 12. Regression—Quasi-Poisson Regression Q. Available from: https://wiki.q-researchsoftware.com/wiki/Regression_-_Quasi-Poisson_Regression
  13. 13. Remington PL, Catlin BB, Gennuso KP. The County Health Rankings: rationale and methods. Population Health Metrics. 2015 Apr 17;13(1). pmid:25931988
  14. 14. Berkowitz SA, Basu S, Meigs JB, Seligman HK. Food Insecurity and Health Care Expenditures in the United States, 2011–2013. Health Services Research. 2017 Jun 13;53(3):1600–1620. pmid:28608473
  15. 15. Ahern M, Brown C, Dukas S. A National Study of the Association Between Food Environments and County-Level Health Outcomes. The Journal of Rural Health. 2011 Apr 27;27(4):367–379. pmid:21967380
  16. 16. The National Ambient Air Quality Standards for Particulate Matter. EPA retains air quality standards for particle pollution (particulate matter): fact sheet. Available from: https://www.epa.gov/sites/default/files/2020-04/documents/fact_sheet_pm_naaqs_proposal.pdf
  17. 17. Danesh Yazdi M, Wang Y, Di Q, Zanobetti A, Schwartz J. Long-term exposure to PM2.5 and ozone and hospital admissions of Medicare participants in the Southeast USA. Environment International. 2019 Sep;130:104879. pmid:31238267
  18. 18. R Core Team. R: The R Project for Statistical Computing. R-project.org. 2022. Available from: https://www.r-project.org/
  19. 19. Hwang M-J, Kim J-H, Koo Y-S, Yun H-Y, Cheong H-K. Impacts of ambient air pollution on glucose metabolism in Korean adults: a Korea National Health and Nutrition Examination Survey study. Environmental Health. 2020 Jun 17;19(1). pmid:32552747
  20. 20. Glycosylated Hemoglobin Test. Winchester Hospital. Available from: https://www.winchesterhospital.org/health-library/article?id=36333
  21. 21. Bowe B, Xie Y, Li T, Yan Y, Xian H, Al-Aly Z. The 2016 global and national burden of diabetes mellitus attributable to PM 2.5 air pollution. The Lancet Planetary Health. 2018 Jul;2(7):e301–12. Available from: https://www.sciencedirect.com/science/article/pii/S2542519618301402?via%3Dihub
  22. 22. Peng C, Bind M-AC, Colicino E, Kloog I, Byun H-M, Cantone L, et al. Particulate Air Pollution and Fasting Blood Glucose in Nondiabetic Individuals: Associations and Epigenetic Mediation in the Normative Aging Study, 2000–2011. Environmental Health Perspectives. 2016 Nov;124(11):1715–1721. pmid:27219535
  23. 23. Ask the doctor: Exercise versus pharmacologic stress testing—Harvard Health. Harvard Health. 2014 Mar 1; Available from: https://www.health.harvard.edu/heart-health/exercise-versus-pharmacologic-stress-testing
  24. 24. Advanced Heart Failure Tests and Diagnosis. Emory Healthcare. Available from: https://www.emoryhealthcare.org/centers-programs/advanced-heart-failure-therapy-center/diagnosis.html
  25. 25. Stress Test. Mayo Clinic. 2018. Available from: https://www.mayoclinic.org/tests-procedures/stress-test/about/pac-20385234
  26. 26. Du Y, Xu X, Chu M, Guo Y, Wang J. Air particulate matter and cardiovascular disease: the epidemiological, biomedical and clinical evidence. Journal of Thoracic Disease. 2016 Jan;8(1):e8–19. pmid:26904258
  27. 27. Fongsodsri K, Chamnanchanunt S, Desakorn V, Thanachartwet V, Sahassananda D, Rojnuckarin P, et al. Particulate Matter 2.5 and Hematological Disorders From Dust to Diseases: A Systematic Review of Available Evidence. Frontiers in Medicine. 2021 Jul 14;8.
  28. 28. Heo J, Adams PJ, Gao HO. Public Health Costs of Primary PM2.5 and Inorganic PM2.5 Precursor Emissions in the United States. Environmental Science & Technology. 2016 May 17;50(11):6061–70. pmid:27153150
  29. 29. CPT Code Lookup, CPT® Codes and Search—Codify by AAPC. www.aapc.com. Available from: https://www.aapc.com/codes/cpt-codes-range/