Skip to main content

Differentiating people who use cannabis heavily through latent class analysis



People who use cannabis daily or near-daily vary considerably in their daily dosage and use frequency, impacting both experienced effects and adverse consequences. This study identified heavy cannabis user groups according to consumption patterns and factors associated with class membership.


We conducted a cross-sectional study of 380 Spanish residents (61.8% male; average age = 30.3 years) who had used cannabis ≥ 3 days/week throughout the past year. Participants were recruited through chain referral and cannabis social clubs. We applied latent class analysis (LCA) to cluster participants according to use intensity. LCA indicators included frequency of weekly cannabis use, joints smoked each day, cannabis dosage, and if cannabis was consumed throughout the day or only at specific times. Associations between class membership and socio-demographics, use patterns, motives, supply sources, adverse outcomes, and use of other substances were measured using ANOVA and chi-squared tests. Multinomial regression identified the factors associated with latent class membership.


Three latent classes (moderately heavy: 21.8%, heavy: 68.2%, very heavy: 10%) had average weekly cannabis intakes of 2.4, 5.5, and 18.3 g, respectively. Very heavy users were older (\({\chi }^{2}\)=17.77, p < 0.01), less educated \({(\chi }^{2}\)=36.80, p < 0.001), and had used cannabis for longer (F = 4.62, p = 0.01). CAST scores (F = 26.51, p < 0.001) increased across the classes. The prevalence of past-month alcohol use was lower among the heaviest users (\({\chi }^{2}\)=5.95, p = 0.05). Cannabis was usually obtained from a club by very heavy users (\({\chi }^{2}\)=20.95, p < 0.001).


People who use cannabis heavily present three groups according to frequency and quantity of cannabis consumption. Use intensity is associated with increased cannabis-related problems. Differences among heavy users must be considered in harm reduction interventions in cannabis clubs and indicated prevention.


Cannabis is the most commonly used illicit drug worldwide [1], and Spain has one of the highest consumption rates in Europe [2]. In 2019, 10.5% of the population (aged 15–64 years) had used cannabis in the past year. In total, 8% had used cannabis in the past 30 days, and 2.9% used it daily [3].

Research has found that daily or near-daily (DND) users (i.e., ≥ 20 days of consumption during the previous month) account for almost 80% of total cannabis consumption [4,5,6]. Their use patterns require more detailed analysis because they consume larger doses and take it several times during the waking day [7,8,9,10]. Therefore, they are more vulnerable to acute and chronic health or psychological harms associated with cannabis use, including impaired psychomotor and cognitive functioning, memory deficits, dependence, respiratory impairments (including bronchitis), onset or amplification of psychosis in predisposed individuals, and driving impairment with risk of traffic injuries [11,12,13,14,15]. Daily use in adolescents and young adults is also associated with early school leaving, increased risk of using other illicit drugs, and cognitive and psychological deterioration [12,13,14,15].

Experts agree that standardized tools for measuring cannabis exposure are necessary to unify the evidence on the antecedents, correlates, and consequences of use [16]. Surveys tend to distinguish people who use cannabis (PWUC) heavily based on frequency of use, and rarely ask about average daily intake. However, DND users consume remarkably diverse amounts of cannabis. For instance, Gamella and Jiménez [17] found that the number of joints smoked monthly by a sample of 81 Spanish daily long-term cannabis users ranged from 10 to 300; similar findings were reported internationally [10, 18]. Therefore, it seems that a more accurate assessment of cannabis use requires assessing both, frequency and quantity. Indeed, Tomko et al. [19] have shown that the sum of the quantity of cannabis used (average grams per administration * number of joints/day) and the frequency of use in the past 30 days significantly improves the prediction of urine cannabinoid level and cannabis-related problems.

This study aimed to analyze groups of heavy users according to their use patterns, and whether their consumption and its consequences varied. We hypothesized that the “heavy use” category comprises a variety of use patterns associated with different vulnerability for cannabis-related problems. We chose latent class analysis (LCA)—a model-based clustering method—to uncover these groups in a sample of heavy users. This method allows to identify different groups, accurately characterizes membership, and predicts which case is likely to belong to each group.

Previous studies have utilized LCA to identify cannabis user types according to several sets of variables. Craft et al. [20] identified seven classes of PWUC characterized by the probability of using different products. Herbal cannabis, sinsemilla, and hashish were associated with increased dependence, whereas the consumption of concentrates was associated with diagnosis of mental health disorders. Davis et al. [21] distinguished four groups in a sample of twins and siblings based on the concurrent and simultaneous use of cannabis with tobacco, alcohol, and other illicit drugs. Simultaneous use of cannabis and tobacco was associated with the most problematic outcomes, including depression, illicit drug use, and cannabis use disorders. Two studies found four and five subgroups, respectively, of persons who use cannabis based on the products they used and their past-month use frequency/intensity [22, 23]. Heaviest users, who consumed plant products, concentrates, or both frequently and spent more time high, were more susceptible to experiment adverse consequences. Another two studies uncovered four and five classes of users, respectively, according to consumption frequency, quantity/intensity, and cannabis-related problems [24, 25]. Both studies identified three groups among heavy users with increasingly more negative consequences.

To conduct our LCA, we included only variables affecting the magnitude of exposure to cannabis and its effects. We excluded consumption method because all our participants principally smoked cannabis in joints. Smoking combusted cannabis is the predominant method among PWUC in Spain. In 2019, 98% of past-month cannabis users reported smoking it in a joint; 87% mixed it with tobacco [3]. We chose for the post-LCA analysis variables describing social, health, or behavioral characteristics or outcomes of cannabis use relevant to public health. We followed the previous literature to select these indicators: sociodemographic factors, settings and sources of cannabis use, other drug use, and health outcomes [26].

