Skip to main content
  • Research article
  • Open access
  • Published:

Clusters of alcohol and drug use and other health-risk behaviors among Thai secondary school students: a latent class analysis



Alcohol and drug use and other health-risk behaviors tend to cluster together among adolescents and contribute a large amount of harm to both themselves and to others. This paper aims to characterize secondary school students based on their clusters of health-risk behaviors and identify factors determining class membership to these behavior-clusters.


Data from a national school survey was used to identify clusters of alcohol and drug use and other health-risk behaviors among secondary school students aged 12–15 years using a latent class regression model. A multinomial logistic regression model was used to identify predictors of the cluster membership.


A total of 25,566 students were included in the analysis, of which 88% were classified as having low-risk behaviors reporting only moderate alcohol use; 11% as having moderate-risk behaviors, such as driving under the influence of alcohol, fighting, carrying a weapon, and alcohol and tobacco use; and 0.6% as having high-risk behaviors, such as use of illicit drugs, particularly kratom and cannabis. Males, older students, those with a poor school performance, not living with parents, drug use by family members and peers, and having a low level of perceived disdain from their friends if they used drugs were significant risk factors for being in the moderate- and high-risk behavior classes.


Alcohol, tobacco and drug use, as well as other health-risk behaviors such as fighting, are clustered in Thai secondary school students. This result highlights the importance of comprehensive prevention and education strategies, particularly for moderate to high-risk groups.

Peer Review reports


Health-risk behaviors consist of many categories: behaviors that contribute to unintentional injuries and violence, sexual behaviors that contribute to unintended pregnancy and sexually transmitted diseases, alcohol, tobacco, and other drug use, unhealthy dietary behaviors and inadequate physical activities [1]. These behaviors are common among adolescents and contribute to a large amount of harm to both themselves and to others with whom they may come into contact [2,3,4].

Generally, alcohol and drug use and other health-risk behaviors tend to be clustered together [5], especially during adolescence when significant physical, sexual and psychological development occurs [3, 6,7,8]. Examples include the co-occurrence of alcohol and tobacco use, early sexual debut, and bullying [2], a clustering of high-fat diet, irregular physical activity and sedentary behaviors [9] and an inter-relationship between alcohol use, marihuana use, sexual risk behaviors and delinquency [10]. Moreover, depressive symptoms were found to occur together with risk behaviors among adolescents and young adults [6]. Such clustering has important implications for research and practice since some behaviors share a similar feature, and as a consequence, changing one behavior affects the prevalence of another [10]. Understanding how young adolescents’ behaviors cluster across individuals in an education setting could help inform how best to adapt prevention programs to more successfully serve different subgroups. Interventions that simultaneously tackled health behaviors which are clustered have been shown to be more effective and less costly [11].

Methods such as cluster analysis or latent class analysis (LCA), the latter aiming to classify people into mutually exclusive groups with which they share similar characteristics, have been commonly used. Examples include studies of health-related behaviors with relevant demographic factors and health-related outcomes [7], polysubstance use [12] and profiles of risk behaviors and mental health concerns among adolescents [13], clustering of smoking, poor nutrition, excess alcohol consumption, and physical inactivity among adults [8]. LCA identifies unobservable subgroups within a population, allowing a better understanding of the impact of an exposure to patterns of multiple risks as well as the antecedents and consequences of complex behaviors so that interventions can be tailored to target the subgroups that will benefit most [14]. Tomczyk and colleagues reported a systematic review of studies using latent class analysis to examine classes of polysubstance use in adolescents (10–19 years), and to describe predictors of polysubstance use. They found that predictors of polysubstance use classes included higher age, higher parental and peer substance use, and poor academic performance. They concluded that latent classes could deliver solid information on polysubstance use in adolescents and help to illustrate differential effects and special groups in prevention and treatment that depend on the actual substance use pattern [12].

Risk behaviors in adolescents are multi-factorial and are associated with individual factors, family characteristics, and social and environmental conditions [15, 16]. Studies have found that alcohol, tobacco and drug use are influenced by factors such as family relationships, peer groups, and school environment [3, 16,17,18]. Based on Family Interaction theory, parenting styles, familial dysfunction, parent-child relations, and parental substance use could contribute to the intergenerational transmission of substance use disorders on children [19,20,21]. Apart from the family, which provides the first context for infants and very young children, as children grow, the wider environment at school and social interaction with friends take on a larger role. The norms created by peers have increasing influence compared to family members on the adolescent’s behavioral development as the individual gets older [22]. Moreover, differences in gender, age, and the family’s socio-economic status were seen for risk-taking behaviors among adolescents, such as alcohol and cannabis consumption, bullying and fighting, early sexual debut, and smoking [2, 23]. Understanding the underlying determinants of health-risk behaviors in young people is essential if interventions are to be developed to prevent these behaviors.

