School climate and academic performance of Italian students: the role of disciplinary behaviour and parental involvement

Educational researchers have increasingly recognised the importance of school climate as a malleable factor for improving academic performance. In this perspective, we exploit the data collected by the Italian Institute for the Evaluation of the Education System (INVALSI) to assess the effect of some school climate related factors on academic performance of tenth-grade Italian students. A Multilevel Bayesian Structural Equation Model (MBSEM) is adopted to highlight the effect of some relevant dimensions of school climate (students’ disciplinary behaviour and parents’ involvement) on academic performance and their role on the relationships between student socioeconomic status and achievement. The main findings show that disciplinary behaviour, on the one hand, directly influences the level of competence of the students, and, on the other hand, it partly mediates the effect of socioeconomic background whereas parents’ involvement does not appear to exert any significant effect on students’ performance.


Introduction
In recent years, there has been a remarkable consensus among educators, policymakers and scholars in the identification of school climate as one of the main levers that can be used to improve national school systems. The first studies describing the positive effects of school climate on the learning process date back to the early 1900s (Perry 1908). Later, Halpin and Croft (1963) started the first empirical analysis of school climate, prompting a rich and broad line of research systemically studying the influence of school climate on student learning and development. Nowadays, most of the existing papers regarding this topic, however, refer to the school system of the United States, where the debate on school climate began over 50 years ago. For example, Sherblom et al. (2006) studied the relationship between aspects of school climate, as reported by students, teachers/staff and parents, and standardised tests in reading and mathematics for third and fourth grade students. Brand et al. (2003) analysed a large-scale sample of over 105,000 students in 188 schools and suggested that school climate dimensions were consistently associated with academic, behavioural, and socioemotional adjustment indices. The results by Stewart (2008) suggest that school climate-in particular, the sense of school cohesion felt by students, teachers, and administrators-is a strong predictor of successful student outcomes. Similarly, the analysis by MacNeil et al. (2009) shows that students achieve higher scores on standardised tests in schools with healthy learning environments. Generally, the learning process can be sustained by a positive and respectful atmosphere (Harris and Chrispeels 2006; Kyriakides and Creemers 2008). Moreover, some longitudinal studies (Hoy et al. 1998) have confirmed that a positive school climate not only contributes to immediate student achievement but also to effects that seem to persist for years.
Another line of research includes several analyses based on large international student assessment data from the OECD Programme for International Student Assessment (PISA) and the Trends in International Mathematics and Science Study (TIMSS) carried out by the International Association for the Evaluation of Educational Achievement (IEA). These data provide opportunities to make comparisons between countries using standardised and high quality information. Additionally, they allow the analysis to be extended not only to school climate but also to factors that characterise the learning processes at the student and school level. Note, however, that exploring specific aspects, such as school climate, requires ad hoc surveys and methodological tools that can hardly be carried out on a large scale. Consequently, most studies based on international databases focus on a few specific dimensions of school climate, such as those related to disciplinary aspects, which are the characteristics most investigated by these surveys.
Despite the strong consensus regarding the importance of a positive climate on educational outcomes, some caution must be adopted in analysing and interpreting its effect. It should be emphasised that school climate can be endogenously linked to other factors, in particular those of a socioeconomic nature, which, in turn, affect school outcomes directly. For example, although there are clear signals of a strong and meaningful relationship between school climate and cognitive outcomes, the interpretation of these relationships should account for the school socioeconomic compo-sition, which is simultaneously related to school climate and academic achievement, as shown by several studies (Ma and Willms 2004;OECD 2010;Rangvid 2007).
Consequently, the analysis of school climate effects must use appropriate statistical methods that can account for the complex structures of relationships that characterise educational processes. From this perspective, our primary research objective is to determine whether some school climate-related aspects can influence the academic outcomes of Italian students by accounting for three main dimensions that characterises the concept of school climate using the data available on Italian students. More precisely, we adopt a Multilevel Bayesian Structural Equation Model (MBSEM) approach to measure three latent multidimensional constructs of school climate through several observed indicators and to assess its influence on academic achievement. The main objective is to assess the causal effect of school climate on students' outcomes to fill a gap in empirical analyses which, with a few exceptions (e.g., Benbenishty et al. 2016), highlight only correlational links. To identify the causal effect from non-experimental cross-sectional data, we include the Socio-Economic Status (SES) of the students as a measured common cause influencing both school climate and academic performance. Conditioning on SES helps in reducing the source of spurious correlation between school climate and performance thus addressing the endogeneity issue in the relationship between climate and performance and recovering the correct causal estimate (Antonakis et al. 2010;Pearl 2012). Using data from the Italian Institute for the Evaluation of the Education System (INVALSI 1 ), this paper contributes to the debate on school climate policies and provides evidence on the role of some aspects related to school climate in the Italian educational context, wherein there is scant literature on this issue.
The paper tries to verify whether there is space to extend the role of school managers to implement effective practices aimed at improving student performance by acting on school climate. In this light, the analysis of school climate is relevant, as school climate can be fully counted among malleable factors, i.e., factors that schools can easily manipulate (Voight et al. 2013). Acting on this factor may result in effective policies that do not need broader system reforms, which can be particularly relevant in school systems, such as the Italian school system, where constraints of autonomy and budgets limit the number of policy levers for academic improvement.
In detail, we aim to analyse (i) the influence of relevant dimensions of school climate on academic performance, namely parents' involvement in interactions with teachers (PI), parents' participation in school activities and organisation (PP) and students' disciplinary behaviour (SDB); (ii) the role of these factors as mediators of the effect of socioeconomic background on performance.
Specifically, our research hypotheses are the following: H1: The three investigated dimensions of school climate positively and significantly affect student achievement.
H2: The socio economic status influences these aspects of school climate, which in turn influence academic achievement, thus pointing to an indirect effect (mediation).
1 INVALSI is a public agency that operates under the supervision of the Italian Ministry of Education. Its task is to promote the improvement of educational attainment using national and international evaluations and to contribute to the development and growth of the Italian educational system. A positive and significant indirect effect of socio-economic background through a specific aspect of school climate would imply that one can use this aspect as a tool to improve academic performance. This paper is structured in the following way. The next section explores the multidimensional concept of school climate and it describes in detail the dimensions that can be analysed in the Italian context on the basis of the available data. The methodology and data are described in Sect. 3. Section 4 reports the empirical results while Sect. 5 presents the conclusions.

