Associations among Drug Acquisition and Use Behaviors, Psychosocial Attributes, and Opioid-Involved Overdoses: A SEM Analysis

Aims This study sought to develop and assess an exploratory model of how demographic and psychosocial attributes, and drug use or acquisition behaviors interact to affect opioid-involved overdoses. Methods We conducted exploratory and confirmatory factor analysis (EFA/CFA) to identify a factor structure for ten drug acquisition and use behaviors. We then evaluated alternative structural equation models incorporating the identified factors, adding demographic and psychosocial attributes as predictors of past-year opioid overdose. We used interview data collected for two studies recruiting opioid-misusing participants receiving services from a community-based syringe service program. The first investigated current attitudes toward drug-checking (N = 150). The second was an RCT assessing a telehealth versus in-person medical appointment for opioid use disorder treatment referral (N = 270). Demographics included gender, age, race/ethnicity, education, and socioeconomic status. Psychosocial measures were homelessness, psychological distress, and trauma. Self-reported drug-related risk behaviors included using alone, having a new supplier, using opioids with benzodiazepines/alcohol, and preferring fentanyl. Past-year opioid-involved overdoses were dichotomized into experiencing none or any. Results The EFA/CFA revealed a two-factor structure with one factor reflecting drug acquisition and the second drug use behaviors. The selected model (CFI = .984, TLI = .981, RMSEA = .024) accounted for 13.1% of overdose probability variance. A latent variable representing psychosocial attributes was indirectly associated with an increase in past-year overdose probability (β=.234, p = .001), as mediated by the EFA/CFA identified latent variables: drug acquisition (β=.683, p < .001) and drug use (β=.567, p = .001). Drug use behaviors (β=.287, p = .04) but not drug acquisition (β=.105, p = .461) also had a significant, positive direct effect on past-year overdose. No demographic attributes were significant direct or indirect overdose predictors. Conclusions Psychosocial attributes, particularly homelessness, increase the probability of an overdose through associations with risky drug acquisition and drug-using behaviors. To increase effectiveness, prevention efforts might address the interacting overdose risks that span multiple functional domains.