As an effect of globalization, Thailand, in recent decades, has seen tremendous changes in people’s lifestyles, family and social relationships and values. Adolescents, even those aged in their early teens, are more independent of their family and as a consequence have become more engaged in health-risk behaviors, including alcohol, tobacco and drug use, premature and unsafe sex, and violent behaviors such as carrying a weapon or fighting [24,25,26]. Our previous school surveys in southern Thailand in 2002–2004 found an increasing trend of substance use among adolescents [27]. In addition, a high prevalence of risk behaviors in the 30 days prior to the survey, such as riding a motorcycle without wearing a helmet (72.7%) or after drinking alcohol (13.1%), carrying a weapon (7.6%) and fighting that involved a physical harm (14.0%) was found in our previous school survey in 2005 [26].

Most previous studies in Thailand have focused on single risk behaviors. This is the first report in Thailand aiming to characterize adolescents based on clusters of their health-risk behaviors and identify individual, family and school environmental factors determining membership of the risk behavior clusters.



The 2009 national school survey was conducted among students in levels 1, 3, and 5 of high schools, which are equivalent to years 7, 9, and 11 of the international education system, and level 2 of vocational schools. A two-stage stratified cluster sample survey was conducted. In Thailand, school administration is divided into three zones in Bangkok, and 12 other educational zones in the rest of the country. In each zone, 40–50% of the provinces were randomly selected, giving 37 provinces throughout the country and three Bangkok zones. In each province, 4–5 schools were randomly selected including one public high school in an urban area, one in a rural area, one private school, one commercial vocational school, and one technical vocational school. Totally 197 schools participated in the study. In each school, 3–5 classes of each level were randomly selected. The sampling technique and data collection method were similar to that described in the national school survey conducted in 2007 [24].

All students in the selected classes were asked to complete a self-administered, voluntary and anonymous questionnaire during their regular class period. Data collection was done between December 2009 and February 2010, the second semester of the school year in Thailand.

In this study, only students from junior high schools (years 7 and 9), representing early teenagers aged between 12 and 15 years, were included. Altogether, 25,566 students from 119 schools were included in this analysis - resulting in a 99% participation rate. About half (54%) were female and most (92%) were Buddhists.

Manifest variables

Twenty manifest variables representing four categories of health-risk behaviors that contribute to the leading causes of death and disability and common among youths [1] were included in the LCA model. They included 1) behaviors that contribute to unintentional injuries and violence in the past 12 months, 2) behaviors that contribute to intentional injuries in the past 12 months, 3) past-month alcohol use, and 3) past-year drug use.

For past-year drug use, students answered whether or not they had ever used any of the following drugs in the past 12 months: anxiolytics/hypnotics, cannabis, cough syrup (a commonly abused substance by young people in Thailand as it contains codeine), cocaine, crystalline amphetamine, ecstasy, heroin, inhalants, ketamine, kratom (Mitragynine speciosa, Kroth.) or kratom cocktail (a mixture of kratom tea, cough syrup and cola), methamphetamine, opium, and tobacco. Examples of common street names of each drug were included in the question. For anxiolytics/hypnotics and cough syrup, the question specified that the use must have been for non-medicinal purposes.

As alcohol consumption was more common than other drugs and quantity of consumption and episodes of binge drinking were considered important in classifying harmful patterns of alcohol consumption, three variables were used to measure alcohol consumption in the past 30 days of the study. These included (1) current alcohol consumption (defined as drinking at least one standard alcoholic drink, not including 1–2 sips for tasting, at least once in the past 30 days prior to the study); (2) heavy drinking (defined as drinking 2 or more standard drinks during a normal drinking occasion) and (3) frequent binge drinking (defined as having at least one episode of drinking 5 or more drinks in a drinking occasion during the past 4 weeks). Behaviors that contributed to unintentional injuries and violence in the past 12 months included driving a motorcycle or car after drinking alcohol; carrying a weapon; and fighting. Behaviors that contributed to intentional injuries in the past 12 months included suicidal risk, defined as having suicidal thoughts, suicidal plans, or suicide attempts in the past 12 months.

All manifest variables were treated as dichotomous. The questionnaire of the national school surveys in Thailand was derived from questions contained in the US Youth Risk Behavior Surveillance System [1], the U.S. Monitoring the Future Project [28], and the National Household Survey of Substance Use in Thailand [29].

Predictors of latent class membership

Individual factors

