Abstract
Purpose
To identify key social and structural factors associated with hazardous alcohol and problematic drug use among sexual and gender minority (SGM) adults.
Methods
We conducted a cross-sectional secondary analysis of 2016 to 2019 data from the Generations Study, a probability sample of 1518 SGM adults. Hazardous alcohol and problematic drug use were measured using the Alcohol Use Disorder Identification Test-Consumption (AUDIT-C) and Drug Use Disorder Identification Test (DUDIT). Candidate predictors included sociodemographic characteristics, minority stressors, adverse childhood experiences, health indicators, and social support. Interpretable machine-learning models with SHapley Additive exPlanations characterized multilevel correlates of 3-level risk categories.
Results
Hazardous alcohol use was linked to past smoking, lower social support, cancer history, housing instability, and identity-related stressors. Problematic drug use was linked to conversion-therapy exposure, housing discrimination, legal encounters, and mental-health indicators; prediction was more accurate for drug- than alcohol-use risk.
Conclusion
Structural vulnerability, discrimination, and identity-related stressors emerged as intervention targets to reduce substance-use disparities in SGM adults.
Plain Language Summary Title
Plain Language Summary
This study looked at what social and life experiences put LGBTQ+ adults in the United States at higher risk for alcohol and drug problems. We analyzed survey data from the Generations Study 2016 to 2019 which included 1518 sexual and gender minority adults (eg, gay, lesbian, transgender, queer, and bisexual).
Alcohol use was measured by how often and how much people drank, while drug use was measured by the frequency and impact of nonalcohol substance use. We also considered many possible influences, such as age, income, childhood experiences, discrimination, health, and support from family and friends.
For alcohol use, the most important factors were past smoking, low social support, cancer, unstable housing, and identity-related stress. For drug use, key factors included exposure to conversion therapy, housing discrimination, legal problems, and mental health concerns such as self-harm. Our models were better at predicting drug use than alcohol use.
These findings show that alcohol and drug use among LGBTQ+ people is shaped by personal behavior and also by stress, discrimination, and unequal access to stable housing, health, and support. Policies and programs to reduce substance use need to go beyond individual choices and address these broader social and structural conditions. Using transparent machine learning methods, we were able to identify complex patterns that point to where interventions and resources are most needed.
Keywords
Introduction
Sexual and gender minority (SGM) individuals, including those who identify as lesbian, gay, bisexual, transgender, queer, nonbinary, and other nonheteronormative identities, experience disproportionately high rates of substance use, including hazardous alcohol and drug consumption, relative to cisgender and heterosexual populations.1–3 These disparities are not solely attributable to individual-level behaviors or choices but are rooted in broader systems of oppression, discrimination, and marginalization that operate across multiple socioecological levels.4–6
The Minority Stress Model7–9 posits that stigma-related stressors (eg, discrimination, victimization, anticipated rejection, and internalized stigma) accumulate over time, contributing to adverse mental health and coping behaviors, including substance use. Intersectionality Theory10,11 adds that multiple marginalized positions (eg, race, socioeconomic status, gender identity) intersect to produce compounded disadvantage.11–14 The Social Determinants of Health (SoDH) framework15–17 further emphasizes how structural conditions (eg, housing instability, discrimination, and criminal-legal exposure) are associated with risk environments beyond individual agency. Together, these lenses underscore how structural inequities and identity-based stressors interact to impact substance-use risk.
However, most existing research has relied on statistical methods that assume linear relationships, independence among variables, and prespecified hypotheses. 18 These methods may fail to capture the complexity of the factors influencing substance use, particularly in populations facing layered social vulnerabilities. There remains a critical need for analytic approaches that can accommodate the high dimensionality and nonlinearity inherent in social and behavioral data.19,20 Alternatively, machine learning (ML) algorithms can model complex interactions among many variables without relying on prior assumptions about their relationships. Additionally, recent advances in interpretable ML, such as SHapley Additive exPlanations (SHAP), 21 allow researchers to retain insight into variable importance and effect directionality, addressing longstanding criticisms of ML's “black box” nature.19,22–24 In the context of SGM health research, ML may uncover previously underrecognized combinations of risk and protective factors spanning individual, interpersonal, and structural domains that influence substance use.
Therefore, we applied interpretable ML techniques that were guided by the Minority Stress Model and the SoDH Framework to: (1) prioritize multilevel correlates of hazardous alcohol use and problematic drug use among SGM adults and (2) translate these findings to actionable screening, clinical workflow, and policy targets. This approach preserves theoretical interpretability (via SHAP value) while accommodating high-dimensional, nonlinear patterns to bridge explanatory theory and decision support. Because alcohol and drug use can elevate HIV risk through decreased PrEP uptake/adherence, condomless sex, and cooccurring mental health burden, identifying upstream social/structural drivers in SGM communities also informs HIV prevention strategies.
Methods
Study Design and Participants
This study employed a cross-sectional secondary data analysis of baseline data (collected 2016-2019) from the Generations Study, which was designed to examine the health, stress, and well-being of SGM adults. 25 The study employed a national probability sample drawn from the Gallup Daily Tracking Survey, which included respondents aged 18 and older from all 50 U.S. states and the District of Columbia. Participants were screened for eligibility using a 2-phase process. In the first phase, individuals were identified via a self-reported question on sexual orientation and gender identity. Then, eligible participants were invited to complete a self-administered survey, either online or by mail. The final dataset included responses from 1518 individuals, with an oversample of Black and Latino participants to ensure adequate representation of these groups. Recruitment procedures ensured diversity in race/ethnicity, geography, and socioeconomic status. Participants provided informed consent and received compensation for their time. 25 This is a cross-sectional secondary analysis that is exploratory in nature and aims to develop and evaluate predictive models for substance-use risk stratification (alcohol and drug use). Reporting follows the STROBE guidelines (Supplemental Table S1).
Outcomes
The primary outcome variables in this study were alcohol use and drug use, assessed using validated instruments. Alcohol use was measured with the Alcohol Use Disorder Identification Test-Consumption (AUDIT-C), a widely used 3-item scale that evaluates hazardous drinking behaviors and potential alcohol-use disorders.26,27 Scores ranged from 0 to 12, with higher scores indicating greater levels of risky alcohol consumption (Cronbach's alpha = .790). For the AUDIT scale, scores of 0 to 3 were classified as low risk, 4 to 6 as moderate risk, and ≥7 as high risk. Drug use was assessed using the Drug Use Disorder Identification Test (DUDIT), an 11-item scale designed to identify individuals with drug-related problems.28,29 DUDIT scores ranged from 0 to 44, with thresholds distinguishing substance abuse from dependence (Cronbach's alpha = .856). For the DUDIT scale, scores of 1 to 10 were classified as low risk, 11 to 15 as moderate risk, and ≥16 as high risk. For both alcohol and drug use, we interpret the 3-level categories (low, moderate, high risk) as reflecting increasing severity of hazardous alcohol use (AUDIT-C) and problematic drug use (DUDIT).
Potential Variables Selected for Model Training
While the outcome variables were the focus of the analysis, variables were selected from a comprehensive dataset spanning multiple domains, including stressors, health outcomes, and healthcare utilization. We selected a series of potential variables informed by the Minority Stress Model, Intersectionality Theory, and the SDoH framework. These variables reflected multiple domains relevant to SGM health, including: (1) sociodemographic factors: age, gender identity, race/ethnicity, education, income; (2) minority stress indicators: enacted discrimination (lifetime and past-year events), anticipated rejection/expectations of rejection, internalized stigma, and identity concealment/outness milestones; we also included conversion-therapy exposure and housing/healthcare discrimination (binary indicators); (3) adverse life experiences: adverse childhood experiences (ACE), family rejection, housing instability, incarceration history; (4) mental and physical health: depressive symptoms, chronic illness (eg, kidney disease), history of self-injury, access to healthcare; and (5) social support and community connection: support from family or friends, identity affirmation (total score; higher = more support), connection to LGBTQ+ communities. We excluded any variables that directly overlapped with items used in the AUDIT-C or DUDIT to avoid artificially inflating associations. To prevent information leakage, we excluded items that directly overlapped with AUDIT-C or DUDIT content. Variables were not prespecified; instead, we used feature-selection algorithms (described below) to identify the most influential predictors of alcohol- and drug-use risk.30,31
Data Preparation
Data preprocessing included excluding variables with more than 20% missing data. Remaining missing values (<20%) were imputed using k-nearest neighbor imputation to preserve sample size while maintaining plausible estimates. This strategy balanced data completeness and model reliability by removing variables with excessive missingness, thereby minimizing bias and improving interpretability.32,33 Questions directly related to alcohol or drug use were excluded from the variables to prevent information leakage. Noninformative columns (eg, study IDs) were removed. Highly correlated variables were screened for multicollinearity, and 1 variable from each correlated pair (Pearson r > .80) was dropped. Categorical variables were binarized, and continuous variables were standardized for comparability.
Statistical Analysis
All analyses were conducted in Python 3.9. The dataset was randomly partitioned into training (80%) and testing (20%) sets, stratified by outcome class (low, moderate, high risk). Within the training set, a 20% validation split was used for hyperparameter tuning and early stopping. All preprocessing procedures (imputation, encoding, scaling, and correlation filtering) were fit on the training data and then applied to validation and test sets to prevent information leakage.
To identify stable and robust predictors, we implemented a feature-selection strategy based on a weighted voting scheme that aggregates feature importance from Elastic Net (EN), Random Forest (RF), and XGBoost (XGB) models. EN provides a linear, sparse representation that handles multicollinearity, whereas RF and XGB capture nonlinearities and higher-order interactions. 34 For each model, feature importance scores were min–max normalized to the 0 to 1 range and combined as Icombo = wenIen + wrfIrf + wxgbIxgb (where Icombo, Ien,Irf, Ixgb are the feature-importance scores from each model; and wen, wrf, wxgb are the nonnegative weights assigned to each model that sum to 1). We evaluated a grid of weight combinations (0.0-0.9 in 0.1 increments) and, for each combination, selected the top 30 features (k = 30), chosen a priori based on the elbow observed in the cumulative importance/predictive plateau. Features were then ranked by selection frequency across weightings and stratified cross-validation folds, and the final feature set was chosen to maximize stability.
Using the selected feature set, we trained 4 multiclass classifiers—XGB, RF, support vector machine, and a shallow feed-forward neural network (NN)—to predict 3-level risk categories (low, moderate, high) for both alcohol (AUDIT-C) and drug (DUDIT) outcomes. For analyses of the full Wave 1 sample (original plus enhancement), we incorporated the Wave 1 probability weight (weight_full) provided by the Generations Study. After preprocessing (imputation, encoding, scaling, correlation filtering), we applied the stability-oriented feature selection procedure described above (weighted voting across EN, RF, and XGB importance scores) and then trained all models using stratified train/validation/test splits to prevent information leakage. The NN comprised 3 fully connected layers with rectified linear unit activations and a SoftMax output layer; to reduce overfitting, we used dropout (rate = 0.30) on hidden layers and early stopping based on validation loss (Adam optimizer, learning rate 1 × 10−3, batch size 64, up to 100 epochs, patience = 10). Model interpretability was assessed using SHAP values, which quantify the direction and relative contribution of each feature to predicted risk. For descriptive comparisons of demographic and key clinical characteristics across alcohol- and drug-use risk categories, we used analysis of variance with appropriate post-hoc tests.
Results
Table 1 shows the distribution of alcohol and drug use across various sociodemographic factors among the study participants. Age was significantly associated with both alcohol and drug use, with individuals in the high-risk drug use group being notably younger (mean age = 33.28 years, standard deviation [SD] = 14.06) compared to those in the low-risk group (mean age = 38.94 years, SD = 15.83, η2 = 0.01), suggesting a small effect. Gender identity had a minor impact on alcohol use (Cramer's V = 0.05) and a slightly larger effect on drug use (Cramer's V = 0.06), with men demonstrating higher rates of drug use compared to women and nonbinary individuals. Similarly, racial and ethnic disparities were evident in substance use patterns. White participants had the highest representation in the high-risk alcohol use group (Cramer's V = 0.03), while Black/African American participants had the lowest representation in the high-risk drug use group (Cramer's V = 0.07). However, effect sizes suggest small to moderate associations. Educational attainment was also associated with substance use, with individuals holding a bachelor's degree or higher being more likely to be in the low-risk categories for both alcohol (Cramer's V = 0.09) and drug use (Cramer's V = 0.09). Employment status had a moderate effect on alcohol use (Cramer's V = 0.13) and a more negligible effect on drug use (Cramer's V = 0.08), as unemployed individuals were overrepresented in high-risk groups for both substances. Census poverty demonstrated a minor impact on drug use (Cramer's V = 0.08), with individuals living in poverty exhibiting higher rates of problematic and hazardous drug use. Self-reported general health status was also linked to substance use, with those reporting excellent or very good health being more likely to be in the low-risk alcohol (Cramer's V = 0.08) and drug (Cramer's V = 0.09) use categories (Table 1).
Demographic Information for Participants.
Notes. &ACE, adverse childhood experience.
*Effect sizes are Cramer's V for chi-squares test, and eta-square for ANOVA.
As shown in Supplemental Tables 2a–2b, post-hoc analyses indicated that hazardous alcohol users and problematic drug users were significantly younger and reported higher ACE scores compared with low- and medium-risk users (p < .01). For categorical variables, pairwise chi-square tests (ie, Holm–Bonferroni adjusted) confirmed that differences in alcohol and drug use risk were most pronounced across education, employment, general health status, and sexual identity categories. In contrast, gender and race/ethnicity differences were minor. Participants with higher education, full-time employment, and very good or excellent health were consistently more likely to fall into the low-risk groups for both alcohol and drug use. In contrast, those reporting economic hardship or poorer health were disproportionately represented among higher-risk users.
Tables 2a and 2b present feature importance rankings and performance metrics from ML models examining alcohol and drug use among SGMs, identifying the most influential features across models. Across all models, social, behavioral, and mental health variables emerged as dominant variables of substance use risk, although the relative influence of specific features varied by model and substance type. For alcohol use, model performance was moderate, with F1 scores ranging from 0.31 to 0.47 and accuracy ranging from 0.37 to 0.44. The XGB model achieved the best balance (F1 = 0.47), while the NN demonstrated consistent performance across both alcohol and drug use tasks (F1 = 0.42 for AUDIT). Lifetime smoking history emerged as the most influential feature in the XGB model, followed by social and identity-related factors such as limited social support, higher frequency of moving, and milestones in sexual identity development (eg, coming-out age to straight friends and age of same-sex relationship). Financial and housing stability indicators, including household income and living arrangements (renting vs owning), were also among the top variables, highlighting protective roles of economic security and residential stability. In addition, recent poor mental or physical health and suicidal ideation contributed substantially to elevated alcohol use risk, underscoring the intersection of psychosocial stress and health vulnerability among SGMs (Table 2a). For drug use, models demonstrated stronger predictive performance, with F1-scores ranging from 0.59 to 0.68 and accuracy values between 0.71 and 0.75, indicating that drug use risk patterns were more distinct and more straightforward to classify than alcohol use patterns. Across models, the most influential features included exposure to conversion therapy, lower perceived social support, and experiences of victimization or discrimination (eg, religious or employment-related). Mental health factors such as poor mental health days, suicidal thoughts, and social well-being ranked prominently, alongside life-course and relational variables, including age of coming out, same-sex relationship age, and family discovery of sexual identity. Health-related factors, such as kidney disease, ACE scores, and STI testing frequency, also contributed to predictive accuracy. The NN model achieved the most stable and generalized performance across both AUDIT and DUDIT outcomes (F1 = 0.42 and 0.68, respectively), reflecting its capacity to capture complex, shared feature representations within the multitask learning framework (Table 2b). Among all classifiers, the NN achieved the most balanced performance across low, moderate, and high-risk categories, with F1-scores that were relatively similar across classes. Supplemental Figures 1-4 displayed SHAP summary plots for the final models, where feature values on the X-axis were standardized (z-scored) rather than raw units (eg, years of age or income in dollars), so that higher values represented relatively higher levels of each characteristic within the sample.
Top 10 Features of the Prediction Model for Alcohol Use Among the LGB Population.
Notes. ACE: adverse childhood experiences.
Top 10 Features of the Prediction Model for Drug Use Among the LGB Population.
Notes. ACE, adverse childhood experiences.
Discussion
This study examined associations between social, psychological, and structural factors and alcohol and drug use among SGM adults using a theoretically grounded, data-driven approach. By applying a combination of ML techniques to a large, nationally representative dataset, we were able to surface both expected and underrecognized factors of substance use across 2 distinct domains: hazardous alcohol use and problematic drug use. While the predictive performance varied by outcomes, our findings reveal meaningful patterns that expand our understanding of substance use risk in SGM groups.
The alcohol use models highlighted factors that resonate with existing literature but also underscore the compounding influence of structural disadvantages and interpersonal stress.35–38 Limited family or social support, particularly in the context of stigma or rejection, further increased the risk of hazardous alcohol use. Notably, housing instability, a widely recognized social determinant of health, was a prominent feature, suggesting that economic precarity and displacement may exacerbate vulnerability to alcohol misuse. The relevance of internalized stigma, identity concealment, and discomfort with specific identity labels suggests that intracommunity dynamics and identity development processes deserve greater attention in alcohol-focused interventions. On the other hand, the alcohol models showed modest discrimination relative to the drug models. This pattern is plausible given that alcohol use is more normative and context-dependent, shaped by social settings, availability/outlet density, and weekday/seasonal routines that are undermeasured in the present dataset.
In contrast, drug use models revealed a more overt pattern of structural and psychological harm. The strong role of conversion therapy exposure aligns with prior findings that such experiences are linked to increased psychological distress, trauma symptoms, and maladaptive coping.39,40 Housing discrimination and criminal legal involvement further exemplify how institutionalized forms of stigma and exclusion contribute to elevated drug use risk. These results are consistent with the Minority Stress Model4,7,8,41 and affirming the SoDH framework.15,16,42 Both frameworks emphasize how material deprivation and systemic bias influence behavioral outcomes.
Our findings carry several implications for intervention and policy. The prioritized social/structural correlates suggest actionable interventions: (1) housing stability supports (eg, rental assistance, eviction protections) in jurisdictions with high housing cost burden; (2) healthcare antidiscrimination enforcement and staff training to reduce stigma in clinical and pharmacy settings; (3) statewide conversion-therapy bans and enforcement; (4) expansion of LGBTQ-affirming primary care and integrated behavioral health; (5) streamlined access to PrEP (same-day starts, confidential or on-site pharmacies, reduced prior authorization); and (6) legal aid and record-expungement programs where criminal-legal exposure emerged as salient.
Our study also highlights the value of ML as a complementary tool for public health equity research.19,20,43,44 Although tree-based methods such as XGB and RF provided strong overall performance, the NN demonstrated the most balanced performance across risk categories. This balanced profile suggests that the NN was less dominated by the majority (low risk) class and retained reasonable discrimination for participants at moderate and high risk of hazardous alcohol use and problematic drug use. This is particularly important in the context of SGM substance use, where the high-risk substance use groups are relatively small but clinically critical, and it suggests that flexible, nonlinear models can improve risk stratification for both lower- and higher-risk SGM adults, rather than primarily optimizing prediction for the majority low-risk group. While the alcohol models showed only modest discrimination, the convergence of salient features across algorithms supports the robustness of the observed patterns. Using SHAP value, we retained transparency about how stressors and supports relate to substance-use classifications, clarifying the direction and relative contribution of multilevel correlates. Because the SHAP plots are based on standardized feature values, the patterns should be interpreted as relative increases or decreases in each predictor (eg, higher vs lower age or income) rather than specific raw units. These models are not intended for diagnosis or deployment; rather, they are hypothesis-generating and designed to inform intervention prioritization and policy strategy. Beyond reaffirming known correlations, interpretable ML highlighted nonlinearities, interactions, and threshold-like effects that conventional linear models may miss. For example, the cooccurrence of stigma and depressive symptoms, and unstable housing with poverty, was linked to disproportionate increases in risk in SHAP interaction plots, even when each factor showed only small bivariate associations. Together, these insights identify that cooccurring conditions that are most strongly linked to hazardous alcohol use or problematic drug use and point to concrete targets for coordinated screening and structural action.
Limitations
Several limitations warrant consideration. First, although we applied for the Generations Study survey weights, residual sampling, and nonresponse bias may persist and could affect generalizability. Second, key variables, including alcohol/drug use and stigmatization, were self-reported, and thus subject to social desirability and recall bias. Third, the cross-sectional design precludes causal inference; results should be interpreted as associations rather than effects. Fourth, bivariate effect sizes were generally small, consistent with complex, multilevel phenomena in which individual variables exhibit modest associations, but combinations of factors carry a stronger signal. Fifth, alcohol consumption may be especially influenced by unmeasured social and environmental norms, which could contribute to the comparatively lower discrimination of alcohol models in this study. Sixth, limited cell sizes for some subgroups (eg, nonbinary participants, Black transgender participants) constrained stratified analyses; future work should oversample underrepresented SGM subgroups to enable precision estimates.
Conclusion
A convergence of identity-based stressors, interpersonal exclusion, and structural inequities shapes high-risk substance use among SGMs.4–6 Our study adds to the growing body of work that calls for comprehensive, trauma-informed, and equity-centered approaches to substance use prevention and care in SGM populations. By combining theory, empirical data, and modern analytic tools, we offer a framework for identifying those most at risk—and the structural levers needed for changes.
Supplemental Material
sj-docx-1-jia-10.1177_23259582251409855 - Supplemental material for Exploring Social and Structural Determinants of Substance Use Among Sexual and Gender Minority Adults in the United States
Supplemental material, sj-docx-1-jia-10.1177_23259582251409855 for Exploring Social and Structural Determinants of Substance Use Among Sexual and Gender Minority Adults in the United States by Chen Zhang, Yao Tang, Wonkyung Kniffen and Yu Liu in Journal of the International Association of Providers of AIDS Care (JIAPAC)
Supplemental Material
sj-docx-2-jia-10.1177_23259582251409855 - Supplemental material for Exploring Social and Structural Determinants of Substance Use Among Sexual and Gender Minority Adults in the United States
Supplemental material, sj-docx-2-jia-10.1177_23259582251409855 for Exploring Social and Structural Determinants of Substance Use Among Sexual and Gender Minority Adults in the United States by Chen Zhang, Yao Tang, Wonkyung Kniffen and Yu Liu in Journal of the International Association of Providers of AIDS Care (JIAPAC)
Supplemental Material
sj-docx-3-jia-10.1177_23259582251409855 - Supplemental material for Exploring Social and Structural Determinants of Substance Use Among Sexual and Gender Minority Adults in the United States
Supplemental material, sj-docx-3-jia-10.1177_23259582251409855 for Exploring Social and Structural Determinants of Substance Use Among Sexual and Gender Minority Adults in the United States by Chen Zhang, Yao Tang, Wonkyung Kniffen and Yu Liu in Journal of the International Association of Providers of AIDS Care (JIAPAC)
Supplemental Material
sj-pdf-4-jia-10.1177_23259582251409855 - Supplemental material for Exploring Social and Structural Determinants of Substance Use Among Sexual and Gender Minority Adults in the United States
Supplemental material, sj-pdf-4-jia-10.1177_23259582251409855 for Exploring Social and Structural Determinants of Substance Use Among Sexual and Gender Minority Adults in the United States by Chen Zhang, Yao Tang, Wonkyung Kniffen and Yu Liu in Journal of the International Association of Providers of AIDS Care (JIAPAC)
Footnotes
Acknowledgments
This publication is funded by the University of Rochester, the Rita and Alex Hillman Foundation, and the generosity of Leonard A. Lauder. The Rita and Alex Hillman Foundation (the “Foundation”) is one of the nation's leading philanthropies dedicated to advancing health equity and improving the health and healthcare of marginalized populations. The Foundation supports innovative, nursing-driven models of care that address the needs of communities experiencing inequity, discrimination, oppression, and indifference. The contents of this manuscript do not necessarily reflect the opinions of the funders.
Ethical Approval and Written Informed Consent Statements
The University of Rochester Research Subjects Review Board reviewed this project (Study ID: STUDY00011070) and determined that it does not constitute human subjects research because it involves secondary analysis of a fully deidentified dataset (ICPSR 37166). As such, additional IRB approval and informed consent were not required for this analysis.
Authorship Contribution
CZ conceived and designed the study, supervised the analytic plan, and drafted the manuscript; YT conducted data management, coding, and ML analyses, and contributed to the interpretation of the findings; WK contributed to study design, interpretation of results in relation to mental health and clinical implications, and critically revised the manuscript for important intellectual content; YL contributed to conceptualization, provided methodological and statistical guidance, and substantially revised the manuscript for clarity and interpretation. All authors reviewed and approved the final version of the manuscript and agree to be accountable for all aspects of the work.
Funding
This study was sponsored by the Research Award of the University of Rochester and the Rita and Alex Hillman Foundation.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Data Availability Statement
All data were downloaded from the data repository from https://www.icpsr.umich.edu/web/DSDR/studies/37166. The current analyses used the publicly available dataset. The original protocol for the study is available: ![]()
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