Background
Opioid-involved overdoses and related fatalities remain persistently high in the United States and Canada despite considerable efforts to increase harm reduction interventions and treatment access (1,2).Per the U.S. National Center for Health Statistics, 85,145 opioid-overdose fatalities are expected for the twelve months ending in June 2023 compared with 82,325 in June 2022, a 3.4% increase (3).Among the more prominent harm reduction efforts undertaken to reduce opioid-involved overdoses and associated fatalities are overdose education and naloxone distribution (OEND); expansion of drug checking services to test illegal drugs through point-of-care testing or by providing users test strips to detect the presence of fentanyl, xylazine, and benzodiazepines; reduction in prescription opioid availability; and the expansion of opioid agonist-and antagonist-based treatments (e.g., buprenorphine, methadone, naltrexone) for opioid use disorder (4)(5)(6)(7)(8)(9).
The continuing year-over-year increases in opioid-related overdoses and fatalities have been attributed largely to the presence of fentanyl and other potent synthetic opioids in the illegal drug supply (2,6,(10)(11)(12)(13).Research has also identi ed a number of risk factors associated with higher rates of opioid-involved overdose rates.For instance, a recent study identi ed differences in overdose rates by biological gender, with males two to three times more likely to experience synthetic opioid-related overdoses and fatalities than females across the lifespan, even after accounting for gender-based differences in misuse and other demographics (14).Speculating on possible reasons for the differential gender-based overdose mortality rates, the authors suggest males might have a greater propensity towards risky drug use behaviors such as "injecting alone, taking larger doses, or using untrusted suppliers".
For a given individual and on any particular drug-use occasion, some or all these factors likely interact in complex ways to affect overdose probability.The nature of how these factors interact to increase overdose risk, however, has largely been unexplored.Studies of multiple risk factors have typically examined them in the context of multiple regression models which, though illuminating, do not analytically consider how independent risk factors might cluster or in uence each other to affect overdose probability (23,24).Thus, there remains a question as to why one opioid user experiences one or more unintentional overdoses over time, whereas another user experiences no or fewer overdoses.This is especially puzzling given the near ubiquity of fentanyl in the illegally manufactured opioid supply and the identi cation of fentanyl as a main contributing factor to drug overdoses unless, possibly, uctuations in fentanyl potency play a larger role than currently understood (2,(29)(30)(31).

The present study
The present study sought to develop and assess an exploratory model of how different behavioral, psychosocial, and contextual/circumstantial risk factors potentially combine to increase opioid-related overdoses.Speci cally, we rst assessed whether items comprising each of these domains clustered into superordinate latent factors using exploratory factor analysis (EFA).Using a second sample of participants, we next assessed this initial factor structure through con rmatory factor analysis (CFA).In the nal step, we used structural equation modeling (SEM) to test whether and how the identi ed latent factors were associated with each other and/or with experiencing an overdose in the past year.
We hypothesized that drug use risk behaviors would directly in uence the probability of an overdose and would themselves be in uenced by psychosocial/demographic factors such as homelessness and gender.We also hypothesized that psychosocial/demographic factors would directly as well as indirectly affect the probability of an overdose, mediated by drug use risk behaviors.

Setting
We used survey interview data collected for two research projects, both of which independently recruited participants from two community-based syringe service program (SSP) sites.The sites are located on Chicago's west and northwest sides in socioeconomically disadvantaged neighborhoods composed largely but not exclusively of racial and ethnic minorities.Among Chicago communities, both neighborhoods have relatively high opioid-involved overdose and fatality rates (32).The SSP sites predominantly serve clients who use opioids as well as other illegally manufactured and sold drugs (e.g., cocaine, benzodiazepines, MDMA).Services provided include but are not limited to syringe exchange, naloxone and fentanyl test strip distribution, drug checking, condom distribution, and medical consultation.
Additional details on both study protocols are available elsewhere (33)(34)(35).The rst study had two components: assessing street-drug users' current attitudes towards drug-checking services and testing the feasibility of using ecological momentary assessment (EMA) to collect data.The second study was a randomized controlled trial (RCT) testing the relative effectiveness of an initial telehealth versus in-person medical appointment on initiation and engagement in medication for opioid use disorder (MOUD) treatment.Participants in both studies provided written informed consent and were compensated with gift cards for their time and participation.

Participants
EFA sample.We used the combined samples from the study on current attitudes towards drug-checking (N = 130) and the feasibility of using EMA to collect drug-checking data from a non-overlapping group of participants (N = 20) to develop a preliminary model of the factor structure of speci c behaviors deemed to increase overdose risk.Both study components used convenience sampling to recruit opioid-using participants at the SSP sites between August 2021 and June 2022.Current opioid use was determined by requiring that all participants provide a small test sample of the drugs they believed to be heroin or another opioid and which they had purchased from an illegal source on the day of recruitment.The provided samples were subsequently tested off-site to assess the actual drug mixtures (36).All participants were also required to provide information on behaviors increasing overdose risk.However, demographic information was not available for twenty participants who opted out of taking a longer survey that included the demographic questions and which was administered to the group participating in the attitudes survey component (34).
CFA sample.The sample used for con rmatory CFA and SEM modeling (N = 274) was obtained from the RCT study.Convenience sampling was also used for this study, with participant recruitment running from August 2020 through July 2022.We recruited participants from among those seeking services at the SSP sites as well as by conducting street outreach in the surrounding communities to recruit current opioid users interested in beginning MOUD treatment.Participants were screened for an OUD at any severity level using the DSM-5 checklist of diagnostic criteria for an opioid use disorder (37).All of the 274 participants met DSM-5 OUD diagnostic criteria, with 269 (98.1%) meeting criteria for a severe OUD.
Owing to missing data on covariates used for the SEM structural analysis step, we excluded four participants (1.4%), yielding a nal analytic sample of 270.

Measures
Drug use behaviors.To the best of our knowledge, there is no existing validated inventory of speci c drug use behaviors that increase overdose risk.We used data collected from the list of drug use behavior questions asked of participants in both studies, developed through a review of the research on drug-use overdose risk behaviors and in consultation with SSP program staff.The ten speci c items were: used in a new location; used with new people; used alone; had a new source/supplier; used a different (more or less) dosage than usual; used multiple drugs; used opioids with benzodiazepine; used opioids with alcohol; used for the rst time in a while; and used a different administration route (e.g., injecting instead of snorting).Participants were presented with this item list and asked to select all that applied to their drug use in the past month.We dichotomized the responses into ten variables by scoring selected risk behaviors as one (1) and those unselected as zero (0).Psychosocial factors.Using the CFA/SEM RCT study sample, we assessed participants for the following psychosocial factors potentially associated with an increased risk for an opioid-related overdose, all of which were coded dichotomously as present (1) or absent (0): severe psychological distress (SPD); homelessness; trauma exposure; and criminal justice involvement.
We assessed for SPD using Kessler's K6-scale, a brief, widely used and validated 6-item scale consisting of questions related to symptoms of anxiety and depression (e.g., How often in the past 30 days did you feel nervous?; How often did you feel hopeless?, etc.).Scores range from 0-24 with scores above 12 indicating SPD and at least moderate functional impairment (38, 39).
Homelessness was assessed with a single question: Thinking about the past year, what kind of place would you say you lived in most of the time throughout the year?Participants who responded they were homeless and lived on the streets most of the time in the past year were counted as experiencing homelessness (1), whereas all others were scored as not homeless (0).
We assessed trauma exposure using two questions: 1) Have you ever experienced violence or trauma in any setting (including community or school violence, domestic violence, etc.) and 2) In the past 30 days, how often have you been hit, kicked, slapped, or otherwise physically hurt?Participants answering yes to both or either question were indicated as having been trauma exposed (1).Those answering no to both questions were coded as non-exposed to trauma (0).
To assess past-year criminal justice involvement, we used responses to the following single question, with those answering yes counted as having criminal justice involvement (1) and those responding no as not criminal justice involved (0): Not counting minor tra c violations, were you arrested, booked, or charged for breaking a law in the past 12 months?Fentanyl preference and drug injection use.As mentioned, a considerable body of research indicates opioid-involved overdoses and fatalities occur owing to fentanyl and other synthetic opioids such as carfentanil mixed with or supplanting heroin and other drugs such as cocaine in the illegally manufactured drug supply (40)(41)(42)(43).Because such a large proportion of our sample had been exposed to fentanyl -baseline urinalysis testing of 240 study participants revealed a majority (N = 225, 93.7%) tested positive for recent fentanyl use -there was not enough variation in fentanyl use to include it as a variable.Instead, to differentiate fentanyl use among participants, we asked how many days in the past month they intentionally used fentanyl, reasoning that those intentionally using fentanyl would be more likely to seek out suppliers thought to provide fentanyl and/or to obtain and use illegally manufactured drugs with higher fentanyl potency.We coded those who said they intentionally used fentanyl on one or more days in the past month as having a fentanyl preference (1) and those who reported no intentional fentanyl use as having no preference for fentanyl (0).
We asked participants at baseline if they had ever injected drugs in their lifetime.Those who responded a rmatively were subsequently presented with a list of ve speci c drugs as well as an option to indicate any drug not listed and asked which drug or drugs they had injected in the past year.Participants who indicated injecting any drug in the past year were counted as being past-year injection drug users (1); those indicating they had never injected drugs or had not injected in the past year were counted as noninjectors (0).
Past-year overdose.The dependent variable in the SEM structural models was any past-year overdose.Participants were rst asked: As of today, how many times in your life would you say you have ever overdosed on drugs?Those who responded they had overdosed one or more times were asked for their most recent overdose and the number of times in the past year they had overdosed.Because this variable was highly right-skewed (skew = 6.43; kurtosis = 64.77),we created a dichotomous indicator representing any past-year overdose (1) versus no past-year overdose (0).Demographics.We obtained the following demographic information from participants at baseline: biological gender at birth (male or female); race (Black/African American; White, Other, Multi-racial); ethnicity (Latino/not Latino); and age in years.Race and ethnicity were coded into a single dichotomous variable represented as non-white (1) and white (0).

