Front. Psychol. Frontiers in Psychology Front. Psychol. 1664-1078 Frontiers Media S.A. 10.3389/fpsyg.2024.1354072 Psychology Original Research The effect of school size and class size on school preparedness Antoniou Faye 1 * Alghamdi Mohammed H. 2 Kawai Kosuke 3 1Department of Educational Studies, National and Kapodistrian University of Athens, Athens, Greece 2Department of Self-Development Skills, King Saud University, Riyadh, Saudi Arabia 3David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, United States

Edited by: Iasonas Lamprianou, University of Cyprus, Cyprus

Reviewed by: Dimitrios Stamovlasis, Aristotle University of Thessaloniki, Greece

Aikaterini Vasiou, University of Crete, Greece

 *Correspondence: Faye Antoniou, fayeantoniou@eds.uoa.gr 26 02 2024 2024 15 1354072 11 12 2023 11 01 2024 Copyright © 2024 Antoniou, Alghamdi and Kawai. 2024 Antoniou, Alghamdi and Kawai

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

The purpose of the present study was to understand students’ school readiness as a function of student and teacher behaviors but also school size and class size using both linear and non-linear analytical approaches. Data came from 21,903 schools distributed across 80 countries as per the 2018 cohort of the PISA database. Results pointed to a preference for the Cusp model in that the relationship between school and class sizes with achievement proved to be best described by the non-linearity of the Cusp catastrophe model. The critical benchmarks were a school size of 801 students and a class size of 27 students for which increases beyond those thresholds were linked to nonlinearity and unpredictability in school readiness. For this reason, we suggest using the cusp catastrophe model from Nonlinear Dynamical Systems Theory (NDST) to understand more fully such complex phenomena.

 cusp catastrophe model NLDPS school readiness PISA 2018 school size class size section-at-acceptance Quantitative Psychology and Measurement The effect of school size and class size on school preparedness

It’s important to note that a school’s size and course offerings greatly affect academic performance. A smaller class size allows the instructor to deliver a more tailored education, which increases the likelihood of meeting each student’s requirements and concerns. This may create a learning environment where students feel recognized and understood, which will boost their engagement and drive (Ajami and Akinyele, 2014). Teachers in smaller courses are also better at using a variety of instructional strategies, accommodating different learning modalities, and creating dynamic and interesting learning environments (Bradley and Taylor, 1998). Since there are fewer students to watch in a classroom, keeping discipline is easier, which may reduce disturbances and improve learning. This is shown by greater results on standardized examinations as well as more positive long-term educational outcomes (Mosteller, 1995; Tseng, 2010; Krassel and Heinesen, 2014; Gereshenson and Langbein, 2015; Lowenthal et al., 2019). According to Blatchford et al. (2011), Nandrup (2016), and Yamamori et al. (2021), the number of pupils in a class has a considerable impact on both the educational experience students have and the academic results they attain.

Class size and the ratio of staff to students are often used as measures to evaluate the quality of higher education (Thom, 1975; Molenaar and Oppenheimer, 1985; Brown, 1995; Bandiera et al., 2010; Martin 2015; Konstantopoulos and Shen, 2023). These studies were conducted by Stewart and Peregoy, 1983; Guastello, 1984, 1987,1992; Bandiera et al. (2010); Konstantopoulos and Shen (2023), and Martin (2015). Some studies suggest that larger class sizes hurt student learning; however, a significant number of studies present findings that are inconclusive or demonstrate a combination of positive and negative effects (Bellante, 1972; Edgell, 1981; Hancock, 1996; Kennedy and Siegfried, 1997; Hill, 1998; Jarvis, 2007; Gleason, 2012; De Paola et al., 2013; Matta et al., 2015; Olson et al., 2011). However, previous research has shown that students tend to have a more positive perception of their learning experience when the number of classes they are required to attend is decreased (Kwan, 1999; Van der Maas and Molennar, 1992; Bedard and Kuhn, 2008; Westerlund, 2008; Mandel and Sussmuth, 2011; Monks and Schmidt, 2011; Benton and Cashin, 2012; Sapelli and Illanes, 2016; Hufford et al., 2003; Stamovlasis, 2006). The use of active learning strategies by teachers in smaller class sizes, in addition to the provision of more individualized attention to students (Lammers and Murphy, 2002; Arias and Walker, 2004; Kokkelenberg et al., 2006; Monks and Schmidt, 2011; Van der Maas et al., 2003), could be the reason for this phenomenon (Lammers and Murphy, 2002; Arias and Walker, 2004; Kokkelenberg et al., 2006; Monks and Schmidt, 2011). However, the amount of material that is currently available about the challenges that are related to the application of active learning approaches in smaller class sizes is quite limited (Wright et al., 2017).

A large, influential study namely, the Student Teacher Achievement Ratio (STAR) study reported significant benefits from class size reductions on students’ achievement, if these reductions take place early with the effects being more pronounced for students from disadvantaged family backgrounds (Word et al., 1990; Mosteller, 1995; Finn and Achilles, 1999; Krueger, 1999; Nye et al., 2000). These findings are backed up by several academic publications, such as those written by Guastello (2001), Word et al. (1990), Mosteller (1995), Finn and Achilles (1999), Krueger (1999), and, Nye et al. (2000), amongst others. Throughout the early years of their schooling, policymakers need to reflect on the most effective way to direct Corporate Social Responsibility (CSR) programs toward children who are starting in life with socioeconomic disadvantages. Policymakers also have the option of choosing to offer funding for CSR programs while at the same time providing local school leaders the liberty to decide how such programs will be implemented. It is essential to take into consideration a cost–benefit analysis of educational policy whenever one is charged with making judgments about the maximum number of students allowed in a given classroom.

