Skip to main content
Advertisement
  • Loading metrics

Social determinants of mortality from COVID-19: A simulation study using NHANES

  • Benjamin Seligman ,

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

    bseligm1@bidmc.harvard.edu

    Affiliations New England Geriatric Research, Education, and Clinical Center, VA Boston Healthcare System, Boston, Massachusetts, United States of America, Harvard Medical School, Boston, Massachusetts, United States of America, Harvard Center for Population and Development Studies, Cambridge, Massachusetts, United States of America

  • Maddalena Ferranna,

    Roles Conceptualization, Writing – review & editing

    Affiliation Department of Global Health and Population, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America

  • David E. Bloom

    Roles Supervision, Writing – original draft, Writing – review & editing

    Affiliation Department of Global Health and Population, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America

Correction

29 Dec 2021: Seligman B, Ferranna M, Bloom DE (2021) Correction: Social determinants of mortality from COVID-19: A simulation study using NHANES. PLOS Medicine 18(12): e1003888. https://doi.org/10.1371/journal.pmed.1003888 View correction

Abstract

Background

The COVID-19 epidemic in the United States is widespread, with more than 200,000 deaths reported as of September 23, 2020. While ecological studies show higher burdens of COVID-19 mortality in areas with higher rates of poverty, little is known about social determinants of COVID-19 mortality at the individual level.

Methods and findings

We estimated the proportions of COVID-19 deaths by age, sex, race/ethnicity, and comorbid conditions using their reported univariate proportions among COVID-19 deaths and correlations among these variables in the general population from the 2017–2018 National Health and Nutrition Examination Survey (NHANES). We used these proportions to randomly sample individuals from NHANES. We analyzed the distributions of COVID-19 deaths by race/ethnicity, income, education level, and veteran status. We analyzed the association of these characteristics with mortality by logistic regression. Summary demographics of deaths include mean age 71.6 years, 45.9% female, and 45.1% non-Hispanic white. We found that disproportionate deaths occurred among individuals with nonwhite race/ethnicity (54.8% of deaths, 95% CI 49.0%–59.6%, p < 0.001), individuals with income below the median (67.5%, 95% CI 63.4%–71.5%, p < 0.001), individuals with less than a high school level of education (25.6%, 95% CI 23.4% –27.9%, p < 0.001), and veterans (19.5%, 95% CI 15.8%–23.4%, p < 0.001). Except for veteran status, these characteristics are significantly associated with COVID-19 mortality in multiple logistic regression. Limitations include the lack of institutionalized people in the sample (e.g., nursing home residents and incarcerated persons), the need to use comorbidity data collected from outside the US, and the assumption of the same correlations among variables for the noninstitutionalized population and COVID-19 decedents.

Conclusions

Substantial inequalities in COVID-19 mortality are likely, with disproportionate burdens falling on those who are of racial/ethnic minorities, are poor, have less education, and are veterans. Healthcare systems must ensure adequate access to these groups. Public health measures should specifically reach these groups, and data on social determinants should be systematically collected from people with COVID-19.

Author summary

Why was this study done?

  • The COVID-19 epidemic in the United States of America (US) has been found to disproportionally affect racial/ethnic minorities and areas with higher concentrations of poverty.
  • Few data on COVID-19 deaths with respect to individual-level social determinants of health other than race/ethnicity are available.

What did the researchers do and find?

  • We used the 2017–2018 National Health and Nutrition Examination Survey (NHANES) and publicly reported data on COVID-19 deaths from the public health agencies of the US, China, the United Kingdom, Spain, Italy, and France to simulate COVID-19 deaths among noninstitutionalized (e.g., not residing in a prison or nursing home) adults aged 20 years or older in the US.
  • We found large social gradients in COVID-19 mortality. Adults from households earning less than the median income made up two-thirds of COVID-19 deaths, while those with less than a high school education accounted for approximately 1 in 4 deaths. Veterans also accounted for nearly 1 in 5 deaths, despite representing less than one-tenth of the population.
  • Our simulation reproduced known racial/ethnic disparities in COVID-19 mortality.

What do these findings mean?

  • The associations between social determinants of health and COVID-19 mortality are similar in scale to those between hypertension and diabetes and COVID-19 mortality.
  • COVID-19 mitigation will require measures to support people from low-income, low-education communities and the healthcare systems that serve them.

Introduction