Relevant dimensions of the school climate: a focus on students' disciplinary behaviour and parental involvement
Despite the broad consensus on the positive effect of the school climate both on academic achievement and healthy behaviour for students (Patton et al. 2006), there is not a universally agreed-upon definition of school climate. The multidimensionality of this concept has led educators and researchers to propose, and use, a wide range of school climate definitions, each focused on specific aspects of this complex construct. Consequently, definitions and measures of school climate have evolved over the past decades. Halpin and Croft (1963) defined school climate in terms of teachers' perceptions of the school's personality, i.e., the teachers' views of school routines, behaviours, and interactions among teachers and administrators. Since this seminal work, the concept of school climate was expanded to include the perceptions of students and other school stakeholders, such as parents and community members (Pyper et al. 1987;Simons-Morton and Crump 2003;Bear et al. 2014;Brand et al. 2008).
Especially in the US, the growing interest in school reforms focused on improving the school climate has led to the development of ad hoc surveys and tools to provide proxy measures of school climate. Since 2004, for example, the California School Climate Survey (CSCS) has been administered to assess school staff perceptions about learning and teaching conditions on many aspects of school climate, such as resource provisions and training, school safety, professional development, student learning environment and student expectations.
Some measures of school climate have been proposed by organizations or groups of researchers; for example, the Comprehensive School Climate Inventory (CSCI) is a scientifically sound and helpful tool suggested by the National School Climate Center (Clifford et al. 2012;Guo et al. 2011) which provides feedback on how students, staff, and parents perceive their school's climate for learning. The National School Climate Center defines the school climate as the quality of school life considering the norms and values, interpersonal relations and social interactions, and organizational process and structures.
Another example of a descriptive measure is the Organizational Climate Index (OCI). This index, proposed by Hoy et al. (2002), analyses four dimensions of school climate: principal leadership, teacher professionalism, achievement press for students to perform academically, and vulnerability to the community. Additionally, Sink and Spencer (2007) suggested the Teacher Version of My Class Inventory, which assesses teachers' perceptions of the classroom climate across five dimensions: (a) overall stu-dent satisfaction with the learning experience, (b) peer relations, (c) difficulty level of classroom materials, (d) student competitiveness, and (e) school counsellor impact on the learning environment.
In conclusion, the concept of school climate has undergone numerous changes over the last few decades, resulting in various definitions affecting how the construct is measured and analysed. Thus, a critical first step for any analysis of school climate starts with defining a construct so that reliable and valid measurement and evaluation can occur (Cohen 2006). In this light, another critical aspect is the availability of surveys, and consequently of data, focused on school climate. Indeed, apart from the U.S. experience, in European and extra-European countries there is a limited availability of data, which does not allow a full description and interpretation of the effect of school climate with consideration of all its many facets. This problem is particularly relevant, especially in Italy, where the school assessment system is relatively young (Castellano and Longobardi, 2020) and many dimensions of educational processes, such as school climate, are not yet adequately investigated.
In this perspective, the availability of data related to the Italian context of this specific topic does not allow to measure and analyse the different dimensions that can be traced back to the multifaceted concept of school climate but only to focus on the role that some relevant aspects of this multidimensional construct play in determining student performance. Specifically, we consider three peculiar dimensions of school climate. The first dimension (PI) concerns the degree of interaction between parents and teachers (e.g. how often parents attend general meeting or parent-teacher conferences), while the second dimension (PP) regards the direct participation of parents in school activities and organization (e.g. whether parents discuss about school funds and/or school structures, contribute to define school programmes or participate in educational activities). The third dimension includes aspects linked to students' disciplinary behaviour (SDB) such as truancy (being late for school and making unexcused absences) and bullying or intimidating behaviour either among students or with respect to teachers.
Although these aspects do not exhaust the wide range of dimensions underlying the concept of school climate, they are certainly relevant because a large literature has analysed their role for improving the student performance.
The strand of literature that investigated the role of parental involvement in determining a positive school climate, and therefore improving the student achievement, is quite broad.
Parental involvement can be interpreted either as home based or school based involvement (Borgonovi and Montt, 2012). The former includes the interactions between parents and their child and it regards helping their child with homework as well as discussing education matters or caring about health and safety of children, etc. The latter refers to the interactions between parents and the school staff by focusing on how often the parents interface with the teachers to follow the school progress of their children, whether they participate in the school governance or in extracurricular/voluntary activities in the school.
In recent years, policy makers have placed great importance on the role of parental involvement as a lever to improve students' school performance and numerous efforts have been taken in many countries to integrate the parent involvement into broader educational initiatives. In particular, in the US experience, parental involvement has been one of the cornerstones of major educational reforms such as President Obama's "Race to the Top" or President Bush's "No Child Left Behind".
Empirical evidence of the relationships between parent involvement and academic achievement is equivocal. Some findings support the claim of a positive correlation between the active participation of parents and their children's academic achievement (Jeynes, 2012;Zellman and Waterman 1998;Burcu and Sungur 2009), other studies highlight the positive effect of parental involvement on social skills (Sheridan et al., 2012), good behaviour (Sheridan et al., 2017) and mental health (Wang and Sheikh-Khalil, 2014). Conversely, many scholars indicates that parent involvement is only modestly associated with student outcomes or has no effect on achievement (e.g. Domina 2005;El Nokali et al., 2010;Epstein 1991;Fan 2001).
Different reasons can explain this apparent inconsistency, for example "…a high level of parental involvement in some school activities, such as volunteering in physical and extracurricular activities, may reflect a lack of school resources" (OECD, 2019) or a high degree of interaction between teachers and parents may occur on the initiative of teachers due to a poor school performance of the child.
As to the student disciplinary behaviour, a series of correlational studies, based on PISA data, have shown that the disciplinary behaviour is positively associated with academic achievement. For example, Shin et al. (2009) suggest that in Japan, Korea and the United States there is a strong correlation between disciplinary climate and mathematics performance. Güzel and Berberoğlu (2005) show that the influence of students' perceptions of classroom disciplinary climate on their reading achievement is positive in Japan, non-significant in Norway, and negative in Brazil. Ma et al. (2013) show that in some Asian countries, disciplinary climate has a positive association with student performance in reading, mathematics, and science literacy. Similarly, Ning et al. (2015) report that the classroom disciplinary climate of schools can explain 11% of the between-school differences in reading achievement over countries.

