Abstract
Background:
Homelessness is an important social determinant of health (SDOH), impacting health outcomes for many medical conditions. Although homelessness is common among people with opioid use disorder (OUD), few studies systematically evaluate homelessness and other SDOH among people enrolled in standard of care treatment for OUD, medication for opioid use disorder (MOUD), or examine whether homelessness affects treatment engagement.
Methods:
Using data from the 2016 to 2018 U.S. Treatment Episode Dataset Discharges (TEDS-D), patient demographic, social, and clinical characteristics were compared between episodes of outpatient MOUD where homelessness was reported at treatment enrollment versus independent housing using pairwise tests adjusted for multiple testing. A logistic regression model examined the relationship between homelessness and treatment length and treatment completion while accounting for covariates.
Results:
There were 188 238 eligible treatment episodes. Homelessness was reported in 17 158 episodes (8.7%). In pairwise analysis, episodes involving homelessness were significantly different from those involving independent living on most demographic, social, and clinical characteristics, with significantly greater social vulnerability in most SDOH variables (P’s < .05). Homelessness was significantly and negatively associated with treatment completion (coefficient = −0.0853, P < 0.001, 95% CI = [−0.114, −0.056], OR = 0.918) and remaining in treatment for greater than 180 days (coefficient = −0.3435, P < 0.001, 95% CI = [−0.371, −0.316], OR = 0.709) after accounting for covariates.
Conclusions:
Patients reporting homelessness at treatment entry in outpatient MOUD in the U.S. represent a clinically distinct and socially vulnerable population from those not reporting homelessness. Homelessness independently predicts poorer engagement in MOUD confirming that homelessness may be an independent predictor for MOUD treatment discontinuation nationally.
Keywords
Introduction
The overdose epidemic is rapidly worsening in the United States. There were more than 105 000 lives lost to overdose in 2021, with a similar number predicted for 2022. 1 Multiple factors contribute to this crisis, including a drug supply that contains increasing quantities of potent synthetic opioids,2,3 but it is in part driven by untreated opioid use disorder (OUD). 4 People experiencing homelessness (PEH) are particularly impacted by the overdose crisis. Overdose is the leading cause of death among PEH, comprising up to one fourth of deaths in this group.5,6 Compared to housed individuals, those experiencing homelessness have a mortality rate 7 to 9 times higher.6,7 As such, addressing OUD and overdose in this high-risk group should be a priority for treatment organizations, public health entities, and governmental agencies.
The relationship between homelessness and substance use disorders (SUDs) is bidirectional. Homelessness is common among people with SUDs and may exacerbate substance use; likewise, SUDs, including OUD, are common among PEH 8 and may contribute to housing instability. 9 Medication for opioid use disorder (MOUD) with methadone, buprenorphine, or naltrexone, is evidence-based treatment for OUD and represents the standard of care. 10 Receipt of MOUD confers many patient benefits, and is the only OUD treatment shown to reduce substance use, infectious disease risk, fatal overdoses, and all-cause mortality.10 -13 MOUD enrollment is also associated with housing-specific improvements, such as increased retention in supportive housing. 14 While many people with OUD receive non-MOUD treatments, 15 PEH with OUD who enter SUD treatment are less likely than their housed counterparts to receive MOUD.16 -18 Expanding engagement in evidence-based OUD treatment is even more important in the face of the aforementioned worsening overdose crisis, which disproportionately impacts PEH.
Despite its prevalence among people with OUD, research on homelessness and other social factors that impact health among this population is scant. Housing is a social determinant of health, evidenced by people experiencing homelessness having higher rates of chronic illness, delayed diagnoses of medical conditions, reduced access to treatment and preventive health, and increased mortality compared to those who are housed.19,20 Social Determinants of Health (SDOH) are individual and systemic factors, such as income, insurance status, or housing, that impact an individual’s physical wellbeing, access to healthcare, and health outcomes.21,22 Homelessness is an important social determinant that impacts health, but it is intertwined with other factors that may also limit healthcare access, like poverty or insurance status. Investigating SDOH among individuals seeking MOUD may reveal potential areas for intervention to increase engagement in MOUD and improve health outcomes and quality of life among people with OUD. For example, one prior study demonstrated that out-of-pocket costs are associated with treatment discontinuation in buprenorphine. 23 However, literature evaluating patient social characteristics in OUD treatment, particularly MOUD, is generally scarce. To our knowledge, there are no published studied examining SDOH based on housing status in MOUD.
MOUD retention rates range widely 24 with most individuals staying in treatment fewer than 6 months.25,26 Few studies examine MOUD retention among PEH. Studies investigating predictors of retention in methadone,27 -29 buprenorphine, 30 and MOUD generally, 25 demonstrate associations between housing and treatment engagement, but few prior studies directly compare MOUD outcomes based on housing status.18,31 -35 A recent study found that homelessness is a significant independent predictor of reduced 12-month retention in low-barrier outpatient methadone treatment, even after adjusting for sociodemographic and clinical characteristics. 31 Studies examining potential associations between housing and engagement in buprenorphine have mixed results.30,33 -35 One prior study utilized the 2018 Treatment Episode Dataset Discharges (TEDS-D), a national database of SUD treatment episodes across the United States in a given year, to examine factors associated with MOUD retention length greater than 6 months and found homelessness was one of several associated factors. 25 A recent study utilizing the 2015 to 2017 TEDS-D database found that, among all outpatient OUD treatment episodes, those experiencing homelessness were less likely to enter MOUD and had reduced treatment completion and 90-day retention in all treatment types (ie, non-pharmacologic and MOUD). 18 However, no published studies have examined whether there are differences in longer-term retention once patients experiencing homelessness enter MOUD.
Prior work suggests that demographic characteristics (eg, age), clinical characteristics (eg, substance use), and treatment-related characteristics (eg, methadone dose) affect retention in MOUD.23,24,36,37 However, there are surprisingly few studies examining whether such characteristics differ in OUD treatment between PEH and their housed counterparts. Two prior studies compared patient demographic and clinical characteristics based on housing status using the TEDS database17,18; Han and colleagues studied all OUD treatment modalities whereas Friesen and Young examined episodes from all outpatient OUD treatment types. They both found significant differences in demographic (eg, sex) and clinical variables (eg, co-occurring psychiatric and substance use disorders) based on housing status, but whether these differences persist for individuals receiving standard of care OUD treatment (MOUD) is unknown.
We used a national database to investigate patient characteristics, social determinants of health, and retention in outpatient MOUD treatment based on housing status. Using data on outpatient MOUD treatment episodes from the 2016 to 2018 Treatment Episode Dataset Discharges (TEDS-D), we compared patient demographics, SDOH, and clinical characteristics based on housing status and examined 2 measures of engagement: 180-day treatment retention of time in treatment and treatment completion. We hypothesized that housing status would be significantly and negatively associated with treatment engagement and that this association would remain after accounting for other covariates.
Materials and Methods
Data Source
Data were retrieved from the 2016 to 2018 Substance Abuse and Mental Health Services Administration (SAMHSA) Treatment Episode Data Set—Discharges (TEDS-D), a national, publicly available database containing information about all SUD treatment episodes that were either publicly funded or conducted at facilities receiving public funding in the U.S. SAMHSA requires all states to report this information for any treatment episodes with providers receiving public funding. The TEDS dataset represents the majority of all SUD treatment episodes conducted in the U.S. in a given year and the TEDS-Discharges (TEDS-D) contains information on both admissions and discharges for these treatment episodes. These data reflect treatment episodes from all states with the exception of Georgia (2016-2018), Oregon (2016-2018), Washington (2018 only), and West Virginia (2016-2018), which were excluded due to a lack of sufficient data. The data reflect treatment episodes rather than individual patients; thus, an individual patient may be included more than once in the dataset if they had multiple treatment episodes. Study investigators contacted the Human Investigations Committee at the Yale School of Medicine to determine whether submission was appropriate for this study, who determined a study utilizing this publicly available dataset did not require review.
Population
We limited our study population to treatment episodes for outpatient MOUD treatment. We restricted the analysis to treatment episodes coded as ambulatory, non-intensive outpatient, then further limited to those where patients reported homelessness or independent living for their living arrangements at admission. We excluded those with dependent living arrangements due to the potential heterogeneity of these settings. We then further restricted the analysis to episodes whose completion reason was listed as: “completed,” “dropped out,” or “terminated by facility.” We excluded episodes whose reason for treatment discontinuation were coded as “transferred to another treatment program or facility,” “incarcerated,” “death,” or “other.” Patient transfers were excluded as the definition or requirements for a transfer may differ across states and individual treatment facilities.
Variables
Treatment Completion
The TEDS database provides data on episode reason for discharge. We considered a patient to have completed treatment if their reason for treatment completion was coded as “Completed.” Participants were considered to have been prematurely terminated from treatment for those whose reasons for discharge were “dropped out,” or “terminated by facility,” representing patient-initiated and program-initiated discharges, respectively.
Treatment Length
Treatment episode length was reported as an interval censored length variable. In the dataset, treatment episodes lasting up to 30 days have the actual number of days in treatment coded, whereas longer lengths of stay are coded as 31 to 45, 46 to 60, 61 to 90, 91 to 120, 121 to 180, 181 to 365, and greater than 365 days. We dichotomized this variable such that episodes were categorized as lasting longer than 180 days (a composite variable composed of the 181 to 365 and greater than 365 days groups) or 180 days or less (a composite variable composed of the remaining 5 groups). This time point was selected to align with other studies of MOUD treatment retention that utilize between 180 and 365 days as a marker of longer-term retention24,25 including one prior study which illustrated significantly decreased retention at 12-months between individuals enrolled in methadone treatment who reported experiencing homelessness at treatment entry compared to those who were living independently. 31
Demographic, Social, and Clinical Characteristics
A full codebook for the TEDS-D dataset is available at https://www.datafiles.samhsa.gov/dataset/teds-d-2018-ds0001-teds-d-2018-ds0001. Race was compressed to American Indian or Alaska Native/Asian, Pacific Islander, or Native Hawaiian/Black/other/2 or more races/white. Ethnicity was compressed to Hispanic/Not Hispanic. Marital status was compressed to Now married/Never married/Previously married. Education was compressed to high school or greater/less than high school.
Statistical Analysis
Housing Status Versus Covariates Pairwise Analysis
To understand associations between housing status and the other covariates, difference in distribution tests were run for each covariate comparing the housed to the unhoused group. Categorical variables were evaluated with a chi-squared test. The ordinal “age” variable was evaluated with a Cochran-Armitage test.38,39 After computing raw P-values for each test, we controlled the False Discovery Rate with the Benjamini-Hochberg (BH) procedure to obtain adjusted P-values. 40
Housing Status Versus Dropout Reason
For the subset of 156 840 treatment episodes where the non-completion reason was listed as “Dropped out” or “Terminated by facility” we tested for an association between non-completion reason and housing status using a chi-squared test.
Retention Regression, Confirmatory Analysis: Housing Status Association With Retention
We tested the hypothesis that the association between housing status and retention would remain after accounting for other covariates by fitting a logistic regression model regressing the retention outcome on 61 covariates. For this analysis we report the p-value for housing status. The logistic regression model was fit using the statsmodels in python. 41 We completed this regression twice to examine 2 outcomes for retention: length of time and reason for treatment termination. For length of time, we compared whether an episode of treatment lasted at greater than 180 days or not. For reason for treatment termination, we compared completed treatment to premature treatment discontinuation, a combination variable containing episodes coded as “dropped out” or “terminated by facility.”
Software and Reproducibility
The statistical analysis was performed using Python and R. A software repository reproducing the entire analysis can be found at https://github.com/idc9/repro_teds_d_retention. In addition to the packages mentioned above, the analysis made use of matplotlib, 42 numpy, 43 pandas, 44 and sklearn. 45
Results
The overall sample contained 188 238 treatment episodes, within which the patients reported independent living in 171 115 episodes and experiencing homelessness in 17 158. Results of the pairwise analysis comparing characteristics of patients for each treatment episode based on housing status are shown in the tables: Table 1 shows demographic characteristics; Table 2 shows social characteristics; and Table 3 shows clinical and treatment-related characteristics.
Demographic Characteristics Based on Housing Status Among 188 238 Treatment Episodes for Outpatient Medication for Opioid Use Disorder Treatment in the United States From 2016 to 2018.
Social Determinants of Health Based on Housing Status Among 188 238 Treatment Episodes for Outpatient Medication for Opioid Use Disorder Treatment in the United States From 2016 to 2018.
Clinical and Treatment-Related Characteristics Based on Housing Status Among 188 238 Treatment Episodes for Medication for Opioid Use Disorder Treatment in the United States From 2016 to 2018.
Abbreviations: EAP, employee assistance program; DUI, driving under the influence; DWI, driving while impaired; HMO, Health Maintenance Organization; AA, alcoholics anonymous; NA, narcotics anonymous.
Baseline Characteristics and Housing Status
All variables except the following were associated with housing status after applying the BH procedure: attendance at substance use self-help groups (AA [“Alcoholics Anonymous”], NA [“Narcotics Anonymous”], other mutual aid) in 30 days prior to admission and discharge. The group experiencing homelessness was more likely to be male, Hispanic, not currently married, without a high school diploma, without employment, not a veteran, and older, and have a recent arrest in the 30 days prior to admission and discharge, have had 1 or more prior treatment episodes, have more frequent substance use at admission and discharge, report first substance use at a younger age, report current intravenous drug use, report co-occurring mental health and substance use disorders, and have no health insurance. There were also significant differences between the primary substance used on admission to treatment, primary payment source for treatment, referral source, region, and racial background between the group of episodes where the patient reported homelessness at treatment admission versus those where the patient was housed.
Housing Status and Reason for Treatment Non-Completion
We examined the relationship between housing status and reason for treatment termination among those treatment episodes that were not coded as “completed.” Housing status was significantly associated with dropout reason (P < .01, df = 1). Compared to those who were housed, individuals experiencing homelessness were significantly less likely to be terminated by the program (15.6% unhoused vs 19.8% housed) and were more likely to experience patient-initiated treatment termination.
Regression Confirmatory Analyses
Housing status was statistically significantly related to retention for both logistic regression analyses. When premature treatment termination was the outcome, the housing status coefficient was estimated to be −0.4283 (P < 0.001, 95% CI = [−0.478, −0.378], odds ratio 0.6516) where the negative sign indicates homelessness has a negative effect on retention. For retention at greater than 180 days the housing status coefficient was estimated to be −0.4 (P < 0.001, 95% CI = [−0.434, −0.367], odds ratio = 0.67).
Discussion
To our knowledge, this study is one of the first investigations of differences in U.S. national MOUD treatment outcomes and social determinants of health based specifically on housing status. There were several main findings. First, the patients reporting homelessness at treatment entry into episodes of outpatient MOUD represent a clinically distinct and more socially vulnerable population from those with independent housing. Second, homelessness was significantly and negatively associated with staying in outpatient MOUD treatment for 180 days or longer, a relationship that remained significant after accounting for covariates. Third, homelessness was significantly negatively associated with treatment completion after adjusting for covariates. Finally, treatment episodes where homelessness was reported at treatment entry were more likely to end in a patient-initiated treatment termination and less likely to experience a program-initiated termination compared to those who were housed.
Our findings reveal a number of demographic differences and social and clinical vulnerabilities among people experiencing homelessness that are distinct from those housed at entry into outpatient MOUD treatment. On nearly every social, demographic, and clinical variable, there were significant differences between the group of episodes representing patients experiencing homelessness and those representing patients who were housed, aligning with national data on homelessness. For example, we found the group experiencing homelessness had significantly more males than the housed group. In comparison to their proportion of the general population, males are overrepresented both among the group of people experiencing homelessness in the U.S. and those with substance use disorders. 46 Similarly, we found a statistically significant association between race and housing status. Looking at Table 1, we show that people experiencing homelessness were empirically more likely to be non-white than the general population, which matches national data on homelessness. 47 Structural factors such as systemic racism may have led to a disproportionately high rate of homelessness among minoritized people in the United States generally and among individuals pursuing outpatient MOUD. 48
Our findings highlight the potentially additive or intersectional nature of social determinants of health for people experiencing homelessness and how they might potentially compound to worsen outcomes in MOUD treatment. Examples from the present study include how, compared to episodes among those who were housed, we found the group of episodes among PEH had significantly higher rates of unemployment, no income, receipt of public assistance, lacking health insurance, and arrests in the 30 days prior to admission and discharge. While the impact of SDOH on morbidity and mortality for many medical conditions49 -51 has been well documented, the literature on the impact of SDOH on MOUD outcomes or interventions is limited.52 -54 One existing study examining SDOH including income inequality and structural racism demonstrated an association with increased rates of injection drug use and worsened related outcomes such as HIV infection and AIDS-related mortality. 55 Similarly, social determinants, including housing status and geography, have been shown to be associated with overdose outcomes. 56 Investigation into interventions that directly address homelessness and other related SDOH is warranted, as are explorations of SDOHs at treatment entry and subsequent treatment planning. Unique to methadone treatment and other outpatient interventions are regulatory requirements of treatment planning at prescribed intervals (eg, 30, 60, 90, and 180 days). This mandatory contact with patients in methadone treatment is unique and could be considered an opportunity to explore and address SDOH.
To our knowledge, this is the first study to directly compare 180-day treatment retention in outpatient MOUD based on patient housing status. One prior study examining factors related to treatment retention in one year of TEDS-D data identified homelessness as one of several risk factors for reduced length of treatment in outpatient MOUD. 25 Our findings build off a recent study that found reduced retention and treatment completion at 90 days in outpatient OUD including MOUD by analyzing a more recent year of TEDS-D data and investigating longer-term treatment engagement, in addition to the expansion of social determinants and patient characteristics discussed above. 18 Our results agree with previous work showing homelessness was an independent risk factor for 12-month treatment discontinuation in low-barrier methadone treatment, where all treatment terminations were patient-initiated 31 to include different types of MOUD programs, such as those dispensing buprenorphine and naltrexone, and include program-initiated discharges. Prior work has indicated homelessness may be a risk factor for decreased retention in MOUD when including housing status as one potential covariate for treatment retention,24,27,29,57 but these studies were in smaller samples and did not directly examine retention outcomes, such as length of time in treatment or reason for treatment termination, based on housing status. In the current study, homelessness was significantly and negatively associated with treatment completion and time in treatment in pairwise analysis and in logistic regression after adjusting for covariates, confirming that homelessness may be an independent predictor for MOUD treatment discontinuation nationally. The fact that episodes involving PEH were more likely than those among people who were housed to experience a patient-initiated, rather than program-initiated, premature termination is notable and requires future investigation into reasons PEH may prematurely terminate treatment. Future work into interventions that improve treatment retention among individuals experiencing homelessness is warranted.
Prior studies have examined housing status within the context of the TEDS dataset. One study examining racial and ethnic disparities in treatment effectiveness of outpatient non-MOUD substance use disorder treatment found that homelessness was associated with decreased improvement during treatment using the 2015 to 2017 dataset. 58 Another study investigating socioeconomic determinants of planned methadone treatment noted that those who were not experiencing homelessness were more likely to receive methadone as part of their treatment 59 and a later study found that experiencing homelessness was associated with decreased odds of receiving MOUD and also lower odds of waiting greater than 1 week to enter treatment. 16 Similarly a study of the use of methadone within “detoxification” settings found those experiencing homelessness were less likely to have methadone planned for their treatment. 15 The present study’s results build off these studies by demonstrating that, in addition to the disparities that may limit treatment entry into MOUD, those individuals experiencing homelessness that are able to enroll into MOUD have reduced retention length and decreased odds of treatment completion. One important area of future direction is to examine how state and city-specific policies impact disparities based on housing status observed in the present study. Future policies to address addiction among PEH will need to intervene on SDOH such as homelessness to have sustained outcomes. 60
Limitations
Our study findings should be considered in light of their limitations. Our definition of housing status relied on treatment centers’ determinations; thus, we were unable to provide a specific definition of homelessness beyond what is reported in the dataset. Our study design utilized housing status at treatment enrollment as the independent variable, limiting examination of longitudinal housing trajectories and their influence on MOUD outcomes. In addition, the TEDS database groups all 3 FDA-approved medications for OUD together, limiting evaluation into outcomes based specifically on treatment type (eg, comparison of methadone to buprenorphine to naltrexone). Finally, while we found many statistically significant differences between the 2 groups, it is possible that not all of these findings are clinically significant given the overall small difference in percentages between the 2 groups (eg, veteran experience, referral source, and prior substance use treatment episodes). Future work should examine how these and our other findings may impact patient experience and outcomes in MOUD.
Conclusion
In a national sample of treatment episodes for outpatient MOUD treatment, the episodes in which patients reported homelessness at treatment entry differed significantly in demographic, social, clinical, and treatment-related characteristics, representing a clinically distinct and more socially vulnerable population. Homelessness was significantly and negatively related to treatment retention in MOUD based on 2 separate metrics, length of time in treatment and treatment completion. Strategies to enhance retention in evidence-based treatment for OUD are needed in the midst of a worsening overdose crisis and may need to include systematic strategies for assessing and improving housing status among people experiencing homelessness enrolled in MOUD treatment.
Footnotes
Acknowledgements
Results from this study were presented as an oral presentation at the Association for Multidisciplinary Education and Research in Substance Use and Addiction (AMERSA) Conference in Boston, Massachusetts on November 10, 2022
Author Contributions
MGG and DTB conceived of the study aims. MGG, DTB, and IDC conceptualized study methodology. IDC conducted the formal analysis, curated the data, and validated the analysis. MGG produced data visualization and oversaw project administration. All authors wrote the manuscript. LMM provided resources for project completion. DTB supervised the project.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This publication was made possible by the James G. Hirsch Endowed Medical Student Research Fellowship and the Yale School of Medicine Medical Student Fellowship (MGG), the National Science Foundation under Award No. 1902440 (IDC), and the National Institutes of Health U01 HL150596-01 and RM1 DA055310 (DTB).
Compliance,Ethical Standards,and Ethical Approval
Study investigators contacted the Human Investigations Committee at the Yale School of Medicine to determine whether submission was appropriate for this study, who determined a study utilizing this publicly available dataset did not require review.
