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

Multilevel modeling of county-level excessive alcohol use, rurality, and COVID-19 case fatality rates in the US

  • George Pro ,

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

    gcpro@uams.edu

    Affiliation Department of Health Behavior and Health Education, Fay W. Boozman College of Public Health, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

  • Paul A. Gilbert,

    Roles Conceptualization, Methodology, Writing – review & editing

    Affiliation Department of Community and Behavioral Health, University of Iowa College of Public Health, Iowa City, Iowa, United States of America

  • Julie A. Baldwin,

    Roles Conceptualization, Methodology, Writing – review & editing

    Affiliation Center for Health Equity Research, Northern Arizona University, Flagstaff, Arizona, United States of America

  • Clare C. Brown,

    Roles Conceptualization, Methodology, Writing – review & editing

    Affiliation Department of Health Policy and Management, Fay W. Boozman College of Public Health, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

  • Sean Young,

    Roles Formal analysis, Methodology, Visualization, Writing – review & editing

    Affiliation Department of Environmental and Occupational Health, Fay W. Boozman College of Public Health, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

  • Nickolas Zaller

    Roles Conceptualization, Methodology, Supervision, Writing – review & editing

    Affiliation Department of Health Behavior and Health Education, Fay W. Boozman College of Public Health, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Abstract

Objective

Reports of disparities in COVID-19 mortality rates are emerging in the public health literature as the pandemic continues to unfold. Alcohol misuse varies across the US and is related to poorer health and comorbidities that likely affect the severity of COVID-19 infection. High levels of pre-pandemic alcohol misuse in some counties may have set the stage for worse COVID-19 outcomes. Furthermore, this relationship may depend on how rural a county is, as access to healthcare in rural communities has lagged behind more urban areas. The objective of this study was to test for associations between county-level COVID-19 mortality, pre-pandemic county-level excessive drinking, and county rurality.

Method

We used national COVID-19 data from the New York Times to calculate county-level case fatality rates (n = 3,039 counties and county equivalents; October 1 –December 31, 2020) and other external county-level data sources for indicators of rurality and health. We used beta regression to model case fatality rates, adjusted for several county-level population characteristics. We included a multilevel component to our model and defined state as a random intercept. Our focal predictor was a single variable representing nine possible combinations of low/mid/high alcohol misuse and low/mid/high rurality.

Results

The median county-level COVID-19 case fatality rate was 1.57%. Compared to counties with low alcohol misuse and low rurality (referent), counties with high levels of alcohol and mid (β = -0.17, p = 0.008) or high levels of rurality (β = -0.24, p<0.001) demonstrated significantly lower case fatality rates.

Conclusions

Our findings highlight the intersecting roles of county-level alcohol consumption, rurality, and COVID-19 mortality.

Introduction

COVID-19, the disease caused by the novel coronavirus SARS-CoV-2, is an ongoing global pandemic. At this writing (January, 2021), the US has seen more cases than any other country, with more than 20 million reported cases and surpassing 350,000 deaths [1]. Co-occurring and pre-existing health conditions exacerbate susceptibility to and the severity of COVID-19 infections [2]. While co-occurring health conditions have received considerable attention, particularly in relation to racial/ethnic disparities [3, 4], there has been less attention paid to behavioral and geographic risk factors associated with COVID-19 outcomes, such as alcohol use and rurality.

Levels of alcohol use and related problems vary across the US. The majority of Americans (73%) report drinking alcohol in the past year, and 13% meet the clinical criteria for a past-year alcohol use disorder [5]. However, alcohol use is not distributed evenly. Previous research has identified distinct regional differences in alcohol consumption [6], and ongoing public health surveillance finds varying levels of hazardous drinking by state [7].

Importantly, moderate to chronic heavy drinking suppresses immune responses and can increase susceptibility to and the seriousness of infectious and respiratory diseases, such as pneumonia, tuberculosis, and acute respiratory stress syndromes [810]. This same alcohol-induced immunosuppression may also increase the severity of COVID-19 infections at a population level. Given the effects of alcohol on infectious disease progression, high levels of county-level alcohol use preceding the COVID-19 pandemic are a particularly concerning population health problem.

Many burdens of heavy alcohol use can be exacerbated by other factors. For example, residence in more rural communities is an important social and environmental determinant of health and is associated with lower life expectancy [11], higher rates of all-cause morbidity and mortality [12], and multiple intersecting barriers in access to healthcare [13]. Roughly 60 million people, or nearly one in five Americans, live in a rural area [14].