Data and procedure
In Italy, education is compulsory from 6 to 16 years of age, and it is divided into three main cycles: (i) pre-primary school for children between 3 and 6 years of age; (ii) the first cycle of education, including primary education (5 years) and lower secondary school (3 years); (iii) the second cycle of education (5 years). The second cycle of education includes three types of upper-secondary school: schools with an academic focus that mainly cover humanities and scientific fields (licei), technical schools and vocational schools.
Each year, all the students of different school grades 2 are subjected to external evaluation by the INVALSI through standardised tests in two domains: reading and mathematics. The main aim of these tests is to help the school principal and the teachers evaluate the performance of their classes and students. Additionally, the INVALSI results help the central government understand the general performance of the school system and make intra-and inter-school comparisons in space (at the national, regional, and provincial level) and time. INVALSI not only evaluates the students' skills but also requires pupils to answer a questionnaire (Student questionnaire) about their family background, which increases the information potential of the data. Although the INVALSI tests are exhaustively administered, for each school level a random sample is extracted with a two-stage method: in the first stage the schools are sampled and in the second stage two whole classes for each school selected at the previous stage are sampled. The purpose of this sampling is to guarantee the reliability of the data collected: in the sample classes there is in fact an external observer with the task of guaranteeing the regularity of the administration of the tests and avoiding that the data are biased by the presence of teacher and/or student cheating (Longobardi et al. 2018). On this sample of schools, INVALSI administers a School questionnaire, completed by the school principal, to collect important information related to various aspects of school life, such as management practices, availability of infrastructure and resources and school climate.
As mentioned in previous sections, the aim of our analysis is to explore the relationships between academic performance (in terms of the students' score in reading and mathematics test, our target variables) and the three dimensions of school climate while accounting for a set of exogenous characteristics; in this light our dataset is the result of a merging that involves (i) data from INVALSI school questionnaire filled by school principals, (ii) the results from INVALSI tests assessing mathematical and reading skills, and (iii) data related to socio economic variables collected by the INVALSI student questionnaire.
Specifically, the analysis was conducted considering the data related to a nationally representative sample of 642 upper secondary schools (tenth grade) with 23,673 students for which the variables obtained from the school principals' questionnaire are available.
A general multilevel SEM approach (MSEM, Muthén and Asparouhov 2008;Preacher et al. 2010) is adopted to account for the data hierarchical structure, students (Level-1 units) nested in schools (Level-2 units). This methodology allows to account for both sources of variability (within-school and between-school in our case) and separate between-school and within-school effects, without introducing bias.
The general MSEM model splits each observed Level-1 variable into two latent (within and between) components (Preacher et al. 2010). In our case, every studentlevel variable is treated as jointly caused by within-and between-school variation: this means that the observations referring to different students in the same school have to be considered as multiple indicators of the latent school-level construct. In addition, at every level multiple indicators can be used to reflect an underlying unob-served construct as in the measurement model of a common SEM approach. In these terms, the general MSEM framework appear as a doubly latent model (Marsh et al. 2012).

Measures
Consistently with the analysed literature, and on the basis of the items of the school questionnaire, we consider three aspects of school climate representing three distinct latent dimensions. The first two refer to the parental involvement: the degree of interaction between parents and teachers (PI); the participation of parents in school activities and organization (PP) while the third dimension is the students' disciplinary behaviour (SDB) 3 . The PI dimension (interactions with teachers) is measured by the school principal's answers to the following items: (i) participation in parent meetings; (ii) voting for the election of school council; (iii) participation in parent-teacher conferences; (iv) meeting the teachers about the progress and behaviours of children. The PP dimension (participation in school activities and organization) is measured by the school principal's answers to the following items: (i) volunteering in building maintenance activities or other outdoor spaces activities; (ii) discussing the academic achievement of all the students; (iii) discussing school fund spending; (iv) discussing the conditions of school structures and buildings; (v) contributing money to good school performance; (vi) volunteering for good school performance; (vii) contributing to the definition of school programmes; (viii) participating actively in educational activities; (ix) discussing the results of the INVALSI assessment tests. These 13 questions were answered on a 4-point Likert-type response option scale (1 = not at all; 2 = little; 3 = to some extent; 4 = a lot). Higher values suggest higher participation of parents in the various aspects of their children's school life.
For student disciplinary behaviour, the responses of the school principal concern the following items: (i) student lateness; (ii) student truancy; (iii) disruptive behaviours in the classroom; (iv) students bullying other students; (v) students intimidating or offending other students; (vi) aggressive behaviours towards teachers or school staff; (vii) students have discontinuous frequency; (viii) students commit vandalism. These items are also measured on a 4-point Likert-type response option scale (1 = not at all; 2 = little; 3 = to some extent; 4 = a lot), where higher values denote a greater frequency of negative student behaviours.
The frequency distributions of the indicators of the three latent constructs are displayed in Table 1. According to the school principal, parents are most involved in meeting the teachers and least involved in contributing to the definition of school programmes and in school building maintenance activities. The most recurrent negative students' behaviours are vandalism and bullying, whereas arriving late is deemed to be the least frequent.  Table 2 includes the descriptive statistics of the target variable (test scores in reading and maths), the student socio-economic status (SES 4 ) and a set of exogenous covariates that typically characterise the Education Production Function (EPF) models (Niederle and Vesterlund 2010;Buchmann and Parrado 2006), namely, gender, grade retention (regular or not) and immigrant status. In addition, four geographical areas 5 are considered to account for the territorial literacy divide (Quintano et al. 2012), i.e. the wide territorial difference in achievement between Northern and Southern students; this gap is confirmed by several studies which show that students at schools in Northern Italy tend to obtain higher achievement scores than their counterparts in other regions, all else being equal (Giambona and Porcu 2018).

Data analysis plan
The first step of our analysis is to test whether each dimension of school climate (PI, PP, SDB 6 ) is accurately measured by the corresponding indicators. In this light, we perform a correlated traits confirmatory factor analysis (CFA) with three latent continuous factors, measured by 21 dependent items. For the correlated traits CFA model, each item is specified to load onto only the factor it is designed to measure. The correlations among the three factors are freely estimated.
The second step consists in the estimation of the multilevel SEM with the aim of addressing our research hypotheses (Fig. 1). Socio-Economic Status (SES) and test 4 The index of Socio Economic Status (SES), provided by INVALSI, is used. This index takes into account parents' occupations and education, along with variables that measure home possession goods (see Campodifiori et al. 2010 for details). 5 The geographical areas are defined as follows: North-West (including Liguria, Lombardia, Piemonte and Valle d'Aosta), North-East (including Emilia Romagna, Friuli-Venezia Giulia, Trentino-Alto Adige and Veneto), Center (including Lazio, Marche, Toscana, and Umbria), South (including Abruzzo, Campania, Molise, Puglia, Basilicata, Calabria, Sardegna and Sicilia). 6 The values of the indicators D16a to D16h have been reversed in such a way that higher values of the transformed indicators (D16ai to D16hi) are associated with a more positive student disciplinary behaviour. scores (Maths Achievement, MA, and Reading Achievement, RA) are student-level (or level-1) variables. It is assumed that they have within-and between-school components. Therefore, they are split into latent within and between components. This means that, in addition to the effect of SES on individual scores within schools, at the school level (or level-2) the mean SES has an effect on the mean scores. Both scores have been included in the model as endogenous variables, and they are allowed to be correlated. At the school level, the effect of SES on the scores is mediated by Parent Involvement (PI), Parent Participation (PP) and Student Disciplinary Behaviour (SDB) that are school-level latent variables with multiple categorical indicators. This is a mediation model with a 1-(2-2-2)-1 design according to Preacher et al. (2010).
To model Level-2 categorical variables, we use a Bayesian estimator with noninformative priors and a Markov chain Monte Carlo (MCMC) algorithm based on the Gibbs sampler. In the case of the estimation of complex multilevel SEMs with categorical indicators, the Bayes estimator provides a valid full-information alternative to the classical weighted least squares mean and variance adjusted (WLSMV) estimation procedures (Asparouhov & Muthén 2010;Muthén & Asparouhov 2012).
Every structural equation includes the above-mentioned control covariates (not shown in Fig. 1).