Individual factors included age, gender, religion, school class and grade point average (GPA). In Thailand, GPA is categorized into four levels of scores ranging from 0 to 4.0, with 4.0 representing the highest grade and 0 representing the lowest. Other individual factors included living with parents or other adult relatives (compared to living alone or with friends), and attitude towards drug use, measured by asking if the students thought that their friends would have a feeling of disdain towards them if they used alcohol, tobacco, and each drug. Each student’s response was coded as 1 for yes and 0 for no and these were summed giving a total score ranging from 0 to 15 where higher scores indicate a higher feeling of disdain. The median was used as a cut-off score to classify the students’ perceived feeling of disdain by their friends towards their own drug use into either high or low level. Perceived roughness of school was measured with seven questions asking if the students perceived that other students in their school had each of the following behavioral problems: physical assault, drug use, suicide, rape or sexual assault, heavy drinking, group fighting, and interpersonal quarreling. Answers to each question ranged from not at all (0), a little (1), moderate (2) and a lot of problems (3). The total score for perceived roughness of the school was derived from the sum of these scores for all seven questions and cut into high and low levels using the median.

Peer and family factors

Peer and family factors included six questions asking if any of the student’s family members or friends currently had alcohol, tobacco or other drug-related problems. Family members listed in the questionnaire included mother, father, step-mother, step-father, grandmother, grandfather, sibling and other close relative.

Statistical analysis

Percentages with standard errors (SE) of alcohol consumption, drug use, and other health risk behaviors were described using sampling weights adjusted for the two-stage stratified sampling design.

Patterns of health-risk behaviors including alcohol, tobacco, and drug use were produced using a latent class regression model. Students were considered to have significant differences in their prior probabilities of latent class membership in the different regions in Thailand [30]. Therefore, region, of which there are five (central, north, northeast, south and Bangkok and its periphery), was included in the latent class regression model as a covariate.

A series of models with increasing number of classes was constructed to determine the one that offered the best fit to the data. The model parameters included class membership probabilities (or class prevalence estimates) and class-specific symptom endorsement probabilities. A heuristic approach was employed for model assessment using the Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), sample size adjusted Bayesian Information Criterion (SSABIC) and the maximum log-likelihood value as measures to determine the model fit which best balanced the number of parameters and information from those parameters as well as the parsimony principle, i.e. fewer parameters estimated [31, 32]. Other considerations of fit were the ability to assign meaningful labels to the classes and having no individual class which contained very few students. Relative entropy was calculated based on the equation given by Gilles Celeux and G. Soromenho [33].

After fitting the models, suitable latent class membership was determined. Then, a multinomial logistic regression model was used to identify independent predictors of latent class membership. The best-fitting model was selected based on the minimum AIC. Relative risk ratios (RRR) and corresponding 95% confidence intervals compared to the baseline class of latent class membership were presented.

Data management, latent class regression, and multinomial logistic regression models were conducted using R with the “poLCA” user-contributed package [30].


Percentages of health risk behaviors

The majority (60.5%) of students (n = 15,470) had no health-risk behaviors during the past year (weighted percent = 59.2%, SE = 0.4%). Table 1 demonstrates that fighting was the most common risk behavior among Thai students accounting for 20.6% (SE = 0.3), followed by alcohol consumption (17.8%; SE = 0.3) and carrying a weapon (14.3%; SE = 0.3). Drugs were not commonly used, the highest prevalence being for tobacco (5.3%), cannabis (1.0%) and cough syrup (0.8%).

Table 1 Distribution of drinking, smoking, drug use and other health-risk behaviors among junior high school students in Thailand, 2009

Latent class regression model

Fit indices for the model on the 20 health-risk behavior variables were calculated with results shown in Table 2. All of the criteria (AIC, BIC, SSABIC, and maximum log-likelihood) suggested that the four-class model fit the data best, followed closely by the three-class model. However, the three-class model had the best balance in the number of parameters and the information from those parameters. The entropy of the three-class model was also higher than that of the four-class model (0.82 vs 0.62), indicating more certainty in the separation of classes.

Table 2 Fit indices for the latent class regression model with region as a covariate for adolescent health-risk behaviors

The three latent classes represented a continuum of levels of different health-risk behaviors. The majority (87.7%) of students were assigned to class 1, which consisted of low probabilities (< 0.2) of alcohol consumption (without binge or heavy drinking), fighting and suicidal behavior. It was labeled as a “low-risk” class with almost no health-risk or drug-related behaviors. About 11.7% of students were assigned to class 2, which was primarily characterized by higher probabilities of fighting, carrying a weapon, driving after drinking alcohol, alcohol consumption and tobacco use. This class was labeled as a “moderate-risk” class due to the students’ low frequency of drug use and high alcohol use and other health-risk behaviors. Finally, a small proportion (0.6%, n = 155) of students were allocated to class 3, which had the highest probability of drug use with alcohol consumption and other health risk behaviors. Consequently, this group was labeled as a “high-risk” class. Figure 1 illustrates the item-response probabilities of health-risk behaviors for the three-class latent membership model grouped by the four health-risk behaviour categories and then sorted by the probability of membership to class 3.