Analyses
We used Stata v.17.1 (44) to generate bivariate descriptive statistics and Mplus v.8.8 for the EFA, CFA, and SEM analyses (45).A series of bivariate statistics for both samples were generated, excluding 20 of the 150 EFA sample participants who completed a shortened interview that did not include demographic information.Signi cant differences were calculated within each sample, disaggregated by past-year overdose status (yes/no).Fisher's Exact Test was used to assess statistical differences among the nominal measures and t-tests for interval/ratio level measures.For the EFA, CFA, and SEM models we used a standard set of t statistics and thresholds to assess the adequacy of model t: comparative t index (CFI) > = .95;Tucker-Lewis index (TLI) > = .95;and root mean square error of approximation (RMSEA) < = .06)(46-48).
An EFA was run using data from ten drug use risk items collected from the rst sample as factor indicators.We used principal component extraction estimated by mean and variance-adjusted weighted least squares (WLSMV) to obtain the optimal factor structure based on chi-square comparisons of model t and retaining factors with eigenvalues > 1.0.Factor extraction was followed by geomin rotation to obtain loadings for each indicator.We next ran a CFA using the RCT-based sample, constraining the model based on the best-tting factor structure from the EFA run on the rst sample.The CFA analysis on the RCT sample also used WLSMV estimation.These analyses yielded a two-factor solution (see results).As one of the factors comprised items related to use of multiple drugs and mixing speci c other drugs such as benzodiazepines and alcohol with opioids, we added fentanyl preference to this factor and reran the CFA model using the RCT sample to con rm the model continued to t the data well.
We then assessed a series of SEM structural models.We began by including psychosocial and demographic variables as well as the two latent drug use risk factors identi ed in the EFA/CFA steps.All SEM models included the dichotomous indicator of past-year overdose as an endogenous variable.Our intent was to identify the best-tting, most parsimonious model that signi cantly predicted past-year overdose.
The rst model (model 1) included all the psychosocial and demographic variables enumerated above and the two identi ed latent risk factors.This model included the psychosocial and demographic variables as individual exogenous predictors of any past-year overdose and incorporated direct as well as indirect pathways mediated by the latent risk factors.For the second model (model 2), we removed the direct pathways from the psychosocial and demographic factors as none were signi cant and included only the indirect pathways mediated by the latent factors.In the next model (model 3), we reduced the number of psychosocial and demographic indicators by removing those that did not have a signi cant mediated path to past-year overdoses.For the fourth model (model 4) we removed the remaining demographic variables and restructured the psychosocial variables by modeling them as a single exogenous latent factor that directly affected the two EFA/CFA-identi ed risk factors.This model included a direct path from the psychosocial latent factor to past-year overdose as well as paths to the drugrelated latent factors to capture indirect effects on past-year overdose.In the nal model (model 5), we removed the direct path between the psychosocial risk factors and past-year overdose.The speci c paths and variables included in each of these models associated model t statistics, and the correlation matrix for the items analyzed for the nal model are provided as supplementary material.