Evaluating the type of relationship between school size, class size, and school outcomes

Past studies have primarily engaged linear modeling to evaluate the role of class size on school outcomes. The idea that the relationship between class size and student, teacher, and school outcomes is linear falls short for the following reasons. Linear models operate under the assumption that the relationship between variables is best depicted by a mathematical straight line. However, human behavior and outcomes, such as learning and teaching, are inherently complex and multi-faceted. A simple linear relationship might not capture all the nuances and intricacies involved.

Empirical evidence has reported both linear and non-linear effects using, for example, quadratic models. The problem with those findings is that they reported both positive and negative nonlinear trajectories that contradict each other (e.g., Foreman-Peck and Foreman-Peck, 2006; Crispin, 2016). It is possible, however, that the impact of class size on outcomes changes when a certain threshold is crossed. For example, reducing a class from 40 to 30 students might have a more significant impact than reducing it from 30 to 20. Such an effect was reported in the Lee and Smith (1997) study as they reported that school sizes ranging in number of students from 600 and 900 were optimal in facilitating reading and math outcomes. Linear models would most likely be ineffective in capturing such non-linear effects. For the above, more elaborate, and complex analytical models that take into account non-linear changes in behavior and operate using multiple predictors who may exert linear and nonlinear effects are needed. Empirical evidence of this effect has been provided by Cobb et al., 1983; Cobb and Zacks, 1985; Bowne et al. (2017) who reported that “both class size and child-teacher ratio showed nonlinear relationships with cognitive and achievement effect sizes” (p. 407) with implications that this relationship may be further complicated by SES associations. In a large-scale study by Cobb, 1981; Lee and Loeb (2000), concerns were raised about the hypothesized relationship between school size and student learning. In the present study, we propose that the relationship between class size and school outcomes becomes non-linear following a crucial threshold in class size, beyond which, student and school outcomes become unpredictable and chaotic (Oliva et al., 1987). This is why the cusp model may provide the most appropriate means to evaluate the proposed relationships.

 Nonlinear dynamical systems theory and the cusp catastrophe

The cusp catastrophe model is essential in nonlinear dynamic systems theory for explaining sudden and discontinuous system state changes generated by continuous variables. The model works across fields with more recent applications in psychology, education, medicine, and public health. It efficiently handles complex linear and nonlinear interactions between independent variables (Chen and Chen, 2017) toward the understanding of behavioral phenomena. As control variables change, the cusp model can decipher abrupt and sudden behavioral changes offering insights into the stability and transitioning of human behavior (Chen et al., 2014). Applications in the field of education include the investigation of cognitive functioning (Stamovlasis, 2011; Tsitsipis et al., 2012), the teaching of physics (Papageorgiou et al., 2010), and chemistry (Stamovlasis et al., 2005). In the examination of stress and trauma the cusp model has contributed significantly to our understanding of changes in human emotions as individuals transition from one psychological state to another (Kira et al., 2019, 2020). Thus, the cusp catastrophe model has recently been popularized in the social sciences to examine scenarios where shuttle changes in a parameter are associated with drastic and dramatic changes in outcome variables The model employs a potential function, f (y; a, b) for a single dependent variable y given linear and nonlinear parameters a and b:

f (y; a, b) = ay + 1/2by2-1/4y4 (Eq. 1).

The phenomenon under study is termed the “catastrophe set” which evaluates outcomes inside the parameter space of coordinates (a, b). As shown in Figure 1, when the “b” terms (bifurcation variables), in our instance school size and class size are at low levels changes in school readiness are expected to be linear and smooth, likely fitting the premises of the linear model (see expected Pattern A to linearity). When, however, increases in school size and/or class size exceed a worrisome, critical level, school readiness oscillates between two behavioral modes, low and high readiness, reaching a state of unpredictability (see Pattern B to non-linearity). Point B in the figure is termed “the cusp point” and reflects unpredictable changes in the outcome variable following midpoint levels in the bifurcation variable. In other words when schools and class sizes are small and move toward medium levels, school readiness is expected to covary in a linear manner; this linear prediction is disturbed following some critical levels of both school size and class size suggesting that increases in these variables are no longer adaptive.

Description of the cusp model with the outcome variable school readiness; it was predicted linearly using the asymmetry terms of student and teacher behaviors and non-linearly using the bifurcation terms of school size and class size (Tasker, 2003).

The purpose of the present study was to understand students’ school readiness as a function of student and teacher behaviors but also school size and class size using both linear and non-linear data analysis procedures to understand more fully such complex relations.

 Method Participants and procedures