Fig. 1
figure 1

Item-response probabilities of health-risk behaviors for a three-class latent membership model. Health-risk behaviors: Violence/unintentional injuries: fig = Fighting; wea = Carrying a weapon; dui = Driving under the influence of alcohol. Intentional injuries: sui = Suicidal intentions. Alcohol consumption: alc = Drinks alcohol; hvy = Heavy drinking; bin = Binge drinking. Drug use: tob. = Tobacco; can = Cannabis; kra. = Kratom/kratom cocktail; met = Methamphetamine; cou = Cough syrup; ice = Crystalline amphetamine; opi = Opium; inh = Inhalants; ecs = Ecstasy; her = Heroin; coc = Cocaine; ket = Ketamine; anx = Anxiolytics

Differences between classes

Table 3 presents a comparison of socio-demographic characteristics and other factors between the three classes. The majority of students in class 1 were female and aged less than 14 years, studying in high school grade 1 (year 7), had higher scores of opposition to drug use and had no peers with drug use problems. In contrast, the majority of students in classes 2 and 3 were male, aged 15 years, were studying in grade 3 (year 9), had a lower level of perceived feeling of disdain by their friends towards their own drug use, and had peers with drug use problems.

Table 3 Distribution of demographic variables by latent class

Multinomial logistic regression

Table 4 shows the results of the multinomial logistic regression analysis using the low-risk class as the baseline outcome group. Males, those aged 15 years and having a lower GPA were more likely to be in the moderate-risk and high-risk behavior classes. Although Muslims were less likely to be in the moderate-class compared to Buddhists, they were more likely to be in the high-risk class. Those who perceived that their friends had less disdain for their own drug use behaviors were more likely to be in the moderate- and high-risk classes.

Table 4 Relative risk ratios (RRR) and 95% confidence intervals (CI) with low-risk behavior class as the baseline outcome group

Compared with students in the low-risk class, those whose family members had tobacco use problems were 1.3 and 1.7 times more likely to be in the moderate- and high-risk classes, respectively, while students whose family members had drug use problems were 2 and 9 times more likely to be in the moderate- and high-risk classes, respectively. Additionally, students whose peers had drug use problems were 1.3 and 3.5 times more likely to be in the moderate- and high-risk classes, respectively, while students who had peers with alcohol use problems were 1.6 and 1.5 times more likely to be in the moderate- and high-risk classes, respectively. Finally, students who had a low level of perceived disdain by their friends if they used drugs were 2.2 and 5.3 times more likely to be in the moderate- and high-risk classes, respectively, while those with a high level of perceived roughness of their own school were 1.2 and 2.0 times more likely to be so.


This study empirically supports the heterogeneity of adolescent health-risk behaviors. Our analyses generated a best-fitting three-class model of adolescent health-risk behaviors in Thai high school students. Class 1 (low-risk behaviors) was characterized by small probabilities of fighting, alcohol and tobacco use without other drug use or other health-risk behaviors while both class 2 (moderate-risk behaviors) and class 3 (high-risk behaviors) were similar in terms of having a higher probability of alcohol and tobacco use and other health-risk behaviors; they could be differentiated by the experience of drug use, which was only present in the high-risk class.

Nearly 90% of students were categorized as having a low level of health-risk behaviors. These students were more likely to be female, younger, and have a high GPA. They also had a smaller percentage of friends who had consumed alcohol, tobacco and other illicit drugs, compared to the other classes. Consistent with research in other countries, this study suggests that most adolescents could go through this transitional period successfully as the highest proportion of adolescents were in this low-risk group [3, 13, 34]. For example, Noel and colleagues found that the majority (80% of a four-class LCA model) of students in secondary schools in New Zealand were ‘healthy’ and did not engage in health-risk behaviors or suffer from emotional health concerns [13] and Ji and colleagues found the highest proportion (52% of a four-class LCA model) of middle school students in rural western China was in the low-risk class with lowest probabilities for all risk behaviors such as binge drinking, self-injurious behaviors, suicide and accidental injury [34]. One explanation of the high proportion of class 1 members in our study may be that these students were still young with about 50% being 13 years and under and studying in year 7 so they had not experienced those risk behaviors yet.

Although a small number of students in our study were in the high-risk class (155 students, 0.6%), they were engaging in multiple risk behaviors, including drug use, especially cannabis, kratom or kratom cocktail, and methamphetamine. Over 75% of these students were in year 9 and aged 14 or 15 years, which, in fact, seems to be too young for such risk behaviors. This warrants concerns of school personnel and parents of the possible harm to the students’ health and social wellbeing. Of specific to Thai adolescents was the high probability of kratom or kratom cocktail use among this group as this drug is highly available in Thailand and perceived as a “soft drug” with little or no harmful effect. The kratom cocktail (dubbed “4 × 100”, or “4 times 100”) is a mixture of four ingredients, namely boiled kratom leaves, soft drink such as coca cola, cough syrup, and either coffee or codeine, and is usually chilled with ice [35,36,37]. It is common among teenagers who drink together in a group or with friends [36, 37]. This issue requires urgent action across multiple sectors, including health, education and social services for screening and early intervention of these high-risk students.