Descriptive statistics
Table 1 shows the demographic, psychosocial, and drug use behaviors, including past year overdoses, for both study samples.Intrasample statistical comparisons between those reporting and those not reporting a past-year overdose yielded signi cant differences.For the EFA-sample, those not experiencing a pastyear overdose had experienced fewer lifetime overdoses (mean = 3.3, sd = 4.6) than those experiencing an overdose (mean = 5.8, sd = 4.9), and were less likely to report using drugs in a new location (46.2% versus 67.3%) or to have used drugs with alcohol (19.2% versus 42.2%).
Insert Table 1 about here Exploratory and rmatory Factor Analyses Table 2. shows the results for the EFA and CFA measurement model.Based on data from the rst sample, EFA yielded a 2-factor model as optimal based on t statistics and scree plot analysis.A chi-square test comparing the 1-and 2-factor models was statistically signi cant (chi-square = 44.3,p < .001)indicating the 2-factor model afforded a better t to the data.The 3-factor model failed to converge, suggesting it was too complex for the data available and hence, we did not consider a 4-factor model.Factor loadings for the rotated 2-factor model solution indicate the rst factor is represented by items related to drug acquisition and context: using in a new location, using with new people, and having a new source/supplier.Items loading more highly on factor 2 re ected drug use risk behaviors and, in particular, polydrug use such as using alcohol and/or benzodiazepines concurrently with opioids as well as using differing amounts, using alone, and using in a different way than usual.

