This study was conducted in Tehran, Iran, in 2020 as per the ethical guidelines approved by the Research Council of National Institutes for Medical Research Development (NIMAD). In this survey, we used data from four stages of the CASPIAN (childhood & adolescence surveillance and prevention of adult non-communicable disease) survey. The CASPIAN study is conducted by the Ministry of Health to assess the health status of children and adolescents (students) and the risk factors of non-communicable diseases between the ages of six and 19 years.
The eligibility criteria were being Iranian, not having any particular diseases, and school attendance at the time of sampling. All students and their parents completed informed consent forms for participating in the study. In this study, a multistage cluster sampling method was applied. The study population was students aged six to 19 years in 31 provinces. The primary sampling unit in all studies were the provinces of Iran that were randomly selected (
18). In each city, students were randomly selected for the study, based on their educational level. After the final sample size estimation, approximately 54,550 students were included in the analysis. Students were from both sexes, urban and rural areas, and all school levels. The ratio of males to females was considered equal.
For collecting demographic, general health, psychiatric problems, and violent behavior data, the Persian versions of standard questionnaires, according to WHO-Global School-Based Student Health Survey (GSHS) were used (
19). These questionnaires were used to measure the number of aggressive behaviors, bullying, and victimization of students during the last 12 months (
9). The validity and reliability of these questionnaires were evaluated by psychiatric and psychology experts after translation. After each pilot study, questions with any difficulty in comprehension were modified. The questionnaire’s internal reliability was assessed and approved (Cronbach’s reliability coefficient > 0.7) (
20,
21). Data collecting was carried out by health experts in the provinces and city health centers. The experts consisting of healthcare professionals and nurses had passed necessary training for the completion of questionnaires (
22). The training process was conducted in a cascading manner. Some staff were initially trained at the Ministry of Health level and were referred to lower and executive levels for training other staff. Staff retraining was also conducted during the CASPIAN study as needed. To measure the engagement in bullying activities, we asked the participants how many times they were involved in bullying others or the number of times they were victimized in the recent year during school hours. The participants were categorized into three groups: (1) never, (2) 1 - 3 times, and (3) equal to or more than four times based on the number of engagements in bullying activities as the one who either bullied or was bullied, based on their roles.
The Persian version of the GSHS questionnaire was also used to screen depression, and anxiety as psychiatric distresses on bullied students and those who were not bullied and in bullies and those who were not bullies. This part of the questionnaire, which has been used in the other previous CASPIAN study, consisted of seven questions that assessed the psychological distress of students in the last six and 12 months (
23).
We applied the multilevel fixed-effect model for analyzing the data. The four survey stages (CASPIAN) were considered as level 1 factors, and age groups were considered as level 2 factors.
We expressed quantitative variables as means (95% CI) and categorized variables as percentages (95% CI). To compare psychiatric problems and violent behavior according to sex and age groups, a Wald chi-square test was applied. We reported all the results acquired based on responses in 6 - 10 and 11 - 19 age groups. In addition to information on those who bullied others and the others who were bullied, we adjusted the analysis based on general health variables. This variable describes some of the health indicators of students, such as physical conditions and so on.
The multilevel fixed-effect model and logistic multivariate regression were applied to adjust the multilevel effects and estimate the odds of anxiety and depression in bullies and bullied students. Estimated odds ratios were reported based on age and sex groups. All assumptions underlying this model were examined. The significant degree of the age-cohort period effect was examined. The collinearity among the covariates was examined using the Spearman correlation test, and no significant collinearity was detected. Given that we used the pairwise analysis method in this project, we did not take any action on missing data. Data were cleaned, coded, and entered into the software. All statistical analyses were performed at a 95% significance level using Stata software version 14 (Stata Corp, College Station, TX, USA).