The SARS-CoV-2 pandemic is a profound challenge to healthcare systems and societies. While no segment of society is unaffected, some groups face disproportionate burdens of illness. Multiple studies have established the increased risks of severe illness and mortality with age and comorbidity [16]. Within the United States of America (US), geographic differences in the prevalence and incidence of COVID-19 morbidity and mortality are known [4]. However, data on social determinants of health and the associated risks of infection with SARS-CoV-2 and death from COVID-19 are less available. US Centers for Disease Control and Prevention (CDC) surveillance has demonstrated the overrepresentation of African Americans among people hospitalized with COVID-19, but more detailed information is lacking in disease statistics [7,8]. Other major sources of COVID-19 statistics likewise lack information on social determinants of health [9,10].

Social determinants of health play an important role in the spread of and mortality from epidemics, from influenza in 1918 to Ebola in 2014 and others, including the first severe acute respiratory syndrome (SARS) epidemic [1116]. They affect susceptibility to acquiring infection, due to differences in social contacts and differences in living circumstances [15,1719]. They affect severity of illness, in part through concentration of comorbidity in susceptible groups [2022]. Finally, they also affect outcomes through differential access to healthcare, which may in turn limit opportunities to identify and contain local outbreaks [17,2325].

Understanding the COVID-19 burden in terms of social determinants of health is important for policymaking and targeting both public health and clinical interventions. Ecological-level data show that areas with higher poverty rates and larger proportions of individuals who identify as racial/ethnic minorities have higher COVID-19 mortality [26]. However, few analyses to date consider social determinants at the individual level. Here, we use simulation to investigate the distribution of COVID-19 mortality with respect to social determinants of health at an individual level.

Methods

Because data on social determinants of health among COVID-19 deaths are limited, we simulated mortality. We started with univariate distributions of COVID-19 deaths by age, sex, race/ethnicity, and comorbid conditions—as reported by multiple public health agencies from the spring to summer of 2020—and the correlations among these variables from the 2017–2018 National Health and Nutrition Examination Survey (NHANES). Using these data, we estimated the joint distribution of deaths by these variables. The probabilities obtained from this joint distribution were used as weights for simulating COVID-19 deaths in the NHANES cohort. This provided a sample with which we could investigate social determinants of health in COVID-19 mortality.

NHANES sample

Data on individuals aged 20 years or older in the general population were taken from the 2017–2018 cycle of NHANES (n = 5,265), a nationally representative study of the health of noninstitutionalized Americans [27]. This study comprised a questionnaire, a physical exam, and selected laboratory studies. Questions asked addressed health status, including comorbidities; demographic information; and social determinants of health, such as income and education.

To code the presence of certain comorbidities we used definitions incorporating both questionnaire and examination or laboratory findings. Hypertension was defined as the presence of either self-reported hypertension in the health questionnaire or an average blood pressure greater than 140/90 mm Hg on examination; we did not incorporate use of antihypertensives due to missingness [28,29]. Diabetes was defined as self-report or a hemoglobin A1C greater than 6.5% [30]. Chronic kidney disease (CKD) was defined as self-report of either kidney failure or dialysis or an estimated glomerular filtration rate less than 60 ml/min/1.73 m2 using CKD-EPI [31,32]. Ischemic heart disease (IHD) was defined as self-report of coronary heart disease, angina, or heart attack. Chronic obstructive pulmonary disease (COPD) was defined as self-report of emphysema, chronic bronchitis, or COPD. Additionally, we created a comorbidity index based on the number of comorbidities each participant had, including the presence of hypertension, diabetes, CKD, IHD, COPD, and self-reported asthma, congestive heart failure, stroke, liver disease, and cancer.

For social determinants of health, we considered race/ethnicity, income, education, and veteran status based on their availability in NHANES. We specifically considered veterans because they are disproportionately older compared with the general population and frequently face service-related health conditions [33,34]. Cases with missing data or with responses of “don’t know” or “refused” were dropped from the final analysis. For race/ethnicity, the categories “Mexican American” and “Other Hispanic” were combined into “Hispanic.” For income, we used household income, limited analysis to responses within specific income ranges, and did not include the answers “under $20,000” and “$20,000 and over.” In the NHANES data, age was top-coded at 80 years, and income was top-coded at $100,000 per year.

Simulation of COVID-19 deaths

