Peer victimization in single‐grade and multigrade classrooms

Abstract Although peer victimization mainly takes place within classrooms, little is known about the impact of the classroom context. To this end, we examined whether single‐grade and multigrade classrooms (referring to classrooms with one and two grades in the same room) differ in victim–bully relationships in a sample of elementary school children (646 students; age 8–12 years; 50% boys). The occurrence of victim–bully relationships was similar in single‐grade and multigrade classrooms formed for administrative reasons, but lower in multigrade classrooms formed for pedagogical reasons. Social network analyses did not provide evidence that peer victimization depended on age differences between children in any of the three classroom contexts. Moreover, in administrative multigrade classrooms, cross‐grade victim–bully relationships were less likely than same‐grade victim–bully relationships. The findings did not indicate that children in administrative multigrade classrooms are better or worse off in terms of victim–bully relationships than are children in single‐grade classrooms.

age differences (Chaux & Castellanos, 2015), such as in multigrade classrooms. Alternatively, from an evolutionary perspective (Ellis et al., 2012), it can be argued that multigrade classrooms encourage prosocial behavior in children, resulting in a lower risk of peer victimization for younger children by older children, particularly in multigrade classrooms formed for pedagogical reasons rather than administrative reasons. In this study, we distinguish between the two classroom types because pedagogical multigrade classrooms aim to stimulate prosocial relations among children by encouraging the provision of help across grades within the same classroom (Gray, 2011;Lillard & Else-Quest, 2006), whereas administrative multigrade classrooms do not have such an explicit goal.
Currently, it is unclear whether multigrade classrooms in comparison to single-grade classrooms are either a risk or protective factor for peer victimization in childhood. We aimed to explore the two competing perspectives by examining victim-bully relationships in single-grade and multigrade classrooms in a sample of elementary school children. Investigating classroom differences is an important direction in bullying research as a better understanding of the role of age differences in the classroom may help to advance prevention or intervention policies.

| A status framework to understand peer victimization
One perspective toward explaining peer victimization is the status framework (Rodkin et al., 2015), which posits that bullying is targeted peer aggression within the context of a relationship of power and abuse. In such relationships of power imbalance, perpetrators are aggressors repeatedly targeting socially vulnerable individuals to maintain or gain a high social status (Salmivalli, 2010;Veenstra, Lindenberg, Munniksma, & Dijkstra, 2010). Targets of peer bullying thus suffer from their disadvantaged and isolated position in the group as they form easy targets because of their low status (Farmer et al., 2011). In this perspective, bullying perpetration is instrumental to acquire social resources and to maintain dominance over marginalized peers (Rodkin et al., 2015), and to demonstrate self-esteem and social skills (Guerra, Williams, & Sadek, 2011). This further enhances the perpetrators' positions in the group and weakens the positions of targets.
Age forms an important aspect through which social dominance hierarchies can be determined in classrooms (Ellis et al., 2012). As children grow older, they become socially, cognitively, and physically more developed (Piaget, 1953). Older children thus have a clear advantage over younger children, to obtain dominant positions in the group. They may use this advantage by targeting younger children to demonstrate social dominance. Younger children compared with older children are more vulnerable for peer victimization (Barker et al., 2008;Chaux, Molano, & Podlesky, 2009;Rivers & Smith, 1994;Scheithauer, Hayer, Petermann, & Jugert, 2006) as they are more likely to have a weaker social position in the peer group (Salmivalli, 2010). This suggests that younger children are an easy target for older children.
It makes sense that an effect of power imbalance through age would be more salient in a context with a large age-range. In contrast to single-grade classrooms where an age difference between the children of 1 year is usually the maximum to be expected, children in multigrade classrooms differ up to the number of grades combined in the classroom. Considering that students spend most of their time in school within the same classroom, this means that in single-grade classrooms, social interactions mostly take place in same-age groups, whereas in multigrade classrooms they also take place among children who clearly differ in age. Researchers have argued that as a consequence of age-mixing, multigrade classrooms may produce peer hierarchies based on students' age and thus, putting younger children at risk for victimization (Kolbert & Crothers, 2003). Hence, in multigrade classrooms older children may target younger children as a means to demonstrate dominance.
Following this line of reasoning, we expected that power imbalance would be more related to age differences in multigrade classrooms compared with single-grade classrooms, and, therefore, in multigrade classrooms we expected to find higher degrees of victimization (H1a), and higher risk of peer victimization for younger children targeted by older children (H2a). In multigrade classrooms, we can separate the age difference in a grade effect and a relative age difference, whereas in single-grade classrooms only the (relative) age difference can be examined.