Among the moderate-risk class students, the probability of alcohol consumption was higher than that in the high-risk class and probability of heavy drinking was similar, although frequent binge drinking was lower. Additionally, the probabilities of tobacco smoking, behaviors related to violence or unintentional injuries, and suicidal behaviors were high in this group. The clustering of alcohol consumption with these behaviors was similar to other studies [8, 13, 34] and emphasized the consequences of alcohol drinking among adolescents such as risky sexual activities [38], school absenteeism [39], and automobile accidents [38, 40]. Drinking and smoking most often start in adolescence. Therefore, strategies to prevent smoking and drinking initiation are needed to improve students’ health behaviors at an early age.

From this study, numerous factors illustrate the complexity of health risk behaviors among students. Males and those with low school performance had a higher risk of health-risk behaviors at an early age. Students who were exposed to social disadvantages, such as living without parents, or those whose families and peers consumed alcohol, smoked tobacco and used drugs, had a higher chance of engaging in health-risk behaviors. These findings are similar to another study among secondary school students where males were more likely to be in the high-risk behavior group, which included behaviors such as unhealthy diet, physical inactivity, unhealthy internet use, accidental injury, tobacco use, self-injurious behavior and suicide risk [34]. Muslims were found more likely to be in the high-risk class but less likely in the low-risk class. This may be because Muslims (80% of the ‘other’ religion in this study) do not drink alcohol, so they were less likely to be classified into the moderate-risk class, a class which was characterized by alcohol consumption risk behaviors.

After adjusting for age and gender, the strongest predictors of high-risk class were drug use problems among family members (RRR = 8.94), followed by a low perceived level of disdain by friends if a student used drugs (RRR = 5.26), and drug use by peers (RR = 3.45). These results are consistent with other studies in which having drug using family members and friends increased the risk of health-risk behaviors among adolescents [41, 42]. Similarly, other studies have found that smoking and alcohol use by friends increased the risk of alcohol use and other health-risk behaviors by the adolescents [43, 44] while perceived peer norms on drug use were influential on alcohol and cannabis use among early adolescents [45]. This suggests that adolescent health-risk behaviors can be conceptualized as a result of both peer influence and selection processes, with peer influences enhancing risk-taking among adolescents [46] and adolescents selecting friends who engage in similar levels of risk behaviors.