Simulating COVID-19 deaths involved 3 steps, with the last 2 repeated in each simulation run. This was necessary given the absence of cross-tabulations of COVID-19 deaths by age, sex, and comorbidity, with only univariate distributions consistently available. In the first step, marginal distributions of age, sex, and comorbidity were taken from multiple public health agencies and used to estimate prior distributions. In the second step, the marginal distributions of age, sex, and comorbidity were randomly drawn from the priors. Using the correlations among these variables from NHANES, we then approximated their joint distribution. In the third step, the joint distribution was used to reweight the NHANES sample to represent COVID-19 deaths. We give more detail on each step in the following.

As joint distributions of characteristics of COVID-19 deaths were unavailable, we first obtained their marginal distributions. We considered the marginal distributions of deaths by age, sex, the absence of comorbidity, and presence of each of hypertension, diabetes, CKD, IHD, COPD, and cancer. These were obtained from the CDC, the United Kingdom’s Office for National Statistics, Santé Publique France, Istituto Superiore di Sanità in Italy, Instituto de Salud Carlos III in Spain, and the China Center for Disease Control and Prevention [3541]; S1 Text provides further details. Comorbidity data from the CDC, the Office for National Statistics, and Santé Publique France were excluded from fitting the model due to their reliance on death certificate data, which underestimate the prevalence of comorbidities compared with reported data from US hospitals and from systematically collected data from Italy, Spain, and China [42,43]. We then fit maximum-likelihood beta or Dirichlet (for age) priors to each marginal distribution. To estimate a joint distribution, race was set as an indicator variable for non-Hispanic white and fixed to the proportion reported by the CDC, as this was the only agency that reported race/ethnicity.

The second step was to approximate the joint distributions of the aforementioned characteristics. To do this, we used marginal distributions for age, sex, and each of the comorbidities drawn from their respective priors, the fixed marginal distribution of non-Hispanic white individuals, and the correlations of these variables from NHANES to estimate a joint distribution using a Gaussian copula [44]. We then assigned the joint probabilities as weights to each NHANES participant such that the participants would, in total, represent 200,000 deaths. We repeated step 2 1,000 times. We checked calibration against CDC-published distributions of deaths by age, by sex, and by race/ethnicity. We report our results and data in terms of gender, which is reported in NHANES. However, data from public health agencies are given in terms of sex.

Analytic approach

With the simulated data we produced distributions across variables of interest. Differences in proportions were assessed by a bootstrap chi-squared test compared to a null distribution of 1,000 replicates randomly sampled from NHANES. To assess the independent contributions of these variables, we also analyzed the data as a case–control study using multivariable logistic regression. In this instance, the NHANES sample was weighted to represent the general population to serve as controls, and to represent COVID-19 deaths to serve as cases. Confidence intervals were bootstrapped from the 1,000 simulation runs, while p-values were bootstrapped through 1,000 replicates of a null. We did not have a prespecified analytic plan.

Neither individually identifiable information nor patient health information was used in this study; all data are publicly available and are described in S1 Text. Analysis was conducted by author BS in R version 4.0.2 using the package GenOrd version 1.4.0 [45,46]. Code used for the simulation is available from https://www.hsph.harvard.edu/pgda/data/.

Results

Table 1 shows the distributions of characteristics of the simulated COVID-19 deaths and the characteristics of their comparison population from NHANES 2017–2018. Fig 1 shows the results of model calibration against age and gender, indicating acceptable calibration of the simulation against the observed age and gender distributions in the US. Clearly, the general pattern of disproportionate deaths among the elderly holds in the simulation, although the simulation indicates a higher proportion of deaths at younger ages than is observed in CDC data. The results for gender show close calibration with observed data and follow the pattern of a greater proportion of deaths among males.

thumbnail
Fig 1. Calibration of simulated distributions of deaths by age and gender against reported distributions from the CDC.

Age (top); gender (bottom). Gray bars and whiskers show mean and 2.5th–97.5th percentile range for simulated proportions, the blue circles show the CDC-reported proportions, and the red triangles show the NHANES-estimated proportions for the general population. CDC, US Centers for Disease Control and Prevention; NHANES, National Health and Nutrition Examination Survey.

https://doi.org/10.1371/journal.pmed.1003490.g001

thumbnail
Table 1. Characteristics of simulated COVID-19 deaths and the weighted NHANES sample.

https://doi.org/10.1371/journal.pmed.1003490.t001