CFA results
We assess the fit of the model through different indices, each measuring the fit in a different way. The chi-square test is significant (χ 2 = 648 with 186 degrees of freedom, p < 0.000), which indicates that the model does not fit the data. However, the use of the chi-square statistic raises several issues in such a study, where the variables are not normally distributed and the large sample size is likely to lead to rejecting a good-fitting model. Together with the chi-square statistic, we comment the Standard- The value of these indices should be compared with a fixed cutoff, derived by Hu and Bentler (1999), that discriminates between a good and a bad model fit. Nevertheless, the fixed cutoffs have to be used with caution as many studies have stressed their lack of generalizability due to the fact that they may depend on some characteristics of the specific study, such as sample size, number of items, number of factors and factor reliability (McNeish and Wolf 2020).
The SRMR value (0.075) falls within the 0.08 threshold, which signals an acceptable fit whereas the RMSEA value (0.062) is a bit over the 0.06 suggested threshold for a good fit. The CFI value is equal to 0.923, with a value of 0.90 or greater suggesting an acceptable fit (Hox and Bechger 1999;Hu and Bentler 1999;Schermelleh-Engel et al. 2003).
As a whole, all the indices seem to point to an adequate model fit.
The standardised loadings (Table 3) are all significant and range from 0.474 for D15_i item ("Contribute money to the good school performance") to 0.850 for D16_ CI item ("Disruptive behaviours in the classroom"). Further, all three latent factors are positively and significantly correlated (p -value < 0.000).
We also calculated the omega coefficient (McDonald, 1999) to measure the internal consistency of the latent constructs (Dunn, et al. 2014;Bandalos 2018). We computed the nonlinear SEM reliability coefficient for ordered categorical items (Green and Yang 2009) 7 . The value of omega index is 0.724 for the PI dimension (95% confidence interval: 0.675-0.768), 0.879 for PP (95% confidence interval: 0.851-0.900), and 0.900 (95% confidence interval: 0.878-0917) for SDB. The lowest reliability (0.724) may cause some concern, considering that more than 3/4 of the variability in the score on the latent factor PI can be attributable to random measurement errors. Nevertheless, values higher than 0.70 are considered to be acceptable by most standards.

Multilevel SEM results
The multilevel SEM with a random intercept was estimated. The results are illustrated in Fig. 2. Solid lines reflect significant paths whereas dashed lines represent non-significant paths. The reported coefficient estimates are in standardised units.
At the student level, the path from Socio-Economic Status to achievement is small and significantly positive, after controlling for gender, migration background and regularity. The standardised coefficient is 0.063 (p-values < 0.000) for reading and 0.057 (p-values < 0.000) for mathematics and the model accounts for 2.9% and 2.5% of the variance of the student level achievement for reading and mathematics, respectively.
At the school level, the model accounts for 25.6%, 9.9% and 23.3% of the variance of the latent constructs (PI, PP and SDB, respectively), whereas it explains 76.8% of Notes: the unstandardized loading of the first indicator of each factor is fixed at 1 to define the metric of the factor the variance of the school-level achievement for reading and 72.4% for mathematics. We found partial support for H1 hypothesis on the direct effects of the latent constructs (PI, PP and SDB) on achievement. Only the path from SDB to reading and mathematics achievement is positive (the standardised coefficient for reading is equal to 0.158, for mathematics is 0.113) and significant (p-value < 0.000), which reveals that a better disciplinary climate in a school is significantly associated with a better school-level performance. Conversely, neither PI nor PP exerts significant effects on students performance, implying that parental involvement does not influence schoollevel outcomes.
H2 hypothesis requires to investigate whether, and to what degree, school-level variability in the latent constructs serves as a mediator of the school-level component of SES on school-level achievement. We found a significant simultaneous role of the multiple mediators in the relationship between SES and achievement. Indeed, in the face of an expected significantly positive direct effect of SES (the standardised path coefficient is equal to 0.580 for reading and 0.400 for mathematics), the results point to a total indirect positive effect through PI, PP and SDB (p-value = 0.001 for reading and p-value = 0.002 for mathematics). The total indirect effect arises from the positive and significant paths from SES to PI, PP and SDB (the corresponding standardised coefficients are 0.464, 0.256 and 0.383, respectively): it represents approximately Reading Achievement Solid lines refer to significant paths/relations, whereas dashed lines refer to non-significant paths/relations Only significant coefficients are reported 10% of the total effect. Nevertheless, we found that SDB only has a significant specific indirect effect (p-value = 0.000).

Conclusions
Empirical analyses, since the work by Coleman et al. (1966), have consistently demonstrated that wealthier students tend to gain better grades and, more generally, that the students' socioeconomic characteristics are strong predictors of academic success.
Since the background factors (e.g., family poverty or immigrant status) are "given" and cannot be controlled by the schools, several studies have identified malleable policy factors to promote academic achievement. Among these malleable school factors, educational researchers have increasingly recognised the importance of school climate on students' development . Indeed, "a large body of evidence connects a positive school climate to improvements in children's learning and in healthy development in schools. A positive school climate is also an essential component within comprehensive school improvement processes." (Berkowitz et al. 2017).
Although, in many countries, especially in the US, policymakers and educational leaders have invested large amounts of resources in developing services to improve school climate, in Europe and, particularly, in Italy, attention to this topic is rather rare. Therefore, we exploit the data from the National Institute for the Evaluation of the Education System (INVALSI) to assess the effect of some relevant aspects of school climate on academic performance of Italian students. Specifically, we analyse three latent constructs of school climate that account for the disciplinary behaviour of the students and two different aspects of parental involvement: the parents-teachers interactions and the parent participation to school activities and organization.
The findings from the estimation of a Multilevel Bayesian Structural Equation Model (MBSEM) suggest that policies aiming to prevent or reduce bullying and truancy behaviours among the students have direct implications for educational achievement; schools with better disciplinary climates score significantly higher when other characteristics (such as regional differences and differences in the composition by gender and migrant status) are considered. Furthermore, programmes that effectively cope with disruptive, intimidating and violent behaviours in the school can also impact educational equity; our findings show that the direct effect of socioeconomic status on school performance slightly decreases due to the mediation of school climate, which occurs, for the most part, through disciplinary behaviour.
These results show that the disciplinary dimension of school climate can be considered an important catalyst for academic success of Italian students.
In Italy, the debate on this topic is still at its initial stages; thus, investing effort and resources to provide focused data and empirical studies is imperative for policy development and school improvement even considering that climate-based policies do not require extensive and difficult structural reforms.
From this perspective, Ansary et al. (2015) describe and discuss some interesting examples (to be adapted and remodelled according to the characteristics of Italian school systems) adopted in Finland (KiVa project), Spain (SAVE Model), the United Kingdom (DFE Sheffield Anti Bullying Project) and Norway (Olweus Bullying Prevention Program).
Regarding the other investigated components of school climate (parent involvement), it seems that more involvement of parents does not have a significant effect on students' achievement. Indeed, the direct effects on students' performance both of the level of interaction between parents and teachers and of the participation of parents in school activities are statistically not significant. At the same time, their role as mediators of the SES level is also not significant. Although this result is consistent with several studies that assign a limited role to the involvement of parents (Domina 2005;Epstein 1991), nevertheless it must be interpreted in the light of the available data that do not allow to investigate in detail this dimension of the school climate. Indeed, we have no information on the extent of parental support at home (homebased involvement) but, at the same time, we are not able to measure and analyse the quality of parental involvement which certainly assumes greater importance than the simple quantity of it as suggested by Pomerantz et al., 2007.