This study aims to contribute to the knowledge of PWUC heavily in Spain, their demographics, use patterns, and experienced consequences. We applied LCA to: [1] identify classes in heavy users according to their cannabis exposure magnitude; [2] examine if there were associations between class membership and other variables (socio-demographics, characteristics of cannabis use, motives, sources of supply, adverse outcomes, and concurrent use of other substances); and [3] identify which of the examined variables remained as predictors of class membership in a multinomial regression model.


Participants and procedure

We recruited 380 individuals who had consumed cannabis ≥ 3 days/week in the past 12 months (average cannabis use days/week = 6.3; 61.8% male; average age = 30.3; age range = 18–76). The Dutch Cannabis Dependence (CanDep) study employed the same eligibility criterion [27]. Participants were recruited through chain referral and cannabis social clubs in Spain from January 2017 to May 2019. Involvement was voluntary and anonymous. Participants were informed about the research aims, the institutions responsible and the funding body. They provided consent to respond to a paper-and-pencil self-administered questionnaire (average completion time: 40 min). In total, 142 participants (37%) were members of cannabis clubs. These clubs are registered non-profit associations that organize the production and distribution of cannabis among their adult members. We invited members of some cannabis clubs to participate with the help of some of the clubs’ staff. They sent us back the completed questionnaires. Other club members appeared in the networks of users found by snowball sampling. Various collaborators help us to follow these networks and to contact participants in different settings.


We selected four variables describing the intensity of cannabis use for LCA:

  • The number of weekly cannabis use-days.

  • The number of joints typically smoked in a use-day was estimated by responding to the following questions: “How many joints did you smoke the last time you smoked alone?” “How many joints did you smoke the last time you smoked with others?” “With how many individuals did you share each joint?” For example, smoking one joint alone and three with another individual gives a total of 2.5 joints. If the last consumption day was atypical, participants responded to the same questions regarding a regular cannabis use day.

  • Cannabis dosage was estimated by self-report of the number of joints made with one gram of cannabis. This method has been previously validated to better estimate average doses per joint at group level compared to other methods such as using a prompt card with real-size pictures of a ruler and different dosages [28]. The estimated median in our sample was 0.28 g per joint, very near the 0.25-gram “standard unit joint” [28].

  • All-day consumption. Finally, individuals had to report if they usually consumed cannabis throughout the waking day or only at specific times (morning, afternoon, evening, before going to bed).

The variables selected for the post-LCA analysis were:

  • Socio-demographics. Age (by cohorts; ordinal variable), sex (male dummy variable), marital status (married or with a partner versus other; categorical variable), education (primary, secondary, vocational training, baccalaureate, or university; ordinal variable), employment (employed, unemployed, or inactive; categorical variable), and monthly income (< 301€, 301–900€, > 901€; ordinal variable).

  • Characteristics of cannabis use. Number of years of cannabis use (continuous). Methods of cannabis consumption in the last month (smoked mixed with tobacco, smoked without tobacco, vaporized, eaten, or drunk; all yes or no). Cannabis use settings in the past month (at home, outdoors [street, square, park, countryside, beach], in a club, in a vehicle; all yes or no).

  • Motives for use. We used the Spanish version of the 25-item Marijuana Motives Measure (MMM) questionnaire [29,30,32]. Participants rated each item on a 5-point Likert scale from 1 (almost never/never) to 5 (almost always/always).

  • Source of supply. Purchased in the illegal market, obtained in a cannabis club, home-grown, obtained for free (all yes or no).

  • Adverse consequences. Experienced cannabis-related problems in the last year (health, psychological, at school or university, family, financial; all yes or no); had received a fine because of cannabis use in the past year (yes or no). Cannabis use disorder (CUD) and cannabis dependence were assessed using the full Spanish version of the Cannabis Abuse Screening Test (CAST) [33,34,35] which assesses the frequency of the following events in the past 12 months: “Have you smoked cannabis before mid-day?”; “Have you smoked cannabis when you were alone?”; “Have you had memory problems when you smoked cannabis?”; “Have friends or family members told you that you should reduce or stop your cannabis use?”; “Have you tried to reduce or stop your cannabis use without succeeding?”; “Have you had problems because of your cannabis use (argument, fight, accident, poor results at school, etc.)?”. All items are answered on a five-point scale from 0 (never) to 4 (very often). Total scores range from 0 to 24. The cut-off points used, according to DSM-IV criteria, were: 9 for CUD and 12 for cannabis dependence (Cronbach’s alpha: 0.75, sensitivities: 0.74 and 0.57, specificities: 0.69 and 0.84) [36].

  • Use of other substances: use of tobacco and alcohol in the last month (both yes or no); use of inhalants, cocaine, amphetamine, ecstasy, LSD, sedative-hypnotics, synthetic cannabinoids, and mushrooms in the past 12 months (all yes or no). A dummy variable was created to classify individuals that had taken one or more non-cannabis illicit substances in the past year.

Statistical analysis