Insert Table 2 about here
The CFA results, based on model 2 and obtained by constraining the measurement model to two factors identi ed by the highest loading items determined by the EFA had a good t: CFI = .972;TLI = .963;RMSEA = .046(95% CI = .020-.069).Because the highest loading items on factor 2 pertained to aspects of drug use per se, we added the variable representing fentanyl preference to this factor and reran the CFA.As the resulting model continued to show a good, even improved t -CFI = .978;TLI = .972;RMSEA = .037(95%CI = .000− .059)-we retained fentanyl preference as an indicator for the drug use factor in the subsequent structural models.

Structural Model
As indicated above, a number of alternative models were considered with details on each provided with the supplementary materials.Here, we focus on the nal model as presented in Fig. 1, which shows the standardized parameter estimates; unstandardized estimates are included with the supplementary materials.This model, in our estimation, best represented the data in terms of both t (CFI = .984,TLI = .981,RMSEA = .024,Chi-square (df=86) = 99.27,p = .155)and parsimony.In this model, the latent variable representing psychosocial risk factors was associated with increases in drug acquisition/contextual risks (=.683***,p = .001)as well as risky drug use behaviors (=.567***,p = .001).As mediated by both latent factors, psychosocial risks were associated with an increase in the probability of a past-year overdose re ected in the total indirect effect (=.234, p = .001).While there was a statistically signi cant association between drug acquisition and drug use risk behaviors (r = .421,p < .001),only drug use risk behaviors (=.287,p = .04)but not drug acquisition/context (=.105,p = .461)had a signi cant and positive direct association with any past-year overdose.In total, the model accounted for thirteen percent (r-square = .131,p = .02) of the variance in the probability of experiencing a past-year overdose.
Insert Fig. 1 about here