Rural counties and towns have faced considerable challenges in dealing with the COVID-19 pandemic, partly attributable to weakened health care infrastructure, health care provider shortages, lower socio-economic status, and higher proportions of aging residents [1517]. In the summer of 2020, the incidence rate of COVID-19 infections began to increase faster in rural areas, following the initial outbreaks in larger and more metropolitan urban centers [18]. Case fatality rates have also been shown to be higher in rural counties [19], rural states [20], and in rural counties with larger proportions of Black and Hispanic populations [21].

Furthermore, alcohol use and misuse vary by rurality. Rural populations are more likely than their urban counterparts to abstain from alcohol, but among non-abstainers, rates of alcohol use disorders (AUD) are higher in more rural areas [22]. Better understanding the association between county-level alcohol use and rurality will inform efforts to mitigate health inequities across regions.

To our knowledge, no studies have investigated the intersection between alcohol misuse, rurality, and COVID-19. This study was designed to measure disparities in COVID-19 infection severity at the county population level, which we operationalized as case fatality rates. We considered two important determinants of population health–excessive alcohol consumption and rurality–and their relationship to case fatality rates at the county level. Although largely an exploratory study, we expected to see higher burdens and more disparities in rural areas.

Materials and methods

Data source and sample

We used the New York Times’ COVID-19 data repository (https://github.com/nytimes/covid-19-data) to calculate case fatality rates in US counties with at least one COVID-19 death reported between October 1, 2020 and December 31, 2020. We restricted our sample to the last quarter of 2020 because this was the most recent data available at the time of analysis, and to highlight relationships during a time of particularly high case and death rates across the US. Case fatality rates were calculated by dividing the total number of deaths by the total number of confirmed cases, indicating the percentage of individuals with COVID-19 who died. Case fatality rates help to illustrate disproportionate mortality burdens between populations or geographies where the case rates may be similar but underlying social, economic, and health conditions differ. Furthermore, case fatality rates allow for more appropriate comparison between counties with substantial variation in population size.

Variables

Our focal predictor was a single combined variable representing levels of pre-pandemic county-level excessive alcohol consumption and rurality. The alcohol component was sourced from the Robert Wood Johnson Foundation’s 2018 County Health Rankings [23] and merged to the analytic dataset by county. County-level excessive alcohol consumption was defined as the percentage of adults who reported exceeding recommended drinking limits, which included: 1) four or more alcoholic drinks consumed on a single occasion in the past month for women and five or more alcoholic drinks on a single occasion in the past month for men (i.e., binge drinking), and/or 2) more than one drink on average per day for women, and more than two drinks on average per day for men (i.e., exceeding recommended daily limits). We used Waldorf and Kim’s Index of Relative Rurality (IRR) to define county-level rurality [24]. The IRR is a threshold-free measurement of rurality at the county level, where greater scores indicate greater rurality. The IRR takes into account population size, population density, remoteness, and built-up infrastructure as a percentage of total land area. Our final analytic variable included nine possible combinations of low/mid/high alcohol consumption and low/mid/high rurality. The low/mid/high cutoffs for both variables were based on tertiles, resulting in nine categorical levels of 1) low alcohol use and low rurality, 2) low alcohol use and mid rurality, 3) low alcohol use and high rurality, 4) mid alcohol use and low rurality, 5) mid alcohol use and mid rurality, 6) mid alcohol use and high rurality, 7) high alcohol use and low rurality, 8) high alcohol use and mid rurality, and 9) high alcohol use and high rurality.

We included several county-level variables in our model that likely confound the relationships between excessive alcohol consumption, rurality, and COVID-19 case fatality rates. All covariates (prevalence of smoking, obesity, diabetes, older age, race/ethnicity, and unemployment) were sourced from the County Health Rankings data. Smoking was defined as the percent of a county’s adults who reported currently smoking tobacco every day or on most days, and who have smoked at least 100 cigarettes in their lifetime. Older age was defined as the percent of a county’s population aged 65 years and older. Obesity was defined as the percentage of a county’s adults, aged 20 years and older, who report a body mass index of 30 kg/m2 or greater. Diabetes was defined as the percentage of a county’s adults, aged 20 years and older, who have been diagnosed with diabetes. Race/ethnicity was represented by the percentage of a county made up of self-identified nonwhite residents. Finally, unemployment was defined as the percentage of a county’s civilian labor force, aged 16 and older, who were unemployed but seeking work during the past week.

