$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
BFEN algorithm performance analysis
The indicators evaluated in the experiments included the following:
Evaluation correction accuracy (ECA) denotes the proportion of deviation items correctly identified and rectified by the model within the teaching assessment results, thereby reflecting the overall effectiveness of the correction mechanism. Higher values indicate superior correction accuracy.
Fairness deviation rate (FDR) measures the magnitude of inter-group score disparities that the model fails to eliminate during the calibration process. It is computed as the ratio of the standard deviation of each sensitive attribute (e.g., gender, region, ethnicity) within the score distribution to the overall standard deviation; lower values imply minimized inter-group variances and a more robust fairness guarantee.
Scenario adaptation rate (SAR) indicates the percentage of correction tasks successfully executed by the model across heterogeneous instructional contexts (e.g., sciences vs. humanities, online vs. offline environments).
Fairness maintenance degree (FMD) is defined as one minus the ratio of the variance of fairness indicators among sensitive attribute groups after correction to that observed before correction, reflecting the capability of the system to preserve its structural fairness during the calibration process.
Disciplinary feature recognition rate (DFRR) calculates the proportion of samples in which core features are correctly identified within the assessment data of each discipline relative to the total assessment sample size, demonstrating the model's capacity to discern and capture domain-specific variations.
Multi-source data fusion efficiency (MDFE) refers to the throughput efficiency of the model during feature alignment, weight allocation, and deviation correction when fusing multi-source heterogeneous evaluation data. This comprehensive metric is defined as the product of the number of assessment samples processed per second (samples/sec) and the compliance rate of correction consistency (at a >95% confidence level), where higher values signify a more efficient and stable fusion pipeline.
Correction result stability (CRS) indicates the consistency of calibration outputs across multiple independent runs, quantified by the reciprocal of the standard deviation of the Euclidean distance between the correction outputs of each run under identical inputs. Higher values signify enhanced system reliability.
Single data correction time (SDCT) represents the average computational latency consumed by the model to complete the calibration of a single assessment sample, measured in milliseconds (ms); lower values indicate a stronger real-time response capability.
Corrected absolute error (CAE) measures the mean absolute error between the calibrated predictions and the ground-truth values, representing the residual deviation of the model relative to the original uncalibrated data. Lower values signify that the calibrated outputs align more closely with true performance levels.
Corrected relative error (CRE) quantifies the average absolute error between the calibrated predictions and ground-truth values expressed as a percentage of the ground truth, with lower values indicating higher relative calibration precision.
Correction accuracy rate (CAR) represents the proportion of samples whose post-calibration absolute error is < 0.5 relative to the total sample size, demonstrating the reliability of the model in satisfying predefined precision constraints. High values substantiate the empirical utility and credibility of the outputs.
Total loss value (TL) represents the comprehensive optimization objective value derived from the weighted summation of individual sub-loss terms upon task completion. Specifically, seven metrics—DFRR, MDFE, CRS, SDCT, CAE, CRE, and CAR—are standardized via Z-score normalization and weighted according to the operational priority of the evaluation tasks. Given the weight vector [0.15, 0.12, 0.1, 0.08, 0.13, 0.12, 0.1], these seven normalized indicator values are aggregated to compute the final loss.
Evaluation feature clustering accuracy (EFCA) evaluates the degree of alignment between the unsupervised clustering configurations generated by the model and the ground-truth categories validated by domain experts.
High-dimensional data processing efficiency (HDPE) measures the number of samples subjected to feature dimensionality reduction and deviation correction per unit time when handling high-dimensional sparse vectors containing ≥50 evaluation features. Measured in samples/second, higher values reflect a more robust capability to process complex, large-scale assessment inputs in real time.
Fairness correction precision (FCP) evaluates the uniformity of the calibration effects across divergent student cohorts. This metric is quantified by calculating the mean of the distribution of the Kolmogorov-Smirnov distance for the corrected absolute errors across the sensitive attribute subgroups; lower values imply a more balanced cross-group performance and stronger fairness assurance.
Cross-scene correction consistency (CSCC) quantifies the distributional offset of calibration results across three typical scenarios—namely, classroom assessment, practical evaluation, and online testing—modeled as the Wasserstein distance ratio.
To verify the performance of the proposed BFEN teaching assessment calibration algorithm, the study compared it with the PRF algorithm, which combined principal component analysis (PCA) and random forest (RF); the EAB algorithm, which combined extreme learning machine (ELM) and adaptive boosting (AdaBoost); and the DBLR algorithm, which combined deep belief network (DBN) and logistic regression (LR). The experiments ran on a system equipped with an Intel Core i9-13900HX processor and NVIDIA RTX 4080 GPU, providing high parallel computing capability and large video memory to efficiently complete feature extraction and calibration computations for large-scale teaching assessment data. The operating system was Windows 11, and Python 3.9 was used for programming. The algorithms were implemented using the Scikit-learn library, deep learning modules were developed via the TensorFlow framework, and Pandas and NumPy libraries were employed for data preprocessing. The development environment used PyCharm Professional 2023.1, and Seaborn was applied for the visualization of calibration results to facilitate intuitive comparison of algorithm performance. To ensure the reliability of the data, the study utilized the Global Education Monitoring Report dataset1 and the Spanish middle school student academic performance dataset2. The dataset contains 12,486 records, covering multimodal teaching evaluation data such as university faculty teaching assessments, student academic performance, student teaching evaluations, and course feedback, encompassing 137 teaching dimension features, including classroom interaction frequency, assignment feedback timeliness, scoring consistency, language inclusivity, and cross-cultural sensitivity. The target variable is the teaching evaluation score, which in this study is mapped to a multidimensional weighted composite value calibrated by experts. The research integrates four types of information sources: student teaching evaluations, peer reviews, supervisory classroom observations, and teaching dossier reviews, with weights of 0.35, 0.25, 0.25, and 0.15, respectively. The training and validation sets are divided by chronological order from the two datasets. The study uses data from September 2024 as the training set and data from October 2024 to June 2025 as the validation set to align with the temporal progression of educational assessments. In the preprocessing stage, a modality-specific strategy is employed, with structured features standardized via Z-score normalization and outliers Winsorized, while textual features are encoded using the BERT-base-chinese model and reduced to 768 dimensions. Adversarial training is applied to enhance robustness against implicit bias expressions and mitigate semantic drift caused by cultural background differences.
To achieve cross-domain training and validation research, the model will be transferred to four heterogeneous teaching scenarios: K-12 basic education, vocational education, continuing education, and international Chinese education for zero or few sample validation. In these four scenarios, the study introduces domain adaptive regularization terms during algorithm training to constrain the consistency of feature distribution of the model in different teaching scenarios. Embedding a lightweight syntax-aware mask mechanism for short sentence noise in K-12 scenarios; To address the ambiguity of professional terminology in vocational education, a dynamic domain dictionary is constructed and integrated with the curriculum outline knowledge graph. Design an age group comparison learning module to address the intergenerational semantic gap in continuing education, aligning semantic anchors for different age groups' evaluation expressions in the latent space. To address the cultural load expression in international Chinese education, cross-language comparison prompts are used for fine-tuning to enhance the robustness of understanding non-literal educational concepts. Research on expert calibration and reliability testing of structured fields in original evaluation records, removing low-reliability entries with Krippendorff's alpha coefficient below 0.75. Double blind manual annotation was implemented on student evaluation texts, with 3 educational measurement experts independently completing deviation-type labeling, achieving Cohen's k = 0.86 consistency in annotation. Finally, the annotated results were cross-validated with the supervision attendance records and teaching file review conclusions to generate the final calibration labels.
The selected datasets originated from distinct populations, measurement systems, and educational backgrounds. This cross-domain validation paradigm rigorously tested the robustness and generalization boundaries of the model within an authentic educational ecosystem, preventing evaluation illusions caused by data homogenization. Both repositories contained substantial volumes of teaching assessment records, offering direct reference value for the deployment of fairness-aware intelligent education algorithms. Both datasets contained substantial teaching assessment data, providing direct reference values for the design of fairness-aware intelligent education algorithms and teaching assessment calibration. The study calibrated 1000 teaching assessment records with the four algorithms and calculated their ECA and Fairness FDR. The results are shown in Figure 8.

