A latent class analysis of patterns of tobacco and cannabis use in Australia and their health-related correlates.

Introduction: The shifting landscape in Australia's tobacco and cannabis policies and emerging new products and modes of administration may increase experimentation and the risks of addiction to these drugs. Methods: We analysed cross‐sectional data from the 2019 National Drug Strategy and Household Survey (n = 22,015) of Australians aged 14 and above. Latent class analysis was used to identify distinct groups based on types of tobacco and cannabis products used. The socio‐demographic, health‐rated correlates and past‐year substance use of each latent class was examined. Results: A four‐class solution was identified: co‐use of tobacco and cannabis (2.4%), cannabis‐only (5.5%), tobacco‐only (8.0%) and non‐user (84.0%). Males (odds ratio [OR] range 1.5–2.9), younger age (OR range 2.4–8.4), moderate to high psychological distress (OR range 1.3–3.0), using illicit substances in the last year (OR range 1.41–22.87) and high risk of alcohol use disorder (OR range 2.0–21.7) were more likely to be in the tobacco/cannabis use classes than non‐users. Within the co‐use class, 78.4% mixed tobacco with cannabis and 89.4% had used alcohol with cannabis at least once. Discussion and Conclusions: Approximately 16% of respondents used tobacco or cannabis, or both substances, and no major distinct subgroups were identified by the use of different product types. Mental health issues and the poly‐substance use were more common in the class who were co‐users of cannabis and tobacco. Existing policies need to minimise cannabis and tobacco‐related harms to reduce the societal burden associated with both substances.

Keywords: cannabis; cigarettes; joints; marijuana; tobacco

Tobacco and cannabis are the most commonly used psychoactive substances after alcohol in Australia [[1]]. These substances have different health impacts and population‐level disease burdens. In 2018, 8.6% of the total disease burden was attributable to tobacco smoking while cannabis use accounted for 0.3% of the total disease burden in Australia [[2]]. A key prevention strategy to reduce the disease burden is to understand the key characteristics of the consumers of single or both substances and the factors related to their patterns of consumption.

Co‐use of tobacco and cannabis is prevalent in many countries including Australia [[3]], with co‐use practices differing by country and region [[5]]. Co‐use of tobacco and cannabis can involve sequential use, where tobacco and cannabis are used separately within the same use episode, or simultaneous use when both substances are mixed or used separately during the same occasion [[6]]. The health effects of co‐using cannabis and tobacco may be more harmful than the use of either substance alone [[7]] from the additive risks of toxicant exposure [[9]].

In Australia, the regulatory environment for tobacco and cannabis is rapidly evolving. To reduce tobacco smoking, Australia has implemented plain packaging, health warnings on tobacco products, tobacco advertising bans [[11]] and annual increments in tobacco excise [[12]]. Many individuals made smoking‐related behaviour changes in response to increasing cigarette prices by quitting smoking, reducing consumption, changing to lower‐priced brands, using loose tobacco or using e‐cigarettes [[13]]. As for cannabis, non‐therapeutic cannabis is regulated as an illicit substance federally but each state sets penalties for its use, possession and cultivation [[14]]. In 2016, the Australian federal government approved the prescription, dispensing, cultivation and manufacture of medicinal cannabis products [[15]]. The Australia Capital Territory implemented new rules in 2020 to legalise the possession of small amounts of cannabis for personal use [[16]]. These rules along with emerging cannabis products and modes of administration [[17]] may increase experimentation and the risk of addiction.

The study aims to use latent class analysis (LCA) to identify distinct groups within the Australian population based on 11 tobacco and cannabis products used. In addition, the socio‐demographic, general health, psychological distress and substance use correlate for each latent class is identified. The main benefit of the LCA approach is that it can group individuals according to the use of a wide range of tobacco and cannabis products, hence, allowing a more nuanced characterisation of how different products might be used in combination [[18]]. Since an individual's pattern of tobacco and cannabis product use can range from not using to using more than one tobacco and/or cannabis product (e.g., using cigarettes only or using joints, edibles and cigarettes, etc.), the LCA approach is preferred over the dichotomisation of (11! = 39,916,800) variables.

Studies examining the patterns of poly‐substance use using LCA have found distinct subpopulations of tobacco and cannabis users that used different products in combination [[19], [21]]. For example, Evans‐Polce and colleagues identified five unique classes of substance use behaviour including a class on 'hookah and marijuana' among college students [[19]]. Based on these studies, it is expected that there will be at least one group who co‐used both cannabis and tobacco products.