Our final analytic sample included 3,039 counties and county equivalents in the US that had at least one death recorded and for which we had complete data for all study variables, (i.e., matching FIPS codes between the New York Times and County Health Rankings datasets), representing 96% of all 3,143 US counties and county equivalents.

Analysis

We used SAS software for all analyses (Version 9.4) [25]. We described the percentage of counties that fall into each level of our categorical predictor, and the median, minimum, and maximum values for all continuous study variables. We reported the IRR scores and case fatality rates for the ten counties with the lowest percentage of excessive drinking and for the ten counties with the highest percentage of excessive drinking. We also developed an original, county-level US map to illustrate the geographic distribution of excessive alcohol consumption and rurality using Mathematica software.

We used beta regression to model case fatality rates. Beta regression is a more appropriate framework than linear regression to model outcomes that are proportions or rates, such as case fatality rates [26]. We considered three separate structures for our model and selected the final model based on how well the model fit the data, as indicated by the lowest Akaike Information Criterion (AIC) value. All candidate models were adjusted for county-level indicators of smoking, obesity, diabetes, age, race/ethnicity, and unemployment. Candidate model #1 used a non-structural, fixed effects design (AIC = -19337). Candidate model #2 included state as a random effect, under the assumption that counties clustered within one state are likely more similar to each other than to counties in other states (AIC = -19845) [27]. Finally, candidate model #3 included state defined as a random effect and accounted for the structural spatial autocorrelation between counties. This model incorporated a Gaussian correlation structure using geographic county centroids (AIC = -19843) [28]. Candidate model #2 was selected as the final analytic model because it demonstrated the lowest AIC value and the best fit to the data. As a multilevel model, we reported the Intraclass Correlation Coefficient (ICC), which is an indicator of how much of the total variation in county case fatality rates is accounted for by the state. Finally, to model our nine-level categorical focal predictor, we included eight dummy variables and used low alcohol use/low rurality as the referent group.

We interpreted study results using a significance threshold of α = 0.05. We also identified no problematic multicollinearity between the study variables, defined as a Pearson correlation coefficient greater than 0.80.

Results

County population characteristics are reported in Table 1. The median county-level COVID-19 case fatality rate was 1.57% (range 0.01% - 16.08%) (Table 2). The high upper limit of this range was occupied by Kenedy County, Texas (16.08%; 184 deaths out of 1,144 reported cases), and the second highest was Franklin County, Massachusetts (11.45%; 7,003 deaths out of 61,156 reported cases). Among all counties, 7.15% (n = 217) were at levels of low alcohol and low rurality, and 10.46% (n = 318) were at levels of high alcohol and high rurality. Eight out of the ten counties with the lowest percentages of excessive drinking were in the US South (Table 3). The ten counties with the lowest alcohol levels had a combined average IRR score of 0.55 (more rural) and an average case fatality rate of 2.31%. Conversely, all ten counties with the highest levels of excessive alcohol use were in Wisconsin, with a combined average IRR score of 0.45 (less rural) and an average case fatality rate of 0.65%. Counties with both high alcohol and high rurality tended to be clustered in the upper Midwest, parts of the mountain West, and in the states of Nevada and Alaska (Fig 1).

thumbnail
Fig 1. US county map of levels of excessive alcohol consumption and rurality.

Original map was created by the study team with Mathematica software.

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

thumbnail
Table 2. COVID-19 case fatality rates by county-level alcohol use and rurality (n = 3,039 counties).

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

thumbnail
Table 3. Characteristics of top 10 and lowest 10 counties by county-level excessive drinking.

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

In our fitted model, compared to counties with low alcohol use/low rurality (referent), counties with mid levels of alcohol use and high levels of rurality (β = -0.14, p = 0.019), as well as counties with high levels of alcohol use and mid (β = -0.17, p = 0.008) and high levels of rurality (β = -0.24, p<0.001), demonstrated significantly lower case fatality rates (Table 4). County prevalence of smoking, obesity, and unemployment was not associated with case fatality rate; however, county prevalence of diabetes, older adults, and non-white residents were each positively associated with greater case fatality rates. The ICC value demonstrated that roughly 4% of the variability of case fatality rates between states was attributed to county clustering within states (p<0.001).