Fig 2 and Table 1 show the distributions of social determinants of health among COVID-19 deaths and in the general population. When we divide the population between non-Hispanic white and all races/ethnicities, the simulation shows that all others are overrepresented among COVID-19 deaths compared with their proportion in the general population (54.8% of deaths, 95% CI 49.0%–59.6%, p < 0.001). When data are broken out by specific race/ethnicity (S1 Fig), we see that the simulation captures some disparities, but understates others, particularly among African Americans.

thumbnail
Fig 2. Simulated proportions of deaths by race/ethnicity, income, veteran status, and education level.

Panels show mean and 2.5th–97.5th percentile range for the simulations (gray bars and whiskers), the CDC-reported proportions (blue circles, only for race/ethnicity), and the NHANES-estimated proportions for the general population (red triangles). AA, associate’s degree; GED, General Educational Development; NHANES, National Health and Nutrition Examination Survey.

https://doi.org/10.1371/journal.pmed.1003490.g002

With respect to income, COVID-19 deaths are disproportionately among middle- and lower-income people. Those making below the median income, which is the $55,000–$64,999 category, make up 67.5% of deaths (95% CI 63.4%–71.5%, p < 0.001). Similarly, disparities exist regarding level of education: Individuals who have less than a high school level of education are overrepresented among COVID-19 deaths (25.6%, 95% CI 23.4%–27.9%, p < 0.001). Disparities with respect to being a veteran are smaller on an absolute scale; however, veterans make up almost 20% of deaths in the simulation (19.5%, 95% CI 15.8%–23.4%, p < 0.001) versus 9% of the population.

Table 2 shows the results of the multivariable logistic regression of COVID-19 mortality against age, gender, and the social determinants of health studied. The coefficients show increasing odds of mortality with age, while female gender is associated with lower odds. These data also show increasing odds of mortality for Hispanics, African Americans, and Asian Americans/Pacific Islanders. For income and education, as the level of either one declines, the odds of mortality rise, except for the “some college” category. In this group, relative to high school graduates, the odds of mortality are increased. Being a veteran was not significantly associated with increased odds of death from COVID-19 independent of the other variables considered.

thumbnail
Table 2. Odds ratios of death from COVID-19 based on age, gender, race/ethnicity, income, education, and veteran status.

https://doi.org/10.1371/journal.pmed.1003490.t002

Discussion

Our simulation provides evidence of the scale of social and economic disparities in the COVID-19 epidemic in the US. Mortality from infection disproportionately strikes individuals from low- or middle-income families, individuals with less education, individuals who are of racial/ethnic minorities, and individuals who have served in the military. The disparities identified here are likely underestimates of their true scope. As the simulation only considers age, gender, race/ethnicity, and comorbidity as drivers of death, direct effects of social determinants of health on vulnerability to infection and mortality are not explicitly modeled. These effects include crowded living conditions, limited access to care, and economic hardship that may force people to continue to risk exposure by working. This analysis also does not consider the low wages of many essential jobs that place workers at increased risk of infection, including healthcare workers such as environmental services staff and patient care attendants [47]. Even with these limitations, however, the associations of COVID-19 mortality with social adversity in our study are comparable to the associations of COVID-19 mortality with diabetes (odds ratio 1.75–1.90) and hypertension (relative risk 2.21) [4851]. Our findings differ from a study that found no difference in income among persons who died or had critical illness during hospitalization for COVID-19 versus those who survived hospitalization without critical illness [52]. This is likely due to our consideration of the general population, versus those who are hospitalized.

Our findings provide further evidence that efforts to reduce COVID-19 mortality should involve prioritizing the needs of disadvantaged communities. This could involve greater assistance for healthcare systems that disproportionately care for low-income or low-education people, such as many public and rural hospital systems, the Veterans Health Administration, and the Indian Health Service. This would also need to involve public health measures, such as paid sick leave, income support, and expansion of health insurance access, to make social distancing more feasible and make care accessible. This may also involve messaging around social distancing and other health behaviors targeted to groups facing social disadvantage. Further, it is crucial that we systematically collect data on social determinants among COVID-19 cases and deaths. These are key data for understanding and controlling the epidemic [53,54].