We applied LCA to cluster participants according to the intensity of their cannabis consumption. LCA is a statistical procedure used to identify latent (unobserved) groups within a sample by sharing some characteristics. The latent groups are inferred from patterns of the observed variables used in the modelling. We conducted different LCAs to identify classes of PWUC heavily based on patterns of cannabis use. We selected four variables to develop the LCA: number of weekly cannabis use-days, number of joints smoked on a typical use-day, number of joints made from one gram of cannabis, and cannabis use throughout the waking day. Then, we explored several solutions, beginning with the most parsimonious (i.e., one class) and increasing the number of latent classes by one to determine the model that featured the best data fit.

For reasons of interpretability and parsimony, we did not use models where the smallest class proportion was < 5.2% [37]. The first model with this characteristic was the 4-class model. To determine the number of classes that have the best fit, we considered the Akaike information criterion (AIC), consistent AIC (cAIC), Bayesian information criterion (BIC), sample-size adjusted BIC (aBIC), approximate weight of evidence criterion (AWE). Lower values indicate better fit and parsimony [38]. The Lo–Mendell–Rubin (LMR) test was also computed to compare k-1 class and k class models: higher p values suggest that the k class model does not fit the data significantly better than a model with one less class. Finally, we estimated the entropy (ranging from 0 to 1) where values > 0.80 indicate a successful classification of individuals into classes.

After deciding the number of latent classes and using the latent class posterior probabilities, individuals were classified in their most likely class. In the next step, ANOVA (Scheffe as post-hoc test) and chi-squared (and Fisher’s exact for groups with small sample size) tests were used to compare continuous and categorical variables, respectively. In the case of multiple comparisons, Holm-Bonferroni correction was performed. Finally, we applied multinomial regression to identify the factors associated with latent class membership; this regression only included statistically significant variables from the previous step. All analyses were conducted in STATA (version 15).


Latent class model

Comparisons of the LCA model-fit statistics suggested that the three-class model provided the best fit (Table 1).

Table 1 Fit indices for one- to three-class models

Figure 1 presents the predicted mean of all variables used in the LCA analysis by class. The three latent classes show a clear division according to cannabis use intensity.

Fig. 1
figure 1

Predicted mean of variables in LCA analysis by classes

  • Class 1 comprised 21.8% of the sample and had the following characteristics. Usage days per week: M = 4.27, SD = 0.77; joints smoked per consumption day: M = 1.9, SD = 1.6; dose of cannabis per joint: M = 264 mg, SD = 94; weekly consumption: M = 2.24 g, SD = 2.02; 20.5% consumed cannabis throughout the waking day. This subgroup was labeled moderately heavy.

  • Class 2 comprised 68.2% of the sample and had the following characteristics. Usage days per week: M = 6.88, SD = 0.33; joints smoked per consumption day: M = 2.6, SD = 1.3; dose of cannabis per joint: M: 305 mg, SD = 154; weekly consumption: M = 5.49 g, SD = 4.68; 33.6% used cannabis throughout the waking day. This subgroup was labeled heavy.

  • Class 3 consisted of 10% of the sample; they were all daily users. Their other characteristics were as follows. Joints smoked per consumption day: M = 7.6, SD = 1.9; dose of cannabis per joint: M = 344 mg, SD = 161; weekly consumption: M = 18.25 g, SD = 8; 89.5% took cannabis throughout the waking day. This subgroup was labeled very heavy.

Bivariate associations


Differences among classes for socio-demographic variables are presented in Table 2. There were significant intergroup differences regarding age—from youngest to oldest: moderately heavy, heavy, and very heavy. Education levels were similar between classes 1 and 2. However, the percentage of respondents with education below university level was significantly higher in the very heavy group. Similarly, unemployment was significantly higher in the very heavy group, and there were more inactive people (mainly students) among the moderately heavy users. The proportions of being employed and monthly income were similar between the three classes. There was no difference in the proportion of being married or having a partner.

Table 2 Latent class membership by socio-demographics

Characteristics of cannabis use

Table 3 panel A shows intergroup differences regarding cannabis use characteristics. The mean number of years between first and current use was 14.2 (SD = 9.1). Class 3 reported the lengthiest period of use (range: 12.2 years [class 1] to 17.5 years [class 3]). However, there was no significant difference regarding the average age of cannabis use onset: 16.1 years (SD = 3.9) overall. The most typical form of cannabis use in the past month was smoking mixed with tobacco in all classes with significant differences (range: 86.7% [class 1] to 100% [class 3]). There were no significant intergroup differences regarding other modes of cannabis administration in the past 30 days: smoked without tobacco, vaped, eaten, or drunk. In the past month, 91.8% had used cannabis at home, 42.1% outdoors (street, square, park, countryside, beach), 39.7% in a cannabis social club, and 31.6% inside a vehicle. Intergroup differences regarding all cannabis use settings in the past month were significant. Class 1 reported the lowest proportion of current use in all the locations, and class 3 the highest, except at home where class 2 was slightly more prevalent. The prevalence of individuals who had consumed cannabis in a club in the last month increased through the classes (range: 19.3% [class 1] to 63.2% [class 3]). Almost half of class 3 had used cannabis in an automobile in the past month. The proportion decreased to 21.7% in class 1

Table 3 Latent class membership by several characteristics

Motives for use