thumbnail
Table 4. Multivariate beta regression modeling case fatality rate (n = 3,039 counties).

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

Discussion

To our knowledge, this is the first study to examine county-level relationships between excessive alcohol consumption before the pandemic, rurality, and COVID-19 mortality. Notably, results indicated a joint contribution of alcohol use patterns and level of rurality. Case fatality rates were generally lowest in more rural counties; compared to low alcohol use/low rurality counties, rates were lower in mid alcohol use/high rurality, high alcohol use/mid rurality and high alcohol use/high rurality counties in particular. Our findings highlight the importance of considering both behavioral and environmental determinants of COVID-19 outcomes.

Our findings did not support our hypothesis of higher burdens and more disparities in rural areas. Instead, we found lower fatality burdens in some mid and high rural counties. In rural settings, other factors besides alcohol use may be more important drivers of COVID-19 fatality, such as a county’s age distribution, racial/ethnic distribution, and unemployment levels. Some rural areas may offer protection against the effect of alcohol on mortality. Transmission may be less efficient than in larger and more concentrated urban populations, and drinking norms in rural areas may include smaller and less frequent group meetings in bars or indoor environments, or higher levels of drinking in isolation. Future studies that incorporate longitudinal changes in alcohol use behaviors could improve our understanding of rural and urban differences in health behavior and substance misuse at multiple time points before, during, and following the pandemic.

Potential limitations

These results should be considered in light of several potential limitations. County-level excessive drinking may have been misestimated as individual survey responses about drinking habits may be subject to social desirability biases [29]. Such underreporting would decrease our ability to detect an association. Second, as a secondary data analysis, some variables were not available, such as asymptomatic COVID-19 cases, COVID-19 testing capacity, relevant county-level comorbidities (i.e., chronic lung disease or asthma, liver disease, serious heart conditions, and other diseases affecting immune response), and county healthcare characteristics (i.e., available ventilators and provider shortage areas). Finally, as an ecological study using county-level data, no individual-level information was available. As such, there is no way to identify the demographic characteristics of the population that made up the cases and deaths. Future studies that evaluate relationships between individual-level alcohol use and COVID-19 mortality outcomes would contribute substantially to the literature around behavioral health and infectious disease outbreaks. As this study was undertaken in the midst of the COVID-19 pandemic, our pattern of results should be confirmed as more data accumulate. Nevertheless, our findings address a gap in knowledge about the association of excessive alcohol consumption, rurality, and COVID-19 outcomes.

Conclusion

Our findings demonstrate lower COVID-19 fatality burdens in some counties, in particular counties with mid/high excessive alcohol use and counties with mid/high levels of rurality. These findings were contrary to our hypothesis and highlight the intersecting roles that excessive alcohol consumption and geography play in baseline risk for COVID-19 outcomes at the community and population levels.