Our approach has several limitations in addition to not explicitly modeling the effects of social determinants of health. We assume that the correlations among age, gender, race/ethnicity, and comorbidity in the population of individuals who have died from COVID-19 are the same as those in the general population. This is not true in the case of the correlation between gender and age, as males are disproportionately represented in all age groups among COVID-19 deaths, whereas females predominate at older ages in the general population. NHANES data do not include nursing home residents, who make up a large fraction of COVID-19 deaths. This likely contributes to differences between the simulated and observed age distribution of deaths. The other model parameters are based on data from multiple countries, which may affect representativeness for the epidemic in the US, especially with regard to comorbidities, where CDC data could not be used for calibration. Finally, we depend on reliable reporting of deaths, and there is concern that not all COVID-19 deaths are being registered [55,56]. We recommend that further work systematically collect and report data on social determinants of health among individuals affected by COVID-19.

As the Introduction notes, social determinants of health have been associated with disease burden in past epidemics. They have potentially facilitated increased transmission as well, suggesting that control may depend in part on addressing the epidemic specifically among individuals who are poor, have less education, or live in poor conditions [12,14]. As more data become available, we may better understand the roles of particular social determinants and be able to design more effective interventions. By protecting the health of the most vulnerable, such measures could mitigate the toll of the COVID-19 pandemic and protect all Americans.

Supporting information

S1 Fig. Simulated proportions of deaths by race/ethnicity.

Figure shows mean and 2.5th–97.5th percentile range for the simulations (gray bars and whiskers), the CDC-reported proportions (blue circles), and the NHANES-estimated proportions for the general population (red triangles). AA, African American; AAPI, Asian American/Pacific Islander; NHW, non-Hispanic white.

https://doi.org/10.1371/journal.pmed.1003490.s001

(TIF)

Acknowledgments

We would like to thank Prof. Marcello Pagano at the Harvard T.H. Chan School of Public Health for his helpful comments on our model.

The work of DEB and MF on this paper was supported by the Value of Vaccination Research Network (VoVRN) through a grant from the Bill & Melinda Gates Foundation (Grant OPP1158136). The content is solely the responsibility of the authors and does not necessarily reflect the views of the VoVRN or the foundation.

This material is the result of work supported with resources and the use of facilities at the VA Boston Healthcare System. The contents do not represent the views of the US Department of Veterans Affairs or the United States government.