Figure 8: Comparison of test results between ECA and FDR. (A) BFEN. (B) PRF. (C) EAB. (D) DBLR. Please click here to view a larger version of this figure.
As shown in Figure 8A, BFEN initially achieved an ECA of 82.47% and an FDR of 5.71% in the assessment calibration task. As the number of calibrated records increased, the ECA rose to 98.75% and stabilized after calibrating 421 records, while the FDR decreased to 2.87% and stabilized after calibrating 514 records. As shown in Figure 8B, the ECA curve of the PRF algorithm gradually increased, while the FDR line slowly decreased. At the end of the calibration task, its ECA value was 92.30%, and PDR value was 6.72%. Figure 8C demonstrates that the ECA curve of the EAB algorithm rose rapidly before flattening, whereas its FDR trajectory displayed severe early-stage fluctuations prior to convergence, terminating with an ECA of 84.64% and an FDR of 8.93%. As delineated in Figure 8D, the DBLR algorithm displayed a sluggish upward ECA trajectory accompanied by marginal fluctuations and limited compression in its FDR line, concluding the calibration task with an ECA of 83.51% and an FDR of 8.42%. These experimental findings confirm that the BFEN algorithm significantly outperforms the baseline approaches in both evaluation correction accuracy and fairness bias control. This superior generalization capability is attributed to the integration of BERT within BFEN, which extracts deep semantic features from unstructured evaluation text to capture hidden bias expressions, while the RL module dynamically optimizes fairness thresholds and correction parameters via a multi-objective reward function to decouple sensitive attributes from the final outputs. The study then tested SAR and FMD for classroom assessment, final exams, practical evaluation, and online assessment, labeling the four scenarios as A, B, C, and D. The results are shown in Figure 9.