The 2019 National Drug Strategy Household Survey (NDSHS) is a multi‐stage, stratified survey that is conducted triennially to assess the attitudes and behaviours in relation to drug use in Australia. The sampling frame includes respondents aged 14 and above in residential households across Australia. A total of 22,015 people (mean age 45.3 years, 42% male) were surveyed in 2019, with a response rate of 49%. Full details of the survey are described elsewhere [[1]].

The four cannabis product variables analysed were: (i) joint; (ii) bong; (iii) edibles; and (iv) vapes. Respondents were asked if they had used cannabis in the last 12 months. If so, they were asked 'how have you used Marijuana/Cannabis?'. Options to this question were: 'Smoked as Joints', 'Smoked from a bong/pipe', 'Inhaled through a vaporising device' and 'By eating it'. Respondents could report more than one form of use. Four dichotomous variables were created for each method of use (yes vs. no for each variable).

The seven tobacco product variables analysed were: (i) manufactured cigarettes; (ii) roll‐your‐own cigarettes; (iii) e‐cigarettes; (iv) cigarillos; (v) cigars; (vi) water pipe tobacco; and (vii) other tobacco products. Respondents were asked 'how often, if at all, do you now smoke manufactured cigarettes, roll‐your‐own cigarettes, cigarillos, cigars, water pipe tobacco, or pipe tobacco?'. Response options to this question were 'Daily', 'At least weekly (but not daily)', 'Less often than weekly' or 'Not at all'. Six dichotomous variables were created for each tobacco product with any endorsement 'daily or at least weekly' or 'less often than weekly' to the frequency question coded as 'yes' to the specific tobacco product. The NDSHS also asked 'Which, if any, of the following products [chewing tobacco, snuff, snus, bidis] have you ever used, and which have you used in the last 12 months?'. Four dichotomous variables were created with endorsement to 'used in last 12 months' coded as 'yes'. Finally, the e‐cigarette question was based on 'How often, do you currently use electronic cigarette?'. Response options to this question were 'Daily', 'At least weekly (but not daily)', 'At least monthly (but not weekly)', 'Less than monthly', 'I used to use them but no longer use', 'I only tried them once or twice' or 'Not at all'. Endorsement to the first three responses was coded 'yes' to e‐cigarette use. Low prevalence (<0.5%) tobacco products (chewing tobacco, snuff, snus, bidis and pipe tobacco) were grouped as 'other tobacco product'.

The following variables were used in the analysis: sex ('male', 'female'), marital status ('never married', 'divorced, separated or widowed', 'currently married including de‐facto'), highest education attainment ('high school or less', 'certificate or diploma', 'bachelor degree or higher'), remoteness of the area of dwelling ('major cities', 'inner regional', 'outer regional, remote or very remote area'), country of birth ('Australia', 'not Australia') and current employment status ('not in labour force', 'unemployed', 'currently employed'). Household annual income were coded by dividing the income distribution into quartiles ('low income', 'low‐average', 'high‐average', 'high income').

Self‐rated general health was measured using the question 'In general, would you say your health is 'excellent', 'very good', 'good', 'fair', 'poor'?'. The Kessler psychological distress scale (K10), a 10‐item questionnaire based on respondents' anxiety and depressive symptoms in the past month, was used to assess the levels of psychological distress—('low', 'moderate', 'high or very high level') [[22]].

Alcohol use status was assessed using the Alcohol Use Disorders Identification Test Consumption short‐form (AUDIT‐C) scale [[23]]. The AUDIT‐C is an alcohol screen that identifies hazardous drinking patterns and a potential alcohol use disorder (AUD). There was a slight difference in how the questions and responses in NDSHS were worded compared to AUDIT‐C. Based on an approximation used by prior studies [[24]], a similar scoring system was derived and NDSHS respondents were classified as 'no or low risk' for AUD (total score <4 for male and <3 for female), 'high risk' for AUD (total score > = 4 for male and > =3 for female). A comparison between AUDIT‐C and NDSHS 2019 questions are provided in Table S1, Supporting Information. Respondents were also asked if they had used any of the following illicit drugs for recreational purposes in the last 12 months: ecstasy, meth/amphetamine, cocaine, hallucinogens, inhalants and opioids. Low prevalence (<2%) drugs were excluded: heroin, gamma‐hydroxybutyrate, synthetic cannabis, novel psychoactive substances, ketamine, tranquilisers, non‐medical steroid use and other opioids (methadone).