References

  1. 1. Guan W, Liang W, Zhao Y, Liang H, Chen Z, Li Y, et al. Comorbidity and its impact on 1590 patients with COVID-19 in China: a nationwide analysis. Eur Respir J. 2020;55(5):2000547. pmid:32217650
  2. 2. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet. 2020;395(10229):1054–62. pmid:32171076
  3. 3. Yang J, Zheng Y, Gou X, Pu K, Chen Z, Guo Q, et al. Prevalence of comorbidities in the novel Wuhan coronavirus (COVID-19) infection: a systematic review and meta-analysis. Int J Infect Dis. 2020;94:91–5. pmid:32173574
  4. 4. Bialek S, Bowen V, Chow N, Curns A, Gierke R, Hall A, et al. Geographic differences in COVID-19 cases, deaths, and incidence—United States, February 12–April 7, 2020. MMWR Morb Mortal Wkly Rep. 2020;69(15):465–71. pmid:32298250
  5. 5. Roxby AC, Greninger AL, Hatfield KM, Lynch JB, Dellit TH, James A, et al. Detection of SARS-CoV-2 among residents and staff members of an independent and assisted living community for older adults—Seattle, Washington, 2020. MMWR Morb Mortal Wkly Rep. 2020;69(14):416–8. pmid:32271726
  6. 6. Chow N, Fleming-Dutra K, Gierke R, Hall A, Hughes M, Pilishvili T, et al. Preliminary estimates of the prevalence of selected underlying health conditions among patients with coronavirus disease 2019—United States, February 12–March 28, 2020. MMWR Morb Mortal Wkly Rep. 2020;69(13):382–6. pmid:32240123
  7. 7. Centers for Disease Control and Prevention. CDC COVID data tracker. Atlanta: Centers for Disease Control and Prevention; 2020 [cited 2020 Apr 24]. Available from: https://www.cdc.gov/coronavirus/2019-ncov/cases-updates/cases-in-us.html.
  8. 8. Garg S, Kim L, Whitaker M, O’Halloran A, Cummings C, Holstein R, et al. Hospitalization rates and characteristics of patients hospitalized with laboratory-confirmed coronavirus disease 2019—COVID-NET, 14 states, March 1–30, 2020. MMWR Morb Mortal Wkly Rep. 2020;69(15):458–64. pmid:32298251
  9. 9. Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect Dis. 2020;20(5):533–4. pmid:32087114
  10. 10. 1Point3Acres. Global COVID-19 tracker & interactive charts. 1Point3Acres; 2020 [cited 2020 Apr 24]. Available from: https://coronavirus.1point3acres.com/en.
  11. 11. Mamelund SE. A socially neutral disease? Individual social class, household wealth and mortality from Spanish influenza in two socially contrasting parishes in Kristiania 1918–19. Soc Sci Med. 2006;62(4):923–40. pmid:16084634
  12. 12. Grantz KH, Rane MS, Salje H, Glass GE, Schachterle SE, Cummings DAT. Disparities in influenza mortality and transmission related to sociodemographic factors within Chicago in the pandemic of 1918. Proc Natl Acad Sci U S A. 2016;113(48):13839–44. pmid:27872284
  13. 13. Mamelund S-E. 1918 pandemic morbidity: the first wave hits the poor, the second wave hits the rich. Influenza Other Respir Viruses. 2018;12(3):307–13. pmid:29356350
  14. 14. Fallah MP, Skrip LA, Gertler S, Yamin D, Galvani AP. Quantifying poverty as a driver of Ebola transmission. PLoS Negl Trop Dis. 2015;9(12):e0004260. pmid:26720278
  15. 15. Fuller TL, Calvet G, Genaro Estevam C, Rafael Angelo J, Abiodun GJ, Halai U-A, et al. Behavioral, climatic, and environmental risk factors for Zika and Chikungunya virus infections in Rio de Janeiro, Brazil, 2015–16. PLoS ONE. 2017;12(11):e0188002. pmid:29145452
  16. 16. Bucchianeri GW. Is SARS a poor man’s disease? Socioeconomic status and risk factors for SARS transmission. Forum Health Econ Policy. 2010;13(2):1–29.
  17. 17. Peters PJ, Pontones P, Hoover KW, Patel MR, Galang RR, Shields J, et al. HIV infection linked to injection use of oxymorphone in Indiana, 2014–2015. N Engl J Med. 2016;375(3):229–39. pmid:27468059
  18. 18. Foster M, Ramachandran S, Myatt K, Donovan D, Bohm S, Fiedler J, et al. Hepatitis A virus outbreaks associated with drug use and homelessness—California, Kentucky, Michigan, and Utah, 2017. MMWR Morb Mortal Wkly Rep. 2018;67(43):1208–10. pmid:30383739
  19. 19. Gardy JL, Johnston JC, Sui SJH, Cook VJ, Shah L, Brodkin E, et al. Whole-genome sequencing and social-network analysis of a tuberculosis outbreak. N Engl J Med. 2011;364(8):730–9. pmid:21345102
  20. 20. Chandrasekhar R, Sloan C, Mitchel E, Ndi D, Alden N, Thomas A, et al. Social determinants of influenza hospitalization in the United States. Influenza Other Respir Viruses. 2017;11(6):479–88. pmid:28872776
  21. 21. Mayoral JM, Alonso J, Garín O, Herrador Z, Astray J, Baricot M, et al. Social factors related to the clinical severity of influenza cases in Spain during the A (H1N1) 2009 virus pandemic. BMC Public Health. 2013;13(1):118. pmid:23391376
  22. 22. Lowcock EC, Rosella LC, Foisy J, McGeer A, Crowcroft N. The social determinants of health and pandemic H1N1 2009 influenza severity. Am J Public Health. 2012;102(8):e51–8. pmid:22698024
  23. 23. Towne SD, Probst JC, Hardin JW, Bell BA, Glover S. Health & access to care among working-age lower income adults in the Great Recession: disparities across race and ethnicity and geospatial factors. Soc Sci Med. 2017;182:30–44. pmid:28411525
  24. 24. Antunes L, Mendonça D, Ribeiro AI, Maringe C, Rachet B. Deprivation-specific life tables using multivariable flexible modelling—trends from 2000–2002 to 2010–2012, Portugal. BMC Public Health. 2019;19(1):276. pmid:30845935
  25. 25. Agabiti N, Pirani M, Schifano P, Cesaroni G, Davoli M, Bisanti L, et al. Income level and chronic ambulatory care sensitive conditions in adults: a multicity population-based study in Italy. BMC Public Health. 2009;9(1):457. pmid:20003336
  26. 26. Wadhera RK, Wadhera P, Gaba P, Figueroa JF, Joynt Maddox KE, Yeh RW, et al. Variation in COVID-19 hospitalizations and deaths across New York City boroughs. JAMA. 2020;323(21):2192–5. pmid:32347898
  27. 27. Centers for Disease Control and Prevention, National Center for Health Statistics. National Health and Nutrition Examination Survey: NHANES 2017–2018. Hyattsville (MD): National Center for Health Statistics; 2020 [cited 2020 Dec 9]. Available from: https://wwwn.cdc.gov/nchs/nhanes/continuousnhanes/default.aspx?BeginYear=2017.
  28. 28. James PA, Oparil S, Carter BL, Cushman WC, Dennison-Himmelfarb C, Handler J, et al. 2014 evidence-based guideline for the management of high blood pressure in adults. JAMA. 2014;311(5):507. pmid:24352797
  29. 29. Crim MT, Yoon SS, Ortiz E, Wall HK, Schober S, Gillespie C, et al. National surveillance definitions for hypertension prevalence and control among adults. Circ Cardiovasc Qual Outcomes. 2012;5(3):343–51. pmid:22550130
  30. 30. American Diabetes Association. 2. Classification and diagnosis of diabetes: standards of medical care in diabetes—2020. Diabetes Care. 2020;43(Suppl 1):S14–31. pmid:31862745
  31. 31. Levey AS, Stevens LA, Schmid CH, Zhang Y, Castro AF, Feldman HI, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–12. pmid:19414839
  32. 32. Murphy D, McCulloch CE, Lin F, Banerjee T, Bragg-Gresham JL, Eberhardt MS, et al. Trends in prevalence of chronic kidney disease in the United States. Ann Intern Med. 2016;165(7):473–81. pmid:27479614
  33. 33. Olenick M, Flowers M, Diaz V. US veterans and their unique issues: enhancing health care professional awareness. Adv Med Educ Pract. 2015;6:635. pmid:26664252
  34. 34. Vespa JE. Those who served: America’s veterans from world war ii to the war on terror. Suitland (MD): US Census Bureau; 2020 [cited 2020 Dec 9]. Available from: https://www.census.gov/content/dam/Census/library/publications/2020/demo/acs-43.pdf.
  35. 35. Centers for Disease Control and Prevention. Daily updates of totals by week and state. Provisional death counts for coronavirus disease 2019 (COVID-19). Atlanta: Centers for Disease Control and Prevention; 2020 [cited 2020 Sep 8]. Available from: https://www.cdc.gov/nchs/nvss/vsrr/COVID19/.
  36. 36. Centers for Disease Control and Prevention. Health disparities: race and Hispanic origin. Provisional death counts for coronavirus disease 2019 (COVID-19). Atlanta: Centers for Disease Control and Prevention; 2020 [cited 2020 Sep 9]. Available from: https://www.cdc.gov/nchs/nvss/vsrr/covid19/health_disparities.htm.
  37. 37. Campbell A, Caul S. Deaths involving COVID-19, England and Wales: deaths occurring in May 2020. Newport: Office for National Statistics; 2020 [cited 2020 Dec 9]. Available from: https://www.ons.gov.uk/peoplepopulationandcommunity/birthsdeathsandmarriages/deaths/bulletins/deathsinvolvingcovid19englandandwales/deathsoccurringinmay2020.
  38. 38. Instituto de Salud Carlos III. [Report no. 32. Situation of COVID-19 in Spain on 21 May 2020.] Madrid: Instituto de Salud Carlos III; 2020 [cited 2020 Dec 9]. Available from: https://www.isciii.es/QueHacemos/Servicios/VigilanciaSaludPublicaRENAVE/EnfermedadesTransmisibles/Paginas/-COVID-19.-Informes-previos.aspx.
  39. 39. SARS-CoV-2 Surveillance Group. Characteristics of COVID-19 patients dying in Italy. Rome: Istituto Superiore di Sanità; 2020 [cited 2020 Dec 9]. Available from: https://www.epicentro.iss.it/en/coronavirus/sars-cov-2-analysis-of-deaths.
  40. 40. Santé Publique France. [COVID-19: epidemiological update for 3 September 2020]. Saint-Maurice: Santé Publique France; 2020 [cited 2020 Sep 10]. Available from: https://www.santepubliquefrance.fr/maladies-et-traumatismes/maladies-et-infections-respiratoires/infection-a-coronavirus/documents/bulletin-national/covid-19-point-epidemiologique-du-3-septembre-2020.
  41. 41. Novel Coronavirus Pneumonia Emergency Response Epidemiology Team. [The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) in China]. Zhonghua Liu Xing Bing Xue Za Zhi. 2020;41(2):145–51. pmid:32064853
  42. 42. Richardson S, Hirsch JS, Narasimhan M, Crawford JM, McGinn T, Davidson KW, et al. Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID-19 in the New York City area. JAMA. 2020;323(20):2052. pmid:32320003
  43. 43. Vahidy FS, Drews AL, Masud FN, Schwartz RL, Askary BB, Boom ML, et al. Characteristics and outcomes of COVID-19 patients during initial peak and resurgence in the Houston metropolitan area. JAMA. 2020;324(10):998–1000. pmid:32789492
  44. 44. Ferrari PA, Barbiero A. Simulating ordinal data. Multivariate Behav Res. 2012;47(4):566–89. pmid:26777670
  45. 45. R Core Team. R: a language and environment for statistical computing. Version 4.0.2. Vienna: R Foundation for Statistical Computing; 2019.
  46. 46. Barbiero A, Ferrari PA. GenOrd: simulation of discrete random variables with given correlation matrix and marginal distributions. Version 1.4.0. Comprehensive R Archive Network; 2015 [cited 2020 Dec 9]. Available from: https://cran.r-project.org/package=GenOrd.
  47. 47. Tomer A, Kane JW. How to protect essential workers during COVID-19. Washington (DC): Brookings; 2020 [cited 2020 Dec 9]. Available from: https://www.brookings.edu/research/how-to-protect-essential-workers-during-covid-19/.
  48. 48. Wu Z, Tang Y, Cheng Q. Diabetes increases the mortality of patients with COVID-19: a meta-analysis. Acta Diabetol. 2020;1:3. pmid:32583078
  49. 49. Gao C, Cai Y, Zhang K, Zhou L, Zhang Y, Zhang X, et al. Association of hypertension and antihypertensive treatment with COVID-19 mortality: a retrospective observational study. Eur Heart J. 2020;41(22):2058–66. pmid:32498076
  50. 50. Huang I, Lim MA, Pranata R. Diabetes mellitus is associated with increased mortality and severity of disease in COVID-19 pneumonia—a systematic review, meta-analysis, and meta-regression: diabetes and COVID-19. Diabetes Metab Syndr Clin Res Rev. 2020;14(4):395–403. pmid:32334395
  51. 51. Pranata R, Lim MA, Huang I, Raharjo SB, Lukito AA. Hypertension is associated with increased mortality and severity of disease in COVID-19 pneumonia: a systematic review, meta-analysis and meta-regression. J Renin Angiotensin Aldosterone Syst. 2020;21(2):1470320320926899. pmid:32408793
  52. 52. McCarthy CP, Murphy S, Jones-O’Connor M, Olshan DS, Khambhati JR, Rehman S, et al. Early clinical and sociodemographic experience with patients hospitalized with COVID-19 at a large American healthcare system. EClinicalMedicine. 2020;26:100504. pmid:32838244
  53. 53. Khalatbari-Soltani S, Cumming RG, Delpierre C, Kelly-Irving M. Importance of collecting data on socioeconomic determinants from the early stage of the COVID-19 outbreak onwards. J Epidemiol Community Health. 2020;84(8):620–3. pmid:32385126
  54. 54. Abrams EM, Szefler SJ. Comment COVID-19 and the impact of social determinants of health. Lancet Respir Med. 2020;8(7):659–61. pmid:32437646
  55. 55. Wu J, McCann A, Katz J, Peltier E. 36,000 missing deaths: tracking the true toll of the coronavirus crisis. The New York Times. 2020 Apr 23.
  56. 56. Brown E, Reinhard B, Thebault R. Which deaths count toward the covid-19 death toll? It depends on the state. The Washington Post. 2020 Apr 16 [cited 2020 Apr 23]. Available from: https://www.washingtonpost.com/investigations/which-deaths-count-toward-the-covid-19-death-toll-it-depends-on-the-state/2020/04/16/bca84ae0-7991-11ea-a130-df573469f094_story.html.