Mathematics interest, self-efficacy, and anxiety predict STEM career choice in emerging adulthood

Introduction
The United Kingdom, like many nations1,2,3, is experiencing an ongoing shortage of people with Science, Technology, Engineering and Mathematics (STEM) expertise, with current estimates indicating a shortfall of over 173,000 workers across the sector4. A national inquiry has highlighted a need to address this skills gap to meet an expected increase in future demand and ensure sustainable economic growth5. Additionally, this issue may be worsened by the continued underrepresentation of women, with 27% of STEM workers being female, compared to 52% of the wider workforce6. Understanding why individuals may or may not choose to pursue a STEM career could be crucial in supporting the development of a globally competitive workforce.
While maths achievement in adolescence has long been demonstrated to predict STEM participation and career choice in adulthood7,8, research has shown that socio-cognitive and emotional factors may be as crucial for STEM career decision-making during this time9,10. Maths anxiety and maths motivational attitudes, such as maths self-efficacy and maths interest, are separately associated with student STEM participation9,11,12,13, although few studies have considered their combined impact on STEM career choice. Here, we tested the relationships between maths anxiety, maths self-efficacy, maths interest and STEM career participation in emerging adulthood using data from a UK-wide cohort.
Social cognitive career theory (SCCT) is a theory of motivation concerned with explaining the socio-cognitive processes behind career-related intentions and choices14. SCCT highlights self-efficacy and interest as key motivational attitudes influencing career-related behaviour14,15. In mathematics, self-efficacy refers to an individual’s perceived competence16, while interest is a motivational construct describing positive affect and orientation towards maths-related activities17. Both have been consistently linked to maths achievement2,3,18 and, in the case of maths self-efficacy, socioeconomic status19,20. In line with SCCT, studies confirm that maths self-efficacy and interest predict STEM participation. For example, a longitudinal study tracking 14,000 US students found that maths self-efficacy predicted enrolment intentions on STEM-related courses, independent of prior maths achievement and socioeconomic status21. Both maths self-efficacy and interest were also found to predict STEM course intentions in an Australian sample after controlling for maths achievement13. However, support for SCCT is not always consistent. Some studies report differential effects for each motivational construct. For example, one longitudinal study found a reciprocal relationship between maths interest and STEM career intentions over time but did not find a relationship between maths self-efficacy and career intentions2. Additionally, current research mainly focuses on career intentions rather than actual career choices and does not sufficiently address the emotional factors that could impact STEM-related behaviour.
Another aspect of mathematics linked to STEM career interest is maths anxiety. Maths motivational attitudes and maths anxiety are related16,22, yet distinct23. Maths anxiety, distinct from general anxiety, is characterised by negative emotional affect and tension at the prospect of, or during, mathematical tasks, and negatively correlates with maths achievement (average r = −0.30)22,23,24. Research also suggests that students from lower SES backgrounds are at increased risk of experiencing maths anxiety24,25. The interplay between maths anxiety and maths motivational attitudes is complex: while high maths anxiety often co-occurs with lower achievement, several studies have shown that motivational attitudes, including self-efficacy and interest, can buffer the adverse effects of maths anxiety on achievement26,27. For example, it has been demonstrated that, in individuals with low overall maths motivation, there is a negative relationship between maths anxiety and maths performance27. Conversely, in individuals with high maths motivation levels, this relationship becomes an inverted U, suggesting that a moderate level of maths anxiety might be performance-enhancing due to heightened alertness and focused attention before it becomes debilitating at higher levels27,28.
Research has also shown that, beyond the influence of previous maths achievement, students with high levels of maths anxiety are more likely to avoid STEM participation and are less inclined to pursue STEM-based careers11,12,29,30. Despite the existing literature on the non-linear joint effects of maths anxiety and maths motivational attitudes on maths achievement, few studies have considered how these factors might jointly influence STEM career choice. One study of high school students found that while maths anxiety did predict students’ self-perceptions of maths ability, maths performance expectancies and value perceptions, it did not predict STEM course enrolment intentions when accounting for these motivational constructs31. The study also found that, of these constructs, only value perceptions of maths predicted enrolment intentions31. This implies that the relationship between maths anxiety and STEM participation may be more complex than has been previously implied. It also suggests that different aspects of maths motivational attitudes may vary in their influence on STEM participation. However, as only course intentions were considered as an outcome, rather than actual STEM participation behaviour, this study’s relevance to understanding why individuals pursue STEM careers in adulthood may be limited. By investigating the joint influence of maths anxiety, maths self-efficacy and maths interest on individuals’ STEM career choices in emerging adulthood, the current investigation aims to address a gap in current research. Sex differences in maths anxiety, maths self-efficacy, maths interest and STEM participation have been consistently demonstrated across adolescence and emerging adulthood. Male participants tend to report higher maths self-efficacy and interest levels than female participants9,13,32,33,34. Conversely, female participants tend to report higher levels of maths anxiety compared to male participants11,12,16,30. They are also less likely to study STEM subjects in higher education12,35,36,37 and commit to STEM-based careers35 despite exhibiting similar levels of previous maths achievement33,38. Research exploring the influence of maths motivational attitudes on STEM participation and career choice reports very few sex differences in the relationship21,35,39. However, one study noted that maths self-efficacy was associated with STEM career interest in females but not males, whereas maths interest was associated for males but not females13. Fewer studies have considered sex differences in the impact of maths anxiety on STEM career choice, however, longitudinal research has shown that, in males, but not females, higher maths anxiety predicts lower future maths self-efficacy37. This finding may indicate that, like maths self-efficacy, the STEM career choices of males may be more negatively influenced by previous maths anxiety than those of females.
To mitigate the STEM skills gap in the UK, it is essential to understand the psychological factors shaping career choices in this sector. Our study employs longitudinal data from a large, nationally representative cohort to examine how maths anxiety, maths self-efficacy and maths interest affect STEM career selection in early adulthood. In this investigation, we aim to address the following questions: (1) What is the joint influence of maths anxiety (with and without controlling for general anxiety), maths self-efficacy and maths interest on STEM career choice in emerging adulthood? (2) Are there sex differences in the joint influence of maths anxiety, maths self-efficacy and maths interest on STEM career choice? Informed by SCCT and previous research work, we expected that maths self-efficacy and maths interest would each positively and independently predict STEM career choice. We also expected that maths anxiety would independently and negatively predict STEM career participation. Given the established associations between these factors and maths achievement, as well as socioeconomic status, we expected that these effects would be partially attenuated when controlling for these variables. Based on prior work considering sex differences, we expected that maths anxiety and maths interest would have stronger associations with STEM career choice in the male cohort than the female cohort and that maths self-efficacy would have a stronger association with STEM career choice in females than males. Our research plan was preregistered in the Open Science Framework (OSF; https://osf.io/y36bg).
Results
Descriptive statistics and mean-level sex differences
The descriptive statistics for continuous variables are presented in Supplementary Table 1. Frequency count and proportion statistics for STEM career choice, a binary variable, are presented in Supplementary Table 2. 2254 participants reported a degree and/or apprenticeship choice. 5.1% (117) reported an apprenticeship choice and 29.4% (662) reported a STEM career choice. Participants who reported a STEM career choice had a significantly higher mean GCSE mathematics achievement (M = 9.5, SD = 1.17) than the whole sample (M = 8.9, SE = 1.44; t(4599) = −24.74, p < 0.001).
Univariate ANOVAs were used to test for mean-level sex differences across all continuous variables, as shown in Supplementary Table 3. Small but significant differences emerged, with sex explaining a small to moderate proportion of the variance across measures (η² = 0.001–0.06). Females had significantly higher levels of maths anxiety (F(948, 541) = 86.68, p < 0.001, η² = 0.06) and general anxiety (F(948, 541) = 57.74, p < 0.001, η² = 0.04), while males showed significantly higher levels of maths self-efficacy (F(1485, 1074) = 149.10, p < 0.001, η² = 0.05) and maths interest (F(1485, 1074) = 39.04, p < 0.001, η² = 0.01). Although males also had significantly higher maths GCSE achievement and SES levels, sex explained 0.1% and 0.007% of the variance, respectively. All variables were therefore adjusted for sex and age using a regression method. The standardised residuals corrected for age and sex were used for all downstream analyses.
Supplementary Table 4 presents contingencies for sex differences in the STEM career choice. 63.9% of the sample were female. 21.9% of female participants chose a STEM career, and 42.5% of male participants chose a STEM career. A Pearson’s Chi-square test and an odds ratio calculation were performed to determine whether the proportion of STEM career choice differed between males and females and the effect size of this difference. The proportion of STEM career choices was found to differ significantly by sex, with males being more than 2.5 times as likely to make a STEM career choice than females, (X2(1, N = 2254) = 105, p < 0.001, OR = 2.63).
The joint influence of maths anxiety, maths self-efficacy, and maths interest on STEM career choice in emerging adulthood
Figure 1 presents Pearson’s correlations between all variables. Supplementary Table 5 reports all correlations with sample sizes and 95% confidence intervals. A negative correlation was observed between maths anxiety and STEM career choice (r = −0.23). A small positive correlation was observed between maths motivation and STEM career choice (r = 0.28). Small associations were observed between all other variables and STEM career choice (r ranging between −0.09 and 0.28). Maths anxiety had a small negative association with maths achievement (r = −0.29) and was moderately correlated to the same extent with both maths self-efficacy and maths interest (r = −0.41). Maths achievement was highly correlated with maths self-efficacy (r = 0.65) and moderately associated with maths interest (r = 0.41).

Zero-order correlations between all variables. MA maths anxiety, MA.GA maths anxiety (independent of general anxiety), MSE maths self-efficacy, MI maths interest, GA general anxiety, GCSE GCSE maths achievement, SES socioeconomic status. *p < 0.05, **p < 0.01, ***p < 0.001.
There was a moderate positive correlation between maths anxiety and general anxiety (r = 0.34). The association between general anxiety and STEM career choice was weak (r = −0.09), and there was no significant difference observed between the correlations for maths anxiety and STEM career choice (r = −0.23, 95% CI [−0.29, −0.17]) and maths anxiety (independent of general anxiety) and STEM career choice (r = −0.21, 95% CI [−0.27, −0.15]). We, therefore, used the maths anxiety variable uncorrected for general anxiety for the regression models.
We conducted three binary logistic regression models to assess the unique contributions of maths anxiety and maths motivation to STEM career choice. In the first model, maths anxiety, maths self-efficacy, and maths interest were jointly assessed. In the second model, maths anxiety, maths self-efficacy and maths interest were adjusted for GCSE maths achievement. In the third model, maths anxiety, maths self-efficacy and maths interest were adjusted for maths achievement and SES. Variance inflation factors (VIF) for each predictor variable in each model were all less than 1.6, meeting the logistic regression assumption of low or no multicollinearity between predictors (see Supplementary Table 6 for full VIF for each predictor variable).
The Nagelkerke pseudo-R-squared values for all models ranged between 0.15 and 0.21. Supplementary Table 7 summarises the findings of each regression model. Figure 2 visualises the relative effects of maths anxiety and maths motivation across all three models. Results from the first model show that maths anxiety (OR = 0.78, 95% CI [0.64, 0.96], p = 0.02) and maths interest (OR = 1.75, 95% CI [1.44, 2.67], p < 0.001) significantly and independently predicted STEM career choice. Maths self-efficacy did not independently predict STEM career choice (OR = 1.25, 95% CI [0.96, 1.62], p = 0.09). Maths interest remained a significant independent predictor of STEM career choice when both maths achievement and SES were adjusted for, although its effects were attenuated (OR = 1.65, 95% CI [1.33, 2.04], p < 0.001). Results from the second model indicate that maths anxiety does not independently predict STEM career choice when the effects of maths achievement are controlled. However, effect sizes in the adjusted models were only slightly attenuated compared to the unadjusted model (ORs ranging between 0.79 and 0.81).

Relative ORs and 95% CIs of maths anxiety, maths self-efficacy and maths interest on STEM career choice across three binary logistic regression models. MA maths anxiety, MSE maths self-efficacy, MI maths interest, GCSE GCSE maths achievement, SES socioeconomic status.
Sex differences in the joint influence of maths anxiety, maths self-efficacy, and maths interest on STEM career choice
Figure 3 presents the sex-stratified Pearson’s correlations between all variables. Supplementary Table 8 presents the sample size and 95% confidence intervals for all correlations. The associations between maths interest and STEM career choice were similar in male and female samples (r = 0.27 and r = 0.30, respectively). The association between maths self-efficacy was greater in the female sample (r = 0.24, 95% CI [0.17, 0.31]) than in the male sample (r = 0.16, 95% CI [0.06, 0.25]), although overlapping confidence interval suggests this difference is not significant. A small negative correlation was observed between maths anxiety and STEM career choice in the female sample (r = −0.21, 95% CI [−0.28, −0.13]). In contrast, the negative correlation was moderate in the male sample (r = −0.32, 95% CI [−0.42, −0.21]), although the overlapping confidence intervals again demonstrate that this difference was not significant. Socioeconomic status (SES) was associated with STEM career choice in the female sample (r = 0.15, 95% CI [0.10, 0.20]) but not in the male sample (r = 0.05, 95% CI [−0.02, 0.11]), although the overlapping confidence intervals show that the effect sizes were not significantly different from each other.

Zero-order correlations between all variables in sex-stratified samples. MA maths anxiety, MA.GA maths anxiety (independent of general anxiety), MSE maths self-efficacy, MI maths interest, GA general anxiety, GCSE maths achievement, SES socioeconomic status. *p < 0.05 **p < 0.01, ***p < 0.001.
To further investigate sex differences in the joint influence of maths anxiety and maths motivation on STEM career choice, we repeated the three binary logistic regression models in sex-stratified samples. As general anxiety was not significantly associated with STEM career choice in either sample, the unadjusted maths anxiety variable was used for both samples. Variance inflation factors (VIF) for predictor variables in each sample were all less than 1.6, meeting the logistic regression assumption of low or no multicollinearity between predictors (see Supplementary Tables 9 and 10 for full VIF for each predictor variable).
Supplementary Table 11 summarises the effects observed for the logistic regression models in the sex-stratified samples. Figure 4 compares the effects of maths anxiety, maths self-efficacy and maths interest on STEM career choice in male and female samples across all three models. The Nagelkerke pseudo-R-squared for models in the female sample ranged between 0.16 and 0.23. In the female sample, maths interest significantly predicted STEM career choice across all models (ORs ranging from 1.63 to 1.77) and maths self-efficacy predicted STEM career choice in the unadjusted model only (OR = 1.77). Maths anxiety did not independently predict STEM career choice in any model performed in the female sample. In the male sample, the Nagelkerke pseudo-R-squared ranged between 0.18 and 0.19. Both maths interest (ORs ranging from 1.71 to 1.88) and maths anxiety (ORs ranging from 0.62 to 0.64) predicted STEM career choice in every model in the male sample. Maths self-efficacy did not predict STEM career choice in any model performed in the male sample.

Comparison of the joint effects of maths anxiety, maths self-efficacy, and maths interest on STEM career choice in sex-stratified logistic regression models. MA maths anxiety, MSE maths self-efficacy, MI maths interest, GCSE GCSE maths achievement, SES socioeconomic status.
Maths anxiety was an independent predictor of STEM career choice in the unadjusted male model (OR = 0.62, 95% CI [0.42, 0.89]) but not in the unadjusted female model (OR = 0.88, 95% CI [0.67, 1.16]). A Wald Chi-square test showed that maths anxiety was not a significantly better predictor of STEM career choice in the male sample than in the female sample (X2 = 2.26, p = 0.13). However, the effect of maths self-efficacy in the unadjusted female model (OR = 1.77, 95% CI [1.20, 2.60]) was larger than the effect observed in the unadjusted male model (OR = 0.98, 95% CI [0.65, 1.50]), and this difference was found to be significant (X2 = 4.07, p = 0.043).
Discussion
Previous research has highlighted the influence of maths anxiety, maths self-efficacy and maths interest on STEM participation9,11,12,13,35 but has not considered the joint effects of these factors on individuals’ STEM career choices. The results of the present study suggest that maths anxiety and maths interest in adolescence independently predict STEM career choices in emerging adulthood. We also showed that while the independent effects of maths interest remained significant when both SES and maths achievement were accounted for, the already weak effect of maths anxiety no longer predicted STEM career choice after accounting for maths achievement. This finding challenges the prevalent narrative that maths anxiety invariably leads to STEM avoidance11,12,40, and underscores the complexity of the pathway from maths anxiety to STEM career choice. However, it is consistent with existing work considering the joint effects of maths anxiety and maths motivational attitudes on general STEM participation31.
Maths self-efficacy at age 16 did not independently predict STEM career choice in our study. This finding contrasts with the expectations of Social Cognitive Career Theory (SCCT), which proposes that self-efficacy in a given domain is a key driver for later career-related actions14. However, we observed a significant and substantial correlation between maths self-efficacy and maths interest (r = 0.52). A longitudinal study of students in primary education has demonstrated that prior maths self-efficacy predicts later maths interest3, suggesting that, consistent with SCCT, self-efficacy may impact later STEM career outcomes indirectly through its influence on developing interest. The same study also showed that prior maths anxiety also predicted later maths interest3, indicating that maths anxiety may also further impact STEM career choice through its impact on maths interest earlier in development. Although our measures of maths self-efficacy, maths interest, and maths anxiety were taken at a single time point, these relationships highlight the need for longitudinal research over adolescence and early adulthood to fully understand how these constructs interact over time to influence STEM career choices in emerging adulthood.
Our results provide some support for the prediction that maths anxiety might be more decisive for STEM career choice in males than in females30, as this study found that maths anxiety was a significant predictor of STEM career choice in the male sample only. This is consistent with a recent study that found that maths anxiety predicted male STEM participation to a greater extent than it did female STEM participation12. Male students have been shown to identify with maths more than females and feel more pressure to perform well41, so high maths anxiety might threaten their identity more, potentially leading to reduced STEM participation. However, we also observed that the difference in the effects of maths anxiety by sex was non-significant, suggesting that the evidence for this sex-specific impact is weak.
We observed that maths interest was a robust predictor of STEM career choice for both males and females, but this investigation also found that maths self-efficacy was a significant predictor of STEM career choice in the female, but not male sample. This is consistent with a previous study, which also found that maths self-efficacy predicted STEM participation in females but not males13 and suggests that self-perceptions of competency are more important for female STEM career participation than male. However, we also found that this effect was attenuated when accounting for maths achievement, showing that the independent association between maths self-efficacy and female STEM participation can be explained by maths achievement levels.
Although only slight sex differences were observed in the joint effects of maths anxiety, maths interest and maths self-efficacy on STEM career choice, analyses of mean-level differences demonstrate that consistent with previous research12,35,36, male participants are more than twice as likely to choose a STEM career than female participants. Sex accounts for only 0.1% of the variance in maths achievement, suggesting that sex differences in maths achievement are unlikely to account for much of this difference. Females had significantly higher maths anxiety, with sex explaining 6% of the variance, in line with prior work11,16,24. Although higher female maths anxiety is thought to reduce STEM participation11,42,43, the weak effects of this construct on STEM career choice demonstrated here, particularly in the female-only sample, suggest that this may not be the case. Conversely, consistent with previous studies13,33,34, maths interest was found to be significantly higher in males, with sex explaining 1% of the variance. Maths interest was also found to be a robust predictor of STEM career choice in both male and female samples, making it a stronger potential candidate than maths anxiety for at least partially explaining the large sex disparity in STEM career choice.
There are some limitations to this study. Robustness checks (see Supplementary Table 12) demonstrated that data missing across all variables were from individuals with a lower mean SES level than those who contributed data. In particular, the maths achievement and STEM career choice variables were shown to have a moderate effect of SES on missing data. This issue is likely to be further affected by the career data availability being limited to degree and apprenticeship decisions, as previous work demonstrates that these intakes tend to favour individuals from higher SES backgrounds in the UK44. While we observed a small effect of SES on STEM career choice, our investigation also highlighted that SES had little attenuating effect on the influence of maths anxiety, maths self-efficacy and maths interest on STEM career choice beyond that of maths achievement. It is possible that this may not be the case for individuals who come from particularly low SES backgrounds. Replicating the results presented in this study in populations with a full range of SES, perhaps by using data for career choices beyond degrees and apprenticeships, would allow for greater generalisability.
A further limitation of this study is that the measures of maths self-efficacy and maths interest were taken two years earlier than the maths anxiety measure. However, longitudinal research has found moderate to high levels of stability in maths interest, maths self-efficacy2,45 and maths anxiety46, suggesting that this two-year time lapse may not significantly impact the findings. Nonetheless, as mentioned earlier, it is crucial to replicate these findings with a more comprehensive understanding of how maths anxiety, maths self-efficacy and maths interest develop over adolescence and into early adulthood.
Finally, while the associations between maths anxiety, maths self-efficacy, maths interest and STEM career choice are evident, this study design does not allow for causal inference. The observed relationships between these variables could be influenced by a third variable that this investigation has not accounted for. It has been suggested that emotional and motivational factors specific to STEM more broadly, rather than maths only (e.g., STEM anxiety and STEM interest and self-efficacy), may also influence the likelihood of STEM career choices21,24. Additionally, the STEM anxiety framework24 highlights how wider social, emotional and motivational factors, such as parental influence, school or teacher influence, financial or career motivation levels, and self-regulation ability, are likely to have important roles to play. Future research should consider whether these proposed factors change the associations between maths anxiety, maths self-efficacy, maths interest and STEM career choice.
In conclusion, our study adds to our understanding of how maths-related motivational and emotional factors in adolescence contribute to STEM career choice in emerging adulthood. It highlights the potential of enhancing maths interest in adolescence over merely reducing maths anxiety and increasing self-efficacy, especially considering the effect of maths interest for both males and females. Our findings suggest developing targeted interventions that foster STEM participation, aligning with the UK’s economic need for increased STEM expertise. The complex interplay of maths anxiety and maths motivational attitudes and their development over time warrants further longitudinal research to accurately inform targeted educational strategies and policies.
Methods
Participants
Participants were from the Twins Early Development Study (TEDS), a longitudinal study of twins born between 1994 and 1996 in England and Wales. TEDS participants were representative of the wider population of England and Wales in terms of ethnicity and socioeconomic status at conception, and although there has been some attrition over time (a 76.5% response rate by early adulthood), the sample remains representative of the birth cohort. More detailed information on the TEDS sample is available elsewhere47. Informed consent was obtained from participants at each wave of data collection, and the study received ethical approval from the King’s College London Ethics Committee (PNM/09/10-104). This investigation used data collected in a subsample of TEDS twins over four waves: first contact (when twins were ~1.5 years old), age ~16 (mean = 16.31, SD = 0.29), age ~18–21 (mean = 20.30, SD = 0.48), and age ~21–24 (mean = 22.28, SD = 0.91). One twin from each pair was randomly selected for analyses to control for non-independence of observations.
At first contact, the parents of twin participants contributed data on their families’ socioeconomic status (SES; N = 7345, 53% female). SES was taken at first contact as this wave offered minimal missing data, and SES has been shown to be developmentally stable, with age-to-age correlations ~0.7048. At age 16, participants contributed data on maths achievement (N = 6409, 52% female) and maths self-efficacy and interest (N = 2558, 58% female). At age 18–21, data on general anxiety and maths anxiety were collected from 2 out of 4 TEDS cohorts (N = 1489, 64% female). At age 21–24, participants contributed data on their choice of degree and/or apprenticeship (N = 2254, 64% female). 94.8% (2137) of individuals reported a degree choice only, and 4.7% (105) reported an apprenticeship choice only. 0.5% (12) reported both a degree and an apprenticeship choice. Sample sizes varied across variables due to participation levels and data collection differences across waves. Individuals with major medical, neurodevelopmental or genetic disorders were excluded from the study sample.
Missing data analysis
Robustness checks indicated MNAR patterns in our data, as individuals with missing data in our longitudinal study had lower socioeconomic status (SES) than those with complete data, consistent with the literature on drop-out predictors49. Independent t-tests across all constructs confirmed these SES differences, as reported in Supplementary Table 12, alongside construct-level missingness percentages. Data contributors had significantly higher SES, with its absence accounting for 7–14% of the variance in anxiety, self-efficacy and interest measures and a substantial 43–44% in maths achievement and STEM career choice. Missing data led to the exclusion of affected individuals from the analysis, with the impact of SES on missingness duly considered in our interpretations.
Key variables
Mathematics self-efficacy and interest
Two scales adapted from the OECD Programme for International Student Assessment measured mathematics interest and self-efficacy. The maths interest scale consisted of three items: ‘I look forward to my mathematics lessons’, ‘I do mathematics because I enjoy it’ and ‘I am interested in the things I learn in mathematics’. The items had to be rated on a four-point scale from 1 (strongly disagree) to 4 (strongly agree). The scale showed good internal validity (α = 0.93). The mathematics self-efficacy scale asked, ‘How confident do you feel having to do the following tasks’ and consisted of eight items rated on a four-point scale from 0 (not at all confident) to 3 (very confident). Items included ‘solving an equation like 2(x + 3) = (x + 3)(x − 3)’ and ‘calculating how much cheaper a TV would be after a 30% discount’. This scale also showed good internal validity (α = 0.90).
Mathematics anxiety
Mathematics anxiety was assessed using an adapted version of the Abbreviated Math Anxiety Scale (AMAS)50. The scale asks participants to rate their anxiety level when facing a range of maths-related situations. Participants rated nine items on a five-point scale, from ‘not nervous at all’ to ‘very nervous’. Items included ‘listening to someone explain a maths formula’, ‘reading a maths book’ and ‘taking an examination in a maths course’. The AMAS has good internal validity (α = 0.94) and test-retest reliability (r = 0.85)50.
General anxiety
The Generalised Anxiety Disorder Scale (GAD-7) was used to collect data on general anxiety51. The scale asks participants to rate seven items on a four-point scale from 1 (not at all) to 4 (nearly every day), ‘how often in the past month have you been bothered by the following problems?’. Items include ‘feeling nervous, anxious or on edge’, ‘being so restless that it is hard to sit still’ and ‘feeling afraid as if something awful might happen’. The GAD-7 has good internal validity (α = 0.89) and moderate test-retest reliability (r = 0.64).
Maths anxiety (independent of general anxiety)
The maths anxiety measure (described above) was residualised for general anxiety (also described above) in the study sample to create a maths anxiety variable that was independent of the effects of general anxiety.
STEM career choice
At age 21, participants self-reported the degree subject they were either studying or had completed or the apprenticeship that they were either actively working in or had completed. 46.5% (2254) of participants at age 21 reported making a degree and/or apprenticeship choice. Degree and apprenticeship career routes were combined in this study, as only 117 (5.2%) of the cohort had chosen an apprenticeship route. Participants’ fields of work and study were classified as STEM or non-STEM using the UK’s Higher Education Statistics Agency (HESA) subject coding system—the Higher Education Classification of Subjects (HECoS)52. STEM career classifications included (but were not limited to) Mathematical Sciences, Computer Sciences and Technology, Medicine and allied professions, and Veterinary Sciences and Agriculture. Non-STEM career classification categories included (but were not limited to) Law, Creative Arts and Design, and Business Administration. STEM career choice was coded as a dichotomous variable where 1 = STEM career choice and 0 = Non-STEM career choice. The 0.5% of participants who had reported both a degree and apprenticeship choice were considered to have made a STEM career choice if one or both choices met the STEM career classification criteria.
Covariates
Maths achievement
Students take General Certificate of Secondary Education (GCSE) exams across all subjects at the end of compulsory education in the United Kingdom, at age 15–16. The mathematics GCSE exam is a mandatory part of the National Curriculum. Grades were obtained from questionnaires completed by parents or the twins themselves or via phone or email contact. Grades were recorded on a scale from 11 (corresponding to A*, the highest possible score) to 4 (corresponding to G, the lowest possible score). For a subsample of twins, self- and parent-reported GCSE grades were verified using data from the National Pupil Database (NPD)53, resulting in a correlation of 0.99 for mathematics54.
Socioeconomic status (SES)
Family SES was assessed at first contact. Mothers and fathers reported their educational level, occupational status and mothers’ age at the birth of their first child. The composite of SES was computed as the mean of the standardised scores of these measures. Socioeconomic status (SES) was taken at first contact as this wave offers minimal missing data, and SES has been shown to be developmentally stable, with age-to-age correlations of ~0.7048.
Statistical analyses
All analyses were conducted in R version 4.2.3. First, descriptive analyses were calculated for all variables. Univariate ANOVAs were used to investigate the mean-level sex differences for all variables. Then, zero-order correlations were carried out in the whole study sample to determine the relationships between all variables. Next, three binary logistic regression models were performed to determine the joint influence of mathematics anxiety (both corrected and uncorrected for general anxiety), maths self-efficacy, maths interest on STEM career choice and to investigate the effects of adjusting for previous maths achievement and socioeconomic status, as illustrated by Fig. 5. All continuous variables were residualised for age and sex, and these were used for all downstream analyses.

The three binary logistic regression models performed in the current study.
The zero-order correlations between all variables and binary logistic regressions outlined in Fig. 5 were then repeated in sex-stratified samples to assess whether the joint influence of maths anxiety, maths self-efficacy, and maths interest on STEM career choice differs between males and females. In the sex-stratified analyses, all continuous variables were residualised for age only before testing.
Responses