Descriptive statistics were computed to estimate the weighted prevalence of tobacco and cannabis products. LCA is a statistical method for identifying unobservable or latent groups of respondents that are similar based on a set of observed categorical variables [[26]]. LCA was used to identify patterns of tobacco and cannabis use by the 11 product type variables (joints, bong, edibles, vapes, manufactured cigarettes, roll‐your‐own cigarettes, e‐cigarettes, cigarillos, cigars, water pipe tobacco and other tobacco products). LCA was conducted using MPLUS version 8.6 [[27]]. Survey weights were applied when estimating the prevalence of each latent class. Missing values in any of the tobacco or cannabis product variables were handled using full information maximum likelihood estimation in MPlus.

We fitted a series of LCA models (2–5 classes). The optimal number of classes was determined using a range of information criteria (Akaike Information Criteria [AIC], Bayesian Information Criteria [BIC], Sample size adjusted BIC) and likelihood ratio test (Vuong‐Lo–Mendell–Rubin adjusted likelihood ratio test) and classification diagnostics (entropy and average posterior probabilities). Lower values in information criterion were preferred as it indicates a better balance between model fit and model parsimony. The Vuong‐Lo–Mendell–Rubin adjusted likelihood ratio test is also used to compare models with k classes versus k‐1 classes. A significant p‐value suggests that the model with k classes fit the data better than models with k‐1 classes [[28]]. Entropy values approaching 1 indicate a clear delineation of classes [[29]], although it was argued that entropy is a poor diagnostic to select the optimal number of class [[30]]. Average posterior probability was estimated for membership in each class, given class assignment. Values that are close to 1 suggest that there is a high degree of certainty about the true class membership for each individual. The unconditional probabilities and BCH weights for all 22,015 respondents were generated in MPlus, and this dataset was merged with the original NDSHS dataset which contains socio‐demographic, health‐rated and substance‐use correlates.

The correlates, patterns of substance use and health factors of the class membership were examined using multivariate multinomial logistic regressions in STATA17.0 [[31]]. All class membership comparisons were weighted by BCH weights. The missing values for the covariates ranged from 0.4% for marital status to 31.0% for household income (see Table S5). Prior to multinomial logistic regression, multiple imputations by chained equations (20 imputations) was performed to impute the missing values by assuming data are missing at random. A complete case analysis was conducted to evaluate the sensitivity of the multiple imputations.

All analyses were weighted to adjust for differential probabilities of selection within households and to adjust for non‐response rates and to match the samples to the population sociodemographic distribution. Significance level, α, was set at 0.0029 (α = 0.05/17) to adjust for multiple comparisons. The codes for the analysis can be found at https://github.com/clim072/patterns-of-tobacco-and-cannabis-use.

Table 1 shows the prevalence of cannabis and tobacco use by product type. In the total sample, joints (9.4%) and manufactured cigarettes (11.2%) were the most common cannabis and tobacco products, respectively. Among current cannabis users, the most prevalent cannabis products were joints (82.3%) while the most prevalent tobacco products were manufactured cigarettes (32.6%). The most prevalent cannabis and tobacco products respectively among current tobacco users were joints (28.1%) and manufactured cigarettes (78.0%).

1 TABLE Prevalence of cannabis and tobacco product use.

1 Abbreviation: CI, confidence interval.

• 2 a Weighted prevalence.
• 3 b Assessed among those with current cannabis use (use in the last year).
• 4 c Assessed among current tobacco smokers (excluding e‐cigarette use) ‐ comprised of daily and non‐daily smokers.
• 5 d Chewing tobacco, pipe tobacco, snuff, snus, bidis.

Table S2 shows the model fit statistics for a 2–5 class solution based on the latent class model on patterns of cannabis and tobacco use. The AIC, BIC, Sample Size Adjusted BIC statistics showed the lowest value for the 5‐class solution but the difference between the 4‐ and 5‐class solution is very small, suggesting both solutions fit the data well. The LRM‐LRT indicated that the 5‐class solution did not fit the data better than the 4‐class solution. Entropy and average posterior probabilities from the 4‐class solution were close to 1, indicating clear class separation. Given that the classes from the 4‐class solution were clearly interpretable and meaningful, the 4‐class solution was chosen as the most optimal.