| An evolutionary model to understand peer victimization
The evolutionary model of risky child/adolescent behavior provides an alternative perspective toward explaining peer victimization in classrooms (Ellis et al., 2012). This model posits that mixed-age settings, rather than age-segregated school and peer environments, are the natural context for child development. The presence of both older and younger children in mixed-age settings provides a natural hierarchy based on age. In this context, both older and younger children settle with their position in the social group, which decreases the tendency to compete for dominance and status by demonstrating aggressiveness.
Evolutionary psychologists argue that older children can serve as positive role models, and that the positive association between status and prosocial behavior reduces the need to gain status through antisocial means (Ellis et al., 2012). The research findings show that when older children are assigned to younger children as caregivers, buddies, or playmates, they tend to behave less aggressive and more prosocial toward younger children and same-age peers in other contexts as well (Ember, 1973;Gray, 2011). In addition, the findings suggest that children in mixed-age school settings interact socially across wide age ranges; older children and younger children associate with each other, and become friends with each other (Miller, 1990;Pratt, 1986). In sum, this study suggests that the presence of younger children in mixed-age settings reduces aggression and promotes nurturance and compassion in children (Gray, 2011). In contrast, age-segregated school and peer environments, such as single-grade classrooms, have been argued to evoke aggression and conflict in children, and, in such a classroom context, children may actively search for dominance (Ellis et al., 2012).
There are two reasons why multigrade classrooms are formed within schools (Veenman, 1995). First, schools deliberately form such classrooms for didactic and pedagogical purposes, for instance to enhance the classroom climate. An example of such schools is Montessori schools, which are characterized by a special set of educational materials, freedom as a student to choose own activities, collaboration between students, and individual and small group instruction in both academic and social skills (Lillard & Else-Quest, 2006). Second, schools tend to form such classrooms for administrative reasons, for example, when dealing with low enrollment and uneven classroom sizes. It is important to make a distinction between these two forms of multigrade classrooms as the basis for their formation may yield different outcomes with regard to the victim-bully relationships. Particularly pedagogical multigrade classrooms are argued and shown to have more positive outcomes because promoting prosocial behavior between the older and younger children in such classrooms is part of the school's educational philosophy (Lillard & Else-Quest, 2006;Moller, Forbes-Jones, & Hightower, 2008). By contrast, teachers in administrative multigrade classrooms were generally found to teach the grades separately (Mulryan-Kyne, 2007;Veenman, 1995; but see Mason & Burns, 1996), which decreases opportunities for prosocial behavior between older and younger children as such multigrade classrooms emphasize individualized work and do not necessarily encourage between-grade interactions (Juvonen, 2018).
Following this line of reasoning, we formulated two additional hypotheses. In pedagogical multigrade classrooms, we expected to find lower degrees of victimization (H1b) and lower risk of victimization for children targeted by older children (H2b).

| The present study
We investigate victim-bully relationships in single-grade and multigrade classrooms. We addressed two main questions: (a) Do multigrade classrooms show different degrees of victimization compared to single-grade classrooms? (b) To what extent does peer victimization depend on age differences between children, and if so, is this different for single-grade and multigrade classrooms? We focused on middle to late childhood, because in that period social dominance hierarchies are established through school bullying (Kolbert & Crothers, 2003). We controlled for sex, because boys are often more dominant and aggressive toward peers than girls (Cook, Williams, Guerra, Kim, & Sadek, 2010).