Figure 9: Comparison of SAR and FMD results in different scenarios. (A) SAR test results. (B) FMD test results. Please click here to view a larger version of this figure.
As shown in Figure 9A, BFEN achieved SAR above 98% across all teaching scenarios, with the highest SAR of 99.01% for classroom tests. The SAR of the comparison algorithm in different scenarios is lower than that of the BFEN algorithm, among which the DBLR algorithm has the lowest SAR for the final assessment scenario among all algorithms, with a SAR of 70.23%. As shown in Figure 9B, the FMD of the BFEN algorithm is still significantly higher than that of the comparison algorithm in different teaching scenarios, with FMDs of 98.47%, 98.05%, 97.81%, and 97.94% in different scenarios, respectively. The FMDs of the benchmark algorithm DBLR are 82.36%, 71.74%, 70.59%, and 69.12%, respectively. The FMDs of the EAB algorithm are 80.79%, 68.21%, 67.65%, and 66.03%, respectively, while the FMDs of the PRF algorithm are 86.42%, 82.17%, 80.98%, and 78.33%, respectively. These results demonstrated that BFEN outperformed comparison algorithms due to GBDT's iterative learning with multiple decision trees, which accurately captured error patterns across scenarios and ensured stable and fair calibration results. The study further evaluated disciplinary feature recognition rate (DFRR) for Chinese, Mathematics, and English using the four algorithms. The results are shown in Figure 10.

Figure 10: Comparison of DFRR test results in different disciplines. (A) Training set. (B) Validation set. Please click here to view a larger version of this figure.
As shown in Figure 10A, BFEN achieved DFRR of 99.23%, 98.94%, and 99.13% for Chinese, Mathematics, and English in the training set. Figure 10B indicated that BFEN also achieved higher DFRR in the validation set compared to other algorithms. Specifically, BFEN's DFRR for Chinese reached 98.41%, 10.8% higher than DBLR's 87.61%; for Mathematics, 97.54%, 9.07% higher than PRF's 88.47%; and for English, 98.13%, 13.34% higher than EAB's 84.79%. These results confirmed that BFEN efficiently adapted to disciplinary assessment calibration by combining BERT's deep capture of semantic differences with GBDT's personalized learning of error patterns. To comprehensively assess BFEN performance, the study evaluated MDFE, CRS, and SDCT for the four algorithms. The results are presented in Table 1.
| Dataset | Algorithm | MDFE (%) | CRS | SDCT (s) |
| Training | BFEN | 98.42 ± 0.28*&@ | 0.975 ± 0.28*&@ | 0.78 ± 0.04*&@ |
| PRF | 83.85 ± 0.41 | 0.779 ± 0.36 | 1.32 ± 0.08 |
| EAB | 85.91 ± 0.50 | 0.806 ± 0.40 | 1.15 ± 0.07 |
| DBLR | 88.76 ± 0.47 | 0.753 ± 0.39 | 1.45 ± 0.08 |
| Validation | BFEN | 97.58 ± 0.25*&@ | 0.968 ± 0.27*&@ | 0.86 ± 0.03*&@ |
| PRF | 81.62 ± 0.32 | 0.687 ± 0.50 | 1.51 ± 0.06 |
| EAB | 84.37 ± 0.40 | 0.749 ± 0.42 | 1.28 ± 0.05 |
| DBLR | 87.53 ± 0.30 | 0.728 ± 0.47 | 1.63 ± 0.07 |
Table 1: Comparison of MDFE, CRS, and SDCT test results. “*” indicates a significant difference (p < 0.05) between BFEN and PRF results. “&” indicates the difference between BFEN and EAB results was significant (p < 0.05). “@” indicates the difference between BFEN and DBLR results is significant (p < 0.05)
Table 1 indicates that BFEN achieved an MDFE of 98.42% in the training set, 14.57% higher than PRF's 83.85%, a CRS of 0.975, 0.222 higher than DBLR's 0.753, and an SDCT of 0.78 s, 0.67 s less than DBLR's 1.45 s. In the validation set, BFEN maintained superior performance, with MDFE, CRS, and SDCT of 97.58%, 0.968, and 0.86 s, respectively, outperforming comparison algorithms. Overall, the results validated the superior comprehensive performance of BFEN in teaching assessment calibration.
GFBFEN teaching assessment calibration model performance validation
Based on the BFEN performance validation, the study further evaluated the performance of the GFBFEN teaching assessment calibration model constructed on BFEN and compared it with calibration models built using PRF, EAB, and DBLR algorithms. To ensure consistent testing conditions and comparability, the global education monitoring report dataset served as the training set, while the university student academic performance dataset served as the validation set. The four models calibrated 500 teaching assessment records, and their CAE and CRE were calculated. The results are shown in Figure 11.

Figure 11: Comparison of CAE and CRE test results. (A) CAE test results. (B) CRE test results. Please click here to view a larger version of this figure.
As shown in Figure 11A, GFBFEN initially achieved a CAE of 0.1279. As calibration proceeded, the CAE decreased to 0.0641 and stabilized after 104 records. The comparison models had higher initial and final CAE than GFBFEN, with EAB having the highest initial and final CAE of 0.1858 and 0.1217, respectively. Figure 11B indicated that GFBFEN's initial and final CRE during the calibration task remained lower than the comparison models, with values of 0.1137 and 0.0912, respectively. These results demonstrated that GFBFEN exhibited strong fitting capability, benefiting from GAT's attention mechanism, which constructed semantic correlation graphs among features, enhanced semantic alignment of multi-source assessment data, and reduced ineffective calibration caused by feature bias. The study then evaluated CAR for student assessment data, teacher assessment data, school assessment data, and online assessment data, labeled as I, II, III, and IV. The results are shown in Figure 12.

Figure 12: Comparison of CAR results of different evaluation data. (A) GFBFEN. (B) PRF. (C) EAB. (D) DBLR. Please click here to view a larger version of this figure.
As shown in Figure 12A, GFBFEN achieved CAR of 99.34%, 98.93%, 99.01%, and 99.12% for student, teacher, school, and online assessment data in the training set. In the validation set, CAR values were 97.89%, 98.56%, 97.57%, and 98.46%, respectively. Figure 12B–D showed that comparison models had lower CAR in both training and validation sets. Among them, EAB performed the worst in the validation set, with CAR of 76.31% for teacher data and 78.29% for online data. These results confirmed that GFBFEN significantly outperformed the comparison models. Next, the study tested the TL in the teaching assessment calibration task to evaluate fitting ability and training stability. The results are shown in Figure 13.

Figure 13: Model fitting ability and training stability test results. (A) GFBFEN. (B) PRF. (C) EAB. (D) DBLR. Please click here to view a larger version of this figure.
As shown in Figure 13A, GFBFEN initially had a TL of 0.1021 in the training set. After 67 iterations, TL decreased to 0.0725 and stabilized. In the validation set, TL decreased from 0.1237 to 0.1018. Figure 13B–D indicated that PRF's final TL in the training set was 0.0824, EAB's final TL in the validation set was 0.1178, and DBLR's TL in both sets was unstable and significantly higher than other models. These results demonstrated that GFBFEN exhibited strong fitting ability and training stability, benefiting from KD-based knowledge transfer, which allowed the model to quickly learn effective features, accelerate loss convergence, and ensure generalization across datasets. To comprehensively validate GFBFEN's core performance, the study tested EFCA, HDPE, FCP, and CSCC for the four models. The results are shown in Figure 14.