Figure 1 shows the item‐response probabilities for each tobacco and cannabis item. Item response probabilities greater than or equal to 0.5 are considered to be high probabilities [[26]]. Class 1 was the co‐use tobacco and cannabis group (2.4%). This group was characterised by high probabilities of using joints (probability, p = 0.96), bongs (p = 0.94), edibles (p = 0.67), manufactured cigarettes (p = 0.84) and roll‐your‐own cigarettes (p = 0.63). Class 2 was the cannabis‐only group (5.5%). Respondents in this group had a high probability of using joints (p = 0.98), bongs (p = 0.86) and edibles (p = 0.67), however, very few vaped cannabis. Class 3 consists of respondents who predominantly used tobacco (8.0%). Respondents in this class had a high probability of using manufactured cigarettes (p = 0.73) and a moderate probability of using roll‐your‐own cigarettes (p = 0.46). Respondents in Class 4 had near zero probabilities of using tobacco and cannabis (84.0%).

Table 2 shows the multinomial logistic regression model using class membership as the outcome and socio‐demographic correlates, general health, psychological distress, past‐year substance use and risk of AUD as the predictors. Non‐user was used as the reference for the outcome while females, persons aged 60+, currently married, high household income, major cities, bachelor or higher education, living outside of Australia, currently employed, excellent general health, low psychological distress, did not endorse past‐year substance use, and no or low risk for AUD were specified as the reference groups for the predictors.

2 TABLE Associations between latent class memberships and socio‐demographics, health‐rated correlates and history of past‐year substance use with non‐use as reference group (Class 4).

• 6 Abbreviations: AUDIT‐C, Alcohol Use Disorders Identification Test Consumption short‐form; CI, confidence interval; OR, odds ratio.
• 7 * Significant at the 0.0029 level, two sided test.

Males, respondents aged less than 60, and those with low household income were more likely to be classified in the co‐use class (Class 1) (p < 0.05). Meanwhile, males, respondents aged less than 60, and never having married or divorced were likely to be assigned to the cannabis‐only class (Class 2) (p < 0.001). Males, respondents aged less than 60, those who were never married or divorced, not having a high household income, not living in major cities and obtained less than a bachelor degree were more likely to be in the tobacco‐only class (Class 3) (p < 0.05).

Respondents reporting fair level of general health were more likely to be in the co‐use class (Class 1) (p < 0.001). Respondents reporting all but excellent level of general health were more likely to be in the tobacco‐only class (Class 3) (p < 0.001). Respondents who reported high to very high levels of psychological distress were likely to be in co‐use and tobacco‐only class (p < 0.001 to p = 0.001) while respondents reported moderate to very high levels of psychological distress were likely to be in cannabis‐only class (p < 0.001).

Respondents that reported past‐year use of ecstasy, methamphetamine, cocaine, hallucinogens or opioids were likely to be in the co‐use class (Class 1) (p < 0.001). Meanwhile, respondents reporting past‐year ecstasy, methamphetamine, cocaine or hallucinogen use were more likely to be assigned to the cannabis‐only class (Class 2) (p < 0.001). Individuals reporting past‐year methamphetamine, cocaine or opioid use were more likely to be in the tobacco‐only class (Class 3) (p < 0.001). Finally, respondents with high risk of AUD were likely to be assigned to all three substance classes (p < 0.001). The prevalence of correlates and other factors for each latent class can be found in Table S3.

Sensitivity analyses were performed by running the multinomial logistic regression on complete cases (listwise deletion of cases) (see Table S4). While the results remained largely similar between complete case analysis and multiple imputation analysis, the results from the complete case analysis had larger standard errors and wider confidence intervals, which may result in less statistical power than multiple imputation results.

Table 3 shows the multinomial logistic regression model using class membership as the outcome. The first column of the table used cannabis class as the reference for the outcome while the second column used tobacco‐only class as the reference for the outcome. The reference groups for the predictors remained the same.

3 TABLE Associations between latent class memberships and correlates, health outcomes and history of past‐year substance use with tobacco and cannabis class as reference group.