All of the aforementioned findings highlight the need for appropriate statistical methods, such as latent class regression, to fully understand alcohol, tobacco, and drug use and other health-risk behaviors and their determinants. Using single, broad measures of health-risk behaviors would result in young people who frequently drink heavily or occasionally use drugs and may have behavioral problems being categorized with those who have a low alcohol intake and no risk behaviors at all. As alcohol and drug use is linked to adverse consequences, our findings emphasize the need to focus on adolescents with alcohol and drug use and experiences of other health-risk behaviors, especially among those whose family members and peers have drug use problems. Thailand’s rapid economic development in the past few decades and the effects of globalization have created many changes in the society and family lives of Thais. More urbanization is evident, and people are under more pressure to compete in an increasingly rapid daily lifestyle. Children are left behind with less parental monitoring. Adolescents are also exposed to a lot of stresses in their lives including peer-pressure and competition at school. It is therefore not surprising to find that our results of clusters of health-risk behaviors and their predictors are similar to what has been found in other countries. Positive youth development programs [47] have been promoted in Thailand recently where children and adolescents are encouraged to engage in productive activities and enhance their strengths. A recent nation-wide school curriculum initiative called “Study less, know more” was launched in 2017 aiming to reduce classroom study time and increase out-of-classroom skills development activities for junior high school students ( Our results could thus serve as a foundation for information on the characteristics of adolescents and thus provide school administrators and teachers with ideas to tailor activities to promote student’s positive strengths through those activities. Our results can also serve as a baseline for future school-based surveys in order to assess the effectiveness of this government initiative.

Our study has strengths in terms of the large sample size representing Thai adolescents of this school age and high response rate. As education is compulsory in Thailand up to year 9, most adolescents of this age group are in school, and therefore our sample well represented young adolescents in the whole country. However, there are a few limitations which should be acknowledged. Firstly, cross-sectional study designs preclude establishment of direction of causation between hypothesized explanatory variables and study outcomes. Secondly, recall bias and social desirability bias may also exist given the 12-month reference periods for risk behaviors and the inquiry about some sensitive questions such as drug use. Thirdly, although four common adolescent health-risk behavior categories were included in the analysis, some were not included, for example unhealthy diet and physical inactivity. Fourthly, although the response rate was very high - the sample included non-absentees and non-dropouts, the biasing effects of school-absenteeism, which is typically highly prevalent among the high-risk-behavior adolescents could not be avoided. Finally, by dichotomizing some of the manifest variables such as the number of standard drinks typically consumed during a drinking occasion and frequency of binge drinking, detailed information may have been missed. However, this did not affect the result of the clustering of risk behaviors using the latent class analysis.


This study supports previous literature in reporting the heterogeneity of adolescents in school and, using a latent class regression analysis, expands our understanding by identifying subgroups of adolescents with multi-dimensional combinations of problem behaviors. With this sample, evidence was found of three classes of students: one having higher drug, alcohol and tobacco use and other health-risk behaviors, one with alcohol and tobacco use and other health-risk behaviors, and one almost without any of these behaviors. Our study in the Thai context of early middle school also confirms previous findings related to predictors of risky behaviors in early adolescents. For example, males, poor academic achievement, and students with families and peers who used drugs were at a higher risk. Accordingly, preventive programs targeting both families and schools are encouraged, as they are able to address all of these problems at the same time. The latent class analysis provides a modern person-centered approach to detect clusters of health-risk behaviors [12], and the results lay a basis for future research and practice in the prevention field of alcohol, drug use and health-risk behaviors among adolescents.



Akaike Information Criterion


Bayesian Information Criterion


relative risk ratio


Grade point average


Latent class analysis


Standard error


Sample size adjusted Bayesian Information Criterion


  1. Centers for Disease Control and Prevention. Youth risk behaviour surveillance system (YRBSS). 2015. Accessed 23 Sept 2018.

  2. Boricic K, Simic S, Eric JM. Demographic and socio-economic factors associated with multiple health risk behaviours among adolescents in Serbia: a cross sectional study. BMC Public Health. 2015;15:157.

    Article  Google Scholar 

  3. Kazmer L, Dzurova D, Csemy L, Spilkova J. Multiple health risk behaviour in Czech adolescents: family, school and geographic factors. Health Place. 2014;29:18–25.

    Article  Google Scholar 

  4. Sychareun V, Thomsen S, Faxelid E. Concurrent multiple health risk behaviors among adolescents in Luangnamtha province, Lao PDR. BMC Public Health. 2011;11:36.

    Article  Google Scholar 

  5. Jessor R. Risk behavior in adolescence: a psychosocial framework for understanding and action. J Adolesc Health. 1991;12(8):597–605.

    Article  CAS  Google Scholar 

  6. Bannink R, Broeren S, Heydelberg J, van’t Klooster E, Raat H. Depressive symptoms and clustering of risk behaviours among adolescents and young adults attending vocational education: a cross-sectional study. BMC Public Health. 2015;15:396.

    Article  Google Scholar 

  7. Busch V, Van Stel HF, Schrijvers AJ, de Leeuw JR. Clustering of health-related behaviors, health outcomes and demographics in Dutch adolescents: a cross-sectional study. BMC Public Health. 2013;13:1118.

    Article  Google Scholar 

  8. Noble N, Paul C, Turon H, Oldmeadow C. Which modifiable health risk behaviours are related? A systematic review of the clustering of smoking, nutrition, alcohol and physical activity ('SNAP') health risk factors. Prev Med. 2015;81:16–41.

    Article  Google Scholar 

  9. Leech RM, McNaughton SA, Timperio A. The clustering of diet, physical activity and sedentary behavior in children and adolescents: a review. Int J Behav Nutr Phys Act. 2014;11(4):4.

    Article  Google Scholar 

  10. Huang DY, Lanza HI, Murphy DA, Hser YI. Parallel development of risk behaviors in adolescence: potential pathways to co-occurrence. Int J Behav Dev. 2012;36(4):247–57.

    Article  Google Scholar 

  11. Prochaska JJ, Prochaska JO. A review of multiple health behavior change interventions for primary prevention. Am J Lifestyle Med. 2011;5(3):208–21.

    Article  Google Scholar 

  12. Tomczyk S, Isensee B, Hanewinkel R. Latent classes of polysubstance use among adolescents-a systematic review. Drug Alcohol Depend. 2016;160:12–29.

    Article  Google Scholar 

  13. Noel H, Denny S, Farrant B, Rossen F, Teevale T, Clark T, et al. Clustering of adolescent health concerns: a latent class analysis of school students in New Zealand. J Paediatr Child Health. 2013;49(11):935–41.

    Article  Google Scholar 

  14. The Methodology Center. Latent class analysis. 2015. Accessed 3 Jan 2016.

  15. de Winter AF, Visser L, Verhulst FC, Vollebergh WA, Reijneveld SA. Longitudinal patterns and predictors of multiple health risk behaviors among adolescents: the TRAILS study. Prev Med. 2016;84:76–82.

    Article  Google Scholar 

  16. Lazzeri G, Azzolini E, Pammolli A, Simi R, Meoni V, Giacchi MV. Factors associated with unhealthy behaviours and health outcomes: a cross-sectional study among Tuscan adolescents (Italy). Int J Equity Health. 2014;13(1):83.

    Article  Google Scholar 

  17. McKay MT. Parental rules, parent and peer attachment, and adolescent drinking behaviors. Subst Use Misuse. 2015;50(2):184–8.

    Article  Google Scholar 

  18. Mrug S, McCay R. Parental and peer disapproval of alcohol use and its relationship to adolescent drinking: age, gender, and racial differences. Psychol Addict Behav. 2013;27(3):604–14.

    Article  Google Scholar 

  19. Straussner SLA, Fewell CH. A review of recent literature on the impact of parental substance use disorders on children and the provision of effective services. Curr Opin Psychiatry. 2018;31(4):363–7.

    Article  Google Scholar 

  20. Mellentin AI, Brink M, Andersen L, Erlangsen A, Stenager E, Bjerregaard LB, et al. The risk of offspring developing substance use disorders when exposed to one versus two parent(s) with alcohol use disorder: a nationwide, register-based cohort study. J Psychiatr Res. 2016;80:52–8.

    Article  Google Scholar 

  21. Brook JS, Brook DW, Gordon AS, Whiteman M, Cohen P. The psychosocial etiology of adolescent drug use: a family interactional approach. Genet Soc Gen Psychol Monogr. 1990;116(2):111–267.

    CAS  PubMed  Google Scholar 

  22. Li Y-H, Mgbere O, Abughosh S, Chen H, Cuccaro P, Essien EJ. Modeling ecodevelopmental context of sexually transmitted disease/HIV risk and protective behaviors among African-American adolescents. HIV AIDS (Auckl). 2017;9:119–35.

    Google Scholar 

  23. Eaton DK, Kann L, Kinchen S, Shanklin S, Flint KH, Hawkins J, et al. Youth risk behavior surveillance - United States, 2011. MMWR Surveill Summ. 2012;61(4):1–162.

    PubMed  Google Scholar 

  24. Assanangkornchai S, Mukthong A, Intanont T. Prevalence and patterns of alcohol consumption and health-risk behaviors among high school students in Thailand. Alcohol Clin Exp Res. 2009;33(12):2037–46.

    Article  Google Scholar 

  25. Saingam D, Assanangkornchai S, Geater AF. Drinking-smoking status and health risk behaviors among high school students in Thailand. J Drug Educ. 2012;42(2):177–93.

    Article  Google Scholar 

  26. Sam-angsri N, Assanangkornchai S, Pattanasattayawong U, Muekthong A. Health-risk behaviors among high-school students in southern Thailand. J Med Assoc Thail. 2010;93(9):1075–83.

    Google Scholar 

  27. Assanangkornchai S, Pattanasattayawong U, Samangsri N, Mukthong A. Substance use among high-school students in southern Thailand: trends over 3 years (2002-2004). Drug Alcohol Depend. 2007;86(2–3):167–74.

    Article  Google Scholar 

  28. Johnston LD, O'Malley PM, Bachman JG, Schulenberg JE. Monitoring the future: National Survey Results on drug use, 1975–2015. Overview, key findings on adolescent drug use. Ann Arbor: Institute for Social Research, The University of Michigan; 2016.

    Google Scholar 

  29. Assanangkornchai S, Sam-Angsri N, Rerngpongpan S, Lertnakorn A. Patterns of alcohol consumption in the Thai population: results of the National Household Survey of 2007. Alcohol Alcohol. 2010;45(3):278–85.

    Article  Google Scholar 

  30. Linzer DA, Lewis J. poLCA: Polytomous Variable Latent Class Analysis R package version 1.1; 2007.

    Google Scholar 

  31. Lanza ST, Rhoades BL. Latent class analysis: an alternative perspective on subgroup analysis in prevention and treatment. Prev Sci. 2013;14(2):157–68.

    Article  Google Scholar 

  32. Gobel K, Scheithauer H, Braker AB, Jonkman H, Soellner R. Substance use patterns among adolescents in Europe: a latent class analysis. Subst Use Misuse. 2016;51(9):1130–8.

    Article  Google Scholar 

  33. Celeux G, Soromenho G. An Entropy criterion for assessing the number of clusters in a mixture model. RR-1874, INRIA. 1993. <inria-00074799>. 2006. Accessed 18 Apr 2017.

  34. Ji Y, Xu H, Zhang Y, Liu Q. Heterogeneity of adolescent health risk behaviors in rural western China: a latent class analysis. PLoS One. 2018;13(6):e0199286.

    Article  Google Scholar 

  35. Chittrakarn S, Penjamras P, Keawpradub N. Quantitative analysis of mitragynine, codeine, caffeine, chlorpheniramine and phenylephrine in a kratom (Mitragyna speciosa Korth.) cocktail using high-performance liquid chromatography. Forensic Sci Int. 2012;217(1–3):81–6.

    Article  CAS  Google Scholar 

  36. Saingam D, Assanangkornchai S, Geater AF, Balthip Q. Pattern and consequences of krathom (Mitragyna speciosa Korth.) use among male villagers in southern Thailand: a qualitative study. Int J Drug Policy. 2013;24(4):351–8.

    Article  Google Scholar 

  37. Tanguay P. Kratom in Thailand 2011. or doi: Accessed October 2018.

  38. Miller JW, Naimi TS, Brewer RD, Jones SE. Binge drinking and associated health risk behaviors among high school students. Pediatrics. 2007;119(1):76–85.

    Article  Google Scholar 

  39. Chou LC, Ho CY, Chen CY, Chen WJ. Truancy and illicit drug use among adolescents surveyed via street outreach. Addict Behav. 2006;31(1):149–54.

    Article  Google Scholar 

  40. Carpenter C, Dobkin C. The effect of alcohol consumption on mortality: regression discontinuity evidence from the minimum drinking age. Am Econ J Appl Econ. 2009;1(1):164–82.

    Article  Google Scholar 

  41. German M, Gonzales NA, Dumka L. Familism values as a protective factor for Mexican-origin adolescents exposed to deviant peers. J Early Adolesc. 2009;29(1):16–42.

    Article  Google Scholar 

  42. Wang C, Hipp JR, Butts CT, Jose R, Lakon CM. Alcohol use among adolescent youth: the role of friendship networks and family factors in multiple school studies. PLoS One. 2015;10(3):e0119965.

    Article  Google Scholar 

  43. Huang GC, Unger JB, Soto D, Fujimoto K, Pentz MA, Jordan-Marsh M, et al. Peer influences: the impact of online and offline friendship networks on adolescent smoking and alcohol use. J Adolesc Health. 2014;54(5):508–14.

    Article  Google Scholar 

  44. Varga S, Piko BF. Being lonely or using substances with friends? A cross-sectional study of Hungarian adolescents’ health risk behaviours. BMC Public Health. 2015;15:1107.

    Article  Google Scholar 

  45. Salvy SJ, Pedersen ER, Miles JN, Tucker JS, D’Amico EJ. Proximal and distal social influence on alcohol consumption and marijuana use among middle school adolescents. Drug Alcohol Depend. 2014;144:93–101.

    Article  Google Scholar 

  46. Cavalca E, Kong G, Liss T, Reynolds EK, Schepis TS, Lejuez CW, et al. A preliminary experimental investigation of peer influence on risk-taking among adolescent smokers and non-smokers. Drug Alcohol Depend. 2013;129(1–2):163–6.

    Article  Google Scholar 

  47. Damon W. What is positive youth development? Ann Am Acad Pol Soc Sci. 2004;591(1):13–24.

    Article  Google Scholar 

Download references


Not applicable.


This 2009 National School Survey was funded by the Centre for Alcohol Studies (CAS), Thai Health Promotion Foundation, with Sawitri Assanangkornchai being the principal investigator. The CAS had no involvement in the study design; in the collection, analysis and interpretation of data; in the writing of the report; nor in the decision to submit the paper for publication. The results presented and their interpretation are the author’s and do not represent the views of CAS or the Thai Health Promotion Foundation.

Availability of data and materials

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

Author information

Authors and Affiliations



Conception or design of the work: SA, JL. Data collection: SA, DS. Data analysis and interpretation: SA, JL, EM, DS. Drafting the article: SA, JL, EM. Critical revision of the article: SA, EM. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Edward McNeil.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee for Research in Human Subjects of the Faculty of Medicine, Prince of Songkla University. All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Students were given information about the study objectives, their rights to participate or not, and the confidentiality of their information before answering the questionnaire. A waiver of the need for a document of consent for minors was approved from the institutional review board. As anonymity was guaranteed, the research involved no more than minimal risk to subjects and therefore the waiver of consent did not affect the rights and welfare of the subjects.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

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 distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Assanangkornchai, S., Li, J., McNeil, E. et al. Clusters of alcohol and drug use and other health-risk behaviors among Thai secondary school students: a latent class analysis. BMC Public Health 18, 1272 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: