Abstract
A growing body of literature has documented the contribution of teacher-student relationship quality to both persistence and reduction in peer aggression incidents in the school context. The research literature indicates that students who are involved in peer aggression also tend to experience lower levels of closeness in their relationships with their teachers. However, these study results have not yet been aggregated, and the size and direction of effects remains unclear. In the present study we quantitatively synthesized 66 individual studies (Nstudents = 352,376) in two meta-analyses by aggregating cross-sectional associations between peer aggression involvement and teacher-student relationship closeness that have been reported in the literature over the last 20 years. A small, negative, and significant association was found between perpetration and victimization and teacher-student relationship closeness, indicating that students who experience greater involvement in peer aggression also have relationships with their teachers that are lacking in closeness. Three moderator analyses were also conducted. No moderating effect was found for school level or measure type; however, a significant moderating effect was found for informant type. The results from the meta-analyses lead to direct recommendations for practice regarding how we can best support students’ psychosocial development in the school context.
A growing body of literature indicates the influential role of the teacher-student relationship (TSR) in the social-emotional well-being of children and youth (Lavy & Naama-Ghanayim, 2020; see Sabol & Pianta, 2012 for a review). It appears that the TSR may be especially important in relation to students’ involvement in peer aggression, including bullying. For the present meta-analyses, the term ‘peer aggression involvement’ includes both perpetration (i.e., committing aggressive acts towards peers) and victimization (i.e., receiving aggressive acts from peers). Continuing research on the associations—both positive and negative—between these variables is critical, given that peer aggression remains a pervasive and persistent global public health problem with roughly one in three children reporting that they experienced peer victimization in the previous month (UNESCO, 2019). It is now widely accepted that involvement in peer aggression has many negative effects for children and youth. Several studies have documented links between perpetration and victimization experiences and various mental health issues (e.g., depression, anxiety), poor physical health, and later substance use (Moore et al., 2018). TSRs have emerged as a possible explanation of both the persistence and the mitigation of involvement in peer aggression in school contexts. By synthesizing the available literature that examines the relationship between involvement in peer aggression and TSR closeness, the present meta-analyses offer robust empirical evidence of the quantitative link between these variables, and it points to possible avenues of remediation to the issue of peer aggression and bullying in schools.
Teacher-Student relationships
The Bioecological Model of development posits that development occurs within a set of interconnected social-cultural environments that reciprocally interact with the child (Bronfenbrenner & Morris, 2006). It has commonly been used to explain the adverse psychosocial consequences on children that are associated with involvement in peer aggression as well as how social contexts, including, for example, relationships that form between students and their teachers, can reduce the incidence and the impacts of peer aggression involvement (Bouchard & Smith, 2017; Wang et al., 2015). The most proximate of the model systems is the microsystem, which comprises the direct interactions that children have with important people in their lives, including family, peers, teachers, and so on. At the next level is the mesosystem, which forms when people from different microsystems containing the developing child interact. For the current study, we situate our analysis of TSRs and involvement in peer aggression at the level of the mesosystem, as classrooms provide numerous opportunities for mesosystems comprising teachers and peers to form and to mutually influence one another.
From this bioecological perspective, social referencing theory is a useful explanation of one potential mechanism by which teachers may influence students’ behaviours in their peer group, including peer aggression involvement (Hughes & Chen, 2011). In brief, social reference theory suggests that students use their teachers as a social reference to inform how students should treat their peers, effectively mirroring to a peer the qualities they perceive in that peer's relationship with the teacher. While some longitudinal research suggests bidirectional links between TSR quality and involvement in peer aggression (e.g., Holfeld & Leadbeater, 2017; Serdiouk et al., 2016), it seems that the teacher is the driving influence in the association between the two constructs. (e.g., Hendrickx et al., 2016). This teacher influence seems to be especially true for conflictual interactions between the teacher and student where students are more likely to dislike their peer if they perceive that their teacher dislikes them.
Attachment Theory (Bowlby, 1969) has also been frequently used to describe how TSRs affect levels of peer aggression involvement. Attachment theory describes the psychosocial development of children within the context of caregiving relationships and offers a useful framework to understanding the quality of relationships that can form between students and teachers. According to Attachment Theory, the dynamics of early relations with caregivers—specifically the levels of sensitivity and attunement they direct toward the child—lead the child to develop an internal working model of relationships (including expectations and core beliefs) that roughly divide between models of secure and insecure attachment. This implies, and research offers some supportive evidence, that the quality of TSRs is related, at least in part, to parental attachment history (Sabol & Pianta, 2012). However, the theory also suggests that non-parental relationships, including those that form between teachers and students, can themselves form into attachment relationships that are at least partially independent from parental attachment histories, and thereby contribute to shaping or revising children's internal relational models (Davis, 2003). Revised models that favour secure attachment in the TSR can then serve to guide behaviour that subsequently contributes to the formation of high-quality peer relationships (Gest & Rodkin, 2011).
Closeness
Generally, TSR quality has been measured on three discrete characteristics: closeness, conflict, and dependency (Sabol & Pianta, 2012). For the purposes of this paper, we exclusively focus on aspects of closeness. TSRs that have high levels of closeness consist of adult involvement that conveys warmth, care, and support for the child (Bouchard & Smith, 2017). Close relationships between teachers and students have been associated with positive outcomes, including increased social competence and academic achievement, as well as lower levels of internalizing behaviour (Pianta & Stuhlman, 2004). Several studies have explored the role of perceived closeness between teachers and students in peer aggression involvement. In younger student populations (i.e., grades 3–5), both cross-sectional and longitudinal research shows that relationships with teachers characterized by closeness are positively associated with peer acceptance and negatively associated with relational and physical aggression and peer victimization (Kremer, 2010; Troop-Gordon & Kopp, 2011). Additionally, teacher-ratings of TSR quality have been found to be predictive of peer-rated student aggression levels in the following academic year (Hughes et al., 1999). Specifically, positive teacher ratings of their relationships with their students in grades 2–4 were associated with reduced rates of peer-rated aggression in the subsequent academic year (e.g., grades 4–5), even when controlling for pre-existing aggression levels. Similar findings have been found in older student populations. In a large-scale cross-sectional study including over 18,000 middle-school students from France, findings revealed that schools who reported more positive TSRs overall also reported lower rates of both physical and relational bullying involvement (Richard et al., 2011).
Within the closeness construct is the concept of social support. Supportive social relationships, including those with teachers, can create a sense of stability and predictability in one's environment and foster a sense of belonging and integration, lending to more frequent positive experiences of self and other for the child (Cohen & Wills, 1985). Supportive TSRs have been found to have a negative correlation with involvement in peer aggression. For example, findings from a study with over 10,000 Israeli students in middle school and high school indicated that rates of peer victimization were lower for students who also reported having a supportive relationship with their teacher (Marachi et al., 2007).
However, it is important to note that close TSRs have not always been found to be significantly associated with peer aggression involvement. In a study investigating the role of TSR quality and its association with peer victimization in fourth and fifth grade students, Elledge et al. (2016) found that close TSRs did not predict a reduction in peer victimization when student social preference (i.e., degree to which student is liked and/or disliked by their peers) was controlled. Additionally, a study completed by Pabian and Vandebosch (2016), with over 2000 early adolescent participants (those in middle and early secondary school), did not find any significant linkages between teacher-student bonding (i.e., a sense of support and connection between teacher and student) and bullying at school. There is some research to suggest that a subset of students who aggress against their peers are experienced as likeable and socially competent, lending to more positive relationships with teachers (Vaillancourt et al., 2003). These particular studies raise questions about the strength and direction of the association between close TSRs and involvement in peer aggression among students.
Moderating variables
Research has documented multifaceted trends in the relationship between peer aggression involvement and TSR closeness as a function of students’ developmental age. For example, a review by Berger (2007) outlines several studies that found physical bullying decreases in high school, but more covert forms of bullying (e.g., verbal, social, and cyberbullying) increase with age, particularly during middle school years (Grade 6–9) (Due et al., 2005). Regarding TSR closeness, literature indicates that younger students tend to view their relationship with their teachers as more positive, compared to older students (Havik, 2017; Serdiouk et al., 2016). It is also well established that as children age their relationships with peers take increasing importance in their lives (Brown & Larson, 2009). Additionally, students in upper grades (e.g., Grade 7 and above) often have several teachers with whom they interact throughout the day rather than a single homeroom teacher. This too, may influence the observed association between involvement in peer aggression and TSR quality. It is for these reasons that school level, a proxy for developmental age, was included as a moderator in our analyses.
Multiple informants are often used in measuring peer aggression and TSR quality (Richard et al., 2011; Runions & Shaw, 2013; Troop-Gordon & Kopp, 2011) given the reporting limitations associated with each type of informant. For example, social desirability may influence self-reports of peer aggression involvement (Volk et al., 2017). It is also known that some systematic reporting differences exist between informants when reporting involvement in peer aggression (Demaray et al., 2013; Stockdale et al., 2002), which is often explained in terms of varying opportunities that different informants have to witness actual peer aggression (Rupp et al., 2018). Given these trends, we decided to include informant type as a moderator in our analysis.
As the final moderator, we will assess the influence of type of measure. Many of the measures used to investigate peer aggression involvement and TSR quality are well-established and extensively validated (Kimberlin & Winterstein, 2008). However, many studies nonetheless included newly created measures that have very limited validity and reliability data compared to the established measures. Consequently, we decided to compare effect sizes based on well-established measures versus novel measures of our key variables.
The current study
The objective of our meta-analyses is to quantitatively synthesize the size and direction of the association between TSR closeness and students’ involvement in peer aggression, namely as perpetrators and as victims. We expected that students who perpetrate against their peers or who are victimized will also have relationships with their teachers that are low in closeness. Peer perpetration and victimization were investigated independently. We divided the analyses in this way since the peer aggression literature routinely and reliably divides children into various roles related to peer aggression involvement (Yang & Salmivalli, 2013). Scholarship also documents unique outcomes related to peer aggression involvement for students who are victimized and students who perpetrate against others (Ttofi & Farrington, 2008). Additionally, we examined the impact of moderating variables (e.g., school level, informant type, and measure type) on the association between peer aggression involvement and TSR closeness. We expected that all three moderator analyses would help explain the heterogeneity across studies. Regarding the school level moderator, we anticipated a stronger association between TSR closeness and peer aggression involvement for younger students compared to older students given that younger students are more dependent on adults in their lives and spend greater amounts of the school day with a single teacher (Wang et al., 2015). Given the limited research indicating how informant or measure type may influence the correlation between involvement in peer aggression and TSR closeness, these two analyses are exploratory in nature.
Methodology
Study selection and screening
A systematic literature search was performed in September 2019 to identify all relevant literature. Four online databases were accessed: ERIC, PsycINFO, Education Source, and ProQuest (for theses and dissertations only). A combination of Medical Subject Heading (MeSH) terms and key words were used to reliably access relevant literature. Some examples of MeSH terms include: relational aggression and teacher student relationship. Some examples of key words used include: peer harass* and teacher adj2 student adj2 relation* (see supplemental material for a full description of the systematic literature search strategy). Backward and forward referencing searches were also completed to ensure all possible studies on the topic were identified. At the end of the search process, there were a total of 816 unique research reports identified for this study. These reports were then uploaded onto an online systematic review manager (i.e., Covidence) for screening.
To be included in the meta-analyses, all studies must have included a quantitative and continuous measure of student peer aggression involvement, either perpetration or victimization, from a student (self), teacher, and/or peer informant that was completed in the school context. Peer aggression involvement, in the context of this study, included aspects of bullying and victimization, as well as other forms of peer aggression (e.g., reactive aggression) that do not necessarily fit within the confines of the traditional bullying definition. We included all forms of peer aggression involvement, including physical, verbal, vandalism, and cyber, that occurred in the context of an interpersonal interaction, and we excluded measures of aggression that targeted objects or animals. An important role identified in the peer aggression literature is the bully-victim role, which comprises students who both aggress against others and are victimized themselves. Due to key differences between students who both perpetrate and are victimized and students who only perpetrate against others or are purely victimized, it would be necessary to examine these roles in separate analyses, however, insufficient number of studies that reported on this specific group of students prevented us from doing so. Additionally, all studies must also have included a quantitative and continuous measure of TSR closeness from either the student (self), teacher, and/or peer informant that was completed in the school context. Measures that assessed positive aspects of TSRs were included (e.g., Aceves et al., 2009; Espelage et al., 2014), such as aspects of support, trust, warmth, and mutual liking. Finally, all studies needed to include individual-level statistics, have study participants who were in grades k-12, and include a correlation coefficient effect size or other statistic that could be transformed into one.
The first author was the primary reviewer during the screening process. When it was not clear whether a study fit within the predetermined inclusion/exclusion criteria, a consensus decision was made following a discussion between the two authors. After all studies were screened and assessed according to the inclusion and exclusion criteria, a total of 66 studies remained. See Figure 1 for a description of the selection and screening process.

PRISMA diagram of study selection process.
Study coding process
See Table 1 for a list of all included studies and their characteristics. The first author coded the entire set of 66 studies. A second coder, a graduate student completing a master's degree in education and who had no knowledge of the study purpose, was trained on how to identify and code each study characteristic. The second coder coded 19 studies (>28%) of the final set of studies to ensure inter-rater reliability. Inter-rater agreement was computed by matching the coding of both raters and calculating the number of agreements against the total number of comparisons. Overall inter-rater agreement was excellent at 97% percent.
Characteristics of studies included in the meta-analyses.
Note. PA = Peer aggression; TSRclos = Teacher-student relationship closeness; P = Peer perpetration; V = Peer victimization.
Aggregation of effect sizes
We aggregated across multiple measures of peer aggression involvement and TSR quality that were reported within a single study. This was required for a total of 35 studies. Aggregation was done by summing the various effects from each measure and dividing it by the number of effect sizes/measures, providing an average effect size across measures used. For the longitudinal studies (total of 16 studies), effect sizes were aggregated if they corresponded with timepoints occurring within the same academic year. When the timepoints of longitudinal studies fell in different academic years (e.g., timepoint 1 = Grade 1, timepoint 2 = Grade 2), only the first timepoint was extracted. This was done to minimize the impact of development on involvement in peer aggression and TSR quality as students age (Havik, 2017). Aggregation for these effect sizes was done by summing the total number of effect sizes occurring in the same academic year and dividing it by the number of effect sizes/time points. One intervention study was included in the final set (Elsaesser et al., 2013), and only the cross-sectional data collected prior to the intervention was used to eliminate effect size bias. A small subset of the unique studies reported on the same sample of students. In these cases, when the data were unique (e.g., one study reported on perpetration while the other reported on victimization; see Troop-Gordon & Kopp, 2011, and Troop-Gordon & Kuntz, 2013), both studies with their respective data were included in the final set. If the data in the two studies were exactly the same, only one study's data were included.
Study characteristics
With our set of 66 research reports we calculated a total of 90 independent effect sizes. Of the 66 included studies, 59 (89.4%) were published studies, and seven (10.6%) were unpublished studies (two master's theses and five doctoral dissertations). Publication dates of the set of studies ranged from 1999–2019. A total of 352, 376 student participants are represented in the present meta-analyses. The mean sample size is 5244 with sample sizes ranging from 89 to 121,311 participants (the largest, a US sample; Doty et al., 2017).
The studies in our dataset represent 20 countries from five continents (see Table 2 for a breakdown of participants across geographical regions). The majority of studies (>80%) reported roughly equal proportions of males and females in their study. Approximately 40% of studies included ethnically heterogeneous samples, and just over 20% of the studies reported ethnically homogeneous samples (i.e., indicated by over 80% of one ethnic background). Roughly 37% of the studies did not provide any ethnicity information about study participants. Students of all ages were included in the present meta-analyses from kindergarten to Grade 12, with 40% of studies including elementary school-aged children, and the remaining 60% including middle/intermediate and secondary school-aged children.
Breakdown of included studies and participants by geographical region.
Three quality assessment criteria, extracted from the Quality Assessment Tool for Studies with Diverse Designs (QATSDD; Sirriyeh et al., 2011) were coded to determine the amount and specificity of information for each indicator of reporting quality. After careful deliberation of other quality assessment tools, the QATSDD was used since it contained items that were highly relevant to the diverse but primarily observational and correlational nature of the included studies. No other quality assessment tools (e.g., Cochrane risk of bias tool, Newcastle-Ottawa Scale) were better suited for our specific set of studies. The three quality assessment items that were coded included clear description of research setting (e.g., target population and clear research questions), clear description of data collection procedure (e.g., when and where data was collected), and statistical assessment of reliability and validity for measurement tools. These three items were selected since they queried critical and over-arching reporting standards observed in high quality studies. Fifty-one studies (77.3%) satisfied all three quality assurance criteria, 13 (19.7%) studies satisfied two out of the three criteria, and only two studies (3.3%) met just one of the quality assurance criteria. The quality assessment item that was most often missing from studies was adequate reliability and validity data for the measures used, with 15.2% of studies missing this criterion. We also examined each study for the presence of validity screening (i.e., the identification and removal of mischievous or careless responders). However, only three out of the total 66 studies included some form of validity screening, thus preventing us from conducting any statistical comparison of studies.
Meta-Analytic procedure
The correlation coefficient (r) is the chosen effect size for our meta-analyses. When a correlation matrix was not provided (N = 16), corresponding authors were emailed to request the bivariate correlations for our variables of interest. Of the remaining studies, those that included odds ratio data were transformed using the Comprehensive Meta-Analysis (CMA) software, version 3. Studies that reported a beta value from a regression analysis (N = 7) were transformed using the method recommended by Peterson and Brown (2005). This method of transforming beta values to correlation coefficients has recently received some criticism (Roth et al., 2018). However, this solution is more favourable than excluding the studies altogether (Borenstein et al. 2011). To account for this recent critique, a moderator analysis was conducted to compare standardized data (beta) transformed to correlation coefficients and direct estimates of Pearson r extracted from studies.
All extracted correlation coefficients were transformed using Fisher's Z to approximate a normal sampling distribution (Borenstein et al., 2011). Effect sizes were weighted using the inverse variance based on the standard error from the study (Card, 2011). Both meta-analyses were conducted with the transformed data, and then effect sizes were transformed back to Pearson r coefficients to make results more readily interpretable. Given that we assumed that differences among study effect sizes were due to a combination of sampling error and other factors, we selected a random effects model. The Der Simonian and Laird estimator was used. Two meta-analyses were conducted: (a) perpetration and TSR closeness and (b) victimization and TSR closeness. Aligned with Cohen's (1992) recommendations for qualitatively assessing sample size values, effect sizes are interpreted as small (r = .10-.29), medium (r = .30-.49), or large (r > .50). CMA, version 3 was used for all computations (Borenstein et al., 2013).
Heterogeneity statistics, the Q statistic and the I2 index, are reported with each of the meta-analyses. The Q statistic is a significance test that tests the assumption of homogeneity of effect sizes (Borenstein et al., 2011). A significant Q statistic indicates that the effect sizes will differ between studies and that sample and methodological characteristics (i.e., potential moderators) play a critical role in the accurate interpretation of meta-analytic results. The I2 index allows for the determination of the amount of heterogeneity among the effect sizes (Card, 2011). Rough estimates of the I2 index are: ∼25% is equal to small heterogeneity beyond sampling error, ∼50% is representative of moderate heterogeneity outside of sampling error, and ∼75% + is considered to be large heterogeneity independent of sampling error (Borenstein et al., 2016).
Three moderator analyses were conducted using dichotomized categorical variables to more accurately understand how the association between peer aggression involvement and TSR quality changes as a function of the level of each moderator. For the moderator analyses, a mixed-effects analysis was selected. When moderator levels included fewer than five studies, a fixed-effects model was used to combine effect sizes (Hedges & Vevea, 1998). Each moderator level included a minimum of 3 effect sizes (k ≥ 3; Borenstein et al., 2011).
School level
School level was entered as a moderator to determine its effect on the association between TSR quality and involvement in peer aggression. School level was divided into two ordinal categories: primary (kindergarten – grade 6) and intermediate/secondary (Grade 7 – Grade 12). We made the decision to categorize the moderator in this way since important developmental changes (e.g., decrease in physical aggression but increase in covert peer aggression; Berger, 2007) and contextual shifts occur around the Grade 6–7 mark as students move to larger schools and interact with greater number of teachers. Additionally, the majority of studies included in the final set of studies reported school or grade level rather than student age.
Informant
Informant type was entered as a moderator to determine what effect informant had on the relationship between TSR closeness and peer aggression involvement. Informants were divided into two categories: self and observer (comprised of peer and teacher reports). Peer and teacher reports were combined because of the limited number of studies in each of these informant categories, and because of the large amount of inter-rater agreement between observer informants identified in the research literature for both peer aggression involvement and TSR quality (Cornell & Brockenbrough, 2004; Li et al., 2012). Two separate moderator analyses were run. First, informant effects for studies that used the same informant across both involvement in peer aggression and TSR measures (self/self or observer/observer) were compared. Second, studies that used the same informant for each construct (self/self or observer/observer) were compared against studies that used different informants for each construct (self/observer), assessing for cross-informant effects. Importantly, only primary school-aged children (kindergarten – grade 6) were included in the informant moderator analyses since no studies in our sample that included students in intermediate or secondary school levels utilized peer or teacher reports for measures of peer aggression involvement.
Measure type
The measure type moderator was divided into two categories: established or novel. Measures were characterized as established if they had been previously subject to validation research and authors reported validity evidence of the measure from prior research. Measures were categorized as novel if there was no available validity evidence and the authors stated that the measure had been created for the specific purposes of their study. Studies were included in this moderator analysis only if they used the same quality of measure for both variables (i.e., established/established or novel/novel for both peer aggression involvement and TSR quality). Studies were excluded if they used one type of measure for one construct (e.g., validated measure), but not the other.
Results
Data screening
A sensitivity analysis was conducted to assess whether any one study had undue influence on the total effect sizes. The analysis indicated that no study was significantly manipulating results for the peer perpetration and peer victimization meta-analyses. No studies were identified as outliers. Despite the large range of effect sizes from the included studies, no study was identified as a true outlier, but was considered within the true range of the studies included.
Publication bias
Publication bias was assessed. Egger's linear regression (Egger et al., 1997) assesses whether or not studies were published due to statistically significant results by assessing the amount of asymmetry when studies are plotted by effect size and sample size (Card, 2011). It is assumed that there is no evidence of publication bias when p > .05. For perpetration and TSR closeness, Egger's linear regression was not significant (p = .69) (see Figure 2 for funnel plot). For victimization and TSR closeness, Egger's linear regression produced a non-significant result (p = .74) (see Figure 3 for funnel plot). Taken together, these results suggest that there are no persuasive indicators of publication bias in these meta-analyses.

Publication bias funnel plot for peer perpetration and teacher-student relationship closeness. Note: CES = combined effect size. The imputed data point represents an adjustment of the combined effect size to account for potentially missing studies using the Trim-and-Fill method.

Publication bias funnel plot for peer victimization and teacher-student relationship closeness. Note: CES = combined effect size.
Perpetration and TSR closeness
Figure 4 displays the distribution of all 42 effect sizes and the summary effect (listed last in the figure) corresponding to peer aggression perpetration and TSR closeness. Effect sizes ranged from -.545 to .100. The overall aggregated effect represents a small significant negative association between perpetration and TSRs characterized by warmth and support [r = -.185 (95% CI. -.235 < r < -.133), Z = -6.886, p < .001]. The results indicate that as students are involved at increasing levels of aggressing against others, they tend to have decreasing closeness within their relationships with teachers. This supports our hypothesis that students who aggress against their peers will also experience reduced feelings of closeness, warmth, and support with their teachers. The test for homogeneity of variance was significant revealing that there was true variability between study effect sizes, beyond sampling error [Q (41) = 1811.805, p < .001]. The magnitude of the heterogeneity outside of sampling error was very large, with an I2 index approaching 100% (I2 index = 97.737).

Forest plot for peer perpetration and teacher-student relationship closeness. Note: Q (41) = 1811.805, p < .001; I2 = 97.737.
Victimization and TSR closeness
The second meta-analysis synthesized 48 effect sizes corresponding to the association between peer victimization and TSR closeness, ranging from -.460 to .120. The distribution of these effect sizes, including the summary effect, can be found in Figure 5. Results of the meta-analysis revealed a small significant negative association between victimization and TSRs characterized by closeness and warmth [r = -.157 (95% CI. -.190 < r < -.123), Z = -8.933, p < .001]. Although the effect is small in scale, these results support our hypothesis, indicating that overall, as victimization increases, closeness in TSRs decreases. The test for homogeneity of variance was significant, indicating that there was true variance among effect sizes, more than what could be attributed to sampling error [Q (47) = 3515.449, p < .001]. The I2 index = 98.663, revealing a very large amount of between-study variance.

Forest plot for peer victimization and teacher-student relationship closeness. Note: Q (47) = 3515.449, p < .001; I2 = 98.663.
Moderator analyses
School level
The results for the school level moderator can be found in Table 3. No significant difference was found in the association between TSR closeness and perpetration (Q (1) = .096, p = .756) between students in primary school level and students in intermediate or secondary school level. The effect of school level on the association between TSR closeness and victimization was likewise non-significant (Q (1) = 1.770, p = .183).
All moderator analysis for each meta-analysis.
Note: P = Perpetration, V = Victimization, TSRclos = Teacher-student relationship closeness. * indicates p-value <.05.
Informant
Same Informant
For perpetration and TSR closeness, there was no significant moderating effect for same informant, (Q (1) = 2.508, p = .113) (see Table 3 for results), indicating no difference in effect size for each informant sub-group. However, a significant moderating effect was found for victimization and TSR closeness (Q (1) = 5.800, p = .016), with stronger study effects attributed to self-reporting informants (r = -.145, p < .001, n = 8) versus observer informants (r = -.090, p < 001, n = 4).
Cross Informant
No cross-informant moderator analysis was conducted for peer perpetration and TSR closeness given the insufficient number of studies in the moderator subgroups. For victimization and TSR closeness, there was a significant moderating effect for cross-informant, (Q (1) = 7.211, p = .007) with stronger associations between peer victimization and TSR closeness when the same informant was used (r = -.138, p < .001, n = 12) than when a cross informant approach was applied (r = -.077, p < .001, n = 3). Results can be found in Table 3.
Measure type
No significant moderating effect was found between established and novel measures for peer perpetration and TSR closeness (Q (1) = .001 p = .981) (see Table 3 for results). For peer victimization and TSR closeness, no moderating effect was found between established and novel measures (Q (1) = 2.228, p = .135).
Computations from standardized coefficients
Results of our analysis of error imparted to the data from transforming standardized betas into correlation coefficients can be found in Table 3. For perpetration and TSR closeness, a significant moderating effect was found (Q (1) = 4.208, p = .040), suggesting that studies with direct correlational estimates were associated with stronger effect sizes (r = -.206, p < .001, n = 36) compared to studies with transformed standardized beta coefficients (r = -.061, p = .354, n = 6). No significant moderating effect was found for victimization and TSR closeness indicating no statistically meaningful differences between studies with direct correlational estimates. Results from this moderator analysis should be interpreted in light of the consideration that the effect sizes were smaller for the studies with beta values because the regression models controlled for covariates.
Discussion
The present study quantitatively synthesizes a substantial body of literature published over the last two decades that has investigated links between peer aggression involvement and TSR closeness. Our results, which aggregate findings of 66 independent studies, arguably provide the clearest indication to date of the overall size and direction of the associations between these variables. A key strength of our meta-analyses is the inclusion of a highly varied set of studies. The 66 studies represent a large and diverse group of students in kindergarten-Grade 12: they numbered more than 350,000; ranged from ages four to 18 years; and originated in 20 countries on five continents. Our results show that there are small but significant negative associations between involvement in peer aggression, including perpetration and victimization, and the closeness that teachers and students experience in their relationship. This means that students who are more highly involved in peer aggression, whether in perpetrator or victim roles, also tend to have relationships with their teachers that lack closeness, which we have conceptualized as consisting of warmth, trust, and support. Studies conducted on the topic have not always found significant and negative correlations between students who are involved in peer aggression and warm and supportive TSRs (e.g., Elledge et al., 2016; Pabian & Vandebosch, 2016). The findings of our meta-analyses clearly indicate that such findings are contrary to overall trends in the literature.
It is important to note that our findings may nonetheless obscure more complex associations between peer aggression involvement, in particular, and certain kinds of individual characteristics. For example, other research has shown that a subset of children who perpetrate bullying against their peers but are not themselves victimized tend to be popular in their peer group and seen by peers as possessing leadership and other social competencies (Vaillancourt et al., 2003). Sutton et al. (1999) found that some children who perpetrate aggression show theory-of-mind skills that are superior to children involved in peer aggression in other roles (e.g., followers, victims) and suggested that these skills for understanding others’ cognitions and emotions can be used in strategic and manipulative ways to perpetrate against others. Another subgroup involved in peer aggression, who we were not able to isolate in our analyses, comprise children in the combined perpetrator-victim group. These children exhibit the highest levels of challenging externalizing behaviours (e.g., aggression, hyperactivity, conduct problems; Kelly et al., 2015; Salmivalli & Nieminen, 2002) which serve to erode closeness and amplify conflict over time in their relationships with teachers (Mejia & Hoglund, 2016; Roorda & Koomen, 2021). In light of this evidence, it is possible that a stronger negative link between peer aggression involvement and TSR closeness may exist for children who act in both perpetrator and victim roles compared to students who exclusively act as perpetrators. In contrast, children who are primarily victimized tend to experience higher levels of internalizing problems relative to the other subgroups of involved children (Schoeler et al., 2018), and these challenges appear not to exert the same strains on relationships with teachers. While internalizing problems seem to have some negative relationship impacts, their impacts are not as strong nor as persistent as the effects of externalizing problems on TSRs (Mejia & Hoglund, 2016). The links we draw here between the aggression involvement groups and TSRs remain speculative, particularly because it is not possible to discern in our data any relationship trends that might apply to students who are involved in peer aggression as both perpetrators and victims. In order to maximize our sample size, we chose to synthesize continuous peer aggression involvement data rather than comparing discrete categories of perpetrator and victim roles, for which there was much less data. This is a limitation of our meta-analyses and points to an important avenue of future research to investigate how exactly various peer aggression involvement roles are linked to TSR closeness.
It is useful to consider how this association between involvement in peer aggression and TSR closeness unfolds in the school context. Longitudinal research findings are mixed about the causal pathways connecting these variables. Some research reveals pathways that are primarily unidirectional, beginning with students’ externalizing behaviours, including aggression, and leading to lower TSR quality (Mejia & Hoglund, 2016; Roorda & Koomen, 2021), while other research indicates bidirectional paths whereby higher peer aggression involvement leads to TSRs low in closeness and, simultaneously, TSRs lacking closeness lead to more peer aggression involvement (Holfeld & Leadbeater, 2017; Pabian & Vandebosch, 2016; Serdiouk et al., 2016). In consideration of these latter findings and the relational systems perspective that informs our own work (e.g., see Watson, 2012), the following scenarios illustrate how these reciprocal interpersonal processes might unfold over time. A teacher who witnesses a student routinely aggress against their peers may begin to dislike the student and withdraw their support and care for them. This may lead the student to feel disconnected from and unsupported by their teacher, which in turn may lead the student to perpetrate against their peers more often. Conversely, a student who continues to be victimized by their peers may perceive their teacher as unhelpful, and unable to provide protection and support, ultimately leading to further disconnection and isolation. It is also prudent to consider that the link between peer aggression involvement and TSR closeness may be due to specific student characteristics that increase risk of involvement in peer aggression and also decrease the likelihood of forming close relationships with their teachers. Smokowski and Kopasz (2005) outline several characteristics often identified in victimized students, including poor social skills and emotional maladjustment.
On a theoretical level, these results highlight the influence that activities in one social context (i.e., the TSR) can have on another social context (i.e., peer interactions) in the school setting and point to the important role of the teacher in the classroom beyond teaching content. Overall, findings from the two meta-analyses emphasize the unique and critical role of the TSR embedded in the social landscape of the classroom, providing yet more evidence for the Bioecological Model of development (Bronfenbrenner & Morris, 2006). This study examined the constructs of peer aggression involvement and TSR closeness at the level of the mesosystem. Mesosystem effects emerge when students learn how to treat their peers by looking to their teacher's interactions with classmates, as described by the social reference theory (Hendrickx et al., 2017). Our quantitative evidence of mesosystem effects, which tend to be more difficult to measure then microsystem effects, is an important strength of the present meta-analyses.
Moderator analyses
To provide a more nuanced understanding of the correlation between peer aggression involvement and TSR closeness, three moderators were investigated, the first of these being school level. The results of the analysis revealed that school level (dichotomized into elementary; k-6 and intermediate/secondary; 7–12) did not significantly moderate the relationship between peer aggression involvement and TSR closeness (see Table 3). It appears that despite clear changes in developmental stage and school context, the association between peer aggression involvement and TSR closeness remains stable. That is, the links between involvement in peer aggression and TSR closeness remain just as important in grade 12 as they do in kindergarten, and all the grades in between. One explanation for this non-significant result is that while the correlation coefficients remain stable, the individual data points may still change. It is entirely possible that two things are happening simultaneously; (1) that both peer aggression involvement rates and TSR closeness ratings are lower for older student compared to younger student populations in our sample of studies – which would align with extant literature (Kennedy, 2021; Lynch & Cicchetti, 1997; Wu & Hughes, 2015), and (2) that peer aggression involvement and TSR closeness ratings are decreasing in parallel. Thus, despite our initial hypothesis, the correlation coefficients for older and younger students would share the same amount of variability, resulting in no statistically meaningful differences. We are not able to confirm this explanation, given the limited availability of individual data points for peer aggression involvement and TSR closeness in our set of studies. It would be informative for future research to explore the association between peer aggression involvement and TSR closeness across student age and school context more deeply.
The informant moderator analyses indicated that informant type (self versus observer, observer including peers or teachers for the purposes of this study) significantly influenced results for one of the two meta-analyses. First, self-reports of peer victimization and TSR closeness yielded stronger effect sizes relative to observer reports. Studies have found that self-reports of involvement in aggressing against others and being victimized share limited agreement with reports from teachers and peers (Cornell & Brockenbrough, 2004). Second, a significant moderating effect of cross informant was found for peer victimization and TSR closeness with larger associations for studies that used the same informant for both involvement in peer aggression and TSR closeness compared to studies that employed a cross-informant approach. This analysis provides quantitative evidence that using the same informant for measures of both peer victimization specifically and TSR closeness may inflate effect sizes. Our results from this cross-informant moderator analysis are aligned with existing literature supporting the use of multiple informants (e.g., Achenbach et al., 1987). Bolstered by these findings and given the discrete differences among informants when assessing both peer aggression involvement and TSR closeness, we suggest that several varied informant perspectives be sought so a more comprehensive understanding can be gathered (Shakoor et al., 2011). When the same informant reports on different constructs on similarly structured scales, biases in the informant can influence responses to both measures in similar ways and cause shared method variance. This can be true of self and other (teacher and peer) informants. Soliciting various informant perspectives removes this methodological bias and potentially provides a more accurate estimate of associations between the two constructs. Importantly, our moderator analyses only included primary school aged children since no studies in our sample that included intermediate/secondary school aged children utilized observer-reports for measures of peer aggression involvement. It would be worthwhile for future research to adopt observer-report measures when investigating involvement of peer aggression and TSR quality in older age groups to further assess the function of age in the association between these two constructs.
While many studies included in our analyses used measures with strong validity evidence (e.g., STRS; Pianta, 2001), several other studies included measures that had been specifically designed for the purposes of the single study, or included a measure that combined items from several different measures. This moderator analysis of measure type suggests that using novel measures does not necessarily result in smaller effect sizes. Thus, attempts to use novel measures should not be entirely dissuaded, particularly if it allows researchers to administer a measure that better suits their study requirements. However, what remains most important, and especially so when utilizing novel measures, is that researchers provide psychometric properties of the included measures whenever possible, to ensure that their chosen measure is indeed assessing the desired construct. The consolidation and validation of measures in these domains will also make future quantitative syntheses more feasible and accurate.
Limitations
An important limitation requiring consideration is the significant heterogeneity of the studies included. Notably, this was an important strength of the paper as it allows the results to be generalized to diverse populations across several different countries and continents. However, the large amount of heterogeneity complicates efforts to explain how and why exactly effect sizes differed from one study to the next. It will be important for future research to explicate the subtleties of the relationship between the two constructs. Additionally, future analyses should assess the reduction of heterogeneity after accounting for particular moderators. While these analyses were ultimately beyond the scope of our present meta-analyses, a closer investigation of the impact of moderators on the link between peer aggression involvement and TSR closeness is recommended to better inform prevention and intervention programs for a variety of students.
While our set of studies was diverse, Africa was not represented in the present analysis and the majority of the included studies were North American, and primarily from the United States. This highlights the need for more geographically diverse research to test the generalizability of our findings to children, youth, and teachers in other parts of the globe.
Additionally, given the cross-sectional nature of our analyses, we are not able to make any conclusions regarding the causal ordering of the variables. It would be worthwhile for future research to focus on quantitatively synthesizing the longitudinal associations between peer aggression involvement and TSR closeness to better understand how the two constructs interact across time.
Implications of study
The results of the present meta-analyses have important implications for both theory and practice. Our results emphasize that child development does not occur in a vacuum; rather, development is influenced heavily by the social-cultural systems that surround the individual child, as posited by the Bioecological Model (Bronfenbrenner & Morris, 2006). Our findings also shed light on the role of the mesosystem within a child's context – providing quantitative evidence of the influence arising from dynamic interactions between two separate micro-systems (i.e., peer relationships and TSRs). Taken together, our meta-analyses further strengthen an ecological understanding of development whereby contextual factors surrounding the child are key indicators of understanding developmental trajectories, both positive and negative, at least as it pertains to involvement in peer aggression.
Importantly, because we used a random effects model for our analyses, the findings from our meta-analyses are generalizable to student populations beyond those included in our final set, which is already large in scale and diverse in nature. As such, our results can be used to inform policy and intervention/prevention programs geared towards reducing peer aggression involvement and improving TSRs for many different groups of students. From a practical perspective, the results support the inclusion of additional training for teachers to promote relationships with their students that are characterized by high levels of closeness. Although findings on causal pathways are mixed, it seems plausible viewed in a systems perspective that a shift in one social context, for example in the quality of TSRs, could lead to changes in the other social context, peer relationships (e.g., Meadows, 2008). In this way, prevention and intervention programs designed to improve levels of closeness in TSRs may cause downstream effects that impact peer aggression involvement rates. One of the core established interventions in this regard is the Banking Time intervention (Williford & Pianta, 2020), which adopts a dyadic and bidirectional approach to improving TSRs that requires the active engagement from both student and teacher. The intervention involves one-on-one sessions in which the teacher intentionally follows the student's lead in activities and interactions with the goal of fostering a sense of closeness and trust. Other interventions have embraced an approach whereby the teacher is primarily responsible for improving levels of closeness in the TSR, which is supported by the social reference theory posed by Hendrickx et al. (2017). Jennings and Greenberg (2009) make a compelling case for the influence of the TSR within the larger classroom network. They state that social and emotional competencies of the teacher (i.e., high levels of social awareness, communication of respect and attunement) are integrally linked to the wellbeing of students and to a well-functioning and healthy classroom environment. Using this conceptual foundation, Smith & Whitley (2022) developed a new intervention (called Teaching with Acceptance and Commitment) aimed at developing teachers’ social-emotional competencies in order to support the development of positive TSRs, although efficacy data on the intervention are still lacking. Other research outlines several specific social-emotional competencies that can nurture close TSRs from a social reference and attachment perspective (Bouchard & Smith, 2017; Yu et al., 2018). Some of these include noticing students’ presence and individual needs both in academic and social contexts, an investment and care for students beyond just broad-level teaching (e.g., taking extra time to make sure student understands content) and emotional support for students (e.g., offering a safe space for students to discuss challenges). While the quality of TSRs is no doubt affected by both child and teacher contributions, interventions have focused mainly on the teacher as the individual most responsible for the relationship quality, given the developmental and professional advantages they have as adults in the relationship. This adult-led approach is promoted by large organizations seeking to reduce acts of peer aggression, such as Promoting Relationships and Ending Violence Network, comprising researchers and educational and community organizations committed to improving the development contexts for children and youth (PREVNet, 2019).
Conclusion
The results of the study illustrate the interconnected social landscape of schools, indicating that students who have relationships with their teachers that are positive, warm, and supportive are less likely to be involved in peer aggression incidents at school, either as a victim or a perpetrator. Our findings shed light on the mutual influence that exists between various social contexts in the school setting, and they highlight the often overlooked importance of teachers beyond just teaching class material and managing classroom behaviour. Findings from our meta-analyses speak to a global need for children and youth of all ages to have teachers that care about them, endeavor to support them, and communicate warmth. In return, children and youth may be more likely to extend the same positive relationship experiences to their peers, potentially decreasing peer aggression involvement overall.
Supplemental Material
sj-docx-1-spi-10.1177_01430343221138038 - Supplemental material for The interconnected school context: Meta-analyses of the associations between peer aggression involvement and teacher-student relationship closeness
Supplemental material, sj-docx-1-spi-10.1177_01430343221138038 for The interconnected school context: Meta-analyses of the associations between peer aggression involvement and teacher-student relationship closeness by Amanda Krause and J. David Smith in School Psychology International
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The study was funded by a Social Sciences and Humanities Research Counsel Canadian Graduate Scholarship – Master's.
Ethics Approval
This was a meta-analysis using already published data. No ethics approval was required.
Supplemental material
Supplemental material for this article is available online.
Author biographies
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