Figure 14: Comprehensive indicator test results of teaching assessment and correction tasks. (A) GFBFEN test results. (B) PRF test results. (C) EAB test results. (D) DBLR test results. Please click here to view a larger version of this figure.
As shown in Figure 14A–D, the GFBFEN model has EFCA values of 99.35% and 98.76% in the training and validation sets, respectively. The EFCA of the PRF model are 88.53% and 87.04%, respectively. In the validation set, the EFCA of the GFBFEN model was 6.59% higher than that of EAB model at 92.17%, and 11.72% higher than that of the DBLR model at 87.04%. In addition, its HDPE, FCP, and CSCC indicators are also comprehensively leading: in the validation set, the HDPE of the GFBFEN model reached 98.05%, the FCP reached 97.93%, and the CSCC reached 0.968, all of which outperformed the PRF, EAB, and DBLR models. Overall, these results validated the superior comprehensive performance of GFBFEN, demonstrating that the model's coordinated optimization of GAT-KD, AM-LSTM, and BFEN effectively improved calibration accuracy, efficiency, and fairness in teaching assessment.
The research gradually constructed FairEduNet, BFEN, FBFEN, and GFBFEN, and the final GFBFEN includes components such as FairEduNet, BERT-RL, AM-LSTM, and GAT-KD. To verify the independent contribution of individual components, research and design ablation experiments are conducted, gradually removing each component and evaluating its impact on overall performance. Comparative indicators include ECA, FDR, SAR, and FMD. The results showed that the ECA value of the FairEduNet component alone was 87.32%, FDR was 12.45%, SAR was 89.67%, and FMD was 11.89%. The ECA of the complete GFBFEN model reached 99.21%, the FDR further dropped to 1.93%, the SAR remained stable at 99.45%, and the FMD was optimized to 7.28%. After removing BERT-RL, ECA decreased to 91.04%, FDR value was 6.23%, SAR value was 92.33%, and FMD increased to 7.85%. The ablation of AM-LSTM increased SAR fluctuations by 14.2%. The removal of GAT-KD led to an increase in FMD to 8.36%, and its suppression mechanism of graph structure bias propagation guided by knowledge distillation was significantly effective in alleviating implicit association bias between populations.
To further test the research methods, the study introduced established fairness standard indicators, namely the Fairness Correction Benchmark (FCB), which covers three types of rigid constraints widely recognized in educational scenarios: (1) The absolute value of the mean difference after correction between groups is ≤ 1.2 points (based on the percentage system); (2) The overlap rate of the corrected confidence intervals for each sensitive attribute subgroup is ≥ 95%; (3) In any single scenario, the joint feasible region of FDR and SAR must satisfy the Pareto frontier constraint (that is, it is impossible to improve SAR without deteriorating FDR, and vice versa). The GFBFEN model was compared with mainstream baseline models, such as EAB, PRF, DBLR, and GBDT, in the experimental evaluation. The experiment was conducted on 421 real assessment data collected from actual teaching scenarios, covering three typical educational scenarios: classroom assessment, practical evaluation, and online testing. The results are shown in Table 2.
| Algorithm | Absolute value of mean score difference | Confidence interval overlap rate(%) | Joint feasible region satisfaction rate(%) | Comprehensive score |
| GFBFEN | 0.87 | 98.3 | 100 | 97.2 |
| EBA | 1.52 | 98.4 | 82.6 | 78.5 |
| PRF | 1.68 | 85.1 | 74.3 | 71.2 |
| DBLR | 1.45 | 87.9 | 79.8 | 75.6 |
| GBDT | 1.33 | 91.2 | 86.4 | 79.9 |
Table 2: Evaluation results of fairness criteria indicators and practical application verification.
As shown in Table 2, GFBFEN achieved full compliance with all four FCB rigid constraints independently, and its comprehensive score significantly outperformed the other baseline models by a margin of +18.7 points, validating the robustness of the proposed multi-stage collaborative correction paradigm. Particularly, in the core indicator reflecting the fairness-efficiency trade-off capability, the joint feasible region coverage rate reached 100%, indicating that the model possesses deployable Pareto-optimal decision-making capabilities in real educational scenarios.
Fairness in teaching evaluation refers to the use of verifiable quantitative constraints as a benchmark while respecting individual differences and educational laws. The calculation logic and threshold setting of equity indicators are based on the theoretical framework of educational equity and empirical data verification standards. They are jointly reviewed by an expert committee composed of five scholars to ensure a high degree of unity between theoretical rigor and adaptability to educational practice. To further validate the adaptability of the model under different fairness standards, the study conducted further validation under multiple standards. Specifically, three fairness standards were selected for verification. Standard 1 is the acceptance rate balance between groups based on Demographic Parity. Standard 2 is the cross-group consistency of true positive rate and false positive rate based on Equalized Odds. Standard 3 is based on Predictive Parity for controlling inter-group bias in prediction accuracy. It was found that GFBFEN consistently met the rigid constraint requirements under all three standards. Standard 1 group acceptance rate deviation ≤1.2%, standard 2 true/false positive rate intergroup difference ≤0.85%, standard 3 prediction accuracy intergroup deviation controlled within ±±0.6%. In contrast, the other models showed a constraint breakthrough of ≥3.7% under at least one criterion, verifying the generalization compatibility of GFBFEN for multivariate fairness paradigms.
DATA AVAILABILITY:
To ensure data reliability, the study used the Global Education Monitoring Report dataset (https://www.education-progress.org/) and the Spanish middle school student academic performance dataset (https://archive.ics.uci.edu/dataset/320/student+performance). The training and validation sets are divided by chronological order from the two datasets. The study uses data from September 2024 as the training set and from October 2024 to June 2025 as the validation set, aligning with the temporal progression of educational assessments.