key: cord-0825546-hia1smo7 authors: Chakrabarti, Suman; Hamlet, Leigh C.; Kaminsky, Jessica; Subramanian, S. V. title: Association of Human Mobility Restrictions and Race/Ethnicity–Based, Sex-Based, and Income-Based Factors With Inequities in Well-being During the COVID-19 Pandemic in the United States date: 2021-04-07 journal: JAMA Netw Open DOI: 10.1001/jamanetworkopen.2021.7373 sha: 7eedfbeb1a0c72ac6185a2d2165a3aca311c6811 doc_id: 825546 cord_uid: hia1smo7 IMPORTANCE: An accurate understanding of the distributional implications of public health policies is critical for ensuring equitable responses to the COVID-19 pandemic and future public health threats. OBJECTIVE: To identify and quantify the association of race/ethnicity–based, sex-based, and income-based inequities of state-specific lockdowns with 6 well-being dimensions in the United States. DESIGN, SETTING, AND PARTICIPANTS: This pooled, repeated cross-sectional study used data from 14 187 762 households who participated in phase 1 of the population-representative US 2020 Household Pulse Survey (HPS). Households were invited to participate by email, text message, and/or telephone as many as 3 times. Data were collected via an online questionnaire from April 23 to July 21, 2020, and participants lived in all 50 US states and the District of Columbia. EXPOSURES: Indicators of race/ethnicity, sex, and income and their intersections. MAIN OUTCOMES AND MEASURES: Unemployment; food insufficiency; mental health problems; no medical care received for health problems; default on last month’s rent or mortgage; and class cancellations with no distance learning. Race/ethnicity, sex, income, and their intersections were used to measure distributional implications across historically marginalized populations; state-specific, time-varying population mobility was used to measure lockdown intensity. Logistic regression models with pooled repeated cross-sections were used to estimate risk of dichotomous outcomes by social group, adjusted for confounding variables. RESULTS: The 1 088 314 respondents (561 570 [51.6%; 95% CI, 51.4%-51.9%] women) were aged 18 to 88 years (mean [SD], 51.55 [15.74] years), and 826 039 (62.8%; 95% CI, 62.5%-63.1%) were non-Hispanic White individuals; 86 958 (12.5%; 95% CI, 12.4%-12.7%), African American individuals; 86 062 (15.2%; 95% CI, 15.0%-15.4%), Hispanic individuals; and 50 227 (5.6%; 95% CI, 5.5%-5.7%), Asian individuals. On average, every 10% reduction in mobility was associated with higher odds of unemployment (odds ratio [OR], 1.3; 95% CI, 1.2-1.4), food insufficiency (OR, 1.1; 95% CI, 1.1-1.2), mental health problems (OR, 1.04; 95% CI, 1.0-1.1), and class cancellations (OR, 1.1; 95% CI, 1.1-1.2). Across most dimensions compared with White men with high income, African American individuals with low income experienced the highest risks (eg, food insufficiency, men: OR, 3.3; 95% CI, 2.8-3.7; mental health problems, women: OR, 1.9; 95% CI, 1.8-2.1; medical care inaccessibility, women: OR, 1.7; 95% CI, 1.6-1.9; unemployment, men: OR, 2.8; 95% CI, 2.5-3.2; rent/mortgage defaults, men: OR, 5.7; 95% CI, 4.7-7.1). Other high-risk groups were Hispanic individuals (eg, unemployment, Hispanic men with low income: OR, 2.9; 95% CI, 2.5-3.4) and women with low income across all races/ethnicities (eg, medical care inaccessibility, non-Hispanic White women: OR, 1.8; 95% CI, 1.7-2.0). CONCLUSIONS AND RELEVANCE: In this cross-sectional study, African American and Hispanic individuals, women, and households with low income had higher odds of experiencing adverse outcomes associated with the COVID-19 pandemic and stay-at-home orders. Blanket public health policies ignoring existing distributions of risk to well-being may be associated with increased race/ethnicity–based, sex-based, and income-based inequities. This supplemental material has been provided by the authors to give readers additional information about their work. Pooling twelve waves of weekly HPS data, we construct a data set that is a repeated cross-section at the individual/household level and a representative panel at the weekstate level. In the HPS, each household is interviewed up to three times. As we treat the data as pooled repeated cross-sections, we do not need to consider the issue of loss to follow-up. We only analyze complete cases, and given that our sample size is large, our estimates are unlikely to be influenced substantially by this decision. Please see the HPS technical documentation for more information on how the Census Bureau handles missing data. 1 We merge these pooled HPS data with the IHME data on mobility, case rates, and death rates by week of survey and state of residence. Thus, in our final data set, all outcomes and covariates from the HPS vary at the individual/household level, while the mobility, case rates, and death rates vary at the state level, operating as ecological exposures (eTable3). For example, all individuals sampled from California in the first week of the HPS are assigned the same level of mobility, and their mobility changes weekly (twelve times corresponding to the twelve weekly surveys of the HPS). For individual/household ℎ in state in week of survey , we estimate multivariable logistic regression models as: There are differences in the levels of measurement of the exposures and outcomes. For food insufficiency, class cancellations, and default on rent/mortgage (eTable3) in which the respondent answers on behalf of the household and him/herself, it is reasonable to assume these outcomes operate at the household and individual levels. For example, we interpret any differences in the class cancellations outcome across sex as being between the children of female versus male respondents. However, given that the question on unemployment asks about the respondent or anyone in the household, and the individual may not be the one who is unemployed, there may be measurement error in the sex gradient when regressing unemployment on sex. As we do not except differential misclassification by sex, any biases for sex differences in unemployment are likely to drive coefficients towards the null. Relaxing the assumption of linearity for mobility dose response In a sensitivity analysis, we relax the assumption of a linear relationship between mobility restriction and outcomes. We recategorize mobility into five dummy variables: 10-19%, 20-29%, 30-39%, 40-49%, and ≥50% restrictions, with 0-9% as the reference category. The range of categories is appropriate given mobility restrictions vary between 0 and 60% in our sample. We fit regression models with these dummy variables instead mobility specification, given that extreme values have higher leverage in the linear specification. However, we interpret these results with caution, given the loss of information by binning the mobility data, which reduces variation and weakens the ability to detect nuanced relationships. Given that the HPS allows individuals to be resampled up to three times during the twelve weeks, we consider the potential issue of autocorrelation, which occurs when observations from the same individual are correlated. While autocorrelation is typically a problem in data sets with long time periods (time points per individual > 50), 3 we check that correlation between repeat observations is not significantly biasing our mobility estimates. We use a Generalized Estimating Equations (GEE) approach that explicitly models intra-individual correlations and corrects for biases stemming from possible autocorrelation or non-independence. 4 Beyond the individual level, another possible source of autocorrelation is due to the fact that our data are structured as repeated cross-sections at the state-level. To gauge the degree of this bias, we collapse our data set to the state-week level and run regressions with state-level mean estimates of the outcomes and covariates (N=612 for 51 states and twelve weeks). We run our main regression models in a panel-data setting with ordinary least squares, accounting for state fixed effects, autocorrelation, and heteroscedasticity. Solid lines represent predicted estimates. Grey shaded areas represent 95% confidence intervals. Mobility restriction represents within-state reductions in mobility considering week-to-week changes in mobility from normal levels. All models control for income, race/ethnicity, age, sex, education, marital status, numbers of individuals in the household, week of survey, state-level heterogeneity, and COVID-19 death and case rates. Estimates are corrected for correlation of observations among the same individual using a Generalized Estimating Equations model with an exchangeable correlation structure. Standard error estimates are robust. Outcome Definition Self-reported causality Unemployment (universe = population over the age of 18 years not retired and willing to work) Respondent or anyone in household (not retired and willing to work) experienced a loss of employment since March 13, 2020 because of (1) sick with coronavirus symptoms or (2) caring for someone with coronavirus symptoms or (3) concerned about getting or spreading the coronavirus or (4) employer experienced a reduction in business (including furlough) due to coronavirus pandemic or (5) laid off due to coronavirus pandemic or (6) employment closed temporarily due to the coronavirus pandemic or (7) employment went out of business due to the coronavirus pandemic? Food insufficiency (universe = all households) In the last 7 days, household sometimes or often did not have enough to eat or household's food security became worse after March 13. Classes cancelled (universe = households with school-aged children) For children in the household, due to coronavirus pandemic classes normally taught in person at school were cancelled without any replacement with online/distance learning or other means. No medical care (universe = all households with a valid response) In the last 4 weeks, respondent needed medical care for something other than coronavirus but did not get it because of the pandemic. Non-self-reported causality Default on rent or mortgage (universe = renters or household with a mortgage) Household did not pay last month's rent or mortgage on time. Mental health problems (universe = all households) In the last 7 days, respondent felt nervous/anxious (score 1(never) -4 (everyday)), could not stop worrying (score 1(never) -4 (everyday)), had little pleasure in doing things (score 1(never) -4 (everyday)), felt depressed/hopeless (score 1(never) -4 (everyday)). Scores from the 4 dimensions were summed and the sum ranged from 4 (good mental health) to 16 (bad mental health)). Respondents with scores ≥12 were considered to be those with mental health problems. Low-inc is the abbreviation for low-income households (< $35,000 per year), and low-mid inc is the abbreviation for lowmiddle-income households ($35,000-$75,000 per year). According the ACS Selected Economic Characteristics Table, 28.1% of households earn less than $35000 and 29.5% earn $35,000-$75,000 per year. b White is the abbreviation for non-Hispanic White respondents, and Black is the abbreviation for African American respondents. The population estimates by race are obtained from the ACS Demographic and Housing Estimates Table. c Eligible population for unemployment is assumed to be the civilian labor force provided in ACS Selected Economic Characteristics The ACS Demographic and Housing Estimates Table provides estimates for the age groups 0-4.9, 5-9.9, 10-14.9, and 15-19.9. We estimate the population of ages 3-19 as: (0-4.9)*0.4 + (5-9.9) + (10-14.9) + (15-19.9). We assume that the populations within one-year age intervals in the category 0-4.9 are approximately equally distributed. 7,8 All models control for income, race/ethnicity, age, sex, education, marital status, numbers of individuals in the household, week of survey, state fixed effects, and COVID-19 weekly death and case rates Design and Operation of the 2020 Household Pulse Survey Initial economic damage from the COVID-19 pandemic in the United States is more widespread across ages and geographies than initial mortality impacts An Effective Approach to the Repeated Cross-Sectional Design Generalized Estimating Equations in Longitudinal Data Analysis: A Review and Recent Developments COVID-19 resources. Institute for Health Metrics and Evaluation The World Bank. Mortality rate, under-5 (per 1,000 live births) -United States | Data The World Bank. Fertility rate, total (births per woman) -United States | Data US Census Bureau. Source of the Data and Accuracy of the Estimates for the 2020 Household Pulse Survey