• 8 Abbreviations: AUDIT‐C, Alcohol Use Disorders Identification Test Consumption short‐form; CI, confidence interval; OR, odds ratio.
• 9 * Significant at the 0.0029 level, two sided test.

Males, those reporting methamphetamine and opioid use in the past year were likely to be assigned to the co‐use class compared to the cannabis‐only class (Class 2) (p < 0.001 to p = 0.002). Those experiencing high to very high psychological distress, reporting cocaine, hallucinogen use or endorse high‐risk AUD were more likely to be in the co‐use class compared to tobacco‐only class (Class 3) (p < 0.001 to p = 0.002).

Table 4 shows among those in Class 1 (co‐use tobacco and cannabis), 78.4% of respondents mixed cannabis and tobacco. They were also likely to use other substances simultaneously in the past year. These included alcohol (89.4%), ecstasy (36.6%), cocaine (30.7%), hallucinogens (25.5%), methamphetamine (22.3%) and opioids (20.6%).

4 TABLE Simultaneous use of other substances among co‐use tobacco and cannabis class in the past year.

• 10 Abbreviation: CI, confidence interval.
• 11 a For non‐medical use.

This study of patterns of tobacco and cannabis use in Australia identified four distinct groups of individuals: co‐use of tobacco and cannabis class (2.4% from the sample), cannabis‐only class (5.5%), tobacco‐only class (8.0%) and a non‐user class (84.0%). We used an LCA approach to explore potential subgroups of individuals who may differ according to the tobacco and cannabis product types used, however, these major distinct subgroups were not identified. Our analysis revealed predominate product types used within each of the latent classes. Our co‐use class was characterised by high probabilities of using joints, bongs, edibles, manufactured cigarettes and roll‐your‐own cigarettes. The cannabis‐only class was characterised by high probabilities of joints, bongs and edible use while the tobacco‐only class was characterised by high probabilities of using manufactured cigarettes and moderate probability of using roll your own. Combustible forms of tobacco and cannabis were common but other routes of administration via inhalation (e.g., vaping) were less common in our study, potentially reflecting the strict regulations imposed on the sales of nicotine and vaping devices [[32]]. The shifting landscape of tobacco and cannabis warrants ongoing surveillance through national surveys, sales data and other forms of surveillance such as online cryptomarkets.

The identified classes in this study are consistent with prior literature in showing an overlap between tobacco and cannabis use [[7], [33]]. Although shared genetic liabilities are a possible explanation for the four distinct classes of tobacco and cannabis, currently there is lack of replication studies to confirm the genes and chromosomal regions related to cannabis [[35]] and tobacco use [[36]]. These genome‐wide association studies also utilised crude measurements on cannabis such as lifetime cannabis use [[37]]. The exact mechanism remains unclear as there are other explanations such as shared risk factors through peer influences, availability of products and regulatory environment.

In terms of socio‐demographics correlates, our study found males, respondents aged less than 60, and having high to very high psychological distress were associated with all tobacco and cannabis classes. The proportion of respondents in the co‐use class was small (2.4%); however, young people (aged 14–29) had an 8‐fold increased odds of being assigned to this class compared to non‐users. The odds of past year use of methamphetamine, opioids, cocaine and hallucinogens were also higher in the co‐use class than tobacco or cannabis‐only classes. Young people are susceptible to risk taking [[39]] and their substance patterns could be influenced by peer influence, stressful events and employment [[7]]. This suggests the need to enhance mental health support among these groups through community campaigns, peer‐support or workplace interventions.

In addition, the present study also showed high rates of combining alcohol (89.4%) with cannabis use among those co‐using both tobacco and cannabis. The high rates of simultaneous use of other substances among co‐users may compromise the success of quit attempts [[7]]. The effects of both drugs may be intensified when combined, and hence increase the chance of negative outcomes, such as drunk driving [[40]]. Physicians should also assess and address alcohol use while treating tobacco and cannabis addictions especially in young people since alcohol use in this age group is extremely common [[1]].