| Sample
Classrooms were drawn from the second wave of the KiVa study at the start of the school year (in October 2012). KiVa is a program aimed to reduce school bullying among children from grades 3-6 in elementary education (8-12 years) in the Netherlands (Huitsing et al., 2019;Kaufman, Kretschmer, Huitsing, & Veenstra, 2018). The 99 participating schools (66 intervention and 33 control schools) contained 25 schools with only single-grade classrooms, 39 schools with only administrative multigrade classrooms, 14 schools with only pedagogical multigrade classrooms, and 21 schools with both single-grade and administrative multigrade classrooms. Only control schools were selected for the analysis to avoid that differences between classrooms were a result of the intervention. Next, we selected schools that had either only single-grade classrooms 3-6 or only multigrade classrooms containing two grades (e.g., 3-4 or 5-6) excluding the 19 schools with other combinations.
For the remaining 14 schools, containing 38 classrooms, we applied an additional set of selection criteria for classroom's eligibility for social network analysis: first, classroom size should be larger than 15; and second, the combinations for sex (boy-boy, girl-girl, boy-girl, and girl-boy) and grade (low-low, high-high, low-high, and high-low) should contain victim-bully relationships for the reference category. The latter criterion was needed for estimating the sex and grade effects comparable across classrooms.
We chose boy-boy (meaning boys nominating other boys as their perpetrator) and low-low (meaning lower grade classmates nominating other lower grade classmates as their perpetrator) as reference.
The two selection criteria resulted in dropping 12 classrooms (see for and was used to categorize the schools into the three categories. Appendix B shows the differences between the reduced sample (the 38 classrooms) and the final network sample (the 26 classrooms). The two samples were comparable in classroom size and percentage of boys, but small differences were observed for average nominations sent for bullying perpetration per student and for the distribution of these nominations. Average degree among pedagogical multigrade classrooms was slightly higher in the network sample compared to the reduced sample (0.5 versus 0.4), whereas the distribution of degrees was slightly more skewed among administrative multigrade classrooms (see Appendix B).
In our subsample of 646 students, 324 (50.2%) students were boys and 322 (49.8%) students were girls. The average age of the sample was 10.2 years (SD = 1.2). Most students (n = 529; 81.9%) were native Dutch; 16.7% were non-Dutch, 1.4% was missing because nine children provided insufficient information about their parent's ethnic background.

| Procedure
Students filled in an internet-based questionnaire in their classroom during regular school hours. The process was administered by the teachers, who were present to answer questions and to assist the students when necessary. Before the data collection, teachers were RAMBARAN ET AL.

| 563
given detailed instructions concerning the administration of questionnaires to students. During the data collection, support was available through phone and email.
At the beginning of the questionnaire, students received information about the goal of the study, and how to fill in the questionnaire. They were told not to talk to each other or to discuss their answers when they filled out the questionnaire or afterward to ensure each other's privacy. It was explained to students that their answers would remain confidential. The teachers ensured that students who could not complete the questionnaire on the day of the data collection participated at another day within a month.
Before the first wave (the pre-assessment) of the KiVa study (and for students who were new in school, after the first wave), schools sent information letters to students' parents. A passive consent procedure allowed students or parents to opt out of students' participation. At the start of data collection (2012), universities in the Netherlands did not require IRB permission for this type of research. All procedures performed were in accordance with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. A few students did not want to participate; also a few parents objected to their child's participation. Accordingly, for the whole sample at the second wave of the KiVa study participation rate was high (96%).

| Dependent variable
Peer victimization (bullying perpetration) was measured with network nominations following a two-stage procedure. Students were first asked to indicate how often they were victimized in general in the previous months (since the summer break), according to Olweus' (1996) self-reported bully/victim items, and, to indicate this for specific form(s) of victimization; physical harm (e.g., kicked), verbal harm (e.g., name-calling), relational harm (e.g., gossiping), and cyber victimization.
When participants indicated that they were victimized by classmates at least "Once or twice" (score 2) on any item, they were presented with a roster showing the names of all classmates, and asked by whom of their classmates they were victimized (referring to "Who starts bullying you?").
Perpetration nominations were measured as either present (1) or absent (0). Students who indicated not being victimized by classmates did not fill out the nomination question. Their "answers" were considered as "structural missing" (no outgoing nomination possible). This means that victims were the only one "allowed" to send a perpetration nomination to classmates, but that everyone could receive a perpetration nomination.
Victim-bully networks were obtained based on all perpetration nominations in a particular classroom (from the target's perspective).

| Independent variables
We included sex (1 = boy) as a control. Students' grade was obtained from the school's office. The multigrade classrooms in our study included two grades in each classroom (3-4 or 5-6); students were categorized as belonging either to the lower grade or to the higher grade in their classroom. Students' birth dates were also obtained from the school's office. Age was derived from birth date and measured in number of months. We computed relative age by subtracting the median age in the classroom (calculated among all students in the classroom) from students' own age. Relative age in multigrade classrooms was obtained by subtracting 12 months from the higher grade students' age, to compensate for a 1-year age difference between the two grades. This way, children's age in singlegrade and multigrade classrooms were made comparable to separate grade differences from (relative) age differences.

| Analytic strategy
The investigation of the victim-bully relationships in single-grade and multigrade classrooms in elementary school was done with two analyses. We tested our first set of hypotheses (H1a and H1b) using nonparametric analysis using Kruskal-Wallis H tests, and our second set of hypotheses (H2a and H2b) using cross-sectional social network analysis using ERGMs (Exponential Random Graph Models; Lusher, Koskinen, & Robins, 2013). The Kruskal-Wallis test is used to test differences in the median scores of tie percentage at the classroom level (average in-degrees) across the three classroom types (single-grade, administrative multigrade, and pedagogical multigrade). ERGMs can be used to analyze cross-sectional social network structures (Robins, 2015), and have been used before to examine victim-bully relationships . This choice allows us to simultaneously investigate individual and dyadic age and grade effects while explicitly taking into account the dependencies in the victim-bully network.
Cross-sectional social network analysis was done in PNet (Wang, Robins, & Pattison, 2009), a software program for the statistical analysis of social network data using ERGMs. The effects were first analyzed per classroom network, and were then meta-analyzed in R with classroom type as "explanatory variable" (referring to a mixedeffects model in metafor, see Viechtbauer, 2010). Each classroom network was estimated with the same model specification. For some classroom networks, however, some parameters were left out, because the accompanying statistics (network configurations) were not there. The usual criteria for convergence (absolute value of t statistics below 0.10 for all parameters; see Wang et al., 2009) was obtained for all classroom networks. Appendix C provides Goodness of Fit (GoF) statistics. As shown, classroom networks adhered to the usual criteria for acceptable GoF statistics (absolute value of t statistics below 2). In four classrooms (7, 8, 12, and 23), the fit was not optimal for the reciprocal age-related effects. Because we already included nonreciprocal age effects, including additional age effects resulted in non-convergence.

| Model specification and effect interpretation
To adequately capture important features of the victim-bully networks, we followed previous research in choosing the structural parameters in the exponential random graph models , using alternating or geometrically weighted versions of the structural statistics. Density indicates the general occurrence of peer victimization ties, comparable with the intercept or grand mean in (generalized) linear models.
The isolates parameter indicates the occurrence of a network configuration where a child does not nominate others as their perpetrators or receive a perpetration nomination from others (noninvolvement), whereas the sinks parameter indicates the occurrence of a configuration where a child receives a perpetration nomination, but does not nominate anyone (referring to perpetrators who are not targets themselves). The multiple two-paths parameter reflects a perpetrator-target, referring to a child receiving a perpetration nomination and nominating a perpetrator. The in-ties spread parameter represents variability in receiving nominations as a perpetrator, whereas the shared in-ties parameter represents variability in sending perpetration nominations by some children attracting more perpetration nominations than others. If a structural parameter is significant, this indicates that the structural configuration occurs more frequently than if perpetration nominations would be formed at random. We included sex as a covariate by specifying dyadic covariates indicating the sex combinations (taking the boy-boy dyad as the reference category).
To test the hypotheses H2a and H2b, individual and dyadic age and grade effects were included in the models. Thus, the effect of age difference is separated in three components that are comparable over the three classroom types, distinguishing differences due to relative age or grade, while controlling for a "main" or general age effect by a receiver effect of the relative age of the child. Note that no grade difference was defined for the single grade classrooms. Table 1 presents the summarized descriptive findings of the three types of classrooms. Appendix D provides information per classroom.

| Descriptive findings
Peer victimization (bullying perpetration) was, on average, almost two times as low in pedagogical multigrade classrooms (av. degree = 0.5) compared with both single-grade classrooms and administrative multigrade classrooms, which were comparable in degrees of peer victimization (av. degree = 0.9).
Table 1 also shows on average higher bullying involvement with boys: they nominated mainly each other as perpetrators, and were also nominated as perpetrators by girls. In multigrade classrooms, bullying perpetration nominations occurred mainly among samegrade students: most students indicated that they were being targeted by a peer in their own grade (see av. densities for the low-low group and the high-high group in Table 1). Appendix D shows the relative age differences between targets' perpetrators and targets' nonperpetrators, showing no clear pattern as children were targeted approximately equally often by older and younger classmates, except for the pedagogical 5-6 grades, where targets more often had older perpetrators.
The Kruskal-Wallis test shows no statistically significant differences in the median scores of tie percentage at the classroom level (average in-degrees) between the three classroom types (singlegrade, administrative multigrade, and pedagogical multigrade). This means that there was no evidence that the three classroom types differ in victimization (χ 2 (2, n = 26) = 4.5, p = .11), with a median rank of 14 for single-grade classrooms, 17 for administrative multigrade classrooms, and 5.8 for pedagogical multigrade classrooms. Table 2 presents the summary of the ERGMs. Appendix E shows the results per classroom. We present our findings leaving out two classrooms that were identified as outliers (see Appendix E).

| ERGM findings
Removing these classrooms from the meta-analysis, led to more reliable results (referring to less variation across classrooms and results were not more "favorable;" compare Table 2 with Appendix F). The first column in Table 2 shows the mean estimates across administrative multigrade classrooms (as reference), the two other columns show the degree to which regular single-grade and pedagogical multigrade classrooms deviate from this.

| Network effects
The results of the basic network effects were comparable across the three classroom types (see Table 2). Taking the degree of peer victimization among boys into account, the negative density effects across all three classroom types indicate the low occurrence of victim-bully relationships. The positive isolates effects indicate that some students were uninvolved in victim-bully relationships: they were neither targets nor perpetrators. The positive sinks effects show that some students, while nominating none as a perpetrator themselves, were considered as perpetrators by classmates. The inties spread effects and shared-in-ties effects, reveal that there is no clear pattern in the occurrence of variability in receiving and sending perpetration nominations. Finally, there was no evidence for the occurrence of bully-victims patterns. Taken together the network effects implied that all three classroom types were comparable in network structure.

| Sex effects
Taking into account all other effects, the ERGMs show on average no clear sex differences in bullying involvement across the three classroom types (see Table 2). Subtle (nonsignificant) differences in the mean estimates were observed for the girl-boy effects (more positive in single-grade classrooms and pedagogical multigrade classrooms) and girl-girl effect (more negative in single-grade classrooms and pedagogical multigrade classrooms). Overall, the results suggest higher involvement in bullying perpetration by boys RAMBARAN ET AL.

| 565
than girls in both single-grade classrooms and pedagogical multigrade classrooms than in administrative multigrade classrooms.

| Age and grade effects
Social network analyses showed no evidence that peer victimization depended on age differences between children in any of the three classroom types: on average, older children did not receive more perpetration nominations from younger children than vice versa (nonsignificant age-difference effects). In addition, there was no evidence that older children received more perpetrator nominations than younger children (nonsignificant age-receiver effects). We checked whether these nonsignificant effects were a result of low sensitivity (by using a stricter cut-off point for peer victimization, referring to being victimized two or three times a month instead of once or two times; see, Solberg & Olweus, 2003), which was not the case (see Appendix G). We also tested for sex differences, by including an interaction term for the sex combinations, but found no significant interaction effects (see Appendix H). Note that age effects were smaller than grade effects in terms of "effect size" even when assessed on a 1-year basis.
Administrative multigrade classrooms showed a negative low-high effect indicating that students in the lower grade were bullied less by students in the higher grade. In contrast, pedagogical multigrade classrooms showed a positive effect. However, this grade dominance effect was driven only by two classrooms (pedagogical multigrade classrooms 22 and 24; see Appendix E). These findings indicate that victim-bully relationships in administrative multigrade classrooms occurred mainly among same-grade students, whereas in pedagogical multigrade classrooms they sometimes also occurred across grades. Av. density high-low -0.04 (0.00-0.11) 0.02 (0.00-0.03) a = Density is the proportion of observed ties present in the network in relation to the total amount of possible ties available in the network. b = Density for sex or grade is based on the proportion of ties present in each group in relation to the total amount of possible ties available in each group.
Taken together the observations, there was no conclusive evidence for the idea that peer victimization depended on (relative) age differences or grade differences between children in classrooms.

| DISCUSSION
We aimed to address the question whether age-differences between classmates can serve as indicators for power imbalance in bullying Notes: + p ≤ .10, *p ≤ .05, **p ≤ .01, ***p ≤ .001. a Significant differences between classrooms. To help interpret the effects in this table, an illustration is provided for each parameter. Victim-bully relationships are represented with directed arrows (referring to targets nominating perpetrators). Students are represented with colored nodes (sex: boys in blue, girls in pink; grade: lower grade students in white, higher grade students in gray; age: older children in orange, younger children in yellow). The parameter statistics of the network effects used in PNet are mentioned in parentheses (short names).
perpetration (Ellis et al., 2012) by examining whether single-grade and multigrade classrooms differ in peer victimization and whether peer victimization depended on (relative) age or grade differences between children. Drawing on a status framework (Rodkin et al., 2015), in multigrade classrooms we expected to find higher peer victimization, especially for younger children targeted by older children. The competing hypotheses, derived from an evolutionary model (Ellis et al., 2012), would suggest otherwise. Peer victimization was treated as a relationship between targets and perpetrators, and therefore analyzed with social network analysis. Our analyses showed no evidence that single-grade and multigrade classrooms differ in degrees of peer victimization nor a systematic indication that peer victimization is based on (relative) age or grade differences. Thus, providing no support for either the statusframework or the evolutionary model.
The analysis showed instead that peer victimization in administrative multigrade classrooms occurred mostly within the same grade, and thus between children whose age differences were relatively small.
An explanation is that it was relatively easy for older children to target younger children and therefore targeting someone younger as they are did not enhance their status. Peers might be more inclined to award perpetrators with status when perpetrators challenge others of equal (or higher) status in the group. Likewise, younger children in the multigrade classroom were probably not in a strong position to take on older perpetrators, and therefore sought out others within their own respective lower grade. This could also explain why we found bullying among same-age or same-grade peers in administrative multigrade classrooms. Another explanation is that teachers in administrative multigrade classrooms generally teach the grades separately (Mulryan-Kyne, 2007;Veenman, 1995), and this decreased opportunities for bullying between cross-age or cross-grade classmates.
The evolutionary model may be relevant for pedagogical multigrade classrooms. We observed differences in degrees of peer victimization between the two types of multigrade classrooms: lower degrees of peer victimization in pedagogical multigrade classrooms compared with administrative multigrade classrooms.
Findings also revealed no clear sex differences in bullying involvement across the three classroom types, after controlling for network structure, and age/grade. This is surprising given evidence on the differences between girls' and boys' friendship networks during the middle to late childhood developmental period (see Rose & Rudolph, 2006), with important differences found in playmate selection, friendship formation (Maccoby, 1998) to control or harm others is peer perceived popularity (Salmivalli, 2010), or more broadly defined status differences, for instance in the form of SES (Chaux & Castellanos, 2015), majority/minority status or LGBT status (Juvonen & Graham, 2014). Clearly, more precision is necessary for the measurement of power imbalance to understand how power imbalance unfolds in victim-bully relationships.

| Implications of the research
Our findings may have implications for school policy about classroom organization. When transitioning to a multigrade classroom organization, parents worry about their child's progress and well-being. Previous research points in the direction that children in administrative multigrade classrooms are not better or worse off in terms of socio-cognitive outcomes (e.g., academic achievement and social adjustment; Mulryan- Kyne, 2005;Veenman, 1995Veenman, , 1996; but see Mason & Burns, 1996). In this study, we extend previous research by showing that there was also no clear indication that multigrade classrooms increased the risk of victimization by classmates compared with single-grade classrooms.
Our findings indicated that multigrade classrooms may be beneficial for peer victimization when they are part of the school's educational philosophy. This finding is also in line with broader theory and research on the positive effects of pedagogical multigrade classrooms (e.g., in Jenaplan or Montessori schools, see Lillard & Else-Quest, 2006;Moller et al., 2008), which stimulate prosocial behavior among children by encouraging the provision of help across the grades within the classroom (Gray, 2011;Lillard & Else-Quest, 2006).
Although based on a few classrooms, it suggests that schools that encourage prosocial relationships among children by encouraging the provision of help across grades within the same classroom (Gray, 2011;Lillard & Else-Quest, 2006), may have less bullying between classmates. A recent study also showed that schools engaged in practices to promote inclusiveness and equity as a social-learning program may foster openness and acceptance of others (Rivas-Drake, Saleem, Schaefer, Medina, & Jagers, 2019). An avenue for further research could be to examine the extent to which single-grade and (administrative and pedagogical) multigrade classrooms operate similarly or differently for positive peer relationships, such as prosocial or helping relationships in the classroom.

| Strengths, limitations, and directions for future research
An important strength of this study is that we investigated the role of (relative) age and grade differences between classmates in victim-bully relationships. Assessing these associations in a sample of single-grade and multigrade classrooms further add to our understanding of classroom differences in elementary education.
By applying social network analysis we did not violate the assumption of independence of observations within classrooms.
To allow for a comparison between different classrooms, our sample was relatively small because we had to exclude many schools with multigrade classrooms with a different grade combination (referring to more than two grades in the same classroom). It is possible that our null-findings are due to low statistical power (related to the low number of classrooms that were included). Future research should validate these findings using a larger sample, preferably with larger within-classroom grade differences.
It should also be pointed out that schools that have different grades within the same classroom as a consequence of grade/ classroom mixing tended to be unstable in terms of student classroom composition over school years. Hence, an obvious next question would be how this grade/classroom mixing affects children in terms of bullying and peer victimization (Rambaran, van Duijn, Dijkstra, & Veenstra, 2019).
We did not measure victimization by peers outside the classroom because the focus of this study was on the classroom.
However, part of the victimization might occur by other peers in school or outside of school. We only measured peer victimization by asking the students "Who starts bullying you?" We had good reason for this because students who initiate the bullying are often considered as "ringleader" perpetrators (Salmivalli, 2010), and thus, they may be more powerful and influential than followers who join in the bullying.

| CONCLUSION
In sum, our findings did not provide evidence that (relative) age differences or grade differences play a role in peer victimization in multigrade classrooms compared with single-grade classrooms.
Apparently, processes of power imbalance in bullying are not (or only weakly) related to students' (relative) age differences or grade differences within the same classroom.