Factor analysis of the five potential motives for cannabis use (MMM) suggested enhancement motives predominated for all classes (M = 3.8, SD 0.7), followed by coping (M = 2.7, SD = 0.9), expansion (M = 2.4, SD = 1.0), and social motives (M = 2.4, SD = 0.9). Conformity motives were rare (M = 1.3, SD = 0.4). Several significant intergroup differences emerged regarding motives for cannabis use. Coping motives significantly increased through the classes (range: 2.4 [class 1] to 2.9 [class 3], and class 2 showed higher expansion motives than the other two classes. (Table 3 panel B).

Supply sources

Table 3 panel C presents the distribution of supply sources among groups. Typically, cannabis was purchased in all classes. Obtaining cannabis for free (or from sharing with others) was more common in class 1. There were substantial significant intergroup differences in the proportion of users who had obtained cannabis from a social club (range: 19.3% [class 1] to 60.5% [class 3]). Only 9% of respondents had grown their cannabis plants. Class 2 included 11.6% cultivators, in contrast with only 2.4% and 5.3% in classes 1 and 3, respectively.

Adverse outcomes

Table 4 panel A provides the bivariate associations between the three subgroups and self-reported adverse outcomes in the past year. Psychological and social troubles were more prevalent than physical problems across subgroups. Intergroup differences emerged concerning the proportion of individuals who had experienced psychological disorders, family disputes, financial problems, or received a fine (higher among very heavy users). There were considerable significant intergroup differences regarding CAST scores. Overall, 59.2% of users reported CUD (CAST ≥ 9) and 36.6% dependence (CAST ≥ 12). The prevalence of CUD (range: 43.4% [class 1] to 89.5% [class 3]) and dependence (range: 21.7% [class 1] to 68.4% [class 3]) increased significantly through the classes.

Table 4 Latent class membership by problems and concurrent use of other substances

Use of other substances

Table 4 panel B shows intergroup differences regarding the concurrent use of other substances. Most respondents had smoked tobacco separately in the past month with significant intergroup differences (range: 62.7% [class 1] to 86.8% [class 3]). The percentage of individuals who had drunk alcohol in the last month was significantly lower in the very heavy group than in the others (range: 73.7% [class 3] to 87.9% [class 2]). More than half had used at least one non-cannabis illicit drug in the past 12 months, most commonly cocaine (range: 23.2% [class 1] to 52.6% [class 3]). Significant intergroup differences emerged regarding the consumption of inhalants, amphetamine, LDS, and synthetic cannabinoids in the last year (higher among very heavy users).

Multivariate associations

We conducted a multinomial regression model with all the variables that were statistically significant in the previous bivariate analyses (Table 5). Almost all the individuals who had consumed cannabis in a club in the past month had also acquired it there. Therefore, we excluded the variable “club as use setting” in the regression model and kept only “club as a supply source.“ In adverse outcomes, we included CAST scores instead of the cut-off points for CUD or dependence. Table 5 shows the resulting model.

Table 5 Multinomial logistic regression

Compared to the moderately heavy group, the heavy group comprised more individuals in the 21–30 years age cohort (b = 1.413, p < 0.05). It was less probable that they obtained cannabis for free (b=-1.018, p < 0.01), and more likely that they acquired cannabis in a club (b = 1.360, p < 0.001). Compared to the moderately heavy group, those in the very heavy group showed a lower proportion of university graduates (b=-3.904, p = 0.01), more years of cannabis use (b = 0.197, p < 0.01), and used it outdoors more (b = 1.371, p = 0.05). It was less probable that they obtained cannabis for free (b = -1.965, p < 0.05), and more probable that they had bought it in a club (b = 3.066, p < 0.001), and had higher CAST scores (b = 0.318, p < 0.001). Compared to the heavy group, subjects in the very heavy group exhibited a lower proportion of university graduates (b=-2.657, p < 0.01) and more years of cannabis use (b = 0.145, p < 0.05). It was more probable that they purchased cannabis in a club (b = 1.706, p < 0.01) and had higher CAST scores (b = 0.223, p < 0.01). It was less probable that they had drunk alcohol in the past month (b=-1.328, p < 0.05). In summary, age, education, number of years of use, using outdoors, obtaining for free, buying cannabis in a club, CAST scores, and alcohol use were associated with cannabis use intensity.


The present study identified three latent classes of PWUC frequently characterized by consumption intensity. Age, education, years of use, and buying cannabis in a club emerged as correlates of class membership. No significant intergroup differences were observed regarding other factors—gender, unemployment, mode of administration, motives, and other substance use (except alcohol). The prevalence of CUD and dependence increased through the classes.

At least two previous LCA-based studies have shown that PWUC are a varied population, with DND users comprising three main groups or categories. Pearson et al. [24] conducted LCA on data collected from a sample of college students who had used cannabis in the past 30 days. They identified four latent classes using three variables related to the intensity of consumption and one more about cannabis-related problems. The largest class consisted of PWUC infrequently; the other three revealed increasingly frequent use and more negative consequences. They concluded that people using cannabis a few times monthly were distinct from DND users. Manning et al. [25] discovered five latent classes among a sample of 374 cannabis-using adults. The three variables used for the LCA were cannabis use frequency, quantity, and problems experienced. Three classes reported more heavy use patterns associated with increased adverse outcomes. Our study confirms that heavy users can be classified into three groups according to cannabis exposure. Besides, it shows that cannabis-related problems increase among DND users in parallel with use intensity. As such, use frequency may not be the key category for distinguishing heavy users. The number of daily doses (joints in our research) and whether users consumed cannabis throughout the waking day resulted in clear, distinct patterns of use and consequences. The heaviest users may be intoxicated throughout the waking day. However, for some DND users, especially if they consume lower doses at specific times, intoxication may not interfere much with their daily life.

In our study, the heaviest users were older, had used cannabis use for more years, and were less educated. Older users may have had more time to develop more heavy use patterns since there were no significant intergroup differences in terms of the age of cannabis use onset. Former research has interpreted early-onset and prolonged cannabis use as predictors of poorer educational outcomes and unemployment [39,40,41,42]. We only observed intergroup differences regarding education. Future studies should examine correlates between the three groups and other demographics such as gender or employment in larger samples of heavy users. For instance, in our research, 79% of the heaviest users were male, but the intergroup gender differences were not significant.

Many studies have reported associations between intensity (frequency and/or quantity) of cannabis use and cannabis-related problems [10, 24, 25, 43,44,45]. Our participants did not widely acknowledge cannabis-related problems. A third reported psychological and financial difficulties, and only 10% health damage. In contrast, the average CAST score was high (M = 10, SD = 4.4), and notable intergroup differences emerged. The proportion of individuals with scores ≥ 9 (indicating CUD) in the very heavy group was more than two times higher than in the moderately heavy group. This difference increased to more than three times when considering scores ≥ 12 (denoting dependence). These findings are consistent with former research, indicating that frequency of consumption is the most significant predictor of CUD, even when controlling for different products and modes of use [22, 46]. Additionally, nearly one in three PWUC frequently develops dependence [47]. Although using multiple products and modes of administration is not common in Spain, some users are under the effects of cannabis most of the day. Consumption intensity is a complex category that must reflect the number of hours under the influence of the substance.

Consistent with previous research [24, 25, 48, 49], enhancement, coping, expansion, and social motives were the most prevalent cannabis use reasons across classes. We did not observe significant intergroup differences regarding reasons for use in the multinomial logistic regression. These findings are consistent with those of Pearson et al. [24].

Concurrent use of other substances, except opioids, was common in our sample. Tobacco and alcohol use in the past month was prevalent across classes, and more than half reported having used at least one illicit drug other than cannabis in the past year (most commonly cocaine, followed by ecstasy and amphetamines). The prevalence of alcohol use in the past month was significantly lower in the very heavy group. Other studies have observed this polydrug use in daily users [50]. According to the literature, tobacco and cannabis seem to be complementary [51]. It is less clear if alcohol is a substitute for cannabis [52]. More research is needed to clarify the relationships between cannabis and other substances.

Spanish cannabis clubs may have become a principal supply source for the heaviest users in Spain. We found significant and broad intergroup differences regarding accessing cannabis clubs to obtain and use cannabis on their premises. Two previous studies of Spanish cannabis clubs [53, 54] found 77% and 68% of their members were daily users, respectively (samples N = 458 and N = 155). Most members were long-term cannabis users, and they did not change their use pattern after joining the club.

Some authors have pointed out that cannabis clubs could play a relevant role in implementing harm reduction practices [55, 56]. The preference of the heaviest users for this source of supply might support this proposal. However, in general, cannabis clubs have to fill some crucial gaps to implement a harm reduction policy: providing information on risks and harms, offering health support services for members, performing lab tests on the cannabis they supply, etc. [57]. Additionally, cannabis clubs must reconsider the maximum quantity of cannabis distributed monthly to each member—currently between 60 and 90 g [56]—, conduct follow-ups with frequent users, advise them to reduce their doses and frequency of use, and help problematic users access treatment and health advisory services. Based on our findings, clinical treatment interventions should also pay special attention to the use patterns of PWUC heavily since they are diverse and related to CUD. Further research will need to identify more correlates of class membership, which will enable more specific interventions for each heavy user class.


The study has several limitations. Firstly, we cannot know the representativeness of the sample. Network sampling is widely used to reduce biases in gathering samples of hard-to-reach and hidden populations, such as the one targeted in this study [58]. This study was not intended to provide prevalence estimates of different user groups that might exist among PWUC but rather to characterize DND cannabis users and identify correlates of class membership. Secondly, all data were self-reported measures, which must be considered when interpreting the results. Although self-report is an accepted method for obtaining population behavior information, individual bias and memory issues can compromise data accuracy [59]. Nonetheless, we have confidence in the validity of our main findings, which are consistent with former studies. We believe the set of questions assessing the quantity of cannabis used is reliable. We have more reservations about the results related to use of other substances. Future research will benefit from combining interview assessments, biological controls of cannabis use, and behavioral tasks to assess more accurate constructs. Thirdly, we could not measure THC content and other cannabinoids in the cannabis herb and resin used by participants, which are paramount to assessing the intensity of consumption and its consequences. However, a previous study [29] did not find large average differences in the potencies of these products in Spain. Nonetheless, future research will need to identify the potency of cannabis products to have better control of the study variables.


This study suggests that PWUC heavily form three well-differentiated classes. Class membership was related to outcomes associated with cannabis use, including increased CUD and dependence. These findings are coherent with former research and highlight the necessity of considering the differences among heavy users to implement harm reduction policies, particularly in cannabis clubs, and clinical treatment of CUD.

Data Availability

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.


  1. UNODC (United Nations Office on Drugs and Crime). World Drug Report 2021. [Internet]. (United Nations publication, Sales No. E.21.XI.8); 2021. Available from:

  2. EMCDDA (European Monitoring Centre for Drugs and Drug Addiction). European drug report 2021. Trends and developments [Internet]. Luxembourg: Publications Office of the European Union. ; 2021. Available from:

  3. OEDA (Observatorio Español de las Drogas y las Adicciones). Informe 2021. Alcohol, tabaco y drogas ilegales en España [Report 2021. Alcohol, tobacco and illegal drugs in Spain] [Internet]. Madrid: Ministerio de Sanidad. Delegación del Gobierno para el Plan Nacional sobre Drogas; 2021. Available from:

  4. Callaghan RC, Sanches M, Benny C, Stockwell T, Sherk A, Kish SJ. Who consumes most of the cannabis in Canada? Profiles of cannabis consumption by quantity. Drug Alcohol Depend [Internet]. 2019;205:107587. Available from:

  5. Caulkins JP, Pardo B, Kilmer B. Intensity of cannabis use: Findings from three online surveys. Int J Drug Policy [Internet]. 2020;79:102740. Available from:

  6. Chan GCK, Hall W. Estimation of the proportion of population cannabis consumption in Australia that is accounted for by daily users using Monte Carlo Simulation. Addiction. 2020;115:1182–6.

    Article  PubMed  Google Scholar 

  7. Burns RM, Caulkins JP, Everingham SS, Kilmer B. Statistics on cannabis users skew perceptions of cannabis use. Front Psychiatry. 2013;4(138):1–10.

    Google Scholar 

  8. Midgette G, Davenport S, Caulkins JP, Kilmer B. What America’s users spend on illegal drugs, 2006–2016 [Internet]. Santa Monica, California: RAND Corporation; 2019. Available from:

  9. Van Laar M, Frijins T, Trautmann F, Lombi L, Cannabis. market: User types, availability and consumption estimates. In: Trautmann F, Kilmer B, Thurnbull P, editors. Further insights into aspects of the illicit EU drugs market [Internet]. Luxembourg: Publications Office of the European Union; 2013. p. 73–182. Available from:

  10. Zeisser C, Thompson K, Stockwell T, Duff C, Chow C, Vallance K, et al. A ‘standard joint’? The role of quantity in predicting cannabis-related problems. Addict Res Theory. 2012;20(1):82–92.

    Article  Google Scholar 

  11. Curran HV, Freeman TP, Mokrysz C, Lewis DA, Morgan C, Parsons JA. Keep off the grass? Cannabis, cognition and addiction. Nat Rev Neurosci. 2016;17(5):293–306.

    Article  CAS  PubMed  Google Scholar 

  12. Hall WD. What has research over the past two decades revealed about the adverse health effects of recreational cannabis use? Addiction. 2015;110(1):19–35.

    Article  PubMed  Google Scholar 

  13. NASEM (National Academies of Sciences, Engineering and Medicine). The health effects of cannabis and cannabinoids: The current state of evidence and recommendations for research [Internet]. Washington, D.C.: The National Academies Press. ; 2017. Available from:

  14. Volkow ND, Swanson JM, Evins AE, DeLisi LE, Meier MH, Gonzalez R, et al. Effects of cannabis use on human behavior, including cognition, motivation, and psychosis: a review. JAMA Psychiatry. 2016;73(3):292–7.

    Article  PubMed  Google Scholar 

  15. WHO (World Health Organization). The health and social effects of nonmedical cannabis use [Internet]. 2016. Available from:

  16. Lorenzetti V, Hindocha C, Petrilli K, Griffiths P, Brown J, Castillo-Carniglia Á, et al. The International Cannabis Toolkit (iCannToolkit): a multidisciplinary expert consensus on minimum standards for measuring cannabis use. Addiction. 2022;117(6):1505–823.

    Article  Google Scholar 

  17. Gamella JF, Jiménez Rodrigo ML. El consumo prolongado de cannabis. Pautas, tendencias y consecuencias [Long-term cannabis use: patterns, trends, and consequences]. Madrid: Fundación de Ayuda Contra la Drogadicción/Junta de Andalucía, Comisionado para la Droga; 2003.

    Google Scholar 

  18. Van der Pol P, Liebregts N, Brunt T, Van Amsterdam J, De Graaf R, Korf D, et al. Cross-sectional and prospective relation of cannabis potency, dosing and smoking behaviour with cannabis dependence: an ecological study. Addiction. 2014;109(7):1101–9.

    Article  PubMed  Google Scholar 

  19. Tomko RL, Baker NL, McClure EA, Sonne SC, McRae-Clark AL, Sherman BJ, et al. Incremental validity of estimated cannabis grams as a predictor of problems and cannabinoid biomarkers: evidence from a clinical trial. Drug Alcohol Depend. 2018;182:1–7.

    Article  CAS  PubMed  Google Scholar 

  20. Craft S, Winstock A, Ferris J, Mackie C, Lynskey MT, Freeman TP. Characterising heterogeneity in the use of different cannabis products: latent class analysis with 55 000 people who use cannabis and associations with severity of cannabis dependence. Psychol Med. 2020;50(14):2364–73.

    Article  PubMed  Google Scholar 

  21. Davis CN, Slutske WS, Martin NG, Agrawal A, Lynskey MT. Identifying subtypes of cannabis users based on simultaneous polysubstance use. Drug Alcohol Depend [Internet]. 2019;205(April):107696. Available from:

  22. Gunn RL, Aston ER, Sokolovsky AW, White HR, Jackson KM. Complex cannabis use patterns: Associations with cannabis consequences and cannabis use disorder symptomatology. Addict Behav [Internet]. 2020;105:106329. Available from:

  23. Krauss MJ, Rajbhandari B, Sowles SJ, Spitznagel EL, Cavazos-Rehg P. A latent class analysis of poly-marijuana use among young adults. Addict Behav. 2017;75:159–65.

    Article  PubMed  Google Scholar 

  24. Pearson MR, Bravo AJ, Conner BT, Anthenien AM, Bravo AJ, Conner BT, et al. Distinguishing subpopulations of marijuana users with latent profile analysis. Drug Alcohol Depend. 2017;172:1–8.

    Article  PubMed  Google Scholar 

  25. Manning K, Garey L, Paulus DJ, Buckner JD, Hogan JBD, Schmidt NB, et al. Typology of cannabis use among adults: a latent class approach to risk and protective factors. Addict Behav. 2019;92:6–13.

    Article  PubMed  Google Scholar 

  26. Fischer B, Rehm J, Irving H, Ialomiteanu A, Fallu J-S, Patra J. Typologies of cannabis users and associated characteristics relevant for public health: a latent class analysis of data from a nationally representative canadian adult survey. Int J Methods Psychiatr Res. 2010;19(2):110–24.

    PubMed  PubMed Central  Google Scholar 

  27. Van der Pol P, Liebregts N, De Graaf R, Korf D, Van den Brink W, Van Laar M. The dutch Cannabis Dependence (CanDep) study on the course of frequent cannabis use and dependence: objectives, methods and sample characteristics. Int J Methods Psychiatr Res. 2011;20(3):169–81.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Van der Pol P, Liebregts N, De Graaf R, Korf DJ, Van den Brink W, Van Laar M. Validation of self-reported cannabis dose and potency: an ecological study. Addiction. 2013;108(10):1801–8.

    Article  PubMed  Google Scholar 

  29. Casajuana Kögel C, Balcells-Olivero MM, López-Pelayo H, Miquel L, Teixidó L, Colom J, et al. The standard joint unit. Drug Alcohol Depend. 2017;176:109–16.

    Article  PubMed  Google Scholar 

  30. Matali Costa J, Simons J, Pardo M, Lleras M, Pérez A, Andión O. Spanish version validation of the Marihuana Motives measure in a drug-consuming adolescent sample. Adicciones. 2018;30(4):282–91.

    Article  PubMed  Google Scholar 

  31. Simons J, Correia CJ, Carey KB, Borsari BE. Validating a five-factor marijuana motives measure: relations with use, problems, and alcohol motives. J Couns Psychol. 1998;45(3):265–73.

    Article  Google Scholar 

  32. Benschop A, Liebregts N, Van der Pol P, Schaap R, Buisman R, Van Laar M, et al. Reliability and validity of the Marijuana Motives measure among young adult frequent cannabis users and associations with cannabis dependence. Addict Behav. 2015;40:91–5.

    Article  PubMed  Google Scholar 

  33. Legleye S, Karila L, Beck F, Reynaud M. Validation of the CAST, a general population Cannabis abuse screening test. J Subst Use. 2007;12(4):233–42.

    Article  Google Scholar 

  34. Legleye S, Guignard R, Richard J, Ludwig K, Pabst A, Beck F. Properties of the Cannabis abuse screening test (CAST) in the general population. Int J Methods Psychiatr Res. 2015;24(2):170–83.

    Article  PubMed  PubMed Central  Google Scholar 

  35. DGPNSD (Delegación del Gobierno para el Plan Nacional sobre Drogas). Consumo problemático de cannabis en estudiantes españoles de 14–18 años: validación de escalas [Problematic cannabis use in Spanish students aged 14–18: validation of scales] [Internet]. 2009. Available from:

  36. Cuenca-Royo AM, Sánchez-Niubó A, Forero CG, Torrens M, Suelves JM, Domingo-Salvany A. Psychometric properties of the CAST and SDS scales in young adult cannabis users. Addict Behav. 2012;37(6):709–15.

    Article  PubMed  Google Scholar 

  37. Weller BE, Bowen NK, Faubert SJ. Latent class analysis: a guide to best practice. J Black Psychol. 2020;46(4):287–311.

    Article  Google Scholar 

  38. Masyn KE. Latent class analysis and finite mixture modeling. In: Little TD, editor. The Oxford handbook of quantitative methods in psychology. New York, NY: Oxford University Press; 2013. pp. 551–611.

    Google Scholar 

  39. Boden JM, Dhakal B, Foulds JA, Horwood LJ. Life-course trajectories of cannabis use: a latent class analysis of a New Zealand birth cohort. Addiction. 2020;115(2):279–90.

    Article  Google Scholar 

  40. Thompson K, Leadbeater B, Ames M, Merrin GJ. Associations between marijuana use trajectories and educational and occupational success in young adulthood. Prev Sci. 2019;20(2):257–69.

    Article  PubMed  Google Scholar 

  41. Zhang C, Brook JS, Leukefeld CG, Brook DW. Trajectories of marijuana use from adolescence to adulthood as predictors of unemployment status in the early forties. Am J Addict. 2016;25(3):203–9.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Beverly HK, Castro Y, Opara I. Age of first marijuana use and its impact on education attainment and employment status. J Drug Issues. 2019;49(2):228–37.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Walden N, Earleywine M. How high: quantity as a predictor of cannabis-related problems. Harm Reduct J. 2008;5:20.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Norberg MM, Mackenzie J, Copeland J. Quantifying cannabis use with the Timeline Followback approach: a psychometric evaluation. Drug Alcohol Depend. 2012;121(3):247–52.

    Article  PubMed  Google Scholar 

  45. Foster KT, Arterberry BJ, Iacono WG, Mcgue M, Hicks BM. Psychosocial functioning among regular cannabis users with and without cannabis use disorder. Psychol Med. 2018;48:1853–61.

    Article  PubMed  Google Scholar 

  46. Swan C, Ferro MA, Thompson K. Does how you use matter? The link between mode of use and cannabis-related risk. Addict Behav [Internet]. 2021;112:106620. Available from:

  47. Van der Pol P, Liebregts N, De Graaf R, Korf DJ, Van den Brink W, Van Laar M. Predicting the transition from frequent cannabis use to cannabis dependence: a three-year prospective study. Drug Alcohol Depend. 2013;133(2):352–9.

    Article  Google Scholar 

  48. Casajuana Kögel C, López-Pelayo H, Oliveras C, Colom J, Gual A, Balcells-Oliveró MM. The relationship between motivations for cannabis consumption and problematic use. Adicciones. 2021;33(1):31–42.

    Article  Google Scholar 

  49. Buckner JD, Zvolensky MJ, Crosby RD, Wonderlich SA, Ecker AH, Richter A. Antecedents and consequences of cannabis use among racially diverse cannabis users: an analysis from ecological momentary assessment. Drug Alcohol Depend. 2015;147:20–5.

    Article  PubMed  Google Scholar 

  50. Hughes JR, Fingar JR, Budney AJ, Naud S, Helzer JE, Callas PW. Marijuana use and intoxication among daily users: an intensive longitudinal study. Addict Behav. 2014;39(10):1464–70.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Agrawal A, Budney AJ, Lynskey MT. The co-occurring use and misuse of cannabis and tobacco: a review. Addiction. 2012;107:1221–33.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Risso C, Boniface S, Subbaraman MS, Englund A. Does cannabis complement or substitute alcohol consumption? A systematic review of human and animal studies. J Psychopharmacol. 2020;34(9):938–54.

    Article  PubMed  Google Scholar 

  53. Parés-Franquero Ò, Jubert-Cortiella X, Olivares-Gálvez S, Díaz-Castellano A, Jiménez-Garrido DF, Bouso JC. Use and habits of the protagonists of the story: Cannabis Social Clubs in Barcelona. J Drug Issues. 2019;49(4):607–24.

    Article  Google Scholar 

  54. Arnoso A, Elgorriaga E. Estudio de las pautas de consumo de cannabis en los clubes sociales de cannabis y evaluación de su eficacia (2016) [Study of cannabis use patterns in cannabis social clubs and evaluation of their efficacy (2016)] [Internet]. San Sebastián: Fundación Renovatio; 2016. Available from:

  55. Belackova V, Tomkova A, Zabransky T. Qualitative research in Spanish Cannabis Social Clubs: “The moment you enter the door, you are minimising the risks. Int J Drug Policy. 2016;34:49–57.

    Article  PubMed  Google Scholar 

  56. Decorte T, Pardal M, Queirolo R, Boidi MF, Sánchez Avilés C, Parés Franquero Ò. Regulating Cannabis Social Clubs: a comparative analysis of legal and self-regulatory practices in Spain, Belgium and Uruguay. Int J Drug Policy. 2017;43:44–56.

    Article  PubMed  Google Scholar 

  57. Obradors-Pineda A, Bouso JC, Parés-Franquero Ò, Romaní JO. Harm reduction and cannabis social clubs: Exploring their true potential. Int J Drug Policy [Internet]. 2021;97:103358. Available from:

  58. Bell DC, Erbaugh EB, Serrano T, Dayton-Shotts CA, Montoya ID. A comparison of network sampling designs for a hidden population of drug users: Random walk vs. respondent-driven sampling. Soc Sci Res. 2017;62:350–61.

    Article  PubMed  Google Scholar 

  59. Rendon A, Livingston M, Suzuki S, Hill W, Walters S. What’s the agreement between self-reported and biochemical verification of drug use? A look at permanent supportive housing residents. Addict Behav. 2017;70:90–6.

    Article  PubMed  Google Scholar 

Download references


We are grateful to all the people who collaborated in conducting the interviews, particularly Verónica Reyero, Sofía España, Agustín Llorca, Alejandro Sáez, Ana Esteban, Joseba Zabala, Otger Amatller, Enriqueta Ochoa, Loreto Belda, Alfons Soler, Begoña Reyero, Cristina Pinet and other members and collaborators of the EVICT project.


This research was supported by the Research Project # 2017 | 52, “Heavy and frequent cannabis use in Spain” (CIRCE), funded by the Government Delegation for National Plan on Drugs, of the Ministry of Health, Social Services and Equality, Spain.

Author information

Authors and Affiliations



A.A-R. organized the research study, conceptualized the research questions, supervised data collection, and drafted the introduction, method, results, and discussion sections. T.G.M. conducted the analyses, prepared the method and results sections, and created the tables and figures. J.F.G. and I.P. assisted with the literature review, drafted the introduction, results, and discussion sections. M.J.D. edited the manuscript. All authors reviewed the manuscript.

Corresponding author

Correspondence to Arturo Alvarez-Roldan.

Ethics declarations

Ethics approval and consent to participate

This study was conducted following the Declaration of Helsinki. Given the nature of the study (no risks of harm, no hazards or discomforts), no specific approval from the ethics committee was needed. According to prevailing Spanish legislation (Organic Law 3/2018), all participants involved in the study received information about the study objectives, the institutions responsible, and the funding body. All participants voluntarily agreed to participate in the study, and informed consent was obtained from all subjects.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Alvarez-Roldan, A., García-Muñoz, T., Gamella, J.F. et al. Differentiating people who use cannabis heavily through latent class analysis. Subst Abuse Treat Prev Policy 18, 31 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Cannabis
  • Marijuana
  • Heavy users
  • Latent class analysis
  • Patterns of use
  • Cannabis clubs