Discussion
Revisiting study hypotheses, the results supported that drug use risk behaviors are in uenced by psychosocial factors such as homelessness and that drug use risks directly in uence the probability of an overdose.The results also supported the hypothesis that psychosocial risk factors indirectly affect the probability of an overdose, mediated by drug use risk behaviors.Conversely, we did not nd that any demographic measure in uenced drug use risk or acquisition, nor were any associated with experiencing an overdose; our assessed models t the data well without including demographic measures.When such measures were included, there was often a decrement in model t among the alternative models considered.We did not have speci c hypotheses about the factor structure of the drug use risk behaviors considered, but the resulting two-factor structure whereby risk items relating to how drugs are acquired tended to cluster on one factor and risk items relating speci cally to how drugs are used clustering on a second made intuitive sense and also t the data well when reassessed using CFA.
Our results suggest psychosocial risks such as being homeless, having serious psychological distress, and trauma exposure, all of which were associated with each other, increased the chances that an opioid misuser will engage in riskier acquisitory/contextual behaviors such as buying drugs from a new dealer or using with new people in a new location.Psychosocial risks were also associated with riskier drug use behavior such as combining opioids with benzodiazepines or alcohol, using alone, or using in a different way than usual.There was, however, no direct association between psychosocial risk and experiencing an overdose.Instead, mediated by riskier drug acquisition and riskier drug-using behaviors, the psychosocial measures were indirectly associated with an increase in overdose probability.The model also indicated that risky drug use behaviors have an additional direct and positive association with an increased probability of past-year overdose.
We might speculate that the joint effect of the psychosocial risks considered as a latent factor represents instability, the effect of which is to narrow a person's options for consistently obtaining illegal opioids from a single supplier and for using with the same group of people in the same place.Homeless persons, who are often involuntarily displaced, are more itinerant than persons with stable housing and consequently are more likely, when acquiring and using drugs, to be among unfamiliar persons and to use in different locations (49).Increased social and environmental instability at a time when local fentanyl concentrations in drug supplies are known to temporally and geographically uctuate (22) could make it more likely a homeless person who uses illegal opioids obtains a "hot dose", with higher-thanexpected fentanyl potency.A recent study of the frequent involuntary displacement of homeless persons who inject drugs found that such displacements contributed to an estimated increase in mortality rates of 15.6-24.4%.(49) Based on our ndings, we expect at least some of the heightened mortality rates in this population is due to fatal overdoses.
In our sample, homelessness was associated with experiencing severe psychological distress and trauma and were more likely to engage in polydrug use, to use different amounts across occasions, and to use alone.Even absent a higher likelihood of obtaining a high-than-expected-potency fentanyl dose, using opioids with benzodiazepines or other depressants such as alcohol also increases the risk of an overdose.(50) In our model, such drug use risk behaviors had a direct and positive association with experiencing an overdose.Riskier drug acquisition and drug use behaviors among those with psychosocial risks combined to further increase the likelihood of an overdose beyond that resulting from riskier drug use behaviors alone.
Somewhat surprisingly, we did not nd signi cant associations among the assessed demographic risk factors.This was especially unexpected given the aforementioned consistent research ndings on higher rates of overdose among males compared with females.One thought is that studies that have identi ed biological gender as associated with higher overdose risk have not examined gender in combination with psychosocial, drug acquisition, and speci c drug use risks.It is possible that differences among males and females on the risk factors assessed in this study could account fully or partially for the gender differences found in prior studies.Unfortunately, we did not have a large enough sample to conduct a multi-group analysis of the model by biological gender to fully test this assertion.Past-year injection drug use was assessed but was also found to not improve model t when added to the drug use latent factor nor was it independently associated with experiencing an overdose.We are not sure of the reasons for this as we expected injection use to be associated with both the drug use latent risk factor as well as with experiencing an overdose.Given that the majority of our sample, about seventy ve percent, were injection users, the variation in injection use might have been too limited to detect an effect.

Limitations
As data for this study were collected in Chicago, Illinois, in the United States and from persons who use illegally manufactured and/or sold opioids, we do not know the extent to which the results apply to other locations and populations also at risk for an opioid-involved overdose such as those misusing prescription opioids.The sample size of 270 participants for the structural model was adequate but minimal and did not afford us the ability to test more complex models with interaction terms or additional latent factors or, as mentioned, to conduct multi-group analyses.All study measures are based on selfreport and could have been subject to over or under-reporting.We dichotomized past-year overdose as the dependent variable, which did not allow us to consider differences among persons who experience fewer or more overdoses over a given time.As noted by Kline, many alternative models can t the data equally or nearly equally as well when using SEM as an analytic technique.(51) Therefore, there is some subjectivity when selecting a "best-tting" model.Although we attempted to examine a number of alternative models before selecting the structural model we believed best represented the data, there are certainly other models that could have been conceived and tested and which might have been stronger statistically and/or conceptually.To remedy these limitations, other investigators might consider testing not only the model we developed but also alternative models with other data sets and samples.

Conclusions
To the of our knowledge, this study represents one of the rst attempts to assess the associations among opioid overdose risks across multiple functional domains: drug acquisition, drug use, and psychosocial.The results show these domains interact with persons who experience homelessness, have severe psychological distress, and who are trauma-exposed more likely to experience an overdose.We have speculated that the likely mechanism underlying the association between unstable social and living circumstances is that they drive acquiring and using drugs from unfamiliar sources and using in unfamiliar circumstances.Such individuals also appear to be more prone toward risky drug use behaviors, particularly polydrug use, that place them at higher risk of an overdose.We hope these ndings, as preliminary as they might be given the exploratory nature of the study, will help shape more comprehensive and more effective harm reduction approaches that reduce opioid-involved overdoses and fatalities, especially among persons who are unstably housed.