- Research article
- Open Access
Individual and contextual factors associated with verbal bullying among Brazilian adolescents
BMC Pediatrics volume 15, Article number: 49 (2015)
Few studies have been carried out in low- middle-income countries assessing contextual characteristics associated with bullying. This study aimed to assess the relative importance of contextual (school and city) and individual-level factors to explain the variance in verbal bullying among a nationally representative sample of Brazilian adolescents.
59,348 students from 1,453 schools and 26 state capitals and the Federal District participated in the National Survey of School Health among 9th Grade Students (PeNSE, 2009). We performed multilevel logistic regression in a three level model (individual, school and city).
The 30-day prevalence of verbal bullying among these students was 14.2%. We found that 1.8% and 0.3% of the total variance in bullying occurred at school-level and city-level, respectively, and 97.9% at individual-level. At city-level, all factors included failed to demonstrate a significant association with bullying (p <0.05) whereas at school-level, private schools presented more bullying than public schools (OR = 1.17, CI 1.04-1.31). At individual-level, male gender, younger age, not living with both parents, exposed to domestic violence, under or overweight were all associated with bullying.
All socioeconomic indicators assessed contributed little to explain the variance in bullying at individual, school or city-level. Population subgroups at risk identified according to their individual profile could be targeted in future interventions in Brazil.
Bullying is a common form of aggression among adolescents worldwide. School-based surveys in high-income countries (HIC) have found average prevalence rates among adolescents ranging from 2% among girls aged 13 years in Italy to 32% within boys aged 11 years in Lithuania . Surveys from low- and middle-income countries (LMIC) also show high prevalence, ranging from 7.8% in Tajikistan to 60.9% in Zambia . In Latin America, relatively high prevalence has been observed in Chile (46.6%) and Venezuela (32.8%) . A recent study from Brazil stratified bullying by frequency, with a prevalence of 5.4% reporting “often and very often being victimized”, and 25.4% “sometimes or rarely being victimized” in the last 30 days . The wide variation in prevalence rates across countries is possibly related to the use of different methodologies and/or cultural or social factors [4,5].
Bullying behavior is of concern to both education and health systems, because of its deleterious consequences to victims. Bullying victims have higher risks of psychosocial problems, like low self-esteem, anxiety, depression  and suicide ; Moreover, victims present poorer school performance, and higher truancy and dropout rates [8,9].
There has been research suggesting that school bullying tends to cluster in classes, schools, cities and countries [10-15]. However, most research conducted so far shows that most of the variance is explained at an individual rather than higher level of aggregation. Notwithstanding this, risk factors should be studied in a comprehensive manner with adequate control of possible confounders  in order to get a broader understanding of the factors associated with bullying.
At the individual-level the following variables have been associated with bullying: male gender , younger age or earlier grades , being overweight or obese [19,20], exposure to physical abuse by adults [21-23], being from families that lack affection and discipline , and living only with one parent . The role of ethnicity and individual socioeconomic status remains unclear [10,25-27].
The association of macroeconomic and social indicators with bullying has also been explored at a higher level (school, city and country). For instance it is argued that income inequality, is responsible for reducing social cohesion across communities and this is considered a preventive factor to reduce violence in general [28,29] including bullying . It is also argued that socio-economic inequalities may contribute to increase abuse inflicted on those with less power .
At school-level, studies in HIC have found that poverty [30-33] or high socioeconomic disparities in schools  had a positive association with bullying. At country-level evidence suggests that income inequality, measured by the GINI index, a commonly used indicator of income inequality, also had a positive association with bullying in HIC [10,34]. Only two studies assessing concomitantly both individual- and contextual-level factors in LMIC have been reported so far [15,35]. However, none of them found an association between income or land ownership inequality or per capita income and bullying.
Violence is another characteristic that has been linked to bullying. The connection between violence and bullying in cities may be explained based on the social learning perspective that postulates that aggression is a behavior that may be learned from the environment in which individuals live. In other words, those who are exposed to violence in their communities learn that aggression, including bullying, may be an acceptable way to achieve goals, consequently developing normative beliefs, which may condone or legitimize aggression [15,36]. A study in Colombia showed that at school-level bullying was associated with higher levels of community violence and stronger beliefs in support of aggression. At city-level, those areas that had experienced more violence from the decades-old-armed conflict a had increased prevalence of bullying in schools, among 5th grade students . On the other hand an international study in 15 LMIC  found no association between homicide rate at country-level and bullying. As the association between social and macroeconomic indicators and bullying has shown disparate findings in these studies in LMIC, more research seems warranted.
Bullying may take different forms, such as physical, verbal, relational/social or electronic (cyberbullying) [37,38]. In Brazil, studies have shown that verbal bullying is the most prevalent form of bullying [39,40], which is similar to the US and UK findings [41,42].
The associations previously described came from studies that analyzed factors associated with bullying as a whole, but may also be important to explore specific predictors of each form of bullying. At the individual-level, some gender differences were found, for example in associations with verbal and physical bullying . Thus, it is possible that school and city-level predictors also vary according to the form of bullying assessed.
Brazil has the largest population in Latin America: 198.7 million people , with about 34.5 million aged 10–19 years and high school attendance, around 84% . Brazil is classified as an upper-middle income country  and in spite of not being considered a poor country it has a high rate of income inequality (GINI index: 0.54) with substantial differences between geographical regions. Also, violence is an important cause of morbidity and mortality in Brazil, with most violent deaths being due to homicide and with an unequal distribution between geographic regions . Social and income inequalities as well as violence are exceedingly common in Brazil and may be associated with bullying. Nonetheless, few studies on the prevalence and risk factors associated with bullying have been conducted in Brazil and all of them explored only individual-level factors [3,39,46,47].
Using a dataset from a Brazilian surveillance, which assessed verbal bullying among adolescents, we hypothesize that: schools and cities may contribute to explain some of the variance in verbal bullying among young people in Brazil. We also intend to explore associations between verbal bullying and individual characteristics, considering that previous research have shown that individual features are more strongly correlated with bullying.
Study population, sampling and data collection
The present study is based on data supplied by the National Survey of School Health among 9th Grade Students - PeNSE, which was carried out between March and July of 2009 .The main objective of this survey was to assess the risk and protective factors for health in adolescents from public and private schools in Brazil. The study population comprised adolescents attending the last year of private and public middle schools (grade 9) from 26 state capitals and the Federal District, which is one of the 27 federative units, in Brazil. In this paper, the Federal district will be analysed as a cluster in the same level as the other 26 state capitals.
The expected age range in 9th grade in Brazil is between 14 and 15 years old , but the real range is wider than this, especially due to students repeating grades or going to school before the age of 6.
The sample was selected from the 2007 School Census database using a complex design that included stratification per cluster and two-stage selection. The sampling strata considered each of the 26 state capitals and the Federal District. The primary sampling units (PSUs) were schools, and the secondary sampling units (SSUs) were school classes. The odds for school selection were proportional to the school size (total number of ninth-year classes), while the classes in each school were chosen by simple random selection. Two classes were selected from the schools with three or more ninth-year classes, and one class was selected from the schools with one or two ninth-year classes. All of the students enrolled in the selected classes were invited to participate in the study .
PeNSE 2009 used a self-reported structured questionnaire available in a Personal Digital Assistantb device, which included the following thematic modules: socio-demographic characteristics, diet, body image, physical activity, smoking, use of alcohol and other drugs, support network (family and friends), sexual behavior, situations at home and school (including bullying), violence and accidents, and anthropometric measures (weight and height). The questionnaire was based on the Global School-based Student Health Survey/World Health Organization-GSHS/WHO and Youth Risk Behavior Surveillance System/Centers for Disease Control and Prevention- YRBSS/CDC.
From the total number of students selected for the sample (n = 63,411), 501 refused to participate, 1,937 did not report their gender and 1,625 did not answer the bullying question and thus were excluded from analysis (response rate 93.6%). The sample was re-weighted (based on gender losses and refusals) to represent students enrolled in the ninth year of schooling who attended regularly. Therefore, the analysis uses data corresponding to 59,348 students from 1,453 schools and 26 cities and the Federal District. Further details on the sampling procedures are available in the PeNSE report .
Assessment of verbal bullying
Verbal victimization was measured using the question “In the past 30 days, how often have you been mocked, teased, called names or intimidated by one of your schoolmates so much that you were hurt/annoyed/upset/offended/ashamed?” Response options were “not at all”, “rarely”, “sometimes”, “most of the time”, and “always”. Students were categorized as victims of bullying if they reported being verbally bullied “sometimes” or more frequently. We assumed that the answer “sometimes” implies more than once, which is equivalent to the frequency of “two or more times within the last month”, a threshold recommended in previous literature [27,50].
Description of independent variables
The following socio-demographic variables were considered in the analysis: gender; age range (≤13, 14, 15, ≥16 years); ethnicity/skin colour (white, black, brown or mixed-race, Asian, native Brazilian Indian); mother’s educational level (incomplete middle school, complete middle school, complete high-school, complete higher education). The socioeconomic-status was assessed as tertiles of the total scores of the Goods and Service index (GS score) . The GS score was built based on the self-report of having: television, refrigerator, stove, microwave, washing machine, landline, mobile phone, DVD player, computer, car, bathroom inside the house and housemaid services. Each item was weighted by the inverse of the frequency of possession or presence in the total study sample. In addition, an imputation procedure was performed to attribute numerical values to missing score items as described previously . The GS score of each student was obtained by adding the weighted scores.
The assessment of weight and height was performed by trained researchers and is described elsewhere . After calculating Body Mass Index (BMI), we used the BMI index for age z-scores as an assessment of nutritional status. The nutritional status categories were: normal weight (≥ −2 z-score ≤ +1), thinness (z-score < −2), overweight (+1 > z-score ≤ +2) and obese (z-score > + 2) according to age and sex following recommendations by the World Health Organization .
Living arrangements were also evaluated and categorized into three groups: living with two parents, living with one parent (mother or father) or other arrangements. There was no differentiation between biological or adoptive parents. Domestic violence against the adolescent was evaluated from the question “in the past 30 days, how many times were you physically assaulted by an adult in your family”, categorized as “none”, “once”, or “two or more times”.
School-level contextual measures
Schools were classified as either public or private. We estimated a school’s socioeconomic level as the mean GS score of the children supplying data from each school. The GS score was distributed into tertiles. We used the coefficient of variation of mean GS score as a measure of a school’s socioeconomic disparity (i.e. a high value indicated a larger variation among students at the school). The coefficient of variation was also divided into tertiles.
Contextual city- level measures
The 2010 per capita income of each city obtained from the Brazilian Institute of Geography and Statistics – IBGE (http://www.ibge.gov.br/home/estatistica/populacao/censo2010/indicadores_sociais_municipais/tabelas_pdf/tab16.pdf), divided in tertiles. It was used as an indicator of overall city income. In addition we obtained GINI indices for each capital, also from IBGE system (2010=http://tabnet.datasus.gov.br/cgi/idb2011/b09capc.htm), which were subsequently divided into tertiles. A GINI index of zero represents perfect income equality and a score of 1 maximum inequality. Violence in the capital was assessed by homicide rate (2009) obtained from the Brazilian Department of Public Health Information  and it was also divided into tertiles.
First the prevalence and distribution of the covariates of interest were analysed incorporating the complex sample design of PeNSE 2009. The relationships between being a victim of verbal bullying and individual-level, school-level and city-level characteristics were assessed by three-level logistic regression models, whereby crude and adjusted odds ratios (ORs) with the corresponding 95% confidence intervals were obtained. The cut-off point for statistical significance was arbitrarily established as p < 0.05.
We performed multilevel analysis at different levels taking into consideration the data hierarchy (students nested in schools and schools nested in cities) to analyze the proportion of the variance explained at different levels using regression models .
Initially we calculated the crude variance coefficients representing between-city and between-school variance in exposure to bullying victimization without considering any individual, school, or city variables (empty model). This model indicates the crude amount of clustering of bullying by school and city. We used the ‘latent variable method’ to obtain intraclass correlation coefficients (ICC), assuming that an individual’s risk of bullying follows a logistic distribution, with individual-level variance equals to 3.29 [56,57].
In a second step, we analyzed the association between each covariate with bullying, without additional adjustments for other variables (Unadjusted model). Then, the influence of all the individual-level covariates on the outcome was examined (model 1). In Model 2 we simultaneously examined the influence of all covariates, student, school and city-level. We also verified the proportional change in variance of bullying across schools and cities for models 1 and 2 .
All the analyses were performed using Stata SE version 13 .
PeNSE 2009 was approved by the National Commission of Research Ethics (ComissãoNacional de ÉticaemPesquisa – Conep), record no. 11.537. It was performed in accordance with the Declaration of Helsinki and all participants gave their informed consent. Database was made publicly available on an IBGE website without any information that could identify subjects.
The descriptive analysis of individual, school and city-level characteristics are presented in Table 1. Students were predominantly within the age range between 14-15years-old (65.39%), white (40.23%) or brown/mixed race (39.11%), and lived with both parents (58.06%). The Goods and Services (GS) Score ranged from 0 to 22, and the mean was 11.92 (SE = 0.07). More than two-tenths of students were overweight and obese; nearly 10% had suffered violence at least once from someone in the family within the last 30 days (Table 1). The prevalence of verbal bullying among students was 14.2% (95% CI: 13.60 -14.75). Most of schools were public (75.95%), the mean of GS score ranged from 5.70 to 21.14 among schools, and the coefficient of variance of GS score ranged from 9.14% to 64.17% among schools.
There was wide variation in the prevalence of bullying among schools, ranging from 1.8% to 40.0%, and between cities ranged from 10.9% (Cuiabá) to 16.3% (Curitiba). In the empty model, the ICC showed that 1.8% and 0.3% of the total variance in verbal bullying occurred at school-level and city-level, respectively, and 97.9% at individual-level. These variances were statistically significant, showing clustering of bullying at these levels. The inclusion of covariates at individual-level in model 1 were able to explain 33.3% of the school-level variance. The inclusion of covariates at school and city-level in model 2 were not able to explain the school and city-level variance (Table 2).
Factors associated with verbal bullying victimization
Contextual-level characteristics (school- and city-level)
In the unadjusted model schools with higher mean of GS scores (third tertile) presented a higher likelihood of victimization and private schools were also more likely to report bullying. In the adjusted model, there were no associations between bullying and the school-level mean of GS scores or the coefficient of variation (representing socio-economic inequalities). However private schools remained associated with more bullying victimization in the fully adjusted model (see Table 3). At the city-level, per capita income was positively associated with bullying only in the unadjusted model. Neither income inequality (GINI index), nor the homicide rate were associated with bullying.
With regards to the association between verbal bullying and individual-level characteristics in adjusted models, there were significant associations for: male gender (OR = 1.39, IC 1.32-1.47), younger age (OR = 1.61, IC1.44-1.80), native ethnic groups (OR = 1.17, CI 1.03 -1.34), not living with both parents (ORone parent = 1.07, CI 1.01-1.13; ORother = 1.14, CI 1.02-1.29), victim of violence within the family at least once in the last month (OR once = 2.09, IC 1.89-2.31; OR twice or more = 3.42, IC 3.08-3.79), overweight (OR = 1.17, IC 1.08-1.26), obesity (OR = 1.41, IC 1.27-1.56) and underweight (OR = 1.50, IC 1.29-1.73) (Table 3).
This is the first large multilevel study in Brazil to investigate the association between verbal bullying and individual and contextual covariates among a representative sample of Brazilian adolescents. The prevalence of self-reported verbal bullying was 14.2%, and most of the variance in bullying occurred at individual-level. Even though, the proportion of the variance explained at school level was fairly small, we found that students from private schools were more likely to report verbal bullying. At city level, all factors failed to demonstrate an association with verbal bullying. At the individual-level, male gender, young, not living with both parents, exposed to domestic violence, underweight, overweight or obese were significantly associated with bullying.
The relative contribution of school and city-level in verbal bullying variance among this large sample of Brazilian students was of small magnitude but in keeping with other studies ranging from 0.6 to 4.0% within schools, and from 1.7 to 9.1% within cities [12,14,15,30,60]. However it is noteworthy that we found a positive association between private schools and bullying. In Brazil there is a distinct socioeconomic difference between those who attend public and private schools, with the lower socioeconomic classes attending public schools. Thus, school status is often considered a socioeconomic indicator in Brazil. However, this association needs to be treated with caution because we did not find an association between socioeconomic level according to the individual GS score or maternal education and bullying. Further research is needed to explore what these complex variables actually represent. Chaux et al.  found similar association in Colombia, but stressed that it is not possible to state whether there is more awareness about bullying in private schools .
The lack of association between income inequality aggregated at school-level with verbal bullying, which had been previously reported to bullying in the literature , might be due to the fact that Brazilian schools are quite homogeneous in this respect.
Despite the fact that we have not found significant association between verbal bullying and income inequality at city-level, this association cannot be ruled out. A limitation of our study is that the indicator of inequality could not be measured at the most appropriate level. Most inequalities in Brazil are more noticeable when comparing large cities and small towns, urban and rural areas and within cities comparing central and peripheral areas but not when comparing state capitals. Nonetheless, we were unable to obtain this data because it was considered that this information could lead to the identification of the schools. It is also possible that the association between inequality and bullying may be more related to other forms of bullying or attributable to individual level factors, such as age and sex, obesity and being victim of domestic violence, rather than contextual variables. Regardless these possibilities, findings from other LMIC have not supported the association between inequality and bullying [15,35].
No association between either poverty or homicide rate and verbal bullying were found in Brazilian data. This suggests that regardless of the theoretical assumptions and findings linking violence and bullying  in Brazilian settings these features do not seem to explain verbal bullying variance among adolescents. The lack of association might be attributable to the measurement of verbal bullying. Perhaps if we had measured bullying as a summary of all its forms (physical, verbal, relational, cyberbullying) these associations could have been found. Despite this, our result is in keeping with those reported by Wilson et al.,  in LMICs, for bullying measured as a whole.
Overall the associations between individual factors and verbal bullying were consistent with the literature, but cannot be extrapolated to other forms of bullying. Boys were more likely to report victimization, as reported by others especially for physical and verbal bullying, which are overt forms of aggression . It is possible that girls may be more subtle in their bullying or less willing to report it than boys . Young students reported victimization more often, something that may reflect the power imbalance when students of different ages attend the same classes. Those adolescents not living with both parents were more likely to report victimization. Growing up in single-parent families seems to be associated with more bullying victimization . Jablonska and Lindberg  stated that social and economic characteristics and physical and mental problems of the parents are more commonly found in single parent households, perhaps explaining, at least in part, the association between single parents and bullying. However, our data do not allow us to examine these causal pathways.
Native Brazilian indian students were more likely to report victimization than white students, and this was the only race/color variable associated with verbal bullying in our sample. Being bullied is often related to being different, thus native individuals, representing a minority group in the Brazilian population (and in the sample of students, 4.1%), may be perceived as different for a number of reasons, including different cultural and social norms .
Individuals who reported violence at home were much more likely to report verbal bullying victimization. Although we cannot establish the causal direction of this association, it has been reported that domestic violence increases anxiety and depressive symptoms, low self-esteem and social isolation, which might interfere with the social skills of these adolescents who could become easy targets for bullies [64,65].
Underweight, overweight and obese adolescents reported more bullying than their eutrophic counterparts, something that had been previously reported in the literature [20,66]. This study, however, uses a more accurate measure of weight and height that were directly assessed and not self-reported by students, as in previous [19,20,67]. It has been reported that those adolescents whose body weight does not fit with what is seen as desirable in a particular setting (ideal body shape) have greater difficulty socializing in school and poorer psychological well-being, something that increases their likelihood to be bullied .
This study has some limitations. First, the measure of bullying based on a single-item is rather limiting but other studies have also had a similar problem. This question seeks to capture those individuals who had experienced negative feelings as a result of bullying. Students classified as victims in this study might be more severe cases of bullying or those who were unable to cope well with the situation. Therefore it is possible that the prevalence of verbal bullying as measure in this study may have been underestimated. Additionally, the associations found or not found, can only be stated for verbal bullying and might differ for other forms of bullying or for bullying measured as a summary of its different forms. However, the lack of these measures in the surveillance used prevented us from exploring other associations.
The second limitation is the use of aggregated measures to assess school characteristics instead of using proper contextual level measures . Unfortunately contextual measures were not available at school level. However, some authors have argued that it is valid to use summarized or aggregated measures to assess contextual characteristics in multilevel analyses . Third, it was not possible to include probability weights in the multilevel regression analysis using Stata, since the full model could not achieve convergence using Gllamm commands. Nevertheless, the sample structure was considered in the analyses using XT commands in Stata, and comparing three-levels models with and without sample weights the results were similar. Finally, when interpreting the findings of this study, the cross-sectional nature of our data did not allow us to make causal inferences.
The major strength of this study is the large random sample, representative of Brazilian schoolchildren of 9th grade and the good response rate. The direct measurement of height and weight of students rather than self-reported measures also reduced potential measurement bias .
This study provided evidence that verbal bullying variance at school and city-level was small, and most of the variance was accounted for at the individual level. This does not mean school or city-level lacks importance, but may suggest that city boundaries comparing only state capitals do not capture differences that shape the relevant environment for bullying occurrence . Those contextual variables analyzed contributed little to explain the variance in verbal bullying, between schools and cities, among 9th grade school students in Brazil. However, the lack of association between inequality at city-level and verbal bullying should be read with caution, since measurement concerns were described.
Those individual variables strongly associated with verbal bullying might allow the identification of possible higher risk sub-groups that could be targeted for interventions in the future and denote the role of family arrangements and violence in verbal bullying victimization. The profile of vulnerable adolescents to verbal bullying in Brazil can guide policy makers, teachers and school staff on bullying prevention.
Finally studies exploring associations between other forms of bullying and other contextual variables, measured at rural/urban or central/peripheral neighborhoods; or comparing state capitals with small towns; including other family and school characteristics, which are more proximal to individuals, may provide a more accurate picture of the possible influence of contextual-level variables on bullying in Brazil.
aDecades-old-armed conflict has happened among left-wing guerrilla groups, right-wing para-militaries and the national-armed forces in Colombia.
Low- and middle-income countries
National Survey of Health among 9th Grade Students
Primary sampling units
Secondary sampling units
- GS score:
Goods and Service scores
Intraclass correlation coefficients
Currie C, Zanotti C, Morgan A, Currie D, Looze M, Roberts C, et al. Social determinants of health and well-being among young people. Health Behaviour in School-aged children (HBSC) study: International Report From the 2009/2010 Survey. In: Health Policy for Children and Adolescents. Copenhagen: WHO Regional Office for Europe; 2012. p. 252.
Fleming LC, Jacobsen KH. Bullying among middle-school students in low and middle income countries. Health Promot Int. 2010;25(1):73–84. Doi: 10.1093/heapro/dap1046. Epub 2009 Nov 1092.
Malta DC, Silva MAI, Mello FCM, Monteiro RA, Sardinha LMV, Crespo C, et al. Bullying nas escolas brasileiras: resultados da Pesquisa Nacional de Saúde do Escolar (PeNSE), 2009. Cien Saude Colet. 2010;15:3065–76.
Smith PK, Cowie H, Olafsson RF, Liefooghe AP, Almeida A, Araki H, et al. Definitions of bullying: a comparison of terms used, and age and gender differences, in a fourteen-country international comparison. Child Dev. 2002;73(4):1119–33.
Liu J, Graves N. Childhood Bullying: A Review of Constructs, Concepts, and Nursing Implications. Public Health Nurs. 2011;28(6):556–68.
Salmon G, James A, Smith DM. Bullying in schools: self reported anxiety, depression, and self esteem in secondary school children. BMJ. 1998;317(7163):924–5.
Klomek AB, Sourander A, Gould M. The Association of Suicide and Bullying in Childhood to Young Adulthood: A Review of Cross-Sectional and Longitudinal Research Findings. Can J Psychiatry. 2010;55(5):282–U221.
Glew GM, Fan MY, Katon W, Rivara FP, Kernic MA. Bullying, psychosocial adjustment, and academic performance in elementary school. Arch Pediatr Adolesc Med. 2005;159(11):1026–31.
Gastic B. School truancy and the disciplinary problems of bullying victims. Educ Rev. 2008;60(4):391–404.
Due P, Merlo J, Harel-Fisch Y, Damsgaard MT, Holstein BE, Hetland J, et al. Socioeconomic inequality in exposure to bullying during adolescence: a comparative, cross-sectional, multilevel study in 35 countries. Am J Public Health. 2009;99(5):907–14.
Karna A, Voeten M, Poskiparta E, Salmivalli C. Vulnerable Children in Varying Classroom Contexts Bystanders’ Behaviors Moderate the Effects of Risk Factors on Victimization. Merrill-Palmer Q J Dev Psychol. 2010;56(3):261–82.
Larochette A-C, Murphy AN, Craig WM. Racial Bullying and Victimization in Canadian School-Aged Children Individual and School Level Effects. Sch Psychol Int. 2010;31(4):389–408.
Vervoort MHM, Scholte RHJ, Overbeek G. Bullying and Victimization Among Adolescents: The Role of Ethnicity and Ethnic Composition of School Class. J Youth Adolesc. 2010;39(1):1–11.
Wei H-S, Williams JH, Chen J-K, Chang H-Y. The effects of individual characteristics, teacher practice, and school organizational factors on students’ bullying: A multilevel analysis of public middle schools in Taiwan. Child Youth Serv Rev. 2010;32(1):137–43.
Chaux E, Molano A, Podlesky P. Socio-Economic, Socio-Political and Socio-Emotional Variables Explaining School Bullying: A Country-Wide Multilevel Analysis. Aggress Behav. 2009;35(6):520–9.
Khoury-Kassabri M, Benbenishty R, Astor RA, Zeira A. The contributions of community, family, and school variables to student victimization. Am J Community Psychol. 2004;34(3–4):187–204.
Nansel TR, Overpeck M, Pilla RS, Ruan WJ, Simons-Morton B, Scheidt P. Bullying behaviors among US youth - Prevalence and association with psychosocial adjustment. JAMA. 2001;285(16):2094–100.
Vieno A, Gini G, Santinello M. Different Forms of Bullying and Their Association to Smoking and Drinking Behavior in Italian Adolescents. J Sch Health. 2011;81(7):393–9.
Janssen I, Craig WM, Boyce WF, Pickett W. Associations between overweight and obesity with bullying behaviors in school-aged children. Pediatrics. 2004;113(5):1187–94.
Puhl RM, Luedicke J, Heuer C. Weight-based victimization toward overweight adolescents: observations and reactions of peers. J Sch Health. 2011;81(11):696–703. doi:610.1111/j.1746-1561.2011.00646.x.
Baldry AC. Bullying in schools and exposure to domestic violence. Child Abuse Negl. 2003;27(7):713–32.
Bowes L, Arseneault L, Maughan B, Taylor A, Caspi A, Moffitt TE. School, Neighborhood, and Family Factors Are Associated With Children’s Bullying Involvement: A Nationally Representative Longitudinal Study. J Am Acad Child Adolesc Psychiatry. 2009;48(5):545–53.
Björkqvist K, Osterman K, Berg P. Higher rates of victimization to physical abuse by adults found among victims of school bullying. Psychol Rep. 2011;109(1):167–8.
Debarbieux E, Blaya C. Violência nas Escolas e Políticas Públicas. Brasília: UNESCO; 2002. p. 268.
Jansen DEMC, Veenstra R, Ormel J, Verhulst FC, Reijneveld SA. Early risk factors for being a bully, victim, or bully/victim in late elementary and early secondary education. The longitudinal TRAILS study. Bmc Public Health. 2011;11:440.
Barboza GE, Schiamberg LB, Oehmke J, Korzeniewski SJ, Post LA, Heraux CG. Individual Characteristics and the Multiple Contexts of Adolescent Bullying: An Ecological Perspective. J Youth Adolesc. 2009;38(1):101–21.
Currie C, Gabhainn SN, Godeau E, Roberts C, Smith R, Currie D, et al. Inequalities in young people’s health. HBSC international report from the 2005/2006 Survey. In: Organization WH, editor. Health Policy for Children and Adolescents, vol. 5. Copenhagen: WHO Regional Office for Europe; 2008. p. 206.
Elgar FJ, Aitken N. Income inequality, trust and homicide in 33 countries. Eur J Public Health. 2011;21(2):241–6. doi: 210.1093/eurpub/ckq1068. Epub 2010 Jun 1094.
Kawachi I, Kennedy BP, Wilkinson RG. Crime: social disorganization and relative deprivation. Soc Sci Med. 1999;48:719–31.
Bradshaw CP, Sawyer AL, O'Brennan LM. A Social Disorganization Perspective on Bullying-Related Attitudes and Behaviors: The Influence of School Context. Am J Community Psychol. 2009;43(3–4):204–20.
Carbone-Lopez K, Esbensen F-A, Brick BT. Correlates and Consequences of Peer Victimization: Gender Differences in Direct and Indirect Forms of Bullying. Youth Violence Juvenile Justice. 2010;8(4):332–50.
Gregory A, Cornell D, Fan X, Sheras P, Shih T-H, Huang F. Authoritative School Discipline: High School Practices Associated With Lower Bullying and Victimization. J Educ Psychol. 2010;102(2):483–96.
Klein J, Cornell D. Is the Link Between Large High Schools and Student Victimization an Illusion? J Educ Psychol. 2010;102(4):933–46.
Elgar FJ, Pickett KE, Pickett W, Craig W, Molcho M, Hurrelmann K, et al. School bullying, homicide and income inequality: a cross-national pooled time series analysis. Int J Public Health. 2013;58(2):237–45.
Wilson ML, Dunlavy AC, Berchtold A. Determinants for Bullying Victimization among 11–16-Year-Olds in 15 Low- and Middle-Income Countries: A Multi-Level Study. Soc Sci. 2013;2:208–20.
Espelage DL, Bosworth K, Simon TR. Examining the social context of bullying behaviors in early adolescence. J Couns Dev. 2000;78(3):326–33.
Olweus D. Bullying at School: What we know and what we can do. Oxford, U.K: Blackwell; 1993.
Smith PK, Mahdavi J, Carvalho M, Fisher S, Russell S, Tippett N. Cyberbullying: its nature and impact in secondary school pupils. J Child Psychol Psychiatry. 2008;49(4):376–85.
de Moura DR, Nova Cruz AC, Quevedo LA. Prevalence and characteristics of school age bullying victims. J De Pediatria. 2011;87(1):19–23.
Santos JA, Cabral-Xavier AF, Paiva SM, Leite-Cavalcanti A. Prevalência e Tipos de Bullying em Escolares Brasileiros de 13 a 17 anos. Rev salud pública. 2014;16(2):173–83.
Wang J, Iannotti RJ, Nansel TR. School Bullying Among US Adolescents: Physical, Verbal, Relational and Cyber. J Adolesc Health. 2009;45(4):368–75.
Rivers I, Smith P. Types of bullying behaviour and their correlates. Aggress Behav. 1994;20:359–68.
WorldBank. 2012. http://data.worldbank.org/country/brazil.
IBGE. Pesquisa Nacional por Amostra de Domicilios. Rio de Janeiro: Brasil; 2008.
Reichenheim ME, de Souza ER, Moraes CL, de Mello Jorge MH, da Silva CM, de Souza Minayo MC. Violence and injuries in Brazil: the effect, progress made, and challenges ahead. Lancet. 2011;377(9781):1962–75. doi:1910.1016/S0140-6736(1911)60053-60056. Epub 62011 May 60059.
Brito CC, Oliveira MT. Bullying and self-esteem in adolescents from public schools. J Pediatr (Rio J). 2013;89(6):601–7. doi:610.1016/j.jped.2013.1004.1001. Epub 2013 Sep 1019.
Rech RR, Halpern R, Tedesco A, Santos DF. Prevalence and characteristics of victims and perpetrators of bullying. J Pediatr (Rio J). 2013;89(2):164–70. doi:110.1016/j.jped.2013.1003.1006.
IBGE. Pesquisa Nacional de Saúde do Escolar 2009. In: Ministério do Planejamento OeGIDdPCdPeIS. Rio de Janeiro: Instituto Brasileiro de Geografia e Estatística; 2009. p. 138.
Ministério da Educação. Ensino Fundamental de nove anos orientações para a inclusão da criança de seis anos de idade. Brasília: Ministério da Educação, Secretaria de Educação Básica; 2007. p. 135 p. : il.
Felix ED, Sharkey JD, Green JG JFM, Tanigawa D. Getting precise and pragmatic about the assesment of bullying: the development of the California Bullying Victimization Scale. Agressive Behav. 2011;37:234–47.
Levy RB, de Castro IRR, Cardoso LO, Tavares LF, Sardinha LMV, Gomes FS, et al. Consumo e comportamento alimentar entre adolescentes brasileiros: Pesquisa Nacional de Saúde do Escolar (PeNSE), 2009. Cien Saude Colet. 2010;15(Supl. 2):3085–97.
Araújo C, Toral N, da Silva A, Velásquez-Melendez G, Dias A. Estado nutricional dos adolescentes e sua relação com variáveis sociodemográficas: Pesquisa Nacional de Saúde do Escolar (PeNSE), 2009. Revista Ciência e Saúde Coletiva. 2010;15(2):3077–84.
Onis M, Onyango A, Borghi E, Siyam A, Nishida C, Siekmann J. Development of a WHO growth reference for school-aged children and adolescents. Bull World Health Organ. 2007;85(9):660–7.
DATASUS. Sistema de Informação sobre Mortalidade. MS/SVS/DASIS. Brasília: MS; 2009.
Snijders TAB, Bosker RJ. Multilevel Analysis: An Introduction to Basic and Advanced Multilevel Modeling. London: Sage Publication; 1999.
Merlo J, Chaix B, Ohlsson H, Beckman A, Johnell K, Hjerpe P, et al. A brief conceptual tutorial of multilevel analysis in social epidemiology: using measures of clustering in multilevel logistic regression to investigate contextual phenomena. J Epidemiol Community Health. 2006;60(4):290–7.
Rasbash J, Steele F, Browne WJ, Goldstein H. A user’s guide to MLwiN. United Kingdom: University of Bristol; 2009.
Merlo J, Yang M, Chaix B, Lynch J, Rastam L. A brief conceptual tutorial on multilevel analysis in social epidemiology: investigating contextual phenomena in different groups of people. J Epidemiol Community Health. 2005;59(9):729–36.
StataCorp. Stata Statistical Software: Release 13. College Station, TX: StataCorp LP; 2013.
Sapouna M. Collective Efficacy in the School Context: Does It Help Explain Victimization and Bullying Among Greek Primary and Secondary School Students? J Interpers Violence. 2010;25(10):1912–27.
Flouri E, Buchanan A. The Role of Mother Involvement and Father Involvement in Adolescent Bullying Behavior. J Interpers Violence. 2003;18(6):634–44.
Jablonska B, Lindberg L. Risk behaviours, victimisation and mental distress among adolescents in different family structures. Soc Psychiatry Psychiatr Epidemiol. 2007;42(8):656–63.
Sweeting H, West P. Being different: correlates of the experience of teasing and bullying at age 11. Res Pap Educ. 2001;16(3):225–46.
Gershoff ET. Corporal Punishment by Parents and Associated Child Behaviors and Experiences: A Meta-Analytic and Theoretical Review. Psychol Bull. 2002;128(4):539–79.
Bordin I, Duarte C, Peres C, Nascimento R, Curto B, Paula C. Severe physical punishment: risk of mental health problems for poor urban children in Brazil. Bull World Health Organ. 2009;87:336–44.
Wang J, Iannotti RJ, Luk JW. Bullying victimization among underweight and overweight U.S. youth: differential associations for boys and girls. J Adolesc Health. 2010;47(1):99–101. doi:110.1016/j.jadohealth.2009.1012.1007. Epub 2010 Feb 1014.
Brixval CS, Rayce SL, Rasmussen M, Holstein BE, Due P. Overweight, body image and bullying--an epidemiological study of 11- to 15-years olds. Eur J Public Health. 2012;22(1):126–30. doi:110.1093/eurpub/ckr1010. Epub 2011 Mar 1097.
Falkner NH, Neumark-Sztainer D, Story M, Jeffery RW, Beuhring T, Resnick MD. Social, educational, and psychological correlates of weight status in adolescents. Obes Res. 2001;9(1):32–42.
Diez-Roux AV. Investigating Neighborhood and Area Effects on Health. Am J Public Health. 2001;91(11):1783–9.
Suzuki E, Yamamoto E, Takao S, Kawachi I, Subramanian SV. Clarifying the Use of Aggregated Exposures in Multilevel Models: Self-Included vs. Self-Excluded Measures. PloS One. 2012;7(12):9.
Sherry B, Jefferds M, Grummer-Strawn LM. Accuracy of Adolescent Self-report of Height and Weight in Assessing Overweight Status: A Literature Review. Arch Pediatr Adolesc Med. 2007;161(12):1154–61.
We would like to thank Márcio Augusto Diniz and Juliana Yukari Kodaira Visconde for providing guidance in the multilevel analysis. This research was supported by the Brazilian National Council of Scientific and Technological Development (CNPq), n°473502/2013-5, awarded to Paulo Rossi Menezes.
The authors declare that they have no competing interests.
CMA participated in conceptualization of the study, data preparation, analysis and interpretation of data, and drafted the manuscript; RBL participated in conceptualization of the study, analysis, interpretation of data and reviewed the manuscript. RA participated in conceptualization of the study, interpretation of data and critically reviewed the manuscript. PRM participated in conceptualization of the study, interpretation of data and reviewed the manuscript. All authors have read and approved the final manuscript.
About this article
Cite this article
Azeredo, C.M., Levy, R.B., Araya, R. et al. Individual and contextual factors associated with verbal bullying among Brazilian adolescents. BMC Pediatr 15, 49 (2015). https://doi.org/10.1186/s12887-015-0367-y
- Associated factor
- Multilevel modelling