Abstract
Deviant peer association has been identified as a risk factor for substance use. Anxiety has been examined as a moderator of this relationship in a limited capacity. Despite this, investigations of specific dimensions of anxiety as moderators of this relationship remain understudied. This study examined three dimensions of anxiety as moderators: physiological anxiety, worry/oversensitivity anxiety, and social concerns/concentration anxiety. The Pathways to Desistance data were used in analyses. Ordered logistic regression models were used to assess relationships of interest for three substance use outcomes: heavy episodic drinking, marijuana use, and cigarette use. Findings indicated that social concerns/concentration anxiety moderated the relationship between deviant peer association and marijuana use, but, contrary to predictions, a protective effect at higher levels of this form of anxiety was observed. Worry/oversensitivity anxiety also appeared to exert a significant and negative direct effect on all three substance use outcomes. Implications are discussed.
Introduction
Deviant peer association has been identified by prior research as a risk factor predicting a range of substance use outcomes (Lee et al., 2017; Van Ryzin et al., 2012; Wojciechowski, 2020). While this is well established, there remain relevant questions regarding how this relationship may be conditioned by other factors. Anxiety is one factor that has been considered in this regard by prior research. Anxiety refers to a condition characterized by persistent unease and worry that may be connected to a discrete upcoming event or happening or may be more general in its presentation also. The condition may be extremely pervasive and may also manifest as a diagnosable mental illness if severe, with a lifetime morbidity risk around 9% in the United States (American Psychiatric Association, 2013). While prior research has indicated that anxiety may condition the relationship between deviant peer association and substance use (Anderson et al., 2011; Marshcall-Lévesque et al., 2014; Rodriguez et al., 2019), there is a dearth of research that has explored how different dimensions of anxiety may impact this relationship when considered in conjunction with one another. This study sought to bridge this gap in the extent literature by examining how three specific dimensions of anxiety may moderate the relationship between deviant peer association and substance use, those dimensions being: physiological anxiety, worry/oversensitivity anxiety, and social concerns/concentration anxiety.
Deviant Peer Association as a Risk Factor for Substance Use
As noted above, deviant peer association has been identified as a robust risk factor predicting substance use (Kim-Spoon et al., 2019; Price et al., 2019; Wojciechowski, 2018). Perhaps the most relevant criminological theoretical perspective aiding to explain this relationship is Akers’ (1973) social learning theory. This framework posits that exposure to individuals who engage in deviant behavior, like substance use, provide models through which individuals learn the motivations, rationalization, and techniques pertaining to engagement in deviant behavior. The theory is general in the sense that prosocial behaviors are transmitted through a similar learning process. Individuals may provide definitions that are favorable or unfavorable toward a specific behavior. Exposure to more definitions that are favorable toward a behavior then should increase the likelihood that an individual will engage in the given behavior themselves. For substance use, this may entail exposure to peers use drugs or alcohol themselves and present models for engagement in the behavior for a given target individual. There also exist other mechanisms by which peers may influence behavior in this regard. Processes pertaining to reinforcement or punishment of a behavior may also increase or decrease the likelihood of future engagement in said behaviors. Reinforcement refers to any action or potential action that may be taken by other individuals to increase the likelihood of future engagement in a behavior. For example, peers may praise or encourage a target individuals first engagement in substance use and this may lead the target individual to feel rewarded for engagement in the behavior, thus, increasing the likelihood that they will do so again in the future. Punishment is the opposite of this process, as it refers to actions or potential actions that may reduce the likelihood of future engagement in a behavior. Prosocial peers may indicate a reticence to maintain relationships with an individual upon discovering that they have been engaging in substance use. The potential for losing these friendship ties then may act as a threatened punishment that should result in lower risk for the individual engaging in substance use in the future. It is through these highlighted mechanisms that peer relationships are relevant for understanding behavior and why deviant peer ties may specifically be important for engaging in antisocial behavior like substance use. Prior research has indicated support for the relevance of these mechanisms for predicting substance use (Pratt et al., 2010; Kruis et al., 2020). For example, Pratt et al., 2010 provide a meta-analysis indicating that all mechanisms of the theory are relevant for predicting behavior, though differential association and modeling/imitation seem most relevant. However, there remains continued need to understand how individual characteristics may amplify or attenuate the salience of deviant peer association for influencing engagement in substance use. One key factor here may be anxiety.
Anxiety Dimensions as Moderators of the Relationship Between Deviant Peer Association and Substance Use
The persistent unease and worry that characterizes the mental condition of anxiety may present an important psychological factor for understanding the relationship between deviant peer association and substance use. It might be expected that adolescents who report greater anxiety may have exacerbated response to potential social punishments from deviant peers. Individuals may perceive that deviant friends may have limited desire to maintain friendships if they themselves to not engage in antisocial behavior like substance use. For individuals high in anxiety, oversensitivity to such concerns may then exacerbate the risk for substance use as they persistently dwell on this worry. While this would make sense, these relationships as they pertain to general anxiety remain understudied. In the case of general anxiety, prior research has even found that anxiety may moderate this relationship in the opposite direction. Rodriguez et al. (2019) found that peer substance use did not influence the substance use behaviors of target adolescents who were high in anxiety. Rather, peer substance use only influenced the substance use of adolescents who reported low anxiety. These relationships also must be understood within the context that some research has suggested that anxiety may also act as a protective factor against substance use (Kaplow et al., 2001; Myers et al., 2003). For example, Kaplow et al. (2001) found that while generalized anxiety was a risk factor for alcohol use, symptoms of separation anxiety was associated with lower risk. This study highlights the fact that, while anxiety may represent a broad range of related symptoms, there do exist distinct dimensions of anxiety and the specific nature of the anxiety may matter for how it impacts the relationships of interest; like that of the relationship between deviant peer association and substance use. It may be that some dimensions of anxiety act as a protective factor against the influence of peers on substance use, whereas other dimensions act in the manner posited above leading to greater risk for substance use. In a situation like this, examining general anxiety as a moderator may be obscuring the moderation effects of more specific dimensions that could act as a moderator in the manner proposed above. This necessitates delineation of these specific dimensions of anxiety in order to determine if certain forms of anxiety protect against peer influence in this regard and if other dimensions exacerbate the effects of peer influence.
Examinations of the roles that specific dimensions of anxiety play for moderating the relationship between deviant peer association and substance use has also been limited. The main branch of research in this area has focused on social anxiety as a moderator specifically. This makes sense considering that social anxiety should directly correspond to the reactions that individuals may have to peers attempting to influence their behavior. If individuals are concerned about how they might interact with peers, this may also lead them to engage in substance use as a means of coping or alleviating perceived social awkwardness. Having greater levels of deviant peer association may then facilitate substance use and make this behavior a realistic option for reducing this form of anxiety during social situations. Findings from studies focused on social anxiety in this role generally indicate that social anxiety does moderate the relationship between deviant peers and substance use in the expected manner, with high levels of both social anxiety and deviant peer association predicting elevated substance use risk (Anderson et al., 2011; Cohen & Prinstein, 2006). For example, Cohen & Prinstein (2006) found that greater social anxiety resulted in greater conformity to peer norms and that this manifested in greater substance use. There may also be gendered effects at play here, as Anderson et al. (2011) found a stronger peer effect on binge drinking specifically for boys who were high in social anxiety. This provides indication that while general anxiety may generate a protective effect against substance use and the influence that peers may exert on substance use, generalized measurement may obscure the increased risk for substance use associated with social anxiety specifically. Past research has also indicated that social anxiety may have a direct effect resulting in increased substance use risk also (Mueller et al., 2021; Villarosa-Hurlocker et al., 2019). For example, Mueller et al., (2021) found that social anxiety increased risk for marijuana use problems and that this relationship was stronger for males. While prior research provides some indication of this, there remains a dearth of research on the other specific forms of anxiety of interest here: worry/oversensitivity anxiety and physiological anxiety.
Worry/oversensitivity anxiety refers to symptoms characterized by vague obsessive concerns about may be nonspecific in nature (e.g., “what if something bad happens?”). This is delineated from social anxiety in the sense that these concerns are not distinctly focused on anxiety driven by concerns regarding interpersonal interaction. Physiological anxiety here refers to the actual physical symptoms that may accompany physical or tangible symptoms of anxiety (e.g., sweating palms, waking up scared, bad dreams, etc.). In the case of worry/oversensitivity anxiety, overwhelming concerns regarding potential consequences or issues stemming from substance use may lead individuals to have lower risk for engagement in substance use behaviors. Such worries may be related to getting into trouble for engagement in illicit behavior, as this may pertain to a range of potential consequences (ex: going to jail, being suspended from school, being grounded by parents, etc.). Because these concerns would not be directly related to the potential for social awkwardness or the desire to feel less inhibition in social situations, this may actually lead individuals to be a protective factor against substance use. In the context of deviant peers, this may also play a protective role, as individuals high in worry/oversensitivity anxiety may be buffered against the effects of deviant peers on their substance use risk. In the case of physiological anxiety, potential moderation effects pertaining to deviant peer association and substance use seem less clear. While individuals may attempt to use substances as a means of coping with these symptoms in a manner consistent with the self-medication hypothesis (Khantzian, 1987, 1997), it is unclear that this would relate to deviant peer association. It would seem to depend on if these physiological symptoms are connected in any way to either of the other forms of anxiety. These physiological symptoms may then be manifestations of either worry/oversensitivity anxiety or social concerns/concentration anxiety. If this is the case, then it would seem that the other forms of anxiety may be more relevant for moderating the relationship between deviant peer association and substance use. Physiological anxiety then should not actually act as a means of moderating this relationship, as this should just be a manifestation of anxiety that corresponds to other types. In this way, it would be expected that physiological anxiety should not significantly moderate the relationship between deviant peer association and substance use.
The moderating role of anxiety for understanding the relationship between deviant peer association and substance use has been examined by prior research. However, findings regarding the nature of this relationship have been mixed (Cohen & Prinstein, 2006; Rodriguez et al., 2019). For example, Cohen and Prinstein (2006) found that youth high in social anxiety were susceptible to influence of both high and low status peers, whereas those low in social anxiety were only susceptible to high status peers; providing indication that high anxiety may increase susceptibility to peer influence. Rodriguez et al. (2019) found the opposite, as youth high in anxiety were not as vulnerable to peer influence on substance use compared to youth with low levels of anxiety. This indicates the need for continued examination of specific anxiety dimensions as moderators of the relationship between deviant peer association and substance use. It may be that general anxiety measures may be too coarse, as worry/oversensitivity anxiety and physiological anxiety may be masking or washing out amplification effects of social concerns/concentration anxiety. The present study sought to add to the extant literature by testing the following hypotheses:
Hypothesis 1: Greater worry/oversensitivity anxiety will predict lower substance use frequency.
Hypothesis 2: Greater social concerns/concentration anxiety will predict greater substance use frequency.
Hypotheses 3: Worry/oversensitivity anxiety will moderate the relationship between deviant peer association and substance use, resulting in a protective effect at greater levels of worry/oversensitivity anxiety.
Hypothesis 4: Social concerns/concentration anxiety will moderate the relationship between deviant peer association and substance use, resulting in an amplification effect at greater levels of social concerns/concentration anxiety.
Hypothesis 5: Physiological anxiety will be unrelated to deviant peer association and substance use both in terms of direct and moderation effects.
Methods
Data
Data utilized in analyses were from the first two waves of the Pathways to Desistance study. This dataset was comprised of the responses of 1354 justice-involved youth who then provided data for the following 7 years after being adjudicated for a serious offense prior to baseline measurements; resulting in a total of 11 waves of data. Serious offenses which met inclusion criteria consisted of all felony offenses, as well as sexual assault and weapons-related misdemeanors. Offenses which qualified participants for inclusion in the study had to have been committed when participants were between the ages of 14 and 17 in order to meet inclusion criteria. The first two waves were utilized because they had the least missing data and provided the most parsimonious means of testing hypotheses while also establishing temporal ordering between the independent and dependent variables as a necessary condition for inferring causality. Participants were recruited from study sites located in Maricopa County, Arizona and Philadelphia, Pennsylvania. The proportion of male drug offenders included in the sample was capped at 15% in order to maintain heterogeneity in these baseline characteristics. Of all qualified youth approached regarding their interest in the study, 20% declined the opportunity to participate. Of the baseline sample, 93.4% were retained and provided data for wave two measurements. There was a 6-month observation period between baseline measurements and wave two measurements.
Data used in this study were collected via participant self-report. Interviews were conducted in locations that were convenient for participants (e.g., libraries, criminal justice facilities, participants’ homes, etc.). The research team provided each participant with a laptop computer during interviews. Participants then manually input responses to verbal prompts from the research team into the laptop. This method of data collection was chosen in order to maximize confidentiality in reporting on sensitive topics.
Measures
Substance Use
The main set of dependent variables examined in this study was focused on the frequency of engagement in several different substance use behaviors at wave two: heavy episodic drinking (HED), marijuana use, and daily cigarette use. 1 HED and marijuana use at wave two were measured using the same nine-point ordinal scale which assessed the general frequency at which participants reported engaging in either behavior during the prior 6 months (0=Not at all; 1=1–2 times; 2=3–5 times; 3=Once a month; 4=2–3 times per month; 5=once per week; 6=2–3 times per week; 7=4–5 times per week; 8=everyday). Daily cigarette use was assessed at wave two using a six-point ordinal scale that asked participants to report about how many cigarettes that they smoked per day during the prior 6 months (0=0 cigarettes; 1=1–4 cigarettes; 2=5–9 cigarettes; 3=10–15 cigarettes; 4=16–20 cigarettes; 5=More than a pack a day). These dependent variables were measured at wave two so that baseline independent variables were assessed prior to substance use, thus, ensuring the temporal ordering necessary for causal inference.
Deviant Peer Influence
One of the key independent variables examined in analyses was deviant peer influence at baseline. This construct was measured using a scale adapted from the Rochester Youth Development Study (Thornberry et al. 1994). This scale was a measure of the general number of peers who participants reported attempt to influence them to engage in antisocial behavior. Seven ordinal items were used to assess how many friends that participants reported having that influenced them to engage in seven different antisocial behaviors. A mean score was then computed from these individual items so that each participant had a single deviant peer association score at baseline. The following items were utilized for the deviant peer association measure: How many of your friends have suggested that you should go out drinking with them?; How many of your friends have suggested or claimed that you have to get drunk to have a good time?; How many of your friends have suggested or claimed that you have to be high on drugs to have a good time?; How many of your friends have suggested that you should sell drugs?; How many of your friends have suggested that you should steal something?; How many of your friends have suggested that you should hit or beat someone up?; How many of your friends have suggested that you should carry a weapon? Cronbach’s alpha indicated good reliability for this measure (.89).
Anxiety
The other set of key independent variables included in analyses focused on different forms of anxiety that participants reported experiencing. Three different forms of anxiety were measured: physiological anxiety, worry/oversensitivity, and social concerns/concentration. The physiological anxiety items focused on somatic manifestations of anxiety. The worry/oversensitivity items focused on internal anxiety about fears related to being hurt or emotionally isolated that may sometimes be vague and difficult to express clearly. Finally, the social concerns/concentration items focused on feelings of anxiety centered distracting thoughts and fears that have a social and/or interpersonal focus. The Revised Children’s Manifest Anxiety Scale was used to assess each type of anxiety at baseline. This consisted of 37 items that asked participants to report whether or not they experienced individual symptoms of anxiety using binary measures. These binary measures allowed participants to report whether or not each symptom was currently descriptive of their own feelings or actions related to the experience of anxiety. Count scores were then computed for each type of anxiety, providing a count of the number of symptoms that participants endorsed for each anxiety dimension. This then resulted in three anxiety scores for each participant at baseline that were included in analyses. The following items were used to assess physiological anxiety: I often have trouble making up my mind; Often I have trouble getting my breath; I get mad easily; It is hard for me to get to sleep at night; Often I feel sick in my stomach; My hands feel sweaty; I am tired a lot; I have bad dreams; I wake up scared some of the time; I wiggle in my seat a lot. The following items were used to assess worry/oversensitivity anxiety: I get nervous when things do not go the right way for me; I worry a lot of the time; I am afraid of a lot of things; I worry about what my parent will say to me; I worry about what other people think about me; My feelings get easily hurt; I worry about what is going to happen; My feelings get easily hurt when I am fussed at; I worry when I go to bed at night; I am nervous; I often worry about something bad happening to me. The following items were used to assess social concerns/concentration anxiety: Others seem to do things easier than I can; I feel that others do not like the way I do things; I feel alone even when there are people with me; Other people are happier than I; I feel someone will tell me I do things the wrong way; It is hard for me to keep my mind on my school work; A lot of people are against me. Cronbach’s alpha scores for these variables indicated adequate reliability for worry/oversensitivity anxiety (.81), but somewhat suboptimum reliability for the physiological anxiety (.66) and social concerns/concentration anxiety measures (.69).
Control Variables
Several control variables were also included in analyses so that risk for bias in estimation may be mitigated. The first of these variables was gender, as prior research has indicated greater differences in substance use risk (McHugh et al., 2018). Gender was assessed at baseline using a binary variable that delineated male and female participants (0=Male; 1=Female).
Race was also included as a control variable, as prior research has indicated that White youth may present greater substance use risk than youth identifying in other Non-White racial groups (Vasilenko et al., 2017). Race was assessed at baseline using a four-category nominal variable with the following racial identification options: Black, Hispanic, White, and Other Race. A series of four dummy variables was then coded, with each dummy variable delineating participants in a single race category from all other participants and was repeated for each race category (e.g., 1=Black; 0=All other participants. The dummy variable corresponding to White participants was then omitted from analyses in order to provide a reference group for interpreting coefficients for the other race variables in comparison with.
Socioeconomic status (SES) was another control variable included in analyses, as prior research has indicated that substance use risk is stratified by SES (Vilsaint et al., 2019). SES was measured at baseline utilizing Hollingshead’s (1957) two-factor index of social position. This instrument measures SES as a weighted score comprised of participants’ parents’ educational attainment and occupational prestige scores. If both parents were available to provide data, then a mean score was computed so that all participants had a single baseline SES score.
Another control variable included in analyses was exposure to violence. This is because past research has indicated that exposure to violence is a risk factor for substance use (Mitchell et al., 2007). This construct was measured at baseline using the Exposure to Violence Inventory (Selner-Ohagan et al., 1998). This scale assessed whether or not participants had ever experienced different forms of direct victimization and/or witnessed violence prior to baseline. A count score was then generated that provided a measure of the number of different types of exposure to violence that each participant experienced prior to baseline.
Impulse control at baseline was also controlled for, as prior research has indicated that low impulse control is a risk factor predicting substance use (Davis et al., 2017). This was measured using the Weinberger Adjustment Inventory (Weinberger et al., 1989). This scale was comprised of a series of ordinal items that asked participants to rate the degree to which a set of statements related to impulse control was true of their own behavior or attitudes (e.g., I say the first thing that comes into my mind without thinking enough about it). Seven of the eight individual items were reverse coded so that higher scores predicted higher impulse control. A mean score was also computed so that all participants had a single impulse control score at baseline. Cronbach’s alpha indicated adequate reliability for this measure (.76).
Age was also included as a control variable, as prior research has indicated that substance use risk may be age-graded (Vasilenko et al., 2017). Age was measured at baseline in single-digit year intervals.
Time spent in secured facilities with no community access was also controlled for (e.g., jails, prisons, psychiatric hospitals), as spending more time without community access during the wave two observation period may have impacted one’s capacity to engage in substance use. This variable was operationalized as a proportion ranging from 0 to 100, with higher scores indicating more time spent in secured facilities during the wave two observation period (e.g., 55 = 55% of the wave two observation period spent in secured facilities).
Wave two observation period length was also controlled for, as longer observation periods offered increased exposure time with which to engage in substance use and this may have impacted rates of use. While all participants’ wave two observation periods were generally 6 months in length, there was some variance in the exact number of days. As such, an indicator of the exact number of days in each participant’s wave two observation period was included in analyses.
Analytic Strategy
Ordered logistic regression was utilized to examine each substance use outcome. This modeling strategy was chosen because of the ordinal measurement of each of the dependent variables. Two models were estimated for each dependent variable. Model 1 for each variable examined the direct effects of the baseline deviant peer association and anxiety measures for predicting substance use. Model 2 for each dependent variable then included an interaction term that modeled the moderating effect of all three anxiety dimensions on the relationship between deviant peer association and substance use. These interaction terms were computed by multiplying each participant’s deviant peer association score by their scores on each anxiety dimension. The deviant peer association and all anxiety measures were also mean-centered for these analyses by subtracting the mean value from individual participants’ values in order to address concerns of multicollinearity in calculating interaction effects. Listwise deletion was used to manage missing data. This missing data management strategy was chosen for several reasons. First, preliminary analyses indicated that there was not a strong correlation between predictors and missingness, with only about 7% of participants from the original sample no longer providing data at wave 2. This 7% missingness was observed fairly evenly across all measures included in analyses, indicating that there was not one specific measure that had significantly more missing data than others. Further, some research has indicated that listwise deletion may actually provide the most accurate estimation of a full dataset for criminal justice populations (Mitchell et al., 2021). Coefficients are described in the form of odds ratios (ORs), indicating the impact that a one-unit increase in an independent variable of interest had on the odds of being in a higher category on the ordinal scale for each substance use outcome. 2
Results
Descriptive statistics for all variables included in analyses are described in Table 1. Table 2 provides a correlation matrix for all variables included in analyses. Table 3 provides regression results pertaining to Model 1 and Model 2 for the binge HED outcome. Table 4 provides these same Model 1 and Model 2 regression results for the marijuana use outcome. Finally, Table 5 provides regression results for Model and Model 2 for the daily cigarette smoking outcome.
Descriptive Statistics.
Correlation Matrix.
=p<.05
Ordered Logistic Regression Covariate Effects on Heavy Episodic Drinking Frequency in Odds Ratios (ORs).
Note. Two-tailed p-value: *p≤05; **p≤.01; ***p≤.001.
Ordered Logistic Regression Covariate Effects on Marijuana Use Frequency in Odds Ratios (ORs).
Note. Two-tailed p-value: *p ≤ 05; **p ≤ .01; ***p ≤ .001.
Ordered Logistic Regression Covariate Effects on Daily Cigarette Use in Odds Ratios (ORs).
Note. Two-tailed p-value: *p ≤ 05; **p ≤ .01; ***p ≤ .001.
Preliminary analyses utilized a Brant test to determine whether the parallel odds assumption was met for all variables included in analyses for each dependent variable. These Brant tests indicated that the overall models did not meet the parallel odds assumption, but that there were also no individual variables in any of the models that violate the parallel odds assumption on their own. Rather, it was the sum of all of the variables included in the model that led to violation of the proportional odds assumption. This indicated that using a generalized ordinal logistic regression model would simply provide the same coefficient estimates as a standard ordinal logistic regression model, since no individual predictor variables will be estimated separately because of their violation of the proportional odds assumption. For this reason, standard ordinal logistic regression was utilized for all analyses as previously planned.
Model 1 results for the HED dependent variable indicated that increased deviant peer association at baseline predicted increased HED frequency at wave two (OR=1.591; Standardized coefficient=5.133; p < .001). Greater oversensitivity/worry anxiety was associated with lower HED frequency also (OR=.862; Standardized coefficient=−4.348; p < .001). Greater physiological anxiety was associated with increased HED frequency (OR = 1.101; Standardized coefficient=.097; p < .029). Social concern/concentration was not significantly associated with HED frequency in this model. Being male, being older, lower impulse control, spending less time in secured facilities, being recruited from Maricopa County, and longer wave two observation periods were also all associated with increased HED frequency in this model. Black participants reported lower HED frequency compared to White participants in this model. Model 2 results indicated that greater deviant peer association still predicted increased HED frequency following inclusion of hypothesized interaction terms (OR=1.650; Standardized coefficient=5.010; p < .001). Greater physiological anxiety predicted increased HED frequency in this model (OR=1.119; Standardized coefficient=.112; p < .013). Greater worry/oversensitivity anxiety predicted lower HED frequency in this model (OR=.868; Standardized coefficient=−4.060; p < .001). Social concern/concentration anxiety did not significantly predicted HED in this model. None of the hypothesized interaction terms were significant in this model either. Being male, lower impulse control, being older, longer wave two observation periods, being recruited from Maricopa County, and spending less time in secured facilities all predicted increased HED frequency also. Black participants reported lower HED frequency than White participants in this model.
Model 1 results for the marijuana use dependent variable indicated that greater deviant peer association predicted increased marijuana use frequency at follow-up (OR=1.633; Standardized coefficient=6.312; p < .001). Increased levels of worry/oversensitivity anxiety predicted decreased marijuana us frequency at follow-up (OR = .873; Standardized coefficient=−4.598; p < .001). Increased levels of social concerns/concentration anxiety predicted increased marijuana use frequency at wave two in this model (OR = 1.168; Standardized coefficient=1.358; p < .001). Physiological anxiety was not a significant predictor of marijuana use frequency in this model. Less time spent in secured facilities and longer wave two observation periods also all predicted increased marijuana use frequency in this model. Model 2 results for the marijuana use dependent variable indicated that greater deviant peer association continued to predict increased marijuana use frequency after inclusion of the interaction terms (OR = 1.722; Standardized coefficient=6.651; p < .001). Worry/oversensitivity anxiety and social concerns/concentration anxiety also continued to exert significant effects on marijuana use frequency in the same directions observed in Model 1 (Worry/oversensitivity OR=.877; Worry/oversensitivity standardized coefficient=−4.384; p < .001; Social concerns/concentration OR=1.182; Social concerns/concentration standardized coefficient=3.814; p < .001). Physiological anxiety remained a nonsignificant predictor of marijuana use frequency in this model. None of the interactions were significant. However, the interaction between deviant peer association and social concerns/concentration anxiety was negative and nearly met the threshold for statistical significance (OR=.918; Standardized coefficient=−1.885; p < .059. Longer wave two observation periods and spending more of the wave two observation period in secured facilities also predicted increased marijuana use frequency in this model.
Model 1 results for the daily cigarette use dependent variable indicated that deviant peer association was a significant predictor of this outcome (OR=1.243; Standardized coefficient=3.082; p < .002). Physiological anxiety was not a significant predictor of daily cigarette use in this model. Higher levels of worry/oversensitivity anxiety were found to predict lower average daily cigarette use in this model (OR=.909; Standardized coefficient=−3.548; p < .001). Higher levels of social concerns/concentration anxiety were found to predict greater average daily cigarette use in this model (OR=1.144; Standardized coefficient=3.439; p < .001). Being female, being older, spending less time in secured facilities during the wave two observation period, being recruited from Philadelphia, and longer wave two observation periods also all predicted increased average daily cigarette use in this model. Black, Hispanic, and Other Race participants all reported lower average daily cigarette use in this model compared to White participants. Model 2 results for this outcome indicated that deviant peer association was a significant predictor of daily cigarette use (OR=1.235; Standardized coefficient=2.934; p < .003). Greater worry/oversensitivity anxiety predicted lower average daily cigarette use (OR=.907; Standardized coefficient=−3.613; p < .001). Greater social concerns/concentration anxiety scores predicted increased cigarette use (OR=1.144; Standardized coefficient=3.416; p < .001). Physiological anxiety was not a significant predictor of average daily cigarette use in this model. None of the hypothesized interactions were significant in this model either. Being female, being older, being recruited from Philadelphia, less time spent in secure facilities, and longer wave two observation periods were also all associated with increased daily cigarette use in this model. Black, Hispanic, and Other Race participants all reported less average daily cigarette consumption relative to White participants in this model.
Sensitivity analyses were estimated to determine the robustness of results. There were concerns that inclusion of all of the anxiety dimensions and their accompanying interaction terms together in the same model may result in collinearity leading to diminished capacity to identify significant results. As such, these sensitivity analyses examined each anxiety dimension as the sole anxiety predictor of each outcome and then estimated the hypothesized interaction as the only interaction in the model also. Findings for physiological anxiety indicated that this dimension was a nonsignificant predictor of HED in both models, inconsistent with Model 2 in the main analyses. However, the interaction with deviant peer association was significant and negative, indicating that the effects of deviant peer association were greatest when physiological anxiety were lower. These sensitivity analysis findings were mirrored for the marijuana use outcome and physiological anxiety. All other effects of physiological anxiety on cigarette use were nonsignificant and consistent with main findings. Social concerns/concentration anxiety significantly interacted with deviant peer association to predict marijuana use, with a negative coefficient which was inconsistent with the main analyses. This same pattern of interaction was observed for worry/oversensitivity anxiety and marijuana use and HED. None of the other findings pertaining to worry/oversensitivity anxiety and social concerns/concentration anxiety differed from the findings observed in the main analyses.
The near significant interaction of deviant peer association with social concerns/concentration anxiety was modeled to provide better understanding of the specific interaction effect. Figure 1 provides visual depiction of these marginal effects. This figure indicates that individuals high in deviant peer association did report greater probability of marijuana use than individuals with low deviant peer association across the lower range of social concerns/concentration anxiety scale. However, the low deviant peer association group overtook the high deviant peer association group in this regard at the higher levels of the anxiety score scale. That said, the 95% confidence intervals cross throughout the entire social concerns/concentration anxiety scale; indicative of the fact that this interaction did not quite reach the threshold for statistical significance in the main analyses. However, this does provide indication of the general direction of this interaction, with greater deviant peer association providing a protective effect on marijuana use at higher levels of the social concerns/concentration anxiety score scale.

Marginal effects of social concerns/concentration anxiety on predicted probability of marijuana use delineated by deviant peer association.
Discussion
Findings from this study indicate the relevance of distinct anxiety dimensions for understanding the relationship between deviant peer association and substance use. A key finding indicated that only social concern/concentration anxiety was significantly moderated this relationship and only for marijuana use. However, this moderation effect was opposite of the manner that was hypothesized, as the greatest effect of deviant peer association on marijuana use frequency appeared to be observed for individuals with lower levels of social concern/concentration anxiety. It should also be noted that this moderation effect did not quite reach the threshold for statistical significance (p < .059). It also cannot be ignored that sensitivity analyses identified similar negative moderation effect for several outcomes and the various forms of anxiety also. It should also be noted that Hypothesis 1 was supported, as greater levels of worry/oversensitivity anxiety robustly predicted lower frequency of use for all three substances. There are several important implications for these findings for public health and criminal justice researchers and professionals concerned with reducing substance use risk among adolescents.
There were several significant direct effects of interest observed here. Greater worry/oversensitivity anxiety was associated with lower substance use frequency. Given that this effect was robust across all three substances, this may indicate that past research that has found anxiety to be a protective factor against substance use may be driven by this specific form of anxiety. Greater levels of social concerns/concentration anxiety were found to predict increased marijuana use only and greater physiological anxiety was found to predict increased HED frequency only. This indicated that, while these effects may be substance-specific, the effects of certain forms of anxiety on substance use may be being masked when employing general measures of anxiety. In terms of treatment implications, this indicates the relevance of screening not only for anxiety, but also for screening to identify the presence of specific dimensions. In the case of justice-involved youth, this may occur at intake into the juvenile/criminal justice system and may result in mandated treatment while under supervision, whether that be community corrections or incarceration. Doing so may allow for both more targeted treatment, but also for identifying which youth may present greater or lesser risk for substance use upon reentry. The effectiveness of such measures, however, remains speculative. Additional research is needed to determine whether or not more specific screening of anxiety presents an improvement in terms of being able to predict justice-involved youth who may be at-risk for substance-related recidivism.
The effect of deviant peer association on marijuana use was moderated by social concerns/concentration anxiety, though the coefficient fell just below the threshold for statistical significance (p < .059). Further, contrary to expectations, this moderation effect was negative in direction. This seems to indicate that greater levels of social concerns/concentration anxiety were protective against the impact on deviant peer association on marijuana use. This finding is inconsistent with prior research that indicates that social anxiety may amplify this relationship (Anderson et al., 2011; Cohen & Prinstein, 2006). This is particularly interesting because this form of anxiety exerted a significant and positive direct effect on marijuana use. One potential explanation for this finding may be that marijuana use may exacerbate these anxiety symptoms when use. One commonly reported effect of marijuana is that it may lead to anxiety felt by the user (Sharpe et al., 2020). It may be that, when used in social situations, marijuana use may actually exacerbate social concerns/concentration anxiety for the justice-involved youth in this sample. If this is the case, then this may lead to infrequent use among youth who have many deviant peers who would influence individuals into using the substance. It may be that the indicated nature of this sample plays a role also, as these social concerns may extend to other interpersonal relationships that these individuals have that may be justice system related. For example, exacerbation of social anxiety that extends to probation officers or law enforcement may lead individuals to be overwhelmed when pressured by peers to use. Because marijuana use was illegal in both study sites during data collection, the increased illicitness of marijuana use compared to the other substances of interest may also explain why this was only observed for this substance. This may then lead these individuals to abstain from use among these peers, thus, leading to diminished use that would be consistent with the moderation effect of interest. In terms of the significant and positive direct effect of social concerns/concentration anxiety on marijuana use, it may be that these anxiety symptoms are not exacerbated if pressure or influence from peers is not present and individuals may use marijuana in a manner that is more consistent with a self-medication role. While this remains speculative, this provides some explanation for why this inconsistent moderation effect was observed here. Additional research is clearly necessary in this area to determine why the effects observed in this study differed from those observed in prior studies regarding the moderating role of social anxiety for understanding the relationship between deviant peer association and substance use. This is further tempered by the fact that this effect was not statistically significant. For this reason, this result should be interpreted with caution.
One final issue observed in this study pertains to the role of the various forms of anxiety as they interacted with deviant peer association to influence the different substance use outcomes in sensitivity analyses when examined as the sole anxiety dimensions in the model. Findings in this regard were mixed then when comparing main and sensitivity analyses. These mixed findings make interpreting the role of these symptoms of anxiety somewhat difficult. It may be that the purpose of this study in trying to parse the effects of individual forms of anxiety net of one another was misplaced then, as the effects may have washed one another out. This indicates the potential need to re-examine these analyses while only utilizing single forms of anxiety as predictors and examining interactive effects in this manner. Future research should seek to test the robustness of these findings in this way.
While the present study provided a novel examination of the role of specific anxiety dimensions for understanding the relationship between deviant peer association and substance use, there remain a number of noteworthy limitations. The first of these limitations pertains to the potential lack of generalizability of these findings. This sample was collected purposively. This nonprobability sampling method could potentially lead to a nonrepresentative sample and biased results. Additionally, this sample was comprised solely of justice-involved youth, making generalizability to the general population of adolescents difficult due to differences between these populations. This indicates the need for replicating these analyses using a probability sample from the general population of adolescents in order to determine the robustness of these findings. Another limitation of this study pertains to the fact that the measures used in analyses were non-diagnostic metrics of anxiety and substance use. While the anxiety symptom count measures provided a useful means of understanding these anxiety dimensions on a continuous scale, their clinical significance is not established here. Similarly, while frequency measures of substance use provide a coarse means of understanding severity of use, these measures do not provide a diagnostic measure of substance use disorder that would actually indicate dysfunction stemming from use that would have psychiatric significance. For these reasons, discussion of clinical implications for treatment is somewhat limited here. This indicates the need for continued research on these topics using clinical measures in order to provide psychiatric professionals and researchers evidence that may inform the design and implementation of psychosocial programming and pharmaceutical avenues for treating psychiatric populations with these issues. Another limitation pertains to the inability to utilize all 11 waves of the Pathways to Desistance data. While all 11 waves were available, the anxiety subscales used in these analyses were only measured at baseline. Because adolescence is a period of the life-course characterized by psychological growth and development, it seemed inapplicable to try to model relationships using a time-invariant measure of anxiety as a key independent variable, as anxiety symptoms likely changed within participants during this time. This informed the approach used in this study that at least allowed for fulfilling the necessary condition of temporal ordering for inferring causality. However, the fact that all 11 waves could not be utilized to take a c=more complex and nuanced approach is indeed a limitation of the study. Future research should seek to extend these findings of this study to take such an approach utilizing richer datasets that allow for more complex longitudinal analyses.
Footnotes
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) received no financial support for the research, authorship, and/or publication of this article.
