- 1Department of Computer Science & Engineering, United International University, Dhaka, Bangladesh
- 2Institute of Natural Sciences, United International University, Dhaka, Bangladesh
- 3Institute for Advanced Research, United International University, Dhaka, Bangladesh
- 4Department of Pediatrics, Shaheed Suhrawardy Medical College, Dhaka, Bangladesh
Introduction: Mathematics anxiety (MA) is a distinct negative emotional state or trait that individuals experience when confronted with mathematical problems in everyday life and academic contexts. This study aims to identify the key predictors of MA among secondary-level students in Bangladesh.
Methods: Utilizing a quantitative cross-sectional research design, data were collected from 486 students across 89 institutions. Later, the data were analyzed using Partial Least Squares Structural Equation Modeling (PLS-SEM).
Results: The findings revealed that math related negative past experiences (β = 0.241, t = 4.914, p < 0.001) and a perceived lack of teacher support (β = 0.234, t = 5.440, p < 0.001) significantly contribute to students’ low self-efficacy in mathematics. This low self-efficacy is further influenced by negative attitudes and test anxiety, ultimately leading to increased MA (β = 0.694, t = 22.695, p < 0.001). Additionally, cognitive challenges, particularly working memory difficulties, directly affect MA (β = 0.110, t = 2.659, p = 0.008). The study also found that negative attitudes (β = 0.347, t = 9.063, p < 0.001) and test anxiety (β = 0.251, t = 5.913, p < 0.001) independently exacerbate MA. Moreover, a lack of motivation in learning mathematics is directly influenced by this elevated level of MA (β = 0.384, t = 9.939, p < 0.001).
Discussion: Taken together, the study proposes several key recommendations and policy implications to inform the development of synchronized policies by educational authorities aimed at combatting, reducing MA among secondary-level students in Bangladesh and similar contexts.
1 Introduction
Globally, anxiety disorders have emerged as the most prevalent and widespread mental health issue. According to the World Health Organization (WHO), approximately 4% of the global population currently experiences at least one anxiety disorder (1). While occasional anxiety is a normal aspect of life, it becomes problematic when it transforms into frequent, persistent worry or fear. There are various forms of anxiety disorders, including generalized anxiety disorder, social anxiety disorder, panic disorder, and phobia-related disorders etc. In educational settings, test anxiety and performance anxiety are particularly common among students. However, a distinct form of anxiety disorder, known as MA, has recently gained increasing focus as a frequent and alarming issue on a global scale in all age groups of students (2, 3). MA has been defined in various ways by researchers, educational psychologists (EPs) and practitioners from different point of view. Consequently, there is no universally accepted or rigid definition for the term MA. However, in most cases, MA is defined as a distinct negative emotional state or trait that individuals experience when confronted with mathematical problems in everyday life and academic contexts (4, 5). Beyond its negative effects on mental health and well-being, MA is well-documented for its negative correlation with mathematical competence, performance and achievement (6, 7), making it a potential barrier to success in science, technology, engineering, and mathematics (STEM) careers (8).
Though MA has short, medium, and long-term negative consequences on students’ self-confidence, math performance, math achievement, and even career choices, the specific factors and other events that trigger to the development of MA remain still unclear (9). However, existing studies broadly categorize these factors into personal, psychological, environmental, cognitive, and pedagogical influences. Among personal factors, gender has received significant attention in relation to MA. Many studies have found that females are more prone to MA than males (10–13), which would suggest that females’ math performance should also be more negatively affected. However, despite these findings, research often reveals little to no gender differences in actual math performance (14–16). Another potential factor that has been investigated under the personal factors is genetics. For instance, Malanchini et al. (17) found that 75% of the genetic variance in MA is tied to math-related attitudes, abilities, and achievement. Wang et al. (18) similarly reported that 40% of the variance in MA is genetically influenced, either directly or through contributions to mathematical ability and general anxiety. In terms of psychological factors, self-efficacy is frequently studied in connection with MA (19). For example, Khasawneh et al. (10) found a significant negative correlation between MA and self-efficacy levels, suggesting that improving students’ self-efficacy could enhance their belief in their capabilities, thereby reducing their anxiety toward math. Similar negative correlation was also observed in the study of Rozgonjuk et al. (20), where the authors suggested that since mathematics self-efficacy plays a significant role in MA, boosting self-efficacy could be an effective strategy for reducing MA. In this context, environmental factors, such as teacher support is crucial. Research has shown that student-perceived math teacher support influences MA through the mediation of the teacher–student relationship and math self-efficacy (21, 22). By fostering positive teacher–student interactions and enhancing self-efficacy, MA can be effectively reduced through increasing mathematical problem-solving ability (23). Conversely, negative classroom experiences, such as adverse teacher behavior or failures in math, are negatively correlated with students’ self-efficacy (3, 24). Similarly, other environmental factors, such as parental stereotypes, culture, social pressure also found influencing effect on MA (2, 3, 25). When examining cognitive factors, researchers often focus on the interaction between MA and working memory (26, 27). According to existing studies, this interaction is reciprocal: higher levels of MA impair working memory (26, 28), leading to more errors and poorer performance in math (29). This poor performance, in turn, further increases MA, creating a vicious cycle (30). Finally, in terms of pedagogical factors, the most frequently documented influence on MA is teaching or instructional strategies. Studies have shown that a rigid curriculum and traditional instructional strategies can significantly contribute to the development of MA (2, 31). Conversely, incorporating technology-enhanced and modernized approaches to mathematics instruction can create an engaging, active learning environment that reduces anxiety and fosters a more positive attitude toward learning math (32).
In educational contexts, it is widely recognized by educational psychologists to practitioners that MA represents a significant barrier to student success across all academic levels. Over the years, numerous studies have investigated the origins of MA and its associated factors, yet few have delved into the intricate causal relationships that play a critical role in predicting MA. Furthermore, most of the existing research has predominantly focused on Western countries, while limited attention given to developing nations such as Bangladesh. To address this gap and support the advancement of domain knowledge, this study aims to identify and analyze the key factors that contribute in shaping MA within the context of Bangladesh, particularly at the secondary education level.
The remaining sections of this study organized as follows:Section 2 briefly discusses the materials and methods. Section 3 presents the findings step-by-step. Section 4 interprets these results in the context of previous research. Section 5 offers recommendations and policy implications. And section 6 highlights the study’s limitations and suggests future research directions.
2 Materials and methods
2.1 Conceptual model and hypotheses
This study introduces and empirically evaluates a comprehensive conceptual model designed to identify the factors of mathematics anxiety among secondary school students in Bangladesh (see Figure 1). The model elucidates the multi-faceted nature of factors contributing to mathematics anxiety, emphasizing both direct and indirect pathways. The core of the model comprises eight constructs: (i) negative past experiences (NPE) in mathematics learning, (ii) lack of perceived teacher support (LPTS), (iii) cognitive challenges (CC), (iv) low self-efficacy (LSE), (v) negative attitude (NA) towards mathematics, (vi) test anxiety (TA), (vii) mathematics anxiety (MA), and (viii) lack of motivation (LM). Four of these constructs are classified as exogenous (independent latent variables), while the remaining four are endogenous (dependent latent variables). The model posits that negative past experiences and a lack of perceived teacher support directly contribute to the development of low self-efficacy in mathematics among students. This low self-efficacy, further influenced by negative attitudes and test anxiety (mediators), is then hypothesized to directly lead to MA. Additionally, the model proposes that cognitive challenges, specifically working memory difficulties, can directly influence mathematics anxiety. Negative attitudes and test anxiety are also theorized to have a direct effect on increasing MA. Finally, the model suggests that a lack of motivation, manifested as math avoidance, is directly influenced by higher levels of MA. The selection of constructs and the direction of their hypothesized relationships are grounded in theoretical foundations established within the relevant academic literatures (2, 10–13, 17–19, 23, 26, 27, 30–32). Finally, according to the conceptual model illustrated in Figure 1, the study hypotheses can be formed as follows:
H1: Negative Past Experiences in Maths Learning are positively associated with Low Self-Efficacy.
H2: Lack of Perceived Teacher Support is positively associated with Low Self-Efficacy.
H3: Low Self-Efficacy is positively associated with Mathematics Anxiety.
H3A: Low Self-Efficacy is positively associated with Negative Attitudes toward mathematics.
H3B: Negative Attitudes mediate the relationship between Low Self-Efficacy and Mathematics Anxiety
H3C: Low Self-Efficacy is positively associated with Test Anxiety.
H3D: Test Anxiety mediates the relationship between Low Self-Efficacy and Mathematics Anxiety.
H4: Cognitive Challenges (working memory difficulties) are positively associated with Mathematics Anxiety.
H5: MA is positively associated with Lack of Motivation (Math Avoidance).
2.2 Study design and settings
This study employed a quantitative cross-sectional research design conducted between April 15, 2024, and May 28, 2024. A total of 486 secondary-level students from 89 institutions participated. The cross-sectional design was selected because it analyzes data from a population at a single point in time, rather than following individuals over time (33, 34). This approach is relatively quick and inexpensive, making it an efficient method for generating hypotheses (35). By capturing a snapshot of math anxiety prevalence and examining causal relationships between factors at a specific point in time, this design was ideal and practical for addressing questions about the current state of the phenomenon studied.
2.3 Study population and sample size
The target population of this study is secondary-level students in Bangladesh. According to the Bangladesh Bureau of Educational Information and Statistics (BANBEIS), as of 2023, there were 8166188 secondary-level students in 18968 educational institutions in the country (36). The sample size for this research consisted of 486 secondary-level students from 89 institutions, where 297 (61.1%) were male and 189 (38.8%) were female. The participants were between the ages of 12 and 17, with an average age of 14.64 years (standard deviation = 1.93). Although the sample size might seem insufficient compared to the target population, the study utilized Partial Least Squares Structural Equation Modeling (PLS-SEM), a methodology widely known for its statistical robustness even with smaller sample sizes (34, 37). The conventional method for sample size estimation in PLS-SEM is the “10 times rule,” which stipulates that the sample size must be at least 10 times the number of indicators used to measure a construct (38, 39). However, this rule has been criticized for potentially leading to overestimation or underestimation of sample size requirements, as it is not model-specific (34, 40, 41). Therefore, this study employed the G*Power software (version 3.1.9.4) to estimate the required sample size. The G*Power software is widely used and recognized as a reliable tool for sample size determination (34, 37). The input parameters used for the calculation included: Effect size f²=0.15 (medium), α err prob=0.05, Power (1-β err prob) =0.95, Number of tested predictors=3, and Total number of predictors =7. The output showed that the required minimum sample size was 119, where the Actual power was 0.9507. On the other hand, the actual sample size of this study was 486 which was sufficient to ensure standard statistical power (see Figure 2).
2.4 Data collection instrument and procedure
The data collection of this study occurred between April 15 and May 28, 2024, involving 486 students from 89 secondary-level institutions in Bangladesh. A questionnaire was developed by quantitative research experts to assess the predicting factors of mathematics anxiety. Comprising 55 closed-ended questions, the questionnaire was divided into two primary sections. The first section gathered demographic information, including age, gender, residence, parental education, and mathematics GPA. The second section focused on influencing factors of mathematics anxiety, further subdivided into sub-sections where students rated agreement with statements on a Likert scale of 1 (strongly disagree) to 5 (strongly agree).
Math anxiety (MA) refers specifically to fear and emotional distress related to performing mathematical tasks, whereas test anxiety (TA) is a broader concept that encompasses anxiety experienced during any test-taking situation, regardless of subject. In this study, we measured these constructs separately using distinct scales. MA was assessed using five self-reported items derived from two validated instruments: the modified Abbreviated Math Anxiety Scale (mAMAS) (42) and the Weighted Scoring Based Rating Scale to Identify the Severity Level of Mathematics Anxiety (WSB-MARS) (5). These instruments cover various subdomains, such as computational anxiety and anxiety specifically related to mathematics tests. TA, on the other hand, was measured using items designed based on the to capture general test anxiety experienced during test situations following Jirjees et al. (43).
The validity and reliability of the questionnaire were ensured by content validity and internal consistency reliability respectively. A panel of five domain experts evaluated the questionnaire using a 5-point rating scale for relevance, representativeness, specificity, and clarity, as outlined in study of Ahmmed et al. (37). Cohen’s kappa (k = .87, p<.005***) indicated substantial inter-rater agreement (44). While Cronbach’s alpha coefficient (α=.83, p<.005***) confirmed robust internal consistency reliability (45). Data was collected through both in-person and online surveys, with appropriate approvals and informed consent obtained from participants.
2.5 Statistical analysis
To explore the data, two stages of analysis were conducted. Initially, a descriptive analysis using IBM SPSS Statistics 26 (46) was performed. Subsequently, Partial Least Squares Structural Equation Modeling (PLS-SEM) was employed using SmartPLS4 (47) to investigate the key “driver” constructs of mathematics anxiety. PLS-SEM, a sophisticated second-generation multivariate analysis method, is a variance-based technique used for estimating both the structural model (inner model) and the measurement model (outer model) (34, 48). This method is frequently used in exploratory research to develop theory and examine complex relationships between multiple latent variables (48–50). The use of PLS-SEM in this study is justified due to its flexibility in exploring and experimenting with various configurations. Additionally, as noted in prior research, PLS-SEM is often more appropriate than methods such as Covariance-based Structural Equation Modeling (CB-SEM) when the objective is to explore the key “driver” constructs. For statistical accuracy, this study utilized percentile bootstrapping with 5000 sub-samples, employing a fixed seed to ensure the reproducibility of results. A significance level of 0.05 was applied, using a two-tailed test with parallel processing to efficiently generate robust standard errors and significance levels for the path coefficients. This approach ensured the reliability of the parameter estimates and their associated significance.
2.6 PLS-SEM parameter settings
To optimize the accuracy and convergence of the model, specific settings in SmartPLS4 were configured. The initial weights for the analysis were set to 1.0, ensuring equal starting values for all indicator variables before optimization. A maximum number of 3000 iterations was allowed to ensure the model had sufficient computational cycles to converge, especially in the case of complex models. The stop criterion was set to 107, ensuring the estimation process would stop once the change in parameter estimates across iterations was below this threshold, indicating precise convergence. All results were standardized, allowing for easy interpretation and comparison of coefficients between constructs within the model. The default Lohmoeller settings were not employed, as custom settings were more appropriate for the data and model characteristics. Additionally, the path weighting scheme was selected, as it is particularly suitable for maximizing explained variance (R²) in endogenous constructs, aligning well with the study’s focus on identifying key drivers of mathematics anxiety.
2.7 Ethical consideration
To ensure the research was conducted in accordance with ethical standards, necessary permissions were obtained from the United International University Ethics Review Board (Ref: IREB/2023/018). Informed consent was secured from all participants, and parental or guardian consent was obtained in line with the legal requirements of our jurisdiction. Prior to data collection, parents or legal guardians were invited to the respective schools for a discussion. During this meeting, the research objectives and the questionnaire were shared with them. Written consent was then obtained from those who willingly agreed to participate. This process ensured that both participants and their guardians were fully informed about the study’s purpose, procedures, and the absence of any potential risks. Furthermore, stringent measures were implemented to protect the privacy and confidentiality of participants’ data.
3 Results
3.1 Demographic information
The study sample consisted of 486 secondary-level students from 89 institutions, comprising 312 (64.2%) males and 174 (35.8%) females. The participants represented a diverse age range, with the majority (89.5%) being between 13 and 15 years old. The mean age was 13.75 years (± 2.72). Regarding socio-economic status, most students (71.6%) came from middle-class families, while 17.7% were from lower-class families, and 10.7% were from upper-class families. Notably, a significant number of students (50.3%) lived with their families, 33.9% lived with friends, and 15.8% resided alone. In terms of math anxiety, 41.4% of the total sample experienced high levels of MA, 47.5% had moderate levels, and 11.1% exhibited low levels of MA. Among male students, 20% had low MA, 51.2% had moderate MA, and 28.7% experienced high MA, with an average MA score of 12.34 (± 4.401). In contrast, female students reported higher anxiety levels: only 6.7% had low MA, while 42.7% had moderate MA, and 50.7% had high MA. The average MA score among females was 15.17 (± 4.864), which was significantly higher than that of their male counterparts. These findings highlight a noticeable gender difference in math anxiety, with female students being more prone to higher anxiety levels.
3.2 Measurement model evaluation
The first step in analyzing PLS-SEM results is to evaluate the reflective measurement model, ensuring the reliability and validity of the construct measures to support their inclusion in the path model (50). According to previous studies (34, 48, 49), the standard assessment criteria for this evaluation include three key aspects: internal consistency reliability, convergent validity, and discriminant validity.
3.2.1 Internal consistency reliability
The concept of internal consistency refers to the degree of similarity or homogeneity among the observed indicator variables. Generally, Cronbach’s alpha is used to assess internal consistency reliability (51). This measure estimates reliability based on the inter-correlations of the observed indicators, assuming all indicators have equal reliability and equal outer loadings on the construct. However, in PLS-SEM, indicators are prioritized based on their individual reliability. Consequently, composite reliability is the recommended criterion for measuring internal consistency in PLS-SEM (48). According to the previous studies, Composite reliability values between 0.60 and 0.70 are generally considered acceptable, while those between 0.70 and 0.90 indicate satisfactory internal consistency in exploratory studies (34, 50). As Table 1 demonstrates, all constructs in our analysis surpassed the recommended threshold of 0.60, confirming the adequate internal consistency of the measurement model.
3.2.2 Convergent validity
Convergent validity refers to the degree to which a measure correlates positively with alternative measures of the same construct. It implies that measures with similar or identical constructs should be substantially related. In PLS-SEM, convergent validity is typically assessed using two criteria: Average Variance Extracted (AVE) and outer loadings. The recommended benchmarks are that outer loadings should be at least 0.708 and AVE values should be at least 0.50 (38). As shown in Table 1, all constructs have AVEs exceeding 0.5, indicating that each construct explains over 50% of the variance in its items. Likewise, all outer loadings met the recommended threshold, except for TA02 (outer loading = 0.480), which was subsequently removed to improve the overall AVE. So, it can be concluded that the measurement model has achieved the required level of convergent validity.
3.2.3 Discriminant validity
Discriminant validity reflects the extent to which a construct is distinguishable from other constructs. The main objective of establishing discriminant validity is to ensure that a reflective construct exhibits stronger correlations with its indicators as compared to other constructs. To assess discriminant validity, scholars commonly use the Heterotrait-monotrait ratio (HTMT) and Fornell-Larcker criterion (52). According to exiting studies, HTMT value should be less than 0.90 to establish discriminant validity between reflective constructs (34, 53). As evident in Table 2, the HTMT of all constructs met the recommended threshold. On the other hand, Fornell-Larcker criterion involves comparing the square root of the average variance extracted (AVE) values of each construct with its correlations with other constructs within the same model (54). Here, as per the existing studies, the square root of each construct’s AVE should be greater than its highest correlation with any other construct (34). As shown in Table 3, the findings support the suggested threshold, thereby confirming the discriminant validity of the reflective construct.
3.3 Structural model evaluation
Once the validity and reliability of the measurement model have been confirmed, the next step in analyzing PLS-SEM results is to evaluate the structural model. The main goal at this stage is to examine the model’s predictive power and the relationships between constructs (37, 50, 53). To assess the structural model, the standard criteria include three key aspects: collinearity assessment, evaluation of structural model path coefficients, and the coefficient of determination (R²).
3.3.1 Collinearity statistics
Since the path coefficients in the structural model are based on OLS regressions, it is crucial to check for collinearity problems to avoid biased regression outcomes. Collinearity is determined by the construct’s variance inflation factor (VIF) value (34, 38). To prevent collinearity issues, previous studies recommend that VIF value should be above 0.20 and below 5 (50, 53). As shown in Table 4, all constructs’ VIF values fall within the recommended range, indicating that the structural model does not have collinearity issues.
3.3.2 Path coefficients (direct effects)
The structured model’s path coefficients were evaluated using bootstrapping with 5000 sub-samples, a two-tailed test, and a significance level of 0.05. As shown in Table 5, the proposed relationships were found to be statistically significant, providing empirical support for the conceptual model. Each path coefficient indicates the strength and significance of the hypothesized relationships between the constructs. For instance, the direct paths from LPTS and NPA to LSE were significant, with coefficients of β = 0.234 (t = 5.440, p < 0.001) and β = 0.241 (t = 4.914, p < 0.001), respectively. This underscores the crucial role of teacher support and past experiences in shaping students’ self-efficacy, supporting H1, and H2. Furthermore, LSE demonstrated substantial direct effects on NA, TA, and MA, with coefficients of β = 0.670 (t = 25.249, p < 0.001), β = 0.706 (t = 31.694, p < 0.001), and β = 0.284 (t = 6.704, p < 0.001), respectively. These findings support hypotheses H3A, H3C, and H3, highlighting LSE’s central role in the model. NA and TA also significantly contributed to MA, with coefficients of β = 0.347 (t = 9.063, p < 0.001) and β = 0.251 (t = 5.913, p < 0.001), respectively, thus supporting hypotheses H3C. Notably, the path from CC to MA was also statistically significant with a coefficient of β = 0.110 (t = 2.659, p = 0.008), suggesting that difficulties in working memory directly contribute to increased MA, supporting H4. Finally, MA also significantly influenced LM, as indicated by a coefficient of β = 0.384 (t = 9.939, p < 0.001), supporting H5. The high t-statistics and low p-values across all paths indicate the model’s validity in explaining the multifaceted contributors to MA among the secondary level students in Bangladesh.
3.3.3 Indirect effects
The structured model also evaluated the specific and total indirect effects (see Table 6), highlighting the mediating roles of NA and TA between LSE and MA. A mediating effect is developed when a third variable or construct (in our case, TA and NA) intervenes between two other related constructs (in our case, LSE and MA). To understand the size of the mediating effect, variance accounted for (VAF) is typically suggested by the scholars (48) that determines the size of the indirect effect in relation to the total effect using the following formula (1):
As per the guideline of Hair et al. (48), if the VAF value is less than 20%, and it can be concluded that (almost) no mediation takes place. On the other hand, if the VAF value is larger than 20% and less than 80%, it can be characterized as partial mediation, meaning, a part of the effect of the independent variable (IV) on the dependent variable (DV) is transmitted through the mediator, but the IV still has a direct effect on the DV even after accounting for the mediator. In contrast, when VAF has very large outcomes of above 80%, it can be assuming a full mediation, meaning, indirect effect of an IV on a DV through a mediator is so strong that the direct effect of the IV on the DV becomes non-significant. In this case, the specific indirect effect of LSE on MA through NA was significant, with a coefficient of β = 0.237 (t = 9.787, p < 0.001) and the VAF was 45.4%, indicating that nearly half of the effect of LSE on MA was transmitted through NA, demonstrating a partial mediation effect (48). Similarly, the indirect effect of LSE on MA through TA was also significant, with a coefficient of β = 0.171 (t = 5.795, p < 0.001). The VAF was 37.4%, suggesting that a significant portion of the effect of LSE on MA is mediated by TA, also indicating partial mediation (48). Furthermore, the total indirect effect of LSE on MA, summing the indirect effects through both NA and TA, was substantial, with a coefficient of β = 0.410 (t = 10.232, p < 0.001). The high t statistic and very low p-value confirm the significance of the total indirect effect, emphasizing the substantial impact of LSE on MA through these mediators. So, the findings support the hypotheses H3B and H3D.
3.3.4 Total effects (sum of direct & indirect effects)
As shown in Table 7, Figure 3, and Figure 4, the total effects analysis revealed that all examined relationships were statistically significant, underscoring the key contributors to mathematics anxiety among secondary-level students in Bangladesh. NPE had a significant total effect on LSE, with a coefficient of β = 0.241 (t = 4.914, p < 0.001), indicating that negative past experiences significantly lowered students’ self-efficacy. Similarly, LPTS significantly influenced LSE, with a coefficient of β = 0.234 (t = 5.440, p < 0.001), suggesting that insufficient teacher support also contributed to reduced self-efficacy. The total effect of LSE on MA was substantial, with a coefficient of β = 0.694 (t = 22.695, p < 0.001), highlighting the pivotal role of self-efficacy in the development of MA. CC, specifically working memory difficulties, were also found to significantly impact MA, with a coefficient of β = 0.110 (t = 2.659, p = 0.008), indicating that cognitive difficulties contribute to increased anxiety. LSE was further found to significantly influence both NA and TA, with coefficients of β = 0.670 (t = 25.249, p < 0.001) and β = 0.706 (t = 31.694, p < 0.001), respectively. This suggests that LSE leads to more NA towards mathematics and higher levels of test anxiety. Both NA and TA then significantly contributed to MA, with coefficients of β = 0.354 (t = 9.063, p < 0.001) and β = 0.243 (t = 5.913, p < 0.001), respectively, underscoring their roles in exacerbating MA. Finally, MA significantly affected LM, with a coefficient of β = 0.384 (t = 9.939, p < 0.001), indicating that higher levels of MA correlate with LM to engage in mathematics-related activities.
3.3.5 Coefficient of determination (R2)
The coefficient of determination (R²) for the endogenous constructs is used to evaluate the model’s explanatory power. R² values range from 0 to 1, with higher values indicating greater explanatory strength. As per Hair et al. (49), R² values of 0.75, 0.50, and 0.25 are considered substantial, moderate, and weak, respectively (34, 48). Table 8 presents the R² values for LM, LSE, MA, NA, and TA as 0.147, 0.155, 0.720, 0.449, and 0.492, respectively, demonstrating the model’s satisfactory explanatory power.
3.3.6 PLSpredict (PLS path model estimations)
The Table 9 presents the PLSpredict path model estimations, utilizing the algorithm developed by Shmueli et al. (55), which evaluates predictive performance through Q²predict, RMSE (Root Mean Square Error), and MAE (Mean Absolute Error) metrics. The constructs assessed include LM, LSE, MA, NA, and TA. The Q²predict values, reflecting predictive accuracy, range from 0.100 for LM to 0.208 for MA, with the MA constructs demonstrating the highest predictive performance. The RMSE values, which measure the variability of prediction errors, indicate that the MA model exhibits superior accuracy with the lowest RMSE of 0.893. The MAE values, representing the average magnitude of prediction errors, identify the LSE model as having the smallest error at 0.737. Overall, the MA model is distinguished by its highest Q²predict and lowest RMSE, while the LSE model is notable for its minimal MAE, suggesting that these constructs provide superior predictive performance according to the PLSpredict algorithm.
4 Discussion
This study aimed to identify and analyze the key factors contributing to MA among secondary-level students in Bangladesh, using a Partial Least Squares Structural Equation Modeling (PLS-SEM) approach. The findings offer valuable insights into the complex interplay of factors that exacerbate MA, reinforcing the results of previous research while contributing new evidence from the context of a developing country. The results highlight a moderate but significant impact of cognitive challenges (CC), particularly working memory difficulties, on MA (β = 0.110, t = 2.659, p = 0.008). Although extensive research has established a negative relationship between MA and working memory (56, 57), no prior study has empirically identified working memory difficulties as a predictor of MA. Therefore, our findings can be interpreted as follows: working memory difficulties contribute to poor mathematical performance, which, in turn, may trigger MA, consistent with the Deficit Theory (58). In this sense, our findings align with existing literature, which suggests that cognitive overload impairs mathematical performance, thereby exacerbating anxiety (26, 28). We also found that both Negative Past Experiences (NPE) and Lack of Perceived Teacher Support (LPTS) had significant positive causal effects on Low Self-Efficacy (LSE), with coefficients of β = 0.241 (t = 4.914, p < 0.001) and β = 0.234 (t = 5.440, p < 0.001), respectively. These results indicate that negative past experiences and insufficient teacher support significantly contribute to reduced student self-efficacy. These findings are consistent with those of Shukla et al. (24) and Dowker et al. (3), who reported that negative classroom experiences, such as adverse teacher behavior or failures in mathematics, are negatively correlated with students’ self-efficacy. Furthermore, the strong positive association between low self-efficacy (LSE) and MA (β = 0.694, t = 22.695, p < 0.001), demonstrating that students who doubt their mathematical abilities are more likely to experience MA. Notably, the total indirect effect of LSE on MA, mediated through both Negative Attitude (NA) towards math and Test Anxiety (TA), was substantial (β = 0.410, t = 10.232, p < 0.001). This finding underscores the significant impact of LSE on MA through these mediators (NA and TA), aligning with studies by Khasawneh et al. (10) and Rozgonjuk et al. (20), which highlight the critical role of self-efficacy in managing mathematical anxiety. Additionally, the study identified that both NA toward mathematics and TA independently have a significant impact on MA (math achievement). Finally, this study also revealed that MA has a causal positive relationship with lack of motivation (LM) in math learning. The connection between math anxiety and decreased motivation is also evident in previous studies (59–62). This is because the anxiety associated with math can lead to avoidance behaviors, reduced confidence in their math abilities, and adversely affects their future math performance (60).
5 Recommendations andpolicy implications
Based on the findings of this study, several key recommendations and policy implications are suggested to address and mitigate MA among secondary-level students in Bangladesh:
A. Enhance Teacher Training: Prioritize the development and implementation of teacher training programs that focus on creating supportive classroom environments and employing teaching methods that build student confidence and self-efficacy in mathematics.
B. Invest in Professional Development: Ensure ongoing professional development for teachers, emphasizing the latest research in educational psychology, teaching strategies, and tools for effectively managing classroom dynamics and student anxiety.
C. Real-World Application: Integrate activities and curricula that highlight the practical applications of mathematics in everyday life. Foster a growth mindset by encouraging students to view mathematical challenges as opportunities for growth.
D. Support for Students with Cognitive Challenges: Provide additional support for students facing cognitive challenges, particularly related to working memory. This could involve personalized tutoring, memory aids, and teaching strategies that break down complex problems. Leveraging educational technologies can support cognitive development and offer interactive learning experiences.
E. Comprehensive Support Systems: Develop integrated support systems including counselors, psychologists, and special educators to address the diverse needs of students.
F. Boost Student Motivation: Employ strategies to increase student motivation and engagement in mathematics, such as gamification, setting achievable goals, and providing timely, constructive feedback.
6 Study limitations and future research directions
Despite potential findings, this study acknowledges several limitations. Firstly, a significant limitation is the exclusion of several potential factors, such as gender, parental attitudes, and instructional strategies, due to the research’s limited scope. Including these variables would have complicated the exploration, as examining all such factors in a single study would be challenging, and incorporating mediator factors would add further complexity to the structural model. Secondly, the use of self-report measurements could introduce potential bias, as this method is often criticized for the difficulty respondents may face in making conscious judgments and the tendency to provide socially desirable answers (63). However, we attempted to mitigate this limitation by ensuring students through informed consent that their responses would remain confidential. Thirdly, in this study, we focused primarily on self-reported cognitive difficulties (working memory difficulties) rather than objective cognitive testing. We recognize that the lack of objective cognitive assessments is a limitation, as it would provide a more comprehensive understanding of the participants’ actual cognitive performance. So, in future, objective cognitive assessment is suggested to gain a more nuanced understanding of the cognitive challenges faced by individuals with MA. Fourthly, there is a scarcity of previous studies in developing countries like Bangladesh that explore the complex causal relationships between factors related to MA. This has made it challenging to adequately compare this study’s findings with other studies in the context of Bangladesh. Fifthly, the study’s quantitative cross-sectional design, with data collected over a short period, limits the ability to fully understand and generalize the relationships between the factors. Finally, this study single-level PLS-SEM due to the nature of the data, which reflects individual-level factors influencing MA. However, a multilevel design, especially considering hierarchical relationships (e.g., student-teacher or school-level effects) could be utilized for more precise finings. So, future research should include longitudinal studies, both quantitative and qualitative with multi-level SEM design, to enhance the understanding and generalizability of these relationships.
7 Conclusion
This study unveiled a complex interplay of factors contributing to mathematics anxiety (MA) among Bangladeshi secondary students. By employing a robust quantitative approach, this study demonstrates that negative past experiences, perceived lack of teacher support, and cognitive challenges, such as working memory difficulties, significantly contribute to MA. These factors, combined with low self-efficacy, negative attitudes, and test anxiety further create a complex web that exacerbates MA and undermines students’ motivation and performance in mathematics. Addressing these issues requires a multifaceted approach, including enhanced teacher training, professional development, and the integration of practical applications of mathematics in curricula. Additionally, supporting students with cognitive challenges and fostering comprehensive support systems are essential for mitigating MA. By implementing these recommendations, educational authorities can better support students, improve their mathematical experiences, and ultimately enhance their future performance in mathematics.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by Institutional Research Ethics Board, United International University (IREB/2023/018). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable data included in this article.
Author contributions
SA: Investigation, Resources, Supervision, Validation, Writing – review & editing. JS: Conceptualization, Investigation, Resources, Software, Supervision, Validation, Writing – review & editing. MT: Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft. KA: Writing – review & editing, Methodology, Supervision, Validation. SI: Writing – review & editing, Data curation, Supervision, Conceptualization, Validation.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study received funding from Institute for Advanced Research (IAR), United International University (Grant Reference: UIU-IAR-02-2022-SE-15). The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.
Acknowledgments
The authors extend their gratitude to Bangladesh Digital Education Research Limited (bdeducation), a prominent R&D-focused educational company in Bangladesh, for their support in data collection.
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.
References
1. Anxiety disorders. World health organization. Retrieved august 7, 2024 (2023). Available online at: https://www.who.int/news-room/fact-sheets/detail/anxiety-disorders. (accessed August 7, 2024).
2. Luttenberger S, Wimmer S, Paechter M. Spotlight on math anxiety. Psychol Res Behav Management Volume. (2018) 11:311–22. doi: 10.2147/prbm.s141421
3. Dowker A, Sarkar A, Looi CY. Mathematics anxiety: what have we learned in 60 years? Front Psychol. (2016) 7:508. doi: 10.3389/fpsyg.2016.00508
4. Cipora K, Santos FH, Kucian K, Dowker A. Mathematics anxiety—where are we and where shall we go? Ann New York Acad Sci. (2022) 1513:10–20. doi: 10.1111/nyas.14770
5. Tamal MA, Akter R, Hossain SA, Rezaul KM. A weighted scoring based rating scale to identify the severity level of mathematics anxiety in students. Int J Interactive Mobile Technol (iJIM). (2021) 15:18. doi: 10.3991/ijim.v15i08.18627
6. Zhang J, Zhao N, Kong QP. The relationship between math anxiety and math performance: A meta-analytic investigation. Front Psychol. (2019) 10. doi: 10.3389/fpsyg.2019.01613
7. Weissgerber SC, Grünberg C, Neufeld L, Steppat T, Reinhard MA. The interplay of math anxiety and math competence for later performance. Soc Psychol Educ. (2022) 25:977–1002. doi: 10.1007/s11218-022-09700-y
8. Saha J, Ahmmed S, Tamal MA. Identifying mathematics anxiety through automated process (IMATAP): A web-based screening tool for teachers. Int J Emerging Technol Learn (iJET). (2023) 18:186–98. doi: 10.3991/ijet.v18i21.40325
9. Szűcs D. The origins of math anxiety and interventions. In: IBE — Science of learning portal (2019). Available at: https://solportal.ibe-unesco.org/articles/the-origins-of-math-anxiety-and-interventions/.
10. Khasawneh E, Gosling C, Williams B. What impact does maths anxiety have on university students? BMC Psychol. (2021) 9. doi: 10.1186/s40359-021-00537-2
11. Vos H, Marinova M, De Léon SC, Sasanguie D, Reynvoet B. Gender differences in young adults’ mathematical performance: Examining the contribution of working memory, math anxiety and gender-related stereotypes. Learn Individ Dif. (2023) 102:102255. doi: 10.1016/j.lindif.2022.102255
12. Szczygiel M. Gender, general anxiety, math anxiety and math achievement in early school-age children. Issues Educ Res. (2020) 30:1126–42. https://www.iier.org.au/iier30/szczygiel.pdf.
13. Wahid SNS, Yusof Y, Nor AHM. Effect of mathematics anxiety on students’ performance in higher education level: A comparative study on gender. AIP Conference Proceedings (2018) 1974:050010. doi: 10.1063/1.5041710
14. Erturan S, Jansen B. An investigation of boys’ and girls’ emotional experience of math, their math performance, and the relation between these variables. Eur J Psychol Educ. (2015) 30:421–35. doi: 10.1007/s10212-015-0248-7
15. Global Education Monitoring Report – Gender Report: Deepening the debate on those still left behind. UNESCO eBooks. Place de Fontenoy, Paris, France: Global Education Monitoring Report Team (2022). doi: 10.54676/rczb6329.
16. Van Mier HI, Schleepen TMJ, Van den Berg FCG. Gender differences regarding the impact of math anxiety on arithmetic performance in second and fourth graders. Front. Psychol. (2019) 9:2690. doi: 10.3389/fpsyg.2018.02690
17. Malanchini M, Rimfeld K, Wang Z, Petrill SA, Tucker-Drob EM, Plomin R, et al. Genetic factors underlie the association between anxiety, attitudes and performance in mathematics. Trans Psychiatry. (2020) 10. doi: 10.1038/s41398-020-0711-3
18. Wang Z, Hart SA, Kovas Y, Lukowski S, Soden B, Thompson LA, et al. Who is afraid of math? Two sources of genetic variance for mathematical anxiety. J Child Psychol Psychiatry. (2014) 55:1056–64. doi: 10.1111/jcpp.12224
19. Hiller SE, Kitsantas A, Cheema JE, Poulou M. Mathematics anxiety and self-efficacy as predictors of mathematics literacy. Int J Math Educ Sci Technol. (2021) 53:2133–51. doi: 10.1080/0020739x.2020.1868589
20. Rozgonjuk D, Kraav T, Mikkor K, Orav-Puurand K, Täht K. Mathematics anxiety among STEM and social sciences students: the roles of mathematics self-efficacy, and deep and surface approach to learning. Int J Stem Educ. (2020) 7. doi: 10.1186/s40594-020-00246-z
21. Wang C, Xu Q, Fei WQ. The effect of student-perceived teacher support on math anxiety: chain mediation of teacher–student relationship and math self-efficacy. Front Psychol. (2024) 15:1333012. doi: 10.3389/fpsyg.2024.1333012
22. Yang Y, Li G, Su Z, Yuan Y. Teacher’s emotional support and math performance: the chain mediating effect of academic self-efficacy and math behavioral engagement. Front Psychol. (2021) 12:651608. doi: 10.3389/fpsyg.2021.651608
23. Zhou D, Du X, Hau KT, Luo H, Feng P, Liu J. Teacher-student relationship and mathematical problem-solving ability: mediating roles of self-efficacy and mathematical anxiety. Educ Psychol. (2019) 40:473–89. doi: 10.1080/01443410.2019.1696947
24. Shukla KD, Kuril S, Chand VS. Does negative teacher behavior influence student self-efficacy and mastery goal orientation? Learn Motivation. (2020) 71:101653. doi: 10.1016/j.lmot.2020.101653
25. Scheibe DA, Was CA, Dunlosky J, Thompson CA. Metacognitive cues, working memory, and math anxiety: the regulated attention in mathematical problem solving (RAMPS) framework. J Intell. (2023) 11(6):117. doi: 10.3390/jintelligence11060117
26. Finell J, Sammallahti E, Korhonen J, Eklöf H, Jonsson B. Working memory and its mediating role on the relationship of math anxiety and math performance: A meta-analysis. Front Psychol. (2022) 12:798090. doi: 10.3389/fpsyg.2021.798090
27. Cuder A, Živković M, Doz E, Pellizzoni S, Passolunghi MC. The relationship between math anxiety and math performance: The moderating role of visuospatial working memory. J Exp Child Psychol. (2023) 233:105688. doi: 10.1016/j.jecp.2023.105688
28. Ashcraft MH, Krause JA. Working memory, math performance, and math anxiety. Psychonomic Bull Rev. (2007) 14:243–8. doi: 10.3758/bf03194059
29. Ashcraft MH, Kirk EP. The relationships among working memory, math anxiety, and performance. J Exp Psychol Gen. (2001) 130:224–37. doi: 10.1037/0096-3445.130.2.224
30. Norton A, Seok Y, Choi-Koh S. Examining mathematics anxiety of undergraduates using a brain-based measurement, EEG. J Behav Brain Sci. (2019) 09:195–209. doi: 10.4236/jbbs.2019.95017
31. Zhang S. Math anxiety: the influence of teaching strategies and teachers’ Attitude. (Rue de Palestro, Paris France: Atlantis Press SARL) (2022). pp. 202–8. doi: 10.2991/978-2-494069-13-8_27.
32. Saha J, Ahmmed S, Ali M, Tamal MA, Rezaul KM. ICT based mathematics skill development program: an initiative to overcome mathematics anxiety. Int J Emerging Technol Learn (iJET). (2020) 15:252. doi: 10.3991/ijet.v15i14.14149
33. Setia MS. Methodology series module 3: Cross-sectional studies. Indian J Dermatol. (2016) 61:261. doi: 10.4103/0019-5154.182410
34. Antora M, Islam MK, Sattar A, Sarker MFH. Factors predicting the Major depressive disorder of engineering students in Bangladesh: a structural equation modeling approach. J Hum Behav Soc Environ. (2024) 34:1–19. doi: 10.1080/10911359.2024.2330461
35. Wang X, Cheng Z. Cross-sectional studies. Chest. (2020) 158:S65–71. doi: 10.1016/j.chest.2020.03.012
36. Bangladesh Education Statistics. Bangladesh bureau of educational information and statistics (BANBEIS) (2023). Available online at: http://www.banbeis.gov.bd/ (Accessed April 3, 2024).
37. Ahmmed S, Saha J, Tamal MA. Effectiveness of need-based teacher’s training program to enhance online teaching quality. Educ Res Int. (2022) 2022:1–13. doi: 10.1155/2022/4118267
38. Ahmmed S, Saha J, Tamal MA. An empirical study for determining the quality indicators for the primary and secondary school of Bangladesh: A structural equation modeling approach. Heliyon. (2022) 8:e10870. doi: 10.1016/j.heliyon.2022.e10870
39. Kock N, Hadaya P. Minimum sample size estimation in PLS-sem: The inverse square root and gamma-exponential methods. Inf Syst J. (2016) 28:227–61. doi: 10.1111/isj.12131
40. Kock N. Should bootstrapping be used in PLS-sem? toward stable P-value calculation methods. J Appl Struct Equation Modeling. (2018) 2:1–12. doi: 10.47263/JASEM.2(1)02
41. Ned K. Minimum sample size estimation in PLS-SEM: an application in tourism and hospitality research. In: Applying partial least squares in tourism and hospitality research. Bingley, United Kingdom: Emerald Group Publishing (2018), 1–16. Available at: https://www.emerald.com/insight/publication/doi/10.1108/9781787566996:~:text=ISBN-,978%2D1%2D78756%2D700%2D9,-eISBN.
42. Carey E, Hill F, Devine A, Szűcs D. The modified abbreviated math anxiety scale: A valid and reliable instrument for use with children. Front Psychol. (2017) 8:11. doi: 10.3389/fpsyg.2017.00011
43. Jirjees F, Odeh M, Al-Haddad A, Ass’ad R, Hassanin Y, Al-Obaidi H, et al. Test anxiety and coping strategies among university students an exploratory study in the UAE. Sci Rep. (2024) 14. doi: 10.1038/s41598-024-59739-4
44. McHugh ML. Interrater reliability: the kappa statistic. Biochemia Med. (2012) 22:276–82. doi: 10.11613/bm.2012.031
45. Tavakol M, Dennick R. Making sense of Cronbach’s alpha. Int J Med Educ. (2011) 2:53–5. doi: 10.5116/ijme.4dfb.8dfd
46. IBM Corp. (Version 26). Armonk, New York: IBM Corp (2019). Available at: https://www.ibm.com/products/spss-statistics.
47. Ringle CM, Wende S, Becker J-M. SmartPLS 4. Oststeinbek: SmartPLS (2022). Available at: https://www.smartpls.com.
48. Hair J, Ringle HMC, Sarstedt M. A primer on partial least squares structural equation modeling (PLS-SEM) (2nd ed.). United States of America: SAGE (2016).
49. Hair J, Risher J, Sarstedt M, Ringle C. When to use and how to report the results of PLS-SEM. Eur Business Rev. (2019) 31:2–24. doi: 10.1108/EBR-11-2018-0203
50. Hair JF, Hult GTM, Ringle CM, Sarstedt M, Danks NP, Ray S. Evaluation of reflective measurement models. In Classroom companion: business. (2021) pp:75–90). doi: 10.1007/978-3-030-80519-7_4
51. Taber KS. The use of cronbach’s alpha when developing and reporting research instruments in science education. Res Sci Educ. (2018) 48:1273–96. doi: 10.1007/s11165-016-9602-2
52. Henseler J, Ringle CM, Sarstedt M. A new criterion for assessing discriminant validity in variance-based structural equation modeling. J Acad Marketing Sci. (2014) 43:115–35. doi: 10.1007/s11747-014-0403-8
53. Hair JF, Hult GTM, Ringle CM, Sarstedt M. A primer on partial least squares structural equation modeling (PLS-SEM) (2nd ed.). Sage: Thousand Oaks (2017).
54. Fornell C, Larcker D. Evaluating structural equation models with unobservable variables and measurement error. J Marketing Res. (1981) 18:39. doi: 10.1177/002224378101800104
55. Shmueli G, Ray S, Estrada JMV, Chatla SB. The elephant in the room: predictive performance of PLS models. J Business Res. (2016) 69:4552–64. doi: 10.1016/j.jbusres.2016.03.049
56. Pellizzoni S, Cargnelutti E, Cuder A, Passolunghi MC. The interplay between math anxiety and working memory on math performance: a longitudinal study. Ann New York Acad Sci. (2021) 1510:132–44. doi: 10.1111/nyas.14722
57. Holmes J, Gathercole SE, Dunning DL. Poor working memory. Adv Child Dev Behav. (2010) 39:1–43. doi: 10.1016/b978-0-12-374748-8.00001-9
58. Carey E, Hill F, Devine A, Szücs D. The chicken or the egg? The direction of the relationship between mathematics anxiety and mathematics performance. Front Psychol. (2016) 6:1987. doi: 10.3389/fpsyg.2015.01987
59. Pollack C, Wilmot D, Centanni TM, Halverson K, Frosch I, D’Mello AM, et al. Anxiety, motivation, and competence in mathematics and reading for children with and without learning difficulties. Front Psychol. (2021) 12:704821. doi: 10.3389/fpsyg.2021.704821
60. Jenifer JB, Rozek CS, Levine SC, Beilock SL. Effort(less) exam preparation: Math anxiety predicts the avoidance of effortful study strategies. J Exp Psychol Gen. (2022) 151:2534–41. doi: 10.1037/xge0001202
61. Zakaria E, Nordin NM. The effects of mathematics anxiety on matriculation students as related to motivation and achievement. Eurasia J Mathematics Sci Technol Educ. (2008) 4:27–30. doi: 10.12973/ejmste/75303
62. Süren N, Kandemir MA. The effects of mathematics anxiety and motivation on students’ Mathematics achievement. Int J Educ Mathematics Sci Technol. (2020) 8:190. doi: 10.46328/ijemst.v8i3.926
Keywords: mathematics anxiety, mathematics performance, factors of mathematics anxiety, partial least squares structural equation modeling (PLS-SEM), adolescents in Bangladesh
Citation: Ahmmed S, Saha J, Tamal MA, Abdullah Al Mamun K and Islam S (2024) Factors predicting the mathematics anxiety of adolescents: a structural equation modeling approach. Front. Psychiatry 15:1484381. doi: 10.3389/fpsyt.2024.1484381
Received: 21 August 2024; Accepted: 07 November 2024;
Published: 05 December 2024.
Edited by:
Takahiro Nemoto, Toho University, JapanReviewed by:
Larry R. Price, Texas State University, United StatesDominika Jarčušková, University of Pavol Jozef Šafárik, Slovakia
Copyright © 2024 Ahmmed, Saha, Tamal, Abdullah Al Mamun and Islam. 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.
*Correspondence: Maruf Ahmed Tamal, dGFtYWwyNS0xMDNAZGl1LmVkdS5iZA==