Data came from 21,903 schools distributed across 80 countries as per the 2018 cohort of the PISA database (OECD, 2019). The unit of analysis was the school. Thus, student and teacher estimates were aggregated per school. Student participants in PISA 2018 had to be in the range of 15 years 3 months and 16 years and 2 months and had to be in grade 7 or above. Procedures regarding ethics and sample selection are described here (https://www.oecd.org/pisa/publications/pisa-2018-results.htm). Data may be accessed at https://www.oecd.org/pisa/data/.

 Measures

All measures were derived from the PISA 2018 most recent cohort.

School readiness

This scale is comprised of 8 items completed by school principals on factors that hinder a school’s capacity to provide instruction. The items use the stem: “Is your school’s capacity to provide instruction hindered by any of the following issues” with the content relating to the shortage of teaching staff, insufficient instructional materials, inadequate physical infrastructure, poor lab equipment, and shortage of ICT resources (Werblow and Duesbery, 2009; Wobmann and West, 2006). Scaling ranged between “not at all” and “a lot” (see Supplementary Appendix A1). Given the high internal consistency reliability of the items with omega being at 0.839, factor scores were estimated using maximum likelihood.

 Student behavior hindering learning

It was assessed using 3 student-reported items that evaluate how often students disrupt lessons (a) with noise and talking, (b) with misbehavior, and (c) by being late or absent (Alexander et al., 1992). The items are scored using a 4-point scaling system anchored between “never” and “almost every day.” Scores were estimated using Weighted Likelihood Estimation (WLE).

 Teacher behavior hindering learning

A 4-item scale was created using student responses to assess how often the teacher (a) explains things in a way that is difficult to understand, (b) does not give enough help when you need it, (c) does not seem to care about whether or not you learn, and (d) does not keep order in the classroom. Items engage a 4-point scaling system from “never” to “almost every day.” Estimated factor scores utilized the WLE estimator.

 Class size and school size

Class size reflected the number of students in the class using a categorically ordered variable with intervals of 5 students. The categories ranged between “fewer than 15 students” to “more than 50 students.” School size was measured as a summative variable expressing the number of students in the school. It is a measure of total enrollment rather than expressing estimates for particular grades, cohorts, genders, or else.

 Statistical data analyses Cusp catastrophe

At present, there exist multiple analytical models that can be utilized for the detection of a cusp catastrophe. Among these models, Cobb’s (1998) methodology and its implementation in R using Grasman’s cusp package (Grasman et al., 2009) are widely recognized and accepted. Additionally, the polynomial regression model proposed by Guastello (2002) and the modifications made by Chen et al. (2020) to Cobb’s method are also noteworthy alternatives although limited by the unavailability of routines in statistical packages. While both Cobb’s and Guastello’s models have been widely used, we consider Cobb’s (1981) method as more closely associated with catastrophe theory whereas Guastello’s methodology is more general and includes various forms of non-linear regression models, such as the quadratic. Thus, we choose the methodology proposed by Cobb (1998), which gained popularity through its implementation in R. To achieve optimality, several conditions must be satisfied by the model. Firstly, the asymmetry and bifurcation parameters should exhibit a significant effect. Secondly, the cusp model should demonstrate superiority over linear and non-linear competing models such as the logistic, the quadratic, and the cusp. Third, a relatively small proportion of observations, approximately 10%, should fall within the bifurcation area. Lastly, there should be evidence indicating the presence of bimodality within the bifurcation area and multimodality elsewhere for the outcome variable. We deferred using the pseudo-R-squared statistic provided by the package given that it can take negative values. Of note here is the sign of the bifurcation term(s) which requires additional elaboration. Assuming that the bifurcation term is scaled so that higher scores are indicative of unpredictability, then a positive coefficient should be observed. A positive coefficient indicates that the relationship between bifurcation and the outcome variable is linear at low levels of the bifurcation term and becomes non-linear later on. Based on that, in the present study, a positive slope in the bifurcation term was desirable. On the contrary, a negative bifurcation term suggests that at low levels of the splitting factor/bifurcation term, the system is chaotic, and as the scores on the splitting factor increase, linearity, and equilibrium are gradually present. This is not the case in the present study for which unpredictability is expected when increases in school and class size move beyond some adaptive level.

 Results Prerequisite assumptions of cusp catastrophe

One of the important assumptions of the cusp catastrophe is that the dependent variable must have more than one mode. In the present study, the latent school readiness variable presented itself with multimodality as shown in Figure 2 satisfying the prerequisites of the cusp model. Further evidence of the multimodality of school readiness is shown in Figure 3 using the mode tree (Minnotte and Scott, 1993). The figure displays the mode locations of the readiness variable for each bandwidth. The figure displays 25 modes suggesting multimodality as does Figure 2.

Plotting the multimodal distribution of the dependent variable school readiness.

Mode locations for each bandwidth of the school readiness factor scores construct. The number of modes reflects those shown in Figure 2.

 Predicting school readiness from student and teacher behaviors linearly and from school and class size nonlinearly

Intercept and slope terms of the cusp model are shown in Table 1 with all terms being significant. School readiness was positively predicted by the linear contribution of student and teacher behaviors (bStudent = 0.196, p < 0.001; bTeacher = 0.264, p < 0.001). Interestingly, both bifurcation terms were also significant signaling non-linearity. Specifically, as school size and class size increase beyond a specific critical threshold, their relationship to school readiness becomes chaotic and unpredictable (bSchool size = 0.001, p < 0.001; bClass size = 0.004, p = 0.027). The critical benchmarks were a school size of 801 students and a class size of 27 students for which increases beyond those thresholds were linked to non-linearity and unpredictability in school readiness. When testing optimal model fit between linear and non-linear models (see Table 2) results pointed to the superiority of the cusp model using the information criteria values of AIC, AICc, and BIC over the competing models (a) linear, (b) quadratic, and (c) logistic, which were consistently lower in the cusp model compared to all other models. Furthermore, a chi-square test contrasted linear and cusp models pointed to a significant misfit of the linear model [χ2(2) = 1732, p < 0.001]. Further evidence for the cusp model’s preference is shown in Figure 4, with multimodal distributions at various areas across the response surface and bimodality within the bifurcation area (bottom right figure). Visually speaking, Figure 5 displays the observations as they oscillate between upper and lower surfaces and within the bifurcation area again fitting the expectations of the cusp catastrophe model. An ancillary to Figure 5, is Figure 6, which displays the relative position of the observations to the upper and lower surfaces using a control plane scatterplot. Observations with darker colors (e.g., purple) are positioned closer to the upper surface and those with lighter colors (e.g., light green) are closer to the lower surface (see Grasman et al., 2009).

Parameter estimates of the cusp model for the prediction of school readiness using a combination of asymmetry (student and teacher behaviors) and bifurcation (school and class size) predictors.

Cusp model intercept and slope terms Unstandardized B S.E. Z-test Value of p
a (Intercept) −0.870 0.027 −31.986 <0.001***
a1 (Student behavior) 0.196 0.012 15.911 <0.001***
a2 (Teacher behavior) 0.264 0.019 13.803 <0.001***
b (Intercept) −0.421 0.077 −5.468 <0.001***
b1 (School size) 0.001 0.001 18.933 <0.001***
b2 (Class size) 0.004 0.001 2.212 0.027*
w (Intercept) −0.496 0.009 −51.496 <0.001***
w (School readiness) 0.770 0.004 173.993 <0.001***

Note: The terms a, b, and w refer to asymmetry, bifurcation, and outcome variables’ intercept terms with those including a subscript the specific slopes for the asymmetry and bifurcation variables. ***p <. 001; **p < 0.01; *p < 0.05.

Comparing linear and cusp models using information criteria.

Models tested Loglikelihood Par AIC AICc BIC
1. Linear −20651.560 6 41315.11 41315.12 41360.87
2. Logistic −20579.590 7 41173.18 41173.18 41226.56
3. Quadratic −20629.389 8 41274.78 41310.37 41335.79
4. Cusp −19784.772 9 39587.55 39587.56 39656.19

Par, number of freely estimated parameters; AIC, akaike criterion; AICc, corrected AIC; BIC, Bayesian information criterion. Preference is given to the BIC. Its sample-adjusted variant is not included because of the large sample size of the PISA 2018 cohort.

Frequency of observations at various places on the response surface as indicated by the shaded areas. The presence of skew and multimodality are apparent in several areas of the lower response surface.

Cusp catastrophe upper and lower surfaces with observations transitioning between the two and within the bifurcation area.

Control plane scatterplot with darker dots showing placement of observations close to the upper surface and light color dots towards the lower surface.

 Discussion

The purpose of the present study was to understand students’ school readiness as a function of student and teacher behaviors but also school size and class size using both linear and non-linear analytical approaches. Results pointed to a preference for the cusp catastrophe model in that the relationship between school and class size with achievement is determined by specific thresholds of these variables.

Past research indicates that the size of a class has a significant influence on the academic achievement of students. The work by Kenayathulla et al. (2019) favored the role and functioning of smaller classrooms and their positive impact on academic achievement. Studies on the Portugal Programme Mais Sucesso Escolar (PMSE) indicate that factors such as class size, composition, and tailored instruction might lead to a decrease in grade repetition and an enhancement in academic achievement (Barata et al., 2015). Nevertheless, these findings have contentious ramifications for educational policy. Other studies on the other hand (e.g., Filges et al., 2018), pointed out that decreasing class sizes has minimal impact and that there are more cost-effective methods for enhancing student achievement. These unequivocal findings demonstrate that class size and outcomes vary based on the circumstances and composition of the student population (Milesi and Gamoran, 2006).

The most important finding of the present study was that preference for the cusp model allowed us to identify important thresholds for which student readiness is no longer predictable. These thresholds were 801 students for school size and 27 students for class size. Interestingly, the estimates for school size agree with earlier suggestions using quadratic models suggesting that between 600 and 900 students is the optimal school size (Lee and Smith, 1997) and also the work of Andrews et al. (2002) who reported dysfunctional schools when exceeding 1,000 students. For class size, earlier work suggested diminishing returns in that reducing class size from 30 to 25 students is more beneficial compared to reducing it from 20 to 15 (Mosteller, 1995; Krueger, 1999). Thus, the currently identified threshold falls within earlier predictions (Word et al., 1990; Finn and Achilles, 1999; Nye et al., 2000).

Study implications for educational policy

The growing body of research suggesting that larger school and class sizes harm student achievement requires a thorough reassessment of educational systems. When class and school sizes become too large, the amount of attention given to each student decreases, which can hinder customized instruction and result in a decrease in academic performance (Blatchford, 2003; Hattie, 2006). This issue highlights the necessity for policymakers to adopt initiatives focused on maintaining or decreasing class and school sizes to cultivate more efficient learning environments. Possible approaches could involve implementing strict class size restrictions, especially in early schooling where personalized attention is vital, and reorganizing bigger educational institutions into smaller learning communities to improve individualization and assistance (Lee and Smith, 1997). Furthermore, it is important to implement laws that provide fair and equal access to small-sized classrooms and schools among all socioeconomic and demographic groups. This will help to resolve any potential inequalities in educational achievements. Furthermore, a transition to smaller educational environments requires corresponding improvements in teacher recruiting, training, and retention methods, guaranteeing that the standard of education remains uncompromised. In conclusion, although there are difficulties in managing and funding efforts to optimize class and school sizes, the possibility of achieving substantial enhancements in student performance makes it an essential area of concentration for educational reform and policy formulation.

 Study limitations and future directions

The variables “school size,” and “classroom size” have long been recognized as significant in research although past studies have presented several methodological and design deficiencies. The existing research on the functioning of school and classroom size has several limitations. Firstly, it lacks generalizability due to the impracticality of randomly assigning students to schools and classes advising caution before generalizing the present findings beyond the specific sample. Secondly, there is a lack of consistent measures, particularly in quantifying the distinction between “large” and “small” schools with proxy measurements including the number of teachers, the number of students, or other ratio variables, which are also challenged by the violation of distributional assumptions. There is also a need to engage multivariate analyses that can consider the variations in sampling across students and schools using both linear and non-linear means. In the present study, the use of international population data using rigorous procedures for sampling and representativeness in each country overcomes one of the major limitations of past studies. Furthermore, the analytical framework utilized here is not without limitations. The cusp catastrophe model has been criticized as overfitting the data, thus, limiting model generality (Poston and Stewart, 1978). Concerns about the reliability of the findings from small samples have also been raised (Zhang, 2016) as well as the model’s generality to real-world phenomena (Schelling, 1973). Others were concerned that the complexity and uncertainty of real-world phenomena cannot be captured by a set of mathematical equations (Borsboom and Cramer, 2013) and specifically one type of asymmetry measured within the cusp catastrophe model (Cramer, 2008).

In the future, we advise the use of the present analytical framework using a per-country analysis as well as the invariance of the findings across important moderating variables such as gender, SES, urbanicity, private or public schooling, and other variables which were found to be important predictors of a school’s climate.

 Conclusion

In this comprehensive examination of the factors influencing students’ school readiness, a significant discovery emerged: the relationship between school and class sizes with achievement is best described by the non-linear complexities of the Cusp catastrophe model. This study, utilizing data from over 21,000 schools across 80 countries, revealed critical thresholds at a school size of 801 students and a class size of 27 students. Beyond these points, size increases are associated with unpredictability and decreased school readiness. This suggests a pronounced shift in the traditional understanding of educational environments, emphasizing the importance of maintaining optimal class and school sizes to ensure effective learning and teaching. The findings underscore the need for policymakers to reconsider current educational structures, advocating for more personalized and manageable learning environments to enhance student achievement. While this study provides a groundbreaking insight into the dynamics of educational settings, its reliance on the cusp catastrophe model and the specific thresholds identified necessitates further investigation and validation to ensure widespread applicability and understanding of its implications in the complex landscape of educational policy and practice.

 Data availability statement

Publicly available datasets were analyzed in this study. This data can be found at: https://www.oecd.org/pisa/data/2018database/.

 Ethics statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent from the patients/ participants or patients/participants legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.

 Author contributions

FA: Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review & editing. MA: Data curation, Funding acquisition, Investigation, Writing – review & editing. KK: Formal analysis, Methodology, Software, Writing – review & editing.

 Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The authors extend their appreciation to the Deputyship for Research and Innovation, “Ministry of Education” in Saudi Arabia for funding this research (IFKSUOR3-420-4).

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyg.2024.1354072/full#supplementary-material

 References Ajami I. R. Akinyele O. B. (2014). Effect of student-teacher ratio on academic achievement of selected secondary school students in Port Harcourt Metropolis, Nigeria. J. Educ. Pract. 5, 100–106. Alexander R. A. Herbert G. R. DeShon R. P. Hanges P. J. (1992). An examination of least-squares regression modeling of catastrophe theory. Psychol. Bull. 111, 366–374. doi: 10.1037/0033-2909.111.2.366 Andrews M. Duncombe W. Yinger J. (2002). Revisiting economies of size in American education: are we any closer to a consensus? Econ. Educ. Rev. 21, 245–262. doi: 10.1016/S0272-7757(01)00006-1 Angrist J. D. Lavy V. (1999). Using Maimonides’ rule to estimate the effect of class size on scholastic achievement. Q. J. Econ. 114, 533–575. doi: 10.1162/003355399556061 Arias J. Walker D. (2004). Additional evidence on the relationship between class size and student performance. J. Eco. Edu. 35, 311–329. Bandiera O. Larcinese V. Rasul I. (2010). Heterogeneous classsize effects: New evidence from a panel of university students. Econ. J. 120, 1365–1398. Bedard K. Kuhn P. (2008). Where class size really matters: Class size and student ratings of instructor effectiveness. Eco. Edu. Rev. 27, 253–265. Bellante D. M. (1972). A summary report on student performance in mass lecture classes of economics. J. Eco. Edu. 4, 53–54. Benton SL Cashin WE (2012) Student ratings of teaching: A summary of research and literature. IDEApaper #50, Manhattan, KS: The Idea Center. Blatchford P. (2003). The class size debate: Is small better? Milton Keynes, England: Open University Press. Blatchford P. Bassett P. Brown P. A. (2011). Examining the effect of class size on classroom engagement and teacher–pupil interaction: differences in relation to pupil prior attainment and primary vs. secondary schools. Learn. Instr. 21, 715–730. doi: 10.1016/j.learninstruc.2011.04.001 Blatchford P. Lai K. C. (2012). “Class size: arguments and evidence” in International encyclopedia of education. eds. McGraw B. Baker E. Peterson P. P.. 3rd ed (Oxford, UK: Elsevier) Borsboom D. Cramer A. O. (2013). Network analysis: an integrative approach to the structure of psychopathology. Annu. Rev. Clin. Psychol. 9, 91–121. doi: 10.1146/annurev-clinpsy-050212-185608 Bowne J. B. Magnuson K. A. Schindler H. S. Duncan G. J. Yoshikawa H. (2017). A Meta-analysis of class sizes and ratios in early childhood education programs: are thresholds of quality associated with greater impacts on cognitive, achievement, and socioemotional outcomes? Educ. Eval. Policy Anal. 39, 407–428. doi: 10.3102/0162373716689489 Bradley S. Taylor J. (1998). The effect of school size on exam performance in secondary schools. Oxf. Bull. Econ. Stat. 60, 291–324. doi: 10.1111/1468-0084.00102 Brown C. (1995). Chaos and catastrophe theories. London: Sage. Chen D. Chen X. (2017). Cusp catastrophe regression and its application in public health and behavioral research. Int. J. Environ. Res. Public Health 14:1220. doi: 10.3390/ijerph14101220, PMID: 29027967 Chen D. Chen X. Lin F. Tang W. Lio Y. Guo Y. (2014). Cusp catastrophe polynomial model: power and sample size estimation. Open J. Stat. 4, 803–813. doi: 10.4236/ojs.2014.410076, PMID: 27158562 Chen X. Wang K. Chen D. G. (2020). “Cusp catastrophe regression analysis of testosterone in bifurcating the age-related changes in PSA, a biomarker for prostate cancer” in Statistical methods for global health and epidemiology: Principles, methods and applications. eds. Chen X. Chen D. G. (Cham, Switzerland: Springer), 353–372. Cobb L. (1981). Parameter estimation for the cusp catastrophe model. Behav. Sci. 26, 75–78. doi: 10.1002/bs.3830260107 Cobb L (1998). An introduction to cusp surface analysis Technical report, Aetheling Consultants, Louisville, CO, USA. Cobb L. Koppstein P. Chen N. H. (1983). Estimation and moment recursion relations for multimodal exponential distributions. J. Am. Stat. Assoc. 78, 124–130. doi: 10.1080/01621459.1983.10477940 Cobb L. Zacks S. (1985). Applications of catastrophe theory for statistical modeling in the biosciences. J. Am. Stat. Assoc. 80, 793–802. doi: 10.1080/01621459.1985.10478184 Cramer A. O. (2008). The cusp catastrophe model of psychosis: testing a bifurcation hypothesis of the structure of psychopathology. J. Abnorm. Psychol. 117, 814–824. Crispin L. M. (2016). School size and student achievement: does one size fit all? East. Econ. J. 42, 630–662. doi: 10.1057/eej.2015.2 De Paola M. Ponzo M. Vincenzo S. (2013). Class size effects on student achievement: Heterogeneity across abilities and fields. Edu. Econ. 21, 135–153. Edgell J. (1981). Effects of class size upon aptitude and attitude of pre-algebra undergraduate students. ERIC Document 203760 Available at: https://eric.ed.gov/?id=ED203760 Ehrenberg R. G. Brewer D. J. Gamoran A. Willms J. D. (2001). Class size and student achievement. Psychol. Sci. Public Interest 2, 1–30. doi: 10.1111/1529-1006.003 Filges T. Sonne-Schmidt Viinholt Nielsen B. (2018). Small class sizes for improving student achievement in primary and secondary schools: A systematic review. Campbell Systematic Reviews 14, 1–107. doi: 10.4073/csr.2018.10 Finn J. D. Achilles C. M. (1999). Tennessee's class size study: findings, implications, misconceptions. Educ. Eval. Policy Anal. 21, 97–109. doi: 10.3102/01623737021002097 Foreman-Peck J. Foreman-Peck L. (2006). Should schools be smaller? The size-performance relationship for welsh schools. Econ. Educ. Rev. 25, 157–171. doi: 10.1016/j.econedurev.2005.01.004 Garrett Z. Newman M. Elbourne D. Bradley S. Noden P. Taylor J. . (2004). “Secondary school size: a systematic review” in Research evidence in education library (London: EPPI-Centre, Social Science Research Unit, Institute of Education, University of London) Gereshenson S. Langbein L. (2015). The effect of primary school size on academic achievement. Educ. Eval. Policy Anal. 37, 135S–155S. doi: 10.3102/0162373715576075 Gleason J. (2012). Using technology-assisted instruction and assessment to reduce the effect of class size on student outcomes in undergraduate mathematics courses. College Teaching 60, 87–94. Grasman R. P. van der Maas H. L. J. Wagenmakers E. J. (2009). Fitting the cusp catastrophe in R: a cusp package primer. J. Stat. Softw. 32, 1–27. doi: 10.18637/jss.v032.i08 Guastello S. (1984). Cusp and butterfly catastrophe modeling of two opponent process models: drug addiction and work performance. Behav. Sci. 29, 258–262. doi: 10.1002/bs.3830290405, PMID: 6517835 Guastello S. J. (1987). A butterfly catastrophe model of motivation in organization: academic performance. J. Appl. Psychol. 72, 165–182. doi: 10.1037/0021-9010.72.1.165 Guastello S. J. (1992). Clash of the paradigms - a critique of an examination of the polynomial regression technique for evaluating catastrophe -theory hypotheses. Psychol. Bull. 111, 375–379. doi: 10.1037/0033-2909.111.2.375 Guastello S. J. (2001). Nonlinear dynamics in psychology. Discret. Dyn. Nat. Soc. 6, 11–29. doi: 10.1155/S1026022601000024 Guastello S. J. (2002). Managing emergent phenomena: Non-linear dynamics in work organizations. Mahwah, NJ: Lawrence. Guclu M. (2002). Crowded classroom problem in elementary schools and the solutions to this problem. Eurasian J. Educ. Res. 9, 52–58. Hancock T. M. (1996). Effects of class size on college student achievement. College Student Journal 30, 479–481. Hattie J. (2006) The paradox of reducing class size and improved learning outcomes. International J. Edu. Res. 43: 387–425. Hill M. C. (1998). Class size and student performance in introductory accounting courses: Further evidence. Issues in Accounting Education 13, 47–64. Howley C. (1996). Compounding disadvantage: the effects of school and district size on student achievement in West Virginia. J. Res. Rural. Educ. 12, 25–32. Howley C. Strange M. Bickel R. (2000). When it comes to schooling…Small works: School size, poverty and student achievement. Knoxville, TN: ERIC Document Reproduction Service. Hufford M. Witkiewitz K. Shields A. Kodya S. Caruso J. (2003). Relapse as a nonlinear dynamic system: application to patients with alcohol use disorders. J. Abnorm. Psychol. 112, 219–227. doi: 10.1037/0021-843x.112.2.219 Jarvis TJ (2007) Class Size and Teacher Effects on Student Achievement and Dropout Rates in University-Level Calculus. Available at: https://www.math.byu.edu/~jarvis/class-size/class-size-preprint.pdf (Accessed December 2014). Jones K. R. Ezeife A. N. (2011). School size as a factor in the academic achievement of elementary school students. Psychology 2, 859–868. doi: 10.4236/psych.2011.28131 Kennedy P. E. Siegfried J. J. (1997). Class size and achievement in introductory economics: evidence from the TUCE III data. Econ.Edu. Rev. 16, 385–394. Kira I. Barger B. Shuwiekh H. Kucharska J. Alhuwailah A. (2019). Cumulative stressors and traumas and suicide: a non-linear cusp dynamic systems model. Psychology 10, 1999–2018. doi: 10.4236/psych.2019.1015128 Kira I. Barger B. Shuwiekh H. Kucharska J. Alhuwailah A. (2020). The threshold non-linear model for the effects of cumulative stressors and traumas: a chained cusp catastrophe analysis. Psychology 11, 385–403. doi: 10.4236/psych.2020.113025 Koc N. Celik B. (2015). The impact of number of students per teacher on student achievement. Proced. Soci. Behav. Sci. 177, 65–70. doi: 10.1016/j.sbspro.2015.02.335 Kokkelenberg EC Dillon M Christy SM (2006) The effects of class size on student grades at a public university. Cornell Higher Education Research Institute Working Paper #88. Available at: http://digitalcommons.ilr.cornell.edu/workingpapers/66/. Konstantopoulos S. Shen T. (2023). Class size and teacher effects on non-cognitive outcomes in grades K-3: a fixed effects analysis of ECLS-K:2011 data. Large-scale Assess. Educ. 11, 1–24. doi: 10.1186/s40536-023-00182-8 Krassel K. F. Heinesen E. (2014). Class-size effects in secondary school. Educ. Econ. 22, 412–426. doi: 10.1080/09645292.2014.902428 Krueger A. (1999). Experimental estimates of education production functions. Q. J. Econ. 114, 497–532. doi: 10.1162/003355399556052 Kwan K. (1999). How fair are student ratings in assessing the teaching performance of university teachers? Assessment and Evaluation in Higher Education 24, 181–195. Lammers W. J. Murphy J. J. (2002). A profile of teaching techniques used in the university classroom: A descriptive profile of a US public university. Active Learning in Higher Education 33, 54–67. Lee V. Loeb S. (2000). School size in Chicago elementary schools: effects on teachers’ attitudes and students’ achievement. Am. Educ. Res. J. 37, 3–31. doi: 10.3102/00028312037001003 Lee V. E. Smith J. B. (1997). High school size: which works best and for whom? Educ. Eval. Policy Anal. 19, 205–227. doi: 10.3102/01623737019003205 Leithwood K. Jantzi D. (2009). A review of empirical evidence about school size effects: a policy perspective. Rev. Educ. Res. 79, 464–490. doi: 10.3102/0034654308326158 Li G. Kharabsheh R. (2022). Dynamic nonlinear system based on complex system theory in the development of vocational education. Applied Math. Nonlinear Sci. 8, 961–968. doi: 10.2478/amns.2022.2.0083 Lowenthal P. R. Nyland R. Jung E. Dunlap J. C. Kepka J. (2019). Does class size matter? An exploration into faculty perceptions of teaching high-enrollment online courses. Am. J. Dist. Educ. 33, 152–168. doi: 10.1080/08923647.2019.1610262 Lytton H. Pyryt M. C. (1998). Predictors of achievement in basic skills: a Canadian effective schools study. Can. J. Educ. 23, 281–301. doi: 10.2307/1585940 Ma X. Klinger D. A. (2000). Hierarchical linear modeling of student and school effects academic achievement. Can. J. Educ. 25, 41–55. doi: 10.2307/1585867 Mandel P. Sussmuth B. (2011). Size matters. The relevance and Hicksian surplus of preferred college class size. Economics of Education Review 30, 1073–1084. Martin J. P. (2015). Moving up in the U.S. News and World Report Rankings. Change: The Magazine of Higher Learning 47, 52–61. Matta B. N. Guzman J. M. Stockly S. K. . (2015). Class size effects on student performance in a Hispanic serving institution. The Review of Black Political Economy 42, 443–457. McGuire K. (1989). School size: the continuing controversy. Educ. Urban Soc. 21, 164–174. doi: 10.1177/0013124589021002005 Meier D. (1996). The big benefits of smallness. Educ. Leadersh. 54, 12–15. Milesi C. Gamoran A. (2006). Effects of Class Size and Instruction on Kindergarten Achievement. Educational Evaluation and Policy Analysis, 28, 287–313. doi: 10.3102/01623737028004287 Minnotte M. C. Scott D. W. (1993). The mode tree: a tool for visualization of nonparametric density features. J. Comput. Graph. Stat. 2, 51–68. doi: 10.2307/1390955 Molenaar P. C. Oppenheimer L. (1985). Dynamic models of development and mechanistic-organismic controversy. New Ideas Psychol. 3, 233–242. doi: 10.1016/0732-118X(85)90017-0 Monks J. Schmidt R. (2011). The impact of class size on outcomes in higher education. Retrieved from. Available at: http://digitalcommons.ilr.cornell.edu/workingpapers/114/ Mosteller F. (1995). The Tennessee study of class size in the early school grades. Futur. Child. 5, 113–127. doi: 10.2307/1602360 Muir E. (2000). Smaller schools: how much more than a fad? American Educator, 24, 40–46. Nandrup A. B. (2016). Do class size effects differ across grades? Educ. Econ. 24, 83–95. doi: 10.1080/09645292.2015.1099616 Newman M. Garrett Z. Elbourne D. Bradley S. Noden P. Taylor J. . (2006). Does secondary school size make a difference? A systematic review. Educ. Res. Rev. 1, 41–60. doi: 10.1016/j.edurev.2006.03.001 Nye B. Hedges L. V. Konstantopoulos S. (2000). The effects of small classes on academic achievement: the results of the Tennessee class size experiment. Am. Educ. Res. J. 37, 123–151. doi: 10.3102/00028312037001123 OECD (2019), PISA 2018 assessment and analytical framework, PISA, OECD Publishing, Paris. Oliva T. Desarbo W. Day D. Jedidi K. (1987). GEMCAT: a general multivariate methodology for estimating catastrophe models. Behav. Sci. 32, 121–137. doi: 10.1002/bs.3830320205 Olson J. C. Cooper S. Lougheed T. (2011). Influences of teaching approaches and class size on undergraduate mathematical learning. PRIMUS 21, 732–751. Papageorgiou G. Stamovlasis D. Johnston P. (2010). Primary teachers’ particle ideas and explanations of physical phenomena: the effect of an in-service training course. Int. J. Sci. Educ. 32, 629–652. doi: 10.1080/09500690902738016 Poston T. Stewart I. N. (1978). Catastrophe theory and its applications. London, UK: Pitman Publishing. Sapelli C. Illanes G. (2016). Class size and teacher effects in higher education. Economics of Education Review 52, 19–28. Scheerens J. Hendricks M. Luyten H. (2014). “School Size Effects: Review and Conceptual Analysis” in School Size Effects Revisited. (Springer, Cham: Springer Briefs in Education), 7–40. Schelling T. C. (1973). Review of structural stability and morphogenesis by René Thom. J. Polit. Econ. 81, 1097–1105. Stamovlasis D. (2006). The nonlinear dynamical hypothesis in science education problem solving: a catastrophe theory approach. Nonlin. Dynam. Psychol. Life Sci. 10, 37–70. Stamovlasis D. (2011). Nonlinear dynamics and neo-piagetian theories in problem solving perspectives on a new epistemology and theory development. Nonlinear Dynamics Psychol. Life Sci. 15, 145–173. PMID: 21382259 Stamovlasis D. Tsaparlis G. Kamilatos C. Zarotiadou E. Papaoikonomou D. (2005). Conceptual understanding versus algorithmic problem solving: further evidence from a National Examination. Chemistry Educ. 6, 104–118. doi: 10.1039/B2RP90001G Stewart I. N. Peregoy P. L. (1983). Catastrophe theory modeling in psychology. Psychol. Bull. 94, 336–362. doi: 10.1037/0033-2909.94.2.336 Tasker M. (2003). Smaller structures in secondary education: A research digest. Bristol: Human Scale Education. Thom R. (1975). Structural stability and morphogenesis. Reading, MA: W.A. Benjamin. Tseng H. (2010). Has the student performance in managerial economics been affected by the class size of principles of microeconomics. J. Econom. Econom. Educ. Res. 11:15. Tsitsipis G. Stamovlasis D. Papageorgiou G. (2012). A probabilistic model for students’ errors and misconceptions on the structure of matter in relation to three cognitive variables. Int. J. Sci. Math Educ. 10, 777–802. doi: 10.1007/s10763-011-9288-x Van der Maas H. L. J. Molenaar P. C. M. van der Pligt J. (2003). Sudden transitions in attitudes. Soci. Meth. Res. 32, 125–152. doi: 10.1177/0049124103253773 Van der Maas H. L. J. Molennar P. C. (1992). Stagewise cognitive development: an application of catastrophe theory. Psychol. Rev. 99, 395–417. doi: 10.1037/0033-295X.99.3.395 Werblow J. Duesbery L. (2009). The impact of high school size on math achievement and dropout rate. High School J. 92, 14–23. doi: 10.1353/hsj.0.0022 Westerlund J. (2008). Class size and student evaluations in Sweden. Education Economics 16, 19–28. Wobmann L. West M. (2006). Class-size effects in school systems around the world: evidence from between-grade variation in TIMSS. Eur. Econ. Rev. 50, 695–736. doi: 10.1016/j.euroecorev.2004.11.005 Word E. Johnston J. Bain H. P. Fulton B. D. (1990). Student/teacher achievement ratio (STAR), Tennessee's K-3 class size study: Final summary report 1985–1990. Nashville, TN: Tennessee State Department of Education. Yamamori K. Tokuoka M. Hagiwara Y. Oouchi Y. Nakamoto K. Isoda T. (2021). Effects of class size and provision of learning goals and feedback on students’ two-year achievement trajectories. The Japanese. J. Educ. Psychol. 69, 297–316. doi: 10.5926/jjep.69.297 Zhang K. (2016). An exploratory statistical cusp catastrophe model.