Our study found that 78% of the respondents in the co‐use class had mixed cannabis with tobacco. The 2018 International Tobacco Control Policy Evaluation Study also showed that co‐users in Australia were more likely to mix tobacco with cannabis [[4]]. Co‐use may also complicate cessation treatment when patients are dependent on both nicotine and cannabis use [[4], [41]]. McClure and colleagues found cannabis use increased by 50% during tobacco cessation quit attempts; similarly, tobacco use increased by 62% during cannabis cessation [[42]]. Evidence on dual interventions (e.g., providing nicotine replacement therapy or medication along with behavioural treatment) to reduce co‐use have been shown to reduce cannabis use but there is no clear effect on tobacco cessation [[43]]. Guidance counsellors, general practitioners, and psychologists could also help the affected individual to understand their substance use behaviour and to introduce coping strategies.

The strengths of the current study were its large sample from a population‐based survey that comprehensively measured a range of substances and relevant socio‐demographic variables and the use of LCA to simultaneously consider how tobacco and cannabis products were used in combination. There were several limitations inherent in the present work. Although the estimates can be broadly generalised to the Australian population, the NDSHS has limitations pertaining to cross‐sectional surveys. The NDSHS do not sample persons in institutional settings such as hospitals, nursing homes, and those lacking permanent addresses. Interviews were not conducted in languages other than English. The surveys also relied on the retrospective recall of respondents' experience with cannabis and tobacco. Our findings are unlikely to be generalisable to high‐risk or specific sub‐populations. The current analysis has only examined the latest NDSHS dataset, future research could consider examining changes in the demographics associated with class membership over time. Although multiple imputation were used, there are limitations because some variables have up to 31% of missing values, and there could be other variables that were not measured by the survey and so we were not able to include a comprehensive list of all possible variables that may help predict the missing values. Nevertheless, multiple imputations can reduce biases from missing data and previous comparative research has shown that it can be robust even in the presence of large percentages of missing values [[44]].

Approximately 16% of the Australian population used either tobacco, cannabis or both substances in 2019. A small proportion co‐used both substances and these were young people with higher levels of psychological distress who are more likely to simultaneously use other illicit substances. Existing policies need to minimise cannabis and tobacco related harms to reduce societal burden associated with both substances.

Carmen C. W. Lim: Conceptualisation, methodology, formal analysis, writing—original draft preparation. Gary C. K. Chan: Conceptualisation, methodology, supervision, writing—review and editing. Janni K. Y. Leung: Conceptualisation, methodology, supervision, writing—review and editing. Shannon Gravely, Coral Gartner, Daniel Stjepanović, Wayne Hall, Jack Y. C. Chung: Supervision, writing—review and editing. Tianze Sun, Vivian Chiu, Jack Y. C. Chung, Roman W. Scheurer: Writing—review and editing. Each author certifies that their contribution to this work meets the standards of the International Committee of Medical Journal Editors.

Open access publishing facilitated by The University of Queensland, as part of the Wiley ‐ The University of Queensland agreement via the Council of Australian University Librarians.

Carmen C. W. Lim is supported by a National Health Medical Research Council (NHMRC) of Australia Postgraduate Scholarship (APP2005317), The University of Queensland Living Stipend and Tuition Scholarship and a National Centre for Youth Substance Use Research top‐up scholarship. Janni K. Y. Leung is supported by an NHMRC Investigator Fellowship (APP2010008) and The University of Queensland development fellowship. Gary C. K. Chan is supported by an NHMRC Investigator Fellowship (APP1176137). The National Centre for Youth Substance Use Research is supported by Commonwealth funding from the Australian Government provided under the Drug and Alcohol Program. The funding bodies had no role in the study design, collection, analysis or interpretation of the data, writing the manuscript, or the decision to submit the paper for publication. The authors would like to acknowledge the Australian Institute of Health and Welfare for providing access to the 2019 National Drug Strategy Household Survey through the Australian Data Archive (ADA Dataverse) and the NDSHS respondents.

The authors declare no known competing financial interest or personal relationships that could influence the work reported in this paper.

This project has been reviewed by the University of Queensland Human Research Ethics Committee and deemed to be exempt from ethics review (ref: 2021/HE001006).

GRAPH: Table S1. A comparison between AUDIT‐C and 2019 NDSHS questions, and the scoring systems used to derive the risk levels of alcohol use status.Table S2. Fit statistics from latent class models on cannabis and tobacco use in the total sample (N = 22,015).Table S3. Sample characteristics by latent class membership.Table S4. Associations between latent class memberships and socio‐demographics, health rated correlates and history of past‐year substance use with non‐user as reference group (Class 4).Table S5. Number and proportion of missing values.