References

  1. 1. Centers for Disease Control and Prevention. Cases in the U.S. Centers for Disease Control and Prevention. 2021;Cases, Data & Surveillance (Updated January 3, 2021). pmid:33474834
  2. 2. Yang J, Zheng Y, Gou X, Pu K, Chen Z, Guo Q, et al. Prevalence of comorbidities and its effects in patients infected with SARS-CoV-2: A systematic review and meta-analysis. International Journal of Infectious Diseases. 2020;94:91–5. pmid:32173574
  3. 3. Yancy C. COVID-19 and African Americans. JAMA. 2020;Published online April 15, 2020. pmid:32293639
  4. 4. Laurencin C, McClinton A. The COVID-19 Pandemic: A call to action toiIdentify and address racial and ethnic disparities. Journal of Racial and Ethnic Health Disparities. 2020;Published online April 18, 2020.
  5. 5. Grant BF, Goldstein RB, Saha TD, Chou SP, Jung J, Zhang H, et al. Epidemiology of DSM-5 Alcohol Use Disorder: Results From the National Epidemiologic Survey on Alcohol and Related Conditions III. JAMA Psychiatry. 2015;72(8):757–66. pmid:26039070
  6. 6. Kerr W. Categorizing US state drinking practices and consumption trends. International Journal of Environmental Research and Public Health. 2010;7(1):269–83. pmid:20195444
  7. 7. Centers for Disease Control and Prevention. Behavioral Risk Factor Surveillance System. Prevalence Data and Analysis Tools: BRFSS Prevalence & Trends Data. 2020.
  8. 8. Pasala S, Barr T, Messaoudi I. Impact of alcohol abuse on the adaptive immune system. Alcohol Research: Current Reviews. 2015;37(2):185–97.
  9. 9. Rehm J. The risks associated with alcohol use and alcoholism. Alcohol Research and Health. 2011;34(2):135–43. pmid:22330211
  10. 10. Sarkar D, Jung M, Wang H. Alcohol and the immune system. Alcohol Research: Current Reviews. 2015;37(2):153–5.
  11. 11. Singh G, Siahpush M. Widening rural-urban disparities in life expectancy, U.S., 1969–2009. American Journal of Preventive Medicine. 2014;46(2):e19–e29. pmid:24439358
  12. 12. Singh G, Siahpush M. Widening rural-urban disparities in all-cause mortality and mortality from major causes of death in the USA, 1969–2009. Journal of Urban Health. 2013;91(2):272–92.
  13. 13. Douthit N, Kiv S, Dwolatzky T, Biswas S. Exposing some important barriers to health care access in the rural USA. Public Health. 2015;129(6):611–20. pmid:26025176
  14. 14. Census Bureau U.S. Rural America. United States Census Bureau. How does the Census Bureau define rural? 2020.
  15. 15. Melvin S, Wiggins C, Burse N, Thompson E, Monger M. The role of public health in COVID-19 emergency response efforts from a rural health perspective. Preventing Chronic Disease. 2020;17(E70):1–6. pmid:32701430
  16. 16. Henning-Smith C. The unique impact of COVID-19 on older adults in rural areas. Journal of Aging & Social Policy. 2020;32:4–5. pmid:32475255
  17. 17. Shah P, Owens J, Franklin J, Mehta A, Heymann W, Sewell W, et al. Demographics, comorbidities and outcomes in hospitalized COVID-19 patients in rural southwest Georgia. Annals of Medicine. 2020;52(7):354–60. pmid:32620056
  18. 18. Centers for Disease Control and Prevention. COVID-19 Incidence, by Urban/Rural Classification—United States, January 22-October 31, 2020. Morbidity and Mortality Weekly Report (MMWR). 2020;69(COVID-19 Stats):1753.
  19. 19. Pro G, Hubach R, Wheeler D, Camplain R, Haberstroh S, Giano Z, et al. Differences in US COVID-19 case rates and case fatality rates across the urban-rural continuum. Rural and Remote Health. 2020;20. pmid:32811154
  20. 20. Ahmed R, Williamson M, Akhter Hamid M, Ashraf N. United States county-level COVID-19 death rates and case fatality rates vary by region and urban status. Healthcare. 2020;8(330).
  21. 21. Cheng K, Sun Y, Monnat S. COVID-19 death rates are higher in rural counties with larger shares of Blacks and Hispanics. The Journal of Rural Health. 2020;36:602–8. pmid:32894612
  22. 22. Borders T, Booth B. Rural, suburban, and urban variations in alcohol consumption in the United States: Findings from the National Epidemiologic Survey on Alcohol and Related Conditions. The Journal of Rural Health. 2007;23(4):314–21. pmid:17868238
  23. 23. Robert Wood Johnson Foundation. County Health Rankings & Roadmaps. Robert Wood Johnson Foundation. 2020;2020 Measures.
  24. 24. Waldorf B, Kim A. Defining and measuring rurality in the US: From typologies to continuous indices. Workshop on Rationalizing Rural Area Classification. 2015;Commissiond Paper: April 16/17, 2015, Keek Center, Washington, D.C.
  25. 25. SAS Institute Inc. Cary, NC. SAS Institute Inc., 2013.
  26. 26. Ferrari S, Cribari-Neto F. Beta regression for modelling rates and proportions. Journal of Applied Statistics. 2004;31(7):799–815.
  27. 27. Ene M, Leighton E, Blue G, Bell B. Multilevel Models for Categorical Data using SAS PROC GLIMMIX: The Basics. SAS Papers. 2015;Paper 3430–2015.
  28. 28. Cressie N. Statistics for Spatial Data. New York, NY: John Wiley & Sons, Inc.; 1991.
  29. 29. Davis C, Thake J, Vilhena N. Social desirability biases in self-reported alcohol consumption and harms. Addictive Behaviors. 2010;35(4):302–11. pmid:19932936