ABSTRACT
Objectives We characterized informally employed US domestic workers’ (DWers) exposure to patterns of workplace hazards, as well as singular hazards, and examined associations with DWers’ work-related and general health.
Methods We analyzed cross-sectional data from the sole nationwide survey of informally employed US DWers with work-related hazards data, conducted in 14 cities (2011-2012; N=2,086). We characterized DWers’ exposures using four approaches: single exposures (n=19 hazards), composite exposure to hazards selected a priori, classification trees, and latent class analysis. We used city fixed effects regression to estimate the risk ratio (RR) of work-related back injury, work-related illness, and fair-to-poor self-rated health associated with exposure as defined by each approach.
Results Across all four approaches—net of individual, household, and occupational characteristics and city fixed effects—exposure to workplace hazards was associated with increased risk of the three health outcomes. For work-related back injury, the estimated RR associated with heavy lifting (the single hazard with the largest RR), exposure to all three hazards selected a priori (did heavy lifting, climbed to clean, worked long hours) versus none, exposure to the two hazards identified by classification trees (heavy lifting, verbally abused) versus “No heavy lifting,” and membership in the most-versus least-exposed latent class were, respectively, 3.4 (95% confidence interval [CI] 2.7 to 4.1); 6.5 (95% CI 4.8 to 8.7); 4.4 (95% CI 3.6 to 5.3), and 6.6 (95% CI 4.6 to 9.4).
Conclusions Measures of joint work-related exposures were more strongly associated than single exposures with informally employed US DWers’ health profiles.
What is already known on this topic Informally employed domestic workers in the US and internationally are frequently exposed to physical and social hazards at work, but only two studies have quantitatively assessed these workers’ exposures to joint patterns of hazards, and neither examined such patterns in relation to health.
What this study adds We characterized informally employed US domestic workers’ exposure to 19 single hazards and to combinations of these hazards, using three distinct approaches: composite exposure to hazards selected a priori, classification trees, and latent class analysis. Across all approaches to defining exposure, domestic workers exposed to worse joint patterns of workplace hazards, as well as to certain single hazards, experienced greater risk of work-related back injury, work-related illness, and fair-to-poor self-rated health.
How this study might affect research, practice, or policy Results underscore the importance of conceptualizing and operationalizing measures that capture domestic workers’ patterns of exposures. Moreover, results support the use of a latent class approach for identifying potential subgroups of workers unduly burdened and—across multiple health metrics—harmed by employer practices.
INTRODUCTION
In 2008, Krieger and colleagues formulated the “inverse hazard law,” stating that “the accumulation of health hazards tends to vary inversely with the power and resources of the populations affected.”[1] Developed in relation to their novel study of the combined health impacts of social and physical workplace hazards in a sample of low-wage US workers, this law highlights how exposures cluster together and jointly impact health, producing social inequities in health within designated occupational categories as well as occupational health inequities.[2-3] The inverse hazard law, thus, adds impetus for research examining the health impacts of— and comparing approaches for characterizing—the joint patterns of exposure experienced by workers.
Despite this, no quantitative studies in the US or internationally have gone beyond examining associations between exposure to single workplace hazards and health among informally employed domestic workers (DWers). As of 2019, at least 800,000 DWers were employed informally in housecleaning, childcare, or adult care work by private households, rather than by third-party agencies or government entities.[4-5] These DWers are subjected to numerous workplace hazards, including wage theft, verbal abuse, racial discrimination, heavy lifting, and more. Although these hazards are not unique to DWers, the overall frequency and social patterning of DWers’ exposures is likely distinct—reflecting their historical and contemporaneous exclusion from laws protecting workers’ rights and the wide range of employer practices and worker social characteristics within this occupational group.[6]
Research examining the health of DWers experiencing combinations of workplace hazards is urgently needed. The sparse literature quantitatively examining workplace exposures and health among DWers in the US and internationally has documented cross-sectional associations between various singular hazards (e.g., lack of rest breaks, cleaning chemicals, physical abuse, late wages) and work-related (e.g., back injury) and general health (e.g., self-rated health).[7-11] Our recent study using a novel latent class approach found that informally employed US DWers experience distinct patterns of workplace hazards, motivating further analyses of the health implications of such exposure patterns.[12] Finally, although a growing body of research indicates the conceptual and empirical importance of examining the health implications of workers’ joint exposures using a range of methods,[1, 13-14] few studies have employed multiple methods simultaneously.[15] Evidence regarding the utility of accessible composite exposure approaches and more complex model-based methods for characterizing workers’ joint exposures, in comparison to a single hazards approach, could inform future research on the health impacts of domestic work (DW) and related occupations. Moreover, evidence regarding which groups of DWers are potentially unduly harmed by employer practices could further inform ongoing work to enforce existing, and create new, DWer protections.
In this study, we used cross-sectional data from the sole nationwide survey of informally employed US DWers with work-related hazards data—conducted in 2011-2012 by the National Domestic Workers Alliance (NDWA), the University of Illinois Chicago Center for Urban Economic Development (UIC CUED), and the DataCenter—to characterize the health of DWers exposed to different combinations of workplace hazards. Informed by ecosocial theory and its emphasis on embodied integration of exposures,[3, 16, 17] we characterized DWers’ exposures using four approaches—single exposures, composite exposure to hazards selected a priori, classification trees, and latent class analysis—and used city fixed effects regression models to examine associations between exposure and work-related and general health.
METHODS
Study population
We analyzed cross-sectional data from the sole nationwide survey of informally employed US DWers with workplace hazards data (N=2,086), conducted from June 2011 to February 2012 in 14 cities (Table 1). NDWA, UIC CUED, and DataCenter investigators selected these cities to “represent every region of the country,”[18] ensure a large sampling frame of domestic workers, and make sure DWers’ community organizations in each city could serve as research partners. City-specific enrollment targets were derived from American Community Survey 2005-2009 five-year estimates (based on occupation, racialized group, and nativity)[19] and used to help ensure the sample reflected each city’s DWer labor force. The UIC Institutional Review Board approved this survey.
Investigators employed participatory methods whereby organizers from 34 NDWA-affiliated community organizations and 190 DWers collaborated on survey design, fielding, and analysis. Because of the secluded, dispersed nature of DW, surveyors recruited participants at public places (i.e., parks, transportation hubs, churches, shopping centers) and used chain-referral of potential participants. Surveys were conducted in nine languages (Cantonese, English, Haitian Creole, Mandarin, Nepali, Polish, Portuguese, Spanish, Tagalog). DWers were eligible to participate if they worked in a private home during the previous week for ≥6 hours as a housecleaning, childcare, or adult care worker; were paid directly by an employing family member; were ≥18 years old; lived in a selected city; and were not members of organizations that advocated for workers’ right—so as to avoid bias potentially related to “participants who had more knowledge about exercising their employment rights.”[18] The Harvard T.H. Chan School of Public Health Institutional Review Board determined our secondary analyses of the NDWA-UIC CUED data was Not Human Subjects Research (IRB21-0855).
Workplace hazards
In the NDWA-UIC CUED survey, participants self-reported exposure (yes/no) to six economic (e.g., paid late), nine social (e.g., verbally abused), and six occupational (e.g., heavy lifting) hazards (n=21 variables) at any DW job in the previous 12 months (Table 1). The option to respond “don’t know” was available only for the “worked with toxic cleaning supplies” variable. We coded as missing the few participants who responded using this option (n=40 [1.9%]). Although the NDWA-UIC CUED survey asked about exposure to workplace sexual harassment and sexual assault, we excluded these variables from the analysis due to severe underreporting,[20,21] thereby reducing the number of exposure variables to 19 (for more details, see Supplemental Methods).
Health outcomes
We used self-reported data on whether DWers reported experiencing back injury (“Back injuries, including pulled back muscles”) and, separately, illness (“Contracted an illness, such as the flu”) while working at any job as a DWer in the previous 12 months. The survey also inquired about participants’ self-rated health in the previous 12 months (1=Excellent, 2=Very good, 3=Good, 4=Fair, 5=Poor), which we categorized as fair-to-poor self-rated health (FPSRH; 1=Fair-to-poor, 0=Excellent, very good, or good). These three outcomes were commonly experienced by DWers in the sample (back injury: 26.8%; illness: 32.6%; FPSRH: 37.1%).
Covariates
The following covariates were included in all health outcome models: racialized group (collected as: “White,” “Latina/o,” “Black,” or “Asian or other,” which participants could, but in no cases did, multiply select), citizenship and immigration status, gender, age, education, main DW occupation, live-in status, years worked as a DWer in the US, hours worked last week for one’s main DW employer, number of DW employers, contract status, number of people relying on the participant’s financial support, income earner status, and household economic insecurity. Covariate selection was informed by our use of ecosocial theory,[3, 16, 17] prior research identifying predictors of workplace exposures and/or health among DWers, observed variation in DWer exposures by covariates, and empirical considerations regarding collinearity and model complexity (for more details, see Supplemental Methods).
Statistical analyses
We used R (Version 4.1.3, R Foundation for Statistical Computing, Vienna, Austria) to conduct all statistical analyses except latent class analysis, which was conducted in LatentGOLD (Version 6.0, Statistical Innovations, Inc., Arlington, VA).[22, 23] All tests for statistical significance were two-sided. We first descriptively quantified the distribution of—and examined relationships between—variables in the dataset, including regarding missingness (e.g., Figure S1).
We then defined exposure to workplace hazards using the four approaches described below and fit models to estimate associations between exposure and each health outcome (Table 2). All health outcome models adjusted for the aforementioned covariates, and used city fixed effects to account for all unobserved city-specific, time-invariant characteristics that may confound exposure-outcome associations.[24] Hausman tests (i.e., tests of the adequacy of city random effects models) indicated the need for city fixed effects to address endogeneity in models for each of the three health outcomes.[25]
Because all three health outcomes were non-rare in the sample, models were fit using log Poisson regression, with parameter estimates interpretable as log risk ratios (RR).[26] All health outcome models employed multiply imputed covariates and health outcome data. We used multiple imputation using chained equations to create 10 imputed datasets under the Missing-At-Random (MAR) assumption.[27] We included in the imputation model all variables to be included in the health outcome models and additional variables that may have otherwise explained nonresponse or variance in the data (for more details, see Supplemental Methods). All 95% confidence intervals are based on cluster robust standard errors, clustered at the city level.
Single exposure
We first estimated the RR of each health outcome associated with exposure to the 19 individual hazards by including them singly in separate models. Because this involved testing multiple dependent hypotheses, we used Benjamini and Yekutieli’s (2001) approach to control the False Discovery Rate.[28]
A priori composite exposure
To examine how exposure to multiple hazards was associated with health—without relying on complex exposure modeling strategies—we selected the hazards we expected would be most strongly associated with each health outcome. For work-related back injury, work-related illness, and FPSRH, respectively, these were: 1) heavy lifting, climbing to clean, and working long hours with no breaks, 2) caring for someone with a contagious illness, and working long hours with no breaks, and 3) working long hours with no breaks, working with toxic cleaning supplies, and being threatened (for more details, see Supplemental Methods). We constructed composite exposure variables reflecting whether DWers were exposed to none, some, or all of these selected hazards. We estimated the RR of each health outcome associated with experiencing some (versus none) or all (versus none) of the selected hazards for that outcome.
Classification trees
We employed classification trees to empirically identify the combinations of workplace exposures most strongly associated with each outcome. Classification trees are a supervised learning method that identifies groups of individuals with similar outcomes based on their patterns of responses to multiple observed variables.[14, 29] This method identifies the single hazard that best predicts the outcome, splits the data into two groups according to exposure to that hazard, and recursively repeats this process with each identified subgroup until splitting by additional variables no longer improves the results (for more details, see Supplemental Methods). Thus, each resulting classification tree was composed of “branches” representing exposures, on the basis of which DWers were grouped into “leaves,” representing the final groupings of DWers. Using these final groupings, we estimated the RR of each outcome associated with experiencing each exposure pattern (versus the least exposed group) for that outcome.
Latent class analysis
We used latent class analysis (LCA) to identify groups of DWers with distinct patterns of exposure to workplace hazards, regardless of the health outcomes they reported. LCA is an unsupervised learning method (i.e., no outcome variables) that identifies groups of individuals, called latent classes, based on the clustering of their responses to multiple observed variables.[30, 31] This method uses the statistical associations between observed variables to learn about the structure of an assumed underlying categorical latent variable with an unknown number of classes. We determined the number of latent classes by comparing the fit, classification certainty, identifiability, and interpretability of multiple fitted LCA models with differing number of latent classes.[12] We initially included all 19 workplace hazard variables in these models. However, violations of the LCA conditional independence assumption—which states that membership in the latent classes explains all associations between variables—can arise due to redundancies between observed variables and lead to the selection of too many latent classes.[32] Thus, we monitored for violations of this assumption, considered combining or dropping redundant hazard variables as needed, and reran the LCA models with this reduced set of variables before deciding how many latent classes to retain. We estimated the RR of each health outcome associated with membership in each latent class (versus the least exposed class). See Wright et al [12] and the Supplemental Methods for more details on our latent class approach.
RESULTS
Descriptive characteristics of the 2,086 informally employed DWers surveyed are presented in Table 1. Regarding single workplace hazards, more than 15 percent of DWers reported being verbally abused or yelled at (15.8%), pressured to work more than their scheduled hours (18.8%), paid late (23.5%), required to work on their knees (35.5%), working long hours with no breaks (35.9%), doing heavy lifting (40.0%), climbing to clean (41.9%), and working with toxic cleaning supplies (46.7%). Exposure to single hazards varied by racialized group, citizenship and immigration status, and main DW occupation. The percent of DWers who reported a work-related back injury, work-related illness, and FPSRH, respectively, equaled 26.8%, 32.6%, and 37.1%.
Results from the classification trees are presented in Figure 1, Figure S2, and Figure S3. Classification trees identified “heavy lifting” and “verbally abused” as the variables most strongly associated with work-related back injury, yielding exposure categories for “No heavy lifting” (n=1,247 [60.0%]), “Heavy lifting, but no verbal abuse” (n=611 [29.4%]), and “Heavy lifting and verbal abuse” (n=222 [10.7%]). For work-related illness, classification trees yielded exposure categories for “Did not work long hours with no breaks” (n=1,338 [64.2%]), “Worked long hours, but did no contagious illness care” (n=588 [28.2%]), and “Worked long hours and did contagious illness care (n=157 [7.5%]). For FPSRH, classification trees identified “Did not climb to clean” (n=1,185 [57.9%]), “Climbed to clean, but was not called insulting names or racial slurs” (n=747 [36.5%]), and “Climbed to clean and was called insulting names or racial slurs” (n=114 [5.6%]). Classification trees yielded similar exposure categories to a priori composite variables (Table 2). For example, a priori selection and classification trees both identified “contagious illness care” and “worked long hours with no breaks” as the most important variables for predicting work-related illness.
Latent class analysis suggested that four latent classes—for which we developed concise summary labels—best represented the patterns of workplace hazard exposure in the data (Figure 1, Table S1). Members of Class 1, which we labeled as “Low hazard domestic work,” had the lowest probability of exposure to all 10 hazards included in the final model (all item response probabilities [IRP], representing the probability of exposure among members of a given class, <0.2) (n=839 [40.2%]). Members of Class 2, labeled “Demanding care work” (n=304 [14.6%]), had moderate probability of exposure to pay violations (IRP=0.42), heavy lifting (IRP=0.45), working long hours with no breaks (IRP=0.57), and contagious illness care (IRP=0.38). Members of Class 3, labeled “Strenuous cleaning work” (n=596 [28.6%]), had high probability of exposure to several cleaning-related occupational hazards (e.g., IRPclimbed to clean=0.87; IRPtoxic cleaning supplies=0.76). Members of Class 4, labeled “Hazardous domestic work” (n=347 [16.6%]), had the highest probability of exposure to all but one hazard (contagious illness care).
Results of city fixed effects models indicated, first, that—across all four approaches to exposure definition—exposure to workplace hazards was associated with the risk of each health outcome, after adjusting for selected covariates (Figure 2). Across all three health outcomes, the single exposure approach produced the smallest (though often still notable and statistically significant, after adjusting for multiple testing [Figure S4]) estimates of the RR, whereas the a priori composite exposure, classification tree, and latent class analysis approaches performed comparably to one another. For example, for work-related back injury, the estimated RR associated with heavy lifting (the single exposure with the largest RR), exposure to all three hazards selected a priori (did heavy lifting, climbed to clean, worked long hours) versus none, exposure to the two hazards identified by classification trees (heavy lifting, verbally abused) versus “No heavy lifting,” and membership in the “hazardous domestic work” latent class versus the least-exposed “low hazard domestic work” class were, respectively, 3.4 (95% CI 2.7 to 4.1); 6.5 (95% CI 4.8 to 8.7); 4.4 (95% CI 3.6 to 5.3), and 6.6 (95% CI 4.6 to 9.4). Additional results from health outcome models are presented in Table S2 (complete results from fully-adjusted models using the a priori, classification tree, and latent class approaches), Table S3 (sensitivity analyses for unmeasured confounding), and Figure S5 (stepwise inclusion of covariates in models). Descriptive figures further characterized the individual, household, and occupational characteristics of DWers in each exposure category, as defined according to each approach (Figure S6).
DISCUSSION
Results from our novel, descriptive study indicate that—net of key individual, household, and occupational characteristics as well as city fixed effects—informally employed US DWers exposed to different combinations of workplace hazards experienced divergent health profiles. Compared to the three approaches used to characterize workers’ joint exposures, the single exposures approach underestimated the risk associated with workers’ lived experiences of multiple hazards. Compared to the a priori composite exposures, the classification tree approach estimated slightly smaller RR, reflecting larger exposure referent groups with greater prevalence of the outcome. Our findings also suggest that an a priori composite exposure approach—when informed explicitly by theory and subject matter expertise—can provide an accessible tool for defining exposure categories that meaningfully capture risk of a given outcome among DWers. Finally, despite being the only approach that did not incorporate information about DWers’ health outcomes in constructing exposure categories, the latent class approach produced RR estimates, across all three outcomes, of similar magnitude to those observed using the a priori and classification tree approaches.
The strengths and limitations of our study must be considered when interpreting results. First, although NDWA-UIC CUED investigators took extensive steps to collect a rich set of covariates and ensure their sample reflected city-specific DW labor force characteristics, we analyzed cross-sectional data from a non-random sample of informally employed DWers. Risk ratios may be underestimated as a result of the healthy worker survivor effect, in which both healthier and more exposed workers are more likely to remain in the active DW labor force and, thus, participate in the survey.[33] Moreover, pre-domestic work employment health status was not measured and may also lead to an underestimation of results (i.e., healthy worker hire effect). Second, although we excluded the only two exposure variables (sexual harassment, sexual assault) we assessed as severely underreported, measurement error of other self-reported hazards, which may vary by participant characteristics, could still bias results (i.e., overestimating the number of latent classes; incorrectly selecting the splitting variables for classification trees; underestimating exposure-outcome associations). Third, results may not be generalizable to contemporary informally employed US DWers. In our recent study using the same data,[12] we partially investigated this concern. We empirically examined—and observing similarities between—the individual, household, and employment characteristics of DWers in this unique 2011-2012 sample relative to those of DWers in the only other nationwide survey of formally and informally employed DWers, which lacked workplace hazards data and was conducted in 2020-2021, solely among Spanish-speaking DWers.[34]
Three lines of evidence lend credibility to our findings. First, the limited literature quantitatively examining exposure to single workplace hazards and health among DWers, in the US and internationally, documents cross-sectional associations of similar magnitude to those observed in this study using the single hazards approach. [7-11] For example, in a 2007 nationally representative sample of formally employed (i.e., agency-employed) home health aides, the odds of having a back injury were 4.64 (95% CI 1.6 to 14.0) times higher among aides who did versus did not report the need for additional ergonomic equipment, net of education, income, racialized group, and job training.[7] Second, a robust qualitative literature on the structure and conditions of DW describes similar patterns of exposures and related health concerns to those identified in our classification tree and latent class approaches.[35-38] For example, in 2007 Hondagneu-Sotelo [35] and in 2002 Romero [36] each reported how cleaning-related hazards cluster together among US DWers primarily responsible for housecleaning, often leading to back pain and injury. Finally, in addition to our recent study that developed the latent class approach further examined here,[12] the only other study we know of that quantitatively analyzed DWers’ joint exposures—which used multiple correspondence and clustering methods but did not examine health outcomes—identified three clusters of workplace abuse among a sample of Portuguese DWers.[39] More broadly, patterns of exposure to multiple hazards such as those observed in our study, and their association with multiple health outcomes, is akin to what other low-wage workers, in the US and internationally, experience at work.[1, 2, 40]
This study addresses an important gap in the literature by providing novel quantitative evidence regarding the work-related and general health of DWers experiencing distinct combinations of workplace hazards and showing the value of using diverse methods to quantify these joint exposure patterns. Taken together, our findings support the use of a latent class approach for identifying potential subgroups of workers unduly burdened and harmed, across multiple health metrics, by employer practices. More broadly, our results underscore the importance of conceptualizing and operationalizing measures that capture the range and patterns of DWers’ exposures. Future research on DWers should analyze joint patterns of exposure, examine the causal role of policy exclusions and employer practices in patterning DWers’ health (including collecting the contemporary, nationally representative, and longitudinal data needed to do so), and investigate both social inequities in health among DWers and also between DWers and more privileged workers.
Data Availability
The 2011-2012 NDWA-UIC CUED survey data underlying this article cannot be shared publicly per the policies of NDWA to protect the privacy of individuals whose information was collected in these surveys. These terms of use are stipulated in the requirements of the authors' Data Use Agreement with the University of Illinois Chicago.
CONTRIBUTORSHIP STATEMENT
Emily Wright initiated the study; led writing of the manuscript, conceptualizing the study design, acquiring the data, conducting the analyses, and interpreting the results; approved submission of the version to be published; and agrees to be accountable for all aspects of the work.
Jarvis T. Chen contributed to: conceptualizing the study, designing the analyses, interpreting the results, and preparing the manuscript; approved submission of the version to be published; and agrees to be accountable for all aspects of the work.
Jason Beckfield contributed to: conceptualizing the study, designing the analyses, interpreting the results, and preparing the manuscript; approved submission of the version to be published; and agrees to be accountable for all aspects of the work.
Nik Theodore led the original study from which the secondary 2011-2012 NDWA-UIC CUED data are derived; assisted with data acquisition; contributed to interpreting the results and preparing the manuscript; approved submission of the version to be published; and agrees to be accountable for all aspects of the work.
Nancy Krieger contributed to: conceptualizing the study, designing the analyses, interpreting the results, and preparing the manuscript; approved submission of the version to be published; and agrees to be accountable for all aspects of the work.
COMPETING INTERESTS
None of the authors have any conflicts of interests to report.
DATA AVAILABILITY STATEMENT
The 2011-2012 NDWA-UIC CUED survey data underlying this article cannot be shared publicly per the policies of NDWA to protect the privacy of individuals whose information was collected in these surveys. These terms of use are stipulated in the requirements of the authors’ Data Use Agreement with the University of Illinois Chicago.
ETHICAL APPROVAL
This study was determined to be Not Human Subjects Research (IRB21-0855) by the Harvard T.H. Chan School of Public Health Institutional Review Board.
FUNDING
The research reported in this article was not supported by any external funding.
ACKNOWLEDGMENTS
The authors acknowledge Paulina López González for her non-author contributions to interpreting the results of the latent class analysis for this manuscript, as well as her contributions (as a co-author) to our recently published manuscript in which we developed the latent class approach employed in this analysis. The authors acknowledge the National Domestic Workers Alliance (NDWA) and the University of Illinois Chicago Center for Urban Economic Development (UIC CUED) for providing the 2011-2012 NDWA-UIC CUED survey data. Nik Theodore thanks the Ford Foundation, Open Society Foundations, and Alexander Soros Foundation for funding the initial study that is the source of the 2011-2012 NDWA-UIC CUED survey data for this secondary analysis.