- 1School of Kinesiology and Health, Capital University of Physical Education and Sports, Beijing, China
- 2The Center of Neuroscience and Sports, Capital University of Physical Education and Sports, Beijing, China
- 3School of Physical Education, Capital University of Physical Education and Sports, Beijing, China
- 4Affiliated Competitive Sports School, Capital University of Physical Education and Sports, Beijing, China
Introduction: Exploring the factors and internal mechanisms influencing school bullying among junior high school students is crucial for preventing and controlling its occurrence. Such research provides a reference for supporting students’ physical and mental health development and promoting the high-quality advancement of education.
Methods: This study based on data from the China Education Panel Survey (2014–2015), this study employed hierarchical linear modeling to examine the individual- and school-level factors associated with bullying in junior high schools. The analysis focused on identifying predictors at both levels and investigating their interactive mechanisms.
Results: At the individual level, gender, academic achievement, health status, disciplinary behaviors, and parental discipline were significant predictors of school bullying. At the school level, school type significantly predicted bullying prevalence. Furthermore, moderation analysis indicated that school type moderated the relationship between parental discipline and bullying. Specifically, the negative association between parental discipline and bullying was stronger in private schools compared to public schools.
Discussion: This study reveals both the direct effects and the interactive mechanisms of factors influencing school bullying among junior high school students. The findings offer valuable insights for developing practical and effective prevention and intervention strategies, particularly by considering the moderating role of school type.
1 Introduction
School bullying is defined as the repeated and prolonged engagement in harmful behaviors by an individual student or a group of students targeting a specific person or group within a school setting (Olweus, 2013). Bullying at school has been extensively examined for its long-term consequences on individuals’ mental and physical health, and several studies have shown that it can be harmful (Hughes et al., 2025). According to Hymel and Swearer, bullying in school settings involves more than just physical assaults; it also includes verbal abuse, social marginalization, and the distinct features of cyberbullying (Hymel and Swearer, 2015; Romera et al., 2022). Technological advancements have been identified as a contributing factor in the increasing prevalence of cyberbullying (Jones et al., 2013). According to the World Health Organization, one in three children worldwide has experienced bullying by peers (Singham et al., 2017). Bullying in educational institutions is a serious problem in China as well. According to research, between 2 and 66% of mainland Chinese adolescents, or the bulk of the adolescent population, experience bullying at school (Han et al., 2017). School bullying demonstrated to be closely associated with a number of factors affecting adolescents, including, but not limited to, academic performance, cognitive development, physical health, mental well-being, and social adaptation (Arseneault et al., 2010; Lenzi et al., 2015; O’Brennan and Furlong, 2010). Research indicates that, in comparison with adolescents who have not been subjected to bullying, those who experience elevated levels of bullying demonstrate diminished social adaptation and academic performance (Arslan, 2022). Furthermore, these individuals are more prone to the development of feelings of inferiority, anxiety, depression, somatization, and behavioral issues. In severe cases, this may result in suicidal behavior (Alikasifoglu et al., 2007; Arslan, 2022; Moore et al., 2017), significantly impacting adolescents’ growth and development.
Individual factors in adolescents are closely associated with the occurrence of school bullying. Among these, gender is a significant factor influencing bullying behavior. The majority of studies indicate that boys experience bullying at a higher frequency than girls (Carrera Fernandez et al., 2013; Feijóo and Rodríguez-Fernández, 2021). Furthermore, there is heterogeneity in the types of bullying experienced by the two groups, boys are more likely to be subjected to physical bullying, while girls are more likely to experience verbal or relational bullying (Carrera Fernandez et al., 2013; García-Fernández et al., 2022). Body size is also an important factor influencing school bullying, with obese children being more likely to experience verbal and relational bullying than children of normal weight (Gong et al., 2020; Morales et al., 2019). Additionally, family economic conditions may also have a certain influence on the occurrence of school bullying. Children from economically advantaged families with higher social status are less likely to become victims of bullying (Elgar et al., 2013). However, parental emotional support can reduce the risk of children from economically disadvantaged families becoming victims of bullying (Jiang and Liang, 2025). Furthermore, parenting styles are also associated with children’s bullying behaviors, harsh parental discipline positively predicted school bullying among Chinese adolescents (Fan et al., 2023; Kokkinos and Panayiotou, 2007). Therefore, to gain a clearer understanding of the risks associated with school bullying among adolescents, it is essential to consider individual-level factors.
Schools are the primary settings where bullying occurs, with most incidents taking place in classrooms, playgrounds, hallways, and restrooms (Francis et al., 2022). Bowes argues that school size is a significant predictor of the occurrence of bullying, as larger schools may underestimate the risk of children being bullied due to teachers’ inability to fully understand each student’s situation (Bowes et al., 2009). School atmosphere also has is also correlated with the occurrence of bullying. A positive school atmosphere reduces the likelihood of bullying, while a negative school atmosphere weakens students’ emotional connection to the school, increasing the likelihood of bullying (Fink et al., 2018). Additionally, the type of school is an important factor influencing bullying. School type is generally determined by its sponsoring institution. Public schools are funded by government departments and tend to have more comprehensive resources and management systems, while non-public schools are funded by individuals or groups and have relatively limited recognition. Research indicates that approximately 7% of the probability of bullying occurring is attributed to differences between schools (Kärnä et al., 2011). Compared to public schools, non-public schools have a higher probability of language and relational bullying occurring (Harris et al., 2019).
However, existing research still has certain limitations. First, most existing studies explore the concepts, phenomena, harms, and related legal systems of school bullying from a theoretical perspective, while relatively few empirical studies. Second, existing studies are mostly based on data from a single region or a few regions, with small sample sizes, resulting in limited representativeness of the findings. Finally, some studies only analyze the causes of school bullying from a single perspective, with few studies considering multiple perspectives comprehensively. Based on this, this study aims to utilize the China Education Panel Survey (CEPS) data, which features rigorous design, a large sample size, and broad regional coverage, to employ Hierarchical Linear Models (HLM) to explore the influencing factors of school bullying perpetration among junior high school students from both the individual and school levels, and to provide valuable references for the prevention and intervention strategies of bullying issues.
2 Methods
2.1 Data sources
The study utilized data from the CEPS (2014–2015), which is China’s first national, longitudinal, and representative large-scale tracking survey targeting secondary school students. The survey employed the Probability Proportional to Size Sampling (PPS), randomly selecting 28 counties, 112 schools, and 438 classes of students nationwide for inclusion in the sample. The initial sample size was 9,920. Missing and outlier values for core variables were removed as needed for the study, resulting in a final effective sample size of 8,183 (Figure 1).
2.2 Variable description
This study examined school bullying behavior as the dependent variable, measured through four CEPS questionnaire items assessing verbal abuse, swearing, quarreling, physical fighting, and bullying of weaker classmates (Nansel et al., 2001; Olweus, 2013), with response coding detailed in Table 1.
Table 1. Definition of the influencing variables of school bullying among junior high school students.
This study uses the factors influencing school bullying as independent variables, including individual-level and school-level variables. Individual-level variables included gender, only child, academic performance, health status, family economic conditions, disciplinary behaviors, and parental discipline. School-level variables included school type, this variable was operationalized as a categorical distinction based on formal governance and funding sources. Public schools were defined as institutions funded and regulated by local or national educational authorities, while non-public schools included private, independent, or parochial institutions funded primarily through tuition, donations, or private grants. No inherent evaluative judgments (e.g., “good” or “bad”) were attached to this classification; it served solely to capture structural differences between school systems. The definitions of each variable and the coding of questionnaire options are shown in Table 1.
2.3 Building a model
The HLM method was originally proposed by Raudenbush and Chan (1993) for the analysis of non-independent data characterized by a multi-level nested structure. In this study, individual students were nested within schools, thus the research model includes two levels: students and schools. In the context of model estimation, the Maximum Likelihood Estimation (MLE) was utilized for the purpose of comparing log-likelihood values and evaluating the overall goodness of fit of the model. The feasibility of utilizing a multilevel model was assessed by calculating intra-class correlation coefficients (ICC), while cross-level interactions were employed to examine the moderating effect of community factors on individual factors. The specific model is as follows:
The initial tier of HLM was the individual level. It is posited that the coefficients at the individual level are capable of explaining the differences in bullying behaviors attributable to student characteristics. In this model, Yij denotes the value of bullying behavior for the i-th student in the j-th school, while β0k represents the value of bullying behavior for students in the k-th school when all variables are set to 0. Xkij signifies the value of the i-th student in the k-th school across the j individual variables, and βkj denotes the estimated coefficient. Finally, εik represents the random error. The specific model is illustrated in Equation 1.
The second level of HLM was the school level. This level provided an explanation of the impact of different school-level characteristics on bullying among junior high school students. In this model, r₀ represents of the constant term at the school level; q is the number of school-level variables; Zqk is the value of the q variables at the kth school; r₀q is the estimated coefficient for Zqk; and μ0k is the random error at the school level, indicating the deviation between the mean of bullying behavior among students at the j-th school and the overall mean of bullying behavior among all students. The specific models are illustrated in Equations 2 and 3.
Finally, substitute the school-level equation into the individual-level equation to obtain the complete HLM model, where μ0k and εik are both assumed to follow a normal distribution and are independent of each other. The specific form is shown in Equation 4.
2.4 Data analysis
The present study utilized SPSS 26.0 to perform data analysis. Firstly, a descriptive statistical analysis was conducted on the selected variables. Secondly, Pearson correlation analysis was employed to explore the correlations among the core variables, thereby providing preliminary evidence for subsequent research. A multilevel regression model incorporating individual and community factors was constructed, based on the HLM analysis method. The moderation effect analysis was conducted using the PROCESS macro (Model 1), systematically examining the independent effects and interaction effects of factors at different levels on school bullying among junior high school students.
3 Results
3.1 The results of descriptive statistics and correlation
Table 2 presents the descriptive statistics and correlation analysis results for each core variable. The results indicate that there is a significant positive correlation between junior high school students’ school bullying and disciplinary violations (r = 0.553, p < 0.001) and school type (r = 0.054, p < 0.001), and a significant negative correlation between school bullying and parental discipline (r = −0.173, p < 0.001).
3.2 Results of HLM
As shown in the results of Model 1 (Table 3), the between-group variance is 0.022 (p < 0.001), indicating that there are significant differences in school-level factors in school bullying behavior. The ICC was calculated as 0.0.022/(0.022 + 0.315) = 0.065, indicating that approximately 0.065 of the variation in bullying behavior among junior high school students is attributable to school-level factors. Based on the criterion that “the ICC is greater than 0.059”(Cohen, 1988), it can be determined that HLM should be used for model construction.
Model 2 was primarily used to examine the independent effects of individual-level variables on school bullying. As shown in the results of Table 4, gender, academic performance, health status, and parental discipline all have a significant negative predictive effect on school bullying. For each additional unit of parental discipline, the incidence of school bullying decreases by 0.066 units. Being an only child and disciplinary violations both positively predicted the occurrence of school bullying. Specifically, for each additional unit of disciplinary violations, the incidence rate of school bullying increased by 0.682 units. However, family economic conditions did not have a significant predictive effect on school bullying.
Model 3 built upon Model 2 by incorporating school-level variables, testing only the independent effects of school-level variables on school bullying without interaction terms. At the school level, the results indicated that school type does not significantly predict school bullying (p > 0.05). Model 4 was the full model, which not only analyzes the independent effects of various influencing factors on school bullying but also examines cross-level interaction effects. The results indicated that there is a significant interaction between school-level variables and individual-level variables related to school bullying among junior high school students. First, the relationship between disciplinary behavior and school bullying differed across public and non-public schools. Specifically, the association was stronger in non-public school compared to public schools. Second, school type significantly moderated the relationship between parental discipline and school bullying, meaning that the negative association between parental discipline and bullying was stronger in non-public schools compared to public schools. Following the approach of Preacher et al. (2006), the simple slope method was used, with M ± 1SD as the standard, to represent high and low values of the moderation effect. The moderation effects were shown in Figures 2, 3.
4 Discussion
This study based on CEPS (2014–2015) data and utilized the HLM to investigate the direct and interactive effects of variables at both individual and school levels on the occurrence of bullying behaviors among junior high school students, aiming to explain the multifaceted factors contributing to school bullying.
4.1 Individual-level variables and bullying among junior school students
Research has found that gender significantly influences school bullying among junior high school students, with boys being considerably more likely to be involved in school bullying than girls. This finding is consistent with Evans et al. (2013). This difference may be related to the socialization differences in gender roles during adolescence, as boys are more inclined to express themselves through overt behaviors and are thus more susceptible to physical bullying such as hitting and slapping (Nansel et al., 2001). Research has found that non-only children are more likely to be bullied. The dispersion of family resources may result in reduced emotional support and material resources for individuals. Neglect in emotional, physical, and familial support affects the development of children’s self-perception (Rees, 2008), making them more prone to involvement in school bullying incidents. Studies have also revealed a significant negative correlation between academic performance and school bullying. Within the context of the Chinese education system, where academic achievement serves as a key evaluation criterion, students with poor academic performance may face risks of social exclusion, discrimination, and ridicule, thereby increasing the likelihood of being bullied (Fu et al., 2013; Shaheen et al., 2018). This finding indicates cultural differences from studies conducted in some Western countries (Jennlfer et al., 2011). In addition, health status has a significant negative effect on school bullying. Students with better health are less likely to become targets of bullying. Conversely, children who are physically weak, timid, or lack self-protection abilities are more vulnerable to being bullied (Smokowski and Kopasz, 2005). Research also shows that disciplinary behaviors significantly affect the occurrence of school bullying. This may be due to the contrast between junior high school students’ pursuit of novelty and excitement and the monotony of school life; some students skip classes or stay away from school in search of new experiences (Zottis et al., 2014), while others resort to school bullying as a means of self-expression, seeking identity recognition, or emotional release.
In addition, research has found that family economic income significantly and positively predicts school bullying. Students from families with better economic income tend to have parents with higher educational attainment, who are more willing to invest time and effort in their children’s education and development, thereby promoting children’s cognitive development (Cave et al., 2022). Furthermore, families with higher economic conditions establish protective mechanisms for their children through increased educational investment and access to quality schools, shielding them from bullying issues. Research also indicates that parental discipline can influence the occurrence of school bullying. Reasonable parental discipline and high levels of support can mitigate the negative impact of low family economic income on children bullying behaviors (Liu et al., 2025), fostering healthy personality development, stronger psychological resilience, and the ability to face life’s challenges with positivity and optimism, thus reducing the likelihood of being bullied (Garaigordobil and Machimbarrena, 2017).
4.2 School-level variables and bullying among junior high school students
Research has found that school type also significantly influences the occurrence of bullying among junior high school students, consistent with the findings of Harris’s study (Harris et al., 2019). Specifically, public schools have the lowest incidence of bullying, and as school type shifts from public to private, the probability of bullying increases (Kelly et al., 2016). In China, most public schools receive government financial support, have high-quality faculty, well-equipped facilities, and strict disciplinary management, thereby better meeting the expectations of parents and society. In contrast, private schools for migrant workers’ children often have limited educational resources, inadequate campus facilities and environments, weak faculty, and high staff turnover, which increases school bullying incidents.
4.3 The moderating effect of school-level variables on bullying in junior schools
Research has found a significant interaction effect between disciplinary behaviors and school type, indicating that school type moderates the relationship between disciplinary behaviors and school bullying. Specifically, the key public schools typically have a long history of operation, are often located in central urban areas, and enjoy high-quality educational resources and well-established management systems, with relatively strict disciplinary standards. Students attending such schools tend to come from families with higher economic income, and these characteristics effectively reduce the likelihood of school bullying (Due et al., 2009; Fu et al., 2013).
The study further indicates that school type significantly moderates the relationship between parental discipline and school bullying. Public schools, with their professional teaching staff and comprehensive support systems, can effectively mitigate the negative effects of inadequate parenting discipline. Ecological Systems Theory posits that interaction among nested environments shape child development. Parental discipline (microsystem) and school type (mesosystem) are interdependent and modulate the protective effect of discipline against bullying. Notably, the buffering effect of parental discipline is more pronounced in public schools with high-quality teachers and standardised management. This is because behavioral norms are aligned between home and school: when parents set clear boundaries, the public school’s unified anti-bullying protocols and structured classroom management reinforce these norms. This consistency eliminates ambiguity in behavioral expectations, thereby reducing the risk of bullying. Additionally, due to their high teaching standards, teachers can promptly identify students’ psychological needs through daily observations and provide necessary emotional support, thereby reducing the risk of such students being subjected to school bullying (Zeng et al., 2023).
4.4 Limitations
This study provides a novel perspective on the factors influencing school bullying among junior high school students based on HLM, but it still has some limitations. Firstly, the phenomenon of school bullying is influenced by a multitude of factors, and the indicators selected in this study are limited. Future research could expand the dimensions of indicator selection for a more comprehensive exploration of the phenomenon. Secondly, the data utilized in this study were collected during the period 2014–2015. It is recommended that future research utilize data from recent years in order to achieve a more precise understanding of the current influencing factors. Third, this study is notable for its exclusive focus on the direct effects and interaction effects of two-level factors on school bullying among junior high school students. Future research could incorporate the family level to explore the roles of three-level factors, thereby further deepening the study. Fourth, our data analysis did not account for the influence of regional hierarchical factors (such as socioeconomic characteristics or cultural norms). Future research focusing on cross-regional comparisons would benefit significantly from adopting a multilevel model framework. This approach would allow researchers to more directly examine how macro-level contextual factors shape relationship patterns at the individual level.
5 Conclusion
This study systematically examined the relevant factors influencing school bullying among junior high school students. At the individual level, gender, academic performance, health status, disciplinary behaviors, and parental discipline all had significant influence school bullying among junior high school students. At the school level, school type had a significant influences school bullying behavior among junior high school students. Furthermore, the school type moderates the relationship between student disciplinary behaviors and bullying among junior high school students, and also moderated the relationship between parental discipline and bullying among junior high school students. This study provided valuable insights into the mechanisms underlying bullying among junior high school students and offered important implications for effective prevention and control of bullying in schools.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics statement
The studies involving humans were approved by the ethics committee of Renmin university of China. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin. Written informed consent was obtained from the individual(s), and minor(s)’ legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
Author contributions
YW: Methodology, Software, Writing – original draft, Writing – review & editing, Formal analysis. XL: Data curation, Investigation, Writing – review & editing, Validation. NS: Formal analysis, Investigation, Methodology, Software, Writing – original draft. CJ: Conceptualization, Funding acquisition, Project administration, Writing – review & editing. PL: Writing – review & editing, Methodology, Supervision.
Funding
The author(s) declared that financial support was not received for this work and/or its publication.
Acknowledgments
We thank the China Education Panel Survey (CEPS) for data access and colleagues at the Capital University of Physical Education and Sports for their feedback. All errors remain our own.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that Generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Alikasifoglu, M., Erginoz, E., Ercan, O., Uysal, O., and Albayrak-Kaymak, D. (2007). Bullying behaviours and psychosocial health: results from a cross-sectional survey among high school students in Istanbul, Turkey. Eur. J. Pediatr. 166, 1253–1260. doi: 10.1007/s00431-006-0411-x,
Arseneault, L., Bowes, L., and Shakoor, S. (2010). Bullying victimization in youths and mental health problems: ‘much ado about nothing’? Psychol. Med. 40, 717–729. doi: 10.1017/S0033291709991383,
Arslan, G. (2022). School bullying and youth internalizing and externalizing behaviors: do school belonging and school achievement matter? Int. J. Ment. Heal. Addict. 20, 2460–2477. doi: 10.1007/s11469-021-00526-x
Bowes, L., Arseneault, L., Maughan, B., Taylor, A., Caspi, A., and Moffitt, T. E. (2009). School, neighborhood, and family factors are associated with children’s bullying involvement: a nationally representative longitudinal study. J. Am. Acad. Child Adolesc. Psychiatry 48, 545–553. doi: 10.1097/CHI.0b013e31819cb017,
Carrera Fernandez, M. V., Lameiras Fernandez, M., Rodriguez Castro, Y., Failde Garrido, J. M., and Calado Otero, M. (2013). Bullying in Spanish secondary schools: gender-based differences. Span. J. Psychol. 16:e21. doi: 10.1017/sjp.2013.37,
Cave, S. N., Wright, M., and von Stumm, S. (2022). Change and stability in the association of parents’ education with children’s intelligence. Intelligence 90:101597. doi: 10.1016/j.intell.2021.101597
Cohen, J. (1988). Statistical power analysis for the behavioral sciences. 2nd Edn. New York: Routledge.
Due, P., Merlo, J., Harel-Fisch, Y., Damsgaard, M. T., scient, M., Holstein, B. E., et al. (2009). Socioeconomic inequality in exposure to bullying during adolescence: a comparative, cross-sectional, multilevel study in 35 countries. Am. J. Public Health 99, 907–914. doi: 10.2105/AJPH.2008.139303,
Elgar, F. J., Pickett, K. E., Pickett, W., Craig, W., Molcho, M., Hurrelmann, K., et al. (2013). School bullying, homicide and income inequality: a cross-national pooled time series analysis. Int. J. Public Health 58, 237–245. doi: 10.1007/s00038-012-0380-y,
Evans, S. E., Steel, A. L., and DiLillo, D. (2013). Child maltreatment severity and adult trauma symptoms: does perceived social support play a buffering role? Child Abuse Negl. 37, 934–943. doi: 10.1016/j.chiabu.2013.03.005,
Fan, H., Xue, L., Xiu, J., Chen, L., and Liu, S. (2023). Harsh parental discipline and school bullying among Chinese adolescents: the role of moral disengagement and deviant peer affiliation. Child Youth Serv. Rev. 145:106767. doi: 10.1016/j.childyouth.2022.106767
Feijóo, S., and Rodríguez-Fernández, R. (2021). A meta-analytical review of gender-based school bullying in Spain. Int. J. Environ. Res. Public Health 18:Article 23. doi: 10.3390/ijerph182312687,
Fink, E., Patalay, P., Sharpe, H., and Wolpert, M. (2018). Child- and school-level predictors of children’s bullying behavior: a multilevel analysis in 648 primary schools. J. Educ. Psychol. 110, 17–26. doi: 10.1037/edu0000204
Francis, J., Strobel, N., Trapp, G., Pearce, N., Vaz, S., Christian, H., et al. (2022). How does the school built environment impact students’ bullying behaviour? A scoping review. Soc. Sci. Med. 314:115451. doi: 10.1016/j.socscimed.2022.115451,
Fu, Q., Land, K. C., and Lamb, V. L. (2013). Bullying victimization, socioeconomic status and behavioral characteristics of 12th graders in the United States, 1989 to 2009: repetitive trends and persistent risk differentials. Child Indic. Res. 6, 1–21. doi: 10.1007/s12187-012-9152-8
Garaigordobil, M., and Machimbarrena, J. M. (2017). Stress, competence, and parental educational styles in victims and aggressors of bullying and cyberbullying. Psicothema 3, 335–340. doi: 10.7334/psicothema2016.258,
García-Fernández, C. M., Moreno-Moya, M., Ortega-Ruiz, R., and Romera, E. M. (2022). Adolescent involvement in Cybergossip: influence on social adjustment, bullying and cyberbullying. Span. J. Psychol. 25:e6. doi: 10.1017/SJP.2022.3,
Gong, Z., Han, Z., Zhang, H., and Zhang, G. (2020). Weight status and school bullying experiences in urban China: the difference between boys and girls. J. Interpers. Viol. 35, 2663–2686. doi: 10.1177/0886260519880170,
Han, Z., Zhang, G., and Zhang, H. (2017). School bullying in urban China: prevalence and correlation with school climate. Int. J. Environ. Res. Public Health 14:1116. doi: 10.3390/ijerph14101116,
Harris, A. B., Bear, G. G., Chen, D., de Macedo Lisboa, C. S., and Holst, B. (2019). Perceptions of bullying victimization: differences between once-retained and multiple-retained students in public and private schools in Brazil. Child Indic. Res. 12, 1677–1696. doi: 10.1007/s12187-018-9604-x
Hughes, A., Grey, E., Haigherty, A., Shepherd, F., Gillison, F., MacArthur, G., et al. (2025). Weight-related bullying in schools: a review of school anti-bullying policies. BMC Public Health 25:2006. doi: 10.1186/s12889-025-23170-9,
Hymel, S., and Swearer, S. M. (2015). Four decades of research on school bullying: an introduction. Am. Psychol. 70, 293–299. doi: 10.1037/a0038928,
Jennlfer, G. G., Dunn, E. C., and Johnson, R. M. (2011). A multilevel investigation of the association between school context and adolescent nonphysical bullying: journal of school violence. J. Sch. Violence 10, 133–149. doi: 10.1080/15388220.2010.539165,
Jiang, T., and Liang, L.-H. (2025). Risk and protective effects of family socioeconomic status and parental emotional support in Asian secondary school students’ bullying victimization. J. Interpers. Viol. doi: 10.1177/08862605251338786,
Jones, L. M., Mitchell, K. J., and Finkelhor, D. (2013). Online harassment in context: trends from three youth internet safety surveys (2000, 2005, 2010). Psychol. Violence 3, 53–69. doi: 10.1037/a0030309
Kärnä, A., Voeten, M., Little, T. D., Poskiparta, E., Alanen, E., and Salmivalli, C. (2011). Going to scale: a nonrandomized nationwide trial of the KiVa antibullying program for grades 1–9. J. Consult. Clin. Psychol. 79, 796–805. doi: 10.1037/a0025740,
Kelly, O., Krishna, A., and Bhabha, J. (2016). Private schooling and gender justice: an empirical snapshot from Rajasthan, India’s largest state. Int. J. Educ. Dev. 46, 175–187. doi: 10.1016/j.ijedudev.2015.10.004
Kokkinos, C. M., and Panayiotou, G. (2007). Parental discipline practices and locus of control: relationship to bullying and victimization experiences of elementary school students. Soc. Psychol. Educ. 10, 281–301. doi: 10.1007/s11218-007-9021-3
Lenzi, M., Furlong, M. J., Dowdy, E., Sharkey, J., Gini, G., and Altoè, G. (2015). The quantity and variety across domains of psychological and social assets associated with school victimization. Psychol. Violence 5, 411–421. doi: 10.1037/a0039696
Liu, H., Lan, Z., Wang, Q., Huang, X., and Zhou, J. (2025). A path to relief from the intertwined effects of school bullying and loneliness: the power of social connectedness and parental support. BMC Public Health 25:1850. doi: 10.1186/s12889-025-23064-w,
Moore, S. E., Norman, R. E., Suetani, S., Thomas, H. J., Sly, P. D., and Scott, J. G. (2017). Consequences of bullying victimization in childhood and adolescence: a systematic review and meta-analysis. World J. Psychiatry 7, 60–76. doi: 10.5498/wjp.v7.i1.60,
Morales, D. X., Grineski Sara, E., and Collins, T. W. (2019). School bullying, body size, and gender: an intersectionality approach to understanding US children’s bullying victimization. Br. J. Sociol. Educ. 40, 1121–1137. doi: 10.1080/01425692.2019.1646115,
Nansel, T. R., Overpeck, M., Pilla, R. S., Ruan, W. J., Simons-Morton, B., and Scheidt, P. (2001). Bullying behaviors among US Youth Prevalence and association with psychosocial adjustment. JAMA 285, 2094–2100. doi: 10.1001/jama.285.16.2094,
O’Brennan, L. M., and Furlong, M. J. (2010). Relations between students’ perceptions of school connectedness and peer victimization. J. Sch. Violence 9, 375–391. doi: 10.1080/15388220.2010.509009
Olweus, D. (2013). School bullying: development and some important challenges. Annu. Rev. Clin. Psychol. 9, 751–780. doi: 10.1146/annurev-clinpsy-050212-185516,
Preacher, K. J., Curran, P. J., and Bauer, D. J. (2006). Computational Tools for Probing Interactions in Multiple Linear Regression, Multilevel Modeling, and Latent Curve Analysis. J. Educ. Behav. Stat, 31, 437–448. doi: 10.3102/10769986031004437
Raudenbush, S. W., and Chan, W. S. (1993). Application of a hierarchical linear model to the study of adolescent deviance in an overlapping cohort design. J. Consult. Clin. Psychol. 61, 941–951. doi: 10.1037/0022-006X.61.6.941,
Rees, C. (2008). The influence of emotional neglect on development. Paediatr. Child Health 18, 527–534. doi: 10.1016/j.paed.2008.09.003
Romera, E. M., Luque-González, R., García-Fernández, C. M., and Ortega-Ruiz, R. (2022). Competencia social y bullying: El papel de la edad y el sexo. Educ. XX1 25, 309–333. doi: 10.5944/educxx1.30461
Shaheen, A. M., Hammad, S., Haourani, E. M., and Nassar, O. S. (2018). Factors affecting Jordanian school adolescents’ experience of being bullied. J. Pediatr. Nurs. 38, e66–e71. doi: 10.1016/j.pedn.2017.09.003,
Singham, T., Viding, E., Schoeler, T., Arseneault, L., Ronald, A., Cecil, C. M., et al. (2017). Concurrent and longitudinal contribution of exposure to bullying in childhood to mental health: the role of vulnerability and resilience. JAMA Psychiatry 74, 1112–1119. doi: 10.1001/jamapsychiatry.2017.2678,
Smokowski, P. R., and Kopasz, K. H. (2005). Bullying in school: an overview of types, effects, family characteristics, and intervention strategies. Child. Sch. 27, 101–110. doi: 10.1093/cs/27.2.101
Zeng, L.-H., Hao, Y., Hong, J.-C., and Ye, J.-N. (2023). The relationship between teacher support and bullying in schools: the mediating role of emotional self-efficacy. Curr. Psychol. 42, 31853–31862. doi: 10.1007/s12144-022-04052-4
Keywords: junior high school students, school bullying, hierarchical linear models, China education panel survey data, school type
Citation: Wang Y, Li X, Song N, Jiang C and Liu P (2026) Influential factors of school bullying among junior high school students: an analysis based on hierarchical linear models. Front. Psychol. 16:1664252. doi: 10.3389/fpsyg.2025.1664252
Edited by:
Daniel H. Robinson, The University of Texas at Arlington College of Education, United StatesReviewed by:
CRISTINA Mª GARCÍA-FERNÁNDEZ, University of Cordoba, SpainDexuan Zhao, Beijing Normal University, China
Zhijun Wu, Tsinghua University, China
Copyright © 2026 Wang, Li, Song, Jiang and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Changhao Jiang, Y2hhbmdoYW9famlhbmdAMTI2LmNvbQ==; Pei Liu, bGl1cGVpQGN1cGVzLmVkdS5jbg==
Xingyu Li3