To detect academic misconduct in students' assessment score trajectories under severe class imbalance. The paper compares tabular learners, gradient-boosting models, and sequence-aware deep networks, and proposes a precision–recall–centric evaluation and deployment protocol (calibration, threshold selection and Recall@Top-k%) tailored to rare-event screening in educational settings.
A cohort of 1,527 students (2021–2024) is modeled using ten algorithms: LR, DT, RF, SVM, MLP, XGBoost, CatBoost, LightGBM, GRU-RNN and 1D-CNN. Features encode sequential score dynamics and metadata. Models are tuned via cross-validation; probabilities are calibrated (Platt/Isotonic); operating thresholds are chosen on validation to maximize minority-class F1 or a cost-sensitive utility. Performance is assessed on a hold-out test set with PR-AUC (headline), F1(+), Recall@Top-k%, ROC-AUC, calibration curves, Brier, and bootstrap CIs.
Sequence-aware models dominate: GRU-RNN and 1D-CNN achieve ROC-AUC ˜0.97–0.98 and the highest F1(+) and Recall@Top-5%. Tabular/boosting baselines show ˜0.90 accuracy yet miss most positives at the default 0.5 threshold, highlighting the necessity of calibration and threshold optimization. With PR-centric selection and tuned operating points, deep temporal models yield strong screening utility for limited human review budgets.
Labels reflect suspected–not adjudicated–cheating, introducing noise. The single-institution cohort may limit external validity; temporal shift across semesters can degrade performance. Future work should include multi-site evaluation, collusion/graph modeling, semi-/weak-supervision for noisy labels and governance topics (fairness audits, drift monitoring and uncertainty reporting).
This study highlights how educational institutions can leverage machine learning for early detection of academic dishonesty based on historical performance data. CNN and RNN models are promising tools for identifying anomalous learning patterns. However, practical deployment requires preprocessing techniques to manage class imbalance and threshold optimization to reduce false negatives. The findings provide a roadmap for building automated cheating detection systems in both online and traditional assessment environments.
By improving the ability to detect cheating, this research contributes to fairer academic environments, upholding educational integrity and credibility of credentials. However, ethical considerations must be taken into account to avoid false accusations and ensure student rights. Human-in-the-loop systems are crucial for verifying algorithmic predictions before disciplinary action, thereby fostering transparency and accountability in automated decision-making processes.
The study unifies a minority-focused evaluation protocol with a comprehensive comparison of tabular, boosting, and sequence-aware models for cheating detection from score trajectories. It demonstrates the decisive value of temporal representation learning and provides a reproducible pipeline and operational metrics that align model performance with real investigative workflows.
1. Introduction
Academic grades constitute primary proxies for student learning, encapsulating the mastery of disciplinary knowledge and skills essential for both academic progression and professional success. Consequently, accurate forecasting of final grade outcomes is pivotal for educational stakeholders, guiding timely interventions and policy formulation. However, the heterogeneity and high dimensionality of longitudinal academic records pose significant challenges to the construction of robust predictive models. In response, machine learning techniques—integrating advanced algorithms with rigorous statistical methodologies—have proven effective at distilling key determinants of scholastic performance, mapping learning behaviors and trajectories, and identifying anomalous patterns associated with academic dishonesty. The fusion of these computational strategies underpins the field of Educational Data Mining (EDM), which systematically exploits data generated within educational systems to derive empirically grounded insights and actionable knowledge for both performance prediction and cheating prevention. In this study, we cast the prediction of student cheating on final assessments as a dual task of classification and anomaly detection using score‐based features. Drawing on EDM and machine learning, we propose an end-to-end framework that transforms historical performance data into predictive models capable of stratifying learners by risk: those likely to excel, those requiring support, and those at elevated risk of failure or dishonest behavior. Continuous outcomes are modeled via regression analysis—most commonly linear regression—where dependent variables (final scores) are estimated as functions of continuous or categorical predictors. Discrete outcomes employ classification algorithms, including support vector machines, backpropagation-trained neural networks, and k-nearest neighbors, to allocate students into risk categories. Collectively, these methods facilitate early detection of potential cheating and inform tailored prevention strategies.
In this work, we frame cheating prediction as (1) score-trajectory classification and (2) rank-based triage. Given pronounced class imbalance, we evaluate models using PR-centric metrics—minority-class F1, PR-AUC, and Recall@Top-k%—and treat ROC-AUC/accuracy as secondary summaries. We calibrate probabilities (Platt/Isotonic) and select thresholds on the validation fold to maximize F1(+) or a cost-sensitive utility. We then confirm generalization on a hold-out test set and report bootstrap 95% CIs. The rest of this paper is organized as follows. Section 2 surveys relevant literature. Section 3 details the proposed methodology. Section 4 outlines the experimental setup and presents a comparative analysis of multiple algorithms. Finally, Section 5 concludes with key findings and future research directions.
2. Related works
The deployment of machine learning for predicting and detecting academic dishonesty has seen rapid growth, with studies spanning anomaly detection and behavioral analytics Elrahman et al. [1] demonstrated that interactive eTextbook engagement data can predict student performance, suggesting that fine-grained interaction logs may also reveal irregular behaviors associated with cheating. Subsequent work [2, 3] applies anomaly-detection techniques to continuous assessment results, showing that deviations from projected score trajectories—rather than simple outliers-can flag potential dishonesty in final examinations. Ensemble methods, particularly Random Forests, have proven effective in modeling cheating behavior: RF classifiers successfully identified reported instances of misconduct among elementary students [4], while XGBoost and related deep-learning models have been leveraged to analyze online exam interactions and detect collusive patterns [5, 6]. Sequential architectures, notably LSTM networks, extend this capability by capturing temporal dependencies in test-score sequences, thereby improving sensitivity to evolving anomalies over time [7]. Beyond post‐hoc detection, machine learning insights have informed proactive interventions, such as optimized seating algorithms that minimize collusion risk based on historical cheating patterns [8]. Collectively, these studies illustrate the versatility of both traditional classifiers and advanced neural models in preserving assessment integrity through the systematic analysis of performance and behavioral data [2, 9].
Machine learning has become a cornerstone for predicting and detecting student cheating from exam data, spanning traditional classifiers, sequential models, deep learning architectures, and even realtime video analysis. The principal methodologies are summarized as follows:
Traditional and Ensemble Classification: Binary classifiers—including generalized linear models (GLM) [7], logistic regression (LR) [10], decision trees (DT), and random forests (RF) [4]—leverage features drawn from student demographics and constructs such as the Fraud Triangle (rationalization, opportunity, pressure). These models typically achieve precision rates of 50–75% for the cheating class. RF frequently identify “opportunity” variables (e.g., unsupervised exam settings) as the strongest predictors, whereas GLM, LR, and DT emphasize“rationalization”factors (e.g., self‐reported justifications for cheating).
Sequential Models with Outlier Detection: Recurrent neural networks (RNNs) treat cheating detection as an anomaly‐detection task on time‐ordered assessment sequences. By training RNNs to forecast final‐exam scores from preceding quizzes and midterms, significant deviations between predicted and observed scores are flagged as potential dishonesty [11]. This temporal framing yields true‐positive rates near 95% and false‐positive rates around 5%, outperforming static outlier‐detection baselines by capturing order‐dependent irregularities in score trajectories.
Deep Learning Variants: Advanced neural architectures—dense feedforward networks (DNN), LSTM, and hybrid DenseLSTM models [12]—have been evaluated on academic datasets. These deep models exploit complex temporal and nonlinear interactions among assessments, with DenseLSTM variants reporting up to 95% accuracy in cheating detection, significantly surpassing simpler classifiers.
Video - based Behavioral Monitoring: For realtime invigilation, object-detection frameworks such as YOLOv8 [3], augmented with attention mechanisms, analyze proctoring video streams to identify cheating behaviors (e.g., whispering, note passing, unauthorized device use). Trained on a bespoke dataset of staged infractions, this approach achieves approximately 82.7% detection accuracy, reducing reliance on constant human oversight. Modern remote-proctoring stacks typically comprise: face/pose tracking; gaze/attention estimation; object detection for handheld/secondary screens; and event logic for rule violations. Many pipelines adopt YOLOv8-based detectors (for devices, notes, mouth/hand occlusions), sometimes complemented by temporal modules (e.g., TCN/GRU) or rule-based smoothing to suppress flicker. Reported benefits include real-time operation and modularity; known challenges include domain shift (camera/lights), privacy, and dataset bias. Our work complements these pipelines by focusing on score-trajectory signals, which are modality-agnostic and auditable for after-the-fact investigations [13–15].
These machine learning approaches provide effective tools for early detection and prevention of academic misconduct based on exam results and student behavior. The combination of temporal modeling of exam scores and behavioral video analysis covers a broad spectrum of cheating detection scenarios, enhancing academic integrity in both traditional and online examination settings (see Table 1).
The summary of machine learning algorithms for predicting cheating from exam results
| Approach | Algorithms/models | Data used | Key strengths | Performance highlights |
|---|---|---|---|---|
| Exam score anomaly detection | RNN + Outlier Detection | Sequential exam scores | Captures temporal patterns, unsupervised | TPR∼95%, FPR∼5% |
| Binary classification | GLM, LR, DT, RF | Demographics + psychological features | Interpretable, uses fraud theory factors | Precision 50–75% for cheating class |
| Deep learning temporal models | DNN, LSTM, DenseLSTM, RNN | Exam score sequences | High accuracy on complex patterns | Up to 95% accuracy |
| Exam score anomaly detection | RNN + Outlier Detection | Sequential exam scores | Captures temporal patterns, unsupervised | TPR∼95%, FPR∼5% |
| Video behavior detection | Improved YOLOv8 + Attention mechanism | Exam room video footage | Real-time detection of cheating actions | ∼82.7% accuracy |
| Approach | Algorithms/models | Data used | Key strengths | Performance highlights |
|---|---|---|---|---|
| Exam score anomaly detection | RNN + Outlier Detection | Sequential exam scores | Captures temporal patterns, unsupervised | TPR∼95%, FPR∼5% |
| Binary classification | GLM, LR, DT, RF | Demographics + psychological features | Interpretable, uses fraud theory factors | Precision 50–75% for cheating class |
| Deep learning temporal models | DNN, LSTM, DenseLSTM, RNN | Exam score sequences | High accuracy on complex patterns | Up to 95% accuracy |
| Exam score anomaly detection | RNN + Outlier Detection | Sequential exam scores | Captures temporal patterns, unsupervised | TPR∼95%, FPR∼5% |
| Video behavior detection | Improved YOLOv8 + Attention mechanism | Exam room video footage | Real-time detection of cheating actions | ∼82.7% accuracy |
Parallel to score-based modeling, several EDM 2025 accepted/demo systems treat the examinee's screen stream as a spatiotemporal signal. A common pattern is a hybrid CNN + RNN (e.g., 2D CNN on frame snippets for spatial cues like unauthorized UI widgets, overlaid with GRU/LSTM on frame indices for temporal dynamics). These systems report that short “burst” kernels detect micro-events (copy/paste dialogs, sudden context switches), while the recurrent stack stabilizes predictions across seconds-long sequences, improving precision at low false-alarm rates. We contextualize our CNN/RNN design choices with this finding, emphasizing kernel sizes for burst-like artifacts and sequence-level aggregation [16, 17]. The 2025 EDM surveys converge on three trends: (1) sequence-centric models for longitudinal student traces; (2) multimodal fusion (scores, clickstreams, proctoring media); and (3) imbalance-first evaluation using PR-AUC and cost-sensitive F-metrics rather than accuracy/ROC alone. The surveys also advocate for threshold selection as a model component, not a post-hoc tweak, aligning with our validation-time F1/utility maximization and probability calibration pipeline [18]. Cheating detection is inherently imbalanced (typically <10% positives). Recent guidance emphasizes: (1) PR-AUC as the headline metric (baseline ≈ prevalence), (2) class-weighted or focal losses for deep models, (3) calibration (Platt/Isotonic) before thresholding, (4) minority-aware minibatching and within-fold SMOTE (training only) to avoid leakage, and (5) reporting Recall@Top-k% risk for operational triage. We adopt this toolbox and treat threshold/operating-point choice as part of model selection, not a cosmetic step [19]. Recent ACM work broadens cheating from individual-level anomalies to group-level collusion. Techniques include graph-based similarity over response vectors and timing, pairwise sequence alignment for suspicious synchrony, and clustering on keystroke/interaction embeddings to identify coordinated behavior. Such approaches caution that “normal-looking” scores may still mask collusion detectable only via relational signals [20–25]. While our current study is individual-centric, we outline collusion modeling as a near-term extension. We (1) compare deep temporal models (CNN/RNN) with strong tabular baselines under strict imbalance protocols; (2) report PR-centric metrics, calibration, and threshold sweeps; and (3) propose pathways to multimodal fusion (with screen-sequence/proctoring cues) and collusion graphs, aligning the score-based approach with 2025 EDM/proctoring practice.
3. Methodology
In this section, we propose a predicted student cheating model was analyzed based on the scores in the learning process of the subjects in the previous semester flows shown on Figure 1.
The diagram starts at the top left, with a database icon labeled “Data Students”, which feeds downward into a rectangular block labeled “Data Preprocessing”. Inside this block, two components are labeled “Feature Encoding” and “Feature Scaling”. A small icon is shown at the bottom left in the Data Preprocessing” box. “Data Preprocessing” branched into two: an arrow leads upward to a database icon labeled “Data Train”, and another arrow leads downward to a database icon labeled “Data Test”. From “Data Train”, an arrow points to a large, rounded rectangle labeled “Modelings”. Inside this modeling block, ten labeled algorithm boxes are arranged in two columns, and five rows, labeled row-wise as “Logistic Regression”, “Support Vector Machines”, “Decision Trees”, “X G Boost Regressor”, “Random Forests”, “Artificial Neural Networks”, “CatBoost”, “Recurrent Neural Network”, “Light G B M”, and “Convolutional Neural Network”. An arrow from the “Modelings” block points downward to a rectangular block labeled “Evaluate, Prediction”. Inside this block, evaluation metrics are labeled on the right as “M S E (Mean Squared Error)”, “R 2 Score”, “R M S E (Root M S E)”, and “M A E (Mean Absolute Error)”. On the left of this block, the word “Cheating?” appears at the bottom alongside a small cartoon-style icon of a person standing next to a checklist and symbols. The “Data Test” database icon leads rightward to the “Evaluate, Prediction” box.The predicted student cheating model was analyzed based on the scores in the learning process
The diagram starts at the top left, with a database icon labeled “Data Students”, which feeds downward into a rectangular block labeled “Data Preprocessing”. Inside this block, two components are labeled “Feature Encoding” and “Feature Scaling”. A small icon is shown at the bottom left in the Data Preprocessing” box. “Data Preprocessing” branched into two: an arrow leads upward to a database icon labeled “Data Train”, and another arrow leads downward to a database icon labeled “Data Test”. From “Data Train”, an arrow points to a large, rounded rectangle labeled “Modelings”. Inside this modeling block, ten labeled algorithm boxes are arranged in two columns, and five rows, labeled row-wise as “Logistic Regression”, “Support Vector Machines”, “Decision Trees”, “X G Boost Regressor”, “Random Forests”, “Artificial Neural Networks”, “CatBoost”, “Recurrent Neural Network”, “Light G B M”, and “Convolutional Neural Network”. An arrow from the “Modelings” block points downward to a rectangular block labeled “Evaluate, Prediction”. Inside this block, evaluation metrics are labeled on the right as “M S E (Mean Squared Error)”, “R 2 Score”, “R M S E (Root M S E)”, and “M A E (Mean Absolute Error)”. On the left of this block, the word “Cheating?” appears at the bottom alongside a small cartoon-style icon of a person standing next to a checklist and symbols. The “Data Test” database icon leads rightward to the “Evaluate, Prediction” box.The predicted student cheating model was analyzed based on the scores in the learning process
3.1 The analytical workflow
Figure 2 summarizes the principal stages of our analytical workflow, encompassing data preprocessing (feature encoding and scaling), data partitioning (training versus testing), model development, and performance evaluation. Each stage is detailed below to elucidate the research methodology.
The five-step diagram is arranged from left to right in overlapping rectangular boxes, connected by right arrows. The box on the left is labeled “Step 1 Data Collection” and lists “Raw exam scores” and “Course metadata”, illustrated with a funnel icon with a pie chart shown on a paper. The second box is labeled “Step 2 Data Preprocessing” and includes “Feature Encoding” and “Feature Scaling”, shown with database and connected gear icons. The third box is labeled “Step 3 Train or Test Split” and lists “Data Train” and “Data Test”, illustrated with folder and database icons. The fourth box is labeled “Step 4 Modeling” and lists algorithms including “Logistic Regression”, “Decision Trees”, “Random Forests”, “Support Vector Machines”, “Neural Networks”, “X G Boost Regressor”, “Cat Boost”, and “Light G B M”, accompanied by a brain and circuit-style icon. The fifth box is labeled “Step 5 Evaluation” and lists “R-squared Score”, “M S E (Mean Squared Error)”, “R M S E (Root M S E)”, and “M A E (Mean Absolute Error)”, shown with a clipboard and checklist icon.Our proposal steps for predicting student cheating model
The five-step diagram is arranged from left to right in overlapping rectangular boxes, connected by right arrows. The box on the left is labeled “Step 1 Data Collection” and lists “Raw exam scores” and “Course metadata”, illustrated with a funnel icon with a pie chart shown on a paper. The second box is labeled “Step 2 Data Preprocessing” and includes “Feature Encoding” and “Feature Scaling”, shown with database and connected gear icons. The third box is labeled “Step 3 Train or Test Split” and lists “Data Train” and “Data Test”, illustrated with folder and database icons. The fourth box is labeled “Step 4 Modeling” and lists algorithms including “Logistic Regression”, “Decision Trees”, “Random Forests”, “Support Vector Machines”, “Neural Networks”, “X G Boost Regressor”, “Cat Boost”, and “Light G B M”, accompanied by a brain and circuit-style icon. The fifth box is labeled “Step 5 Evaluation” and lists “R-squared Score”, “M S E (Mean Squared Error)”, “R M S E (Root M S E)”, and “M A E (Mean Absolute Error)”, shown with a clipboard and checklist icon.Our proposal steps for predicting student cheating model
Step 1. Data Collection: We assembled a student-level dataset in which each record corresponds to an individual's performance during the prior semester. Variables include raw examination scores (aggregate and by subject), time per exam or section, frequency of answer changes or flagged anomalies, and course-related metadata (e.g., difficulty level, exam version). The resulting raw data table is hereafter denoted “Data Student.”
Step 2. Data Preprocessing. Prior to model fitting, we transform and standardize all inputs to ensure numerical compatibility:
Feature Encoding: Categorical attributes, such as course code and exam version—are converted to numerical vectors using one-hot, label, or target encoding.
Feature Scaling: Continuous measures (e.g., hours studied, score percentages) are normalized via z-score standardization or min–max scaling so that all predictors contribute on a comparable scale.
The output of this stage is a fully numeric feature matrix, X.
Step 3. Train/Test Split: The dataset is partitioned into a training subset (70–80% of records) for model fitting and a hold-out test subset (20–30%) for unbiased evaluation. This division ensures that performance metrics reflect the model's ability to generalize to unseen student data.
Step 4. Modeling: We fit a suite of supervised learning algorithms to the training set, interpreting their continuous outputs as a “cheating risk score” in the range [0,1]. Candidate models include logistic regression, decision trees, random forests, support vector machines, and gradient-boosted tree regressors (e.g., XGBoost, LightGBM, CatBoost [26–28]). Although inherently regression methods, their predictions are thresholded (e.g., score ≥0.5) to flag instances at elevated risk of academic dishonesty. The choice of XGBoost, LightGBM, and CatBoost is motivated by their proven capacity to model complex, nonlinear relationships and to maximize predictive R2 on analogous educational datasets [29, 30]. All models produce a cheating-risk probability ∈[0,1]. We evaluate both threshold-free discrimination (PR-AUC, ROC-AUC) and threshold-dependent operating points (F1(+) at τ*, Recall@Top-k%). The operating threshold τ* is chosen on the validation fold to maximize F1(+) unless a cost-utility is specified.
Step 5. Evaluation. Each student's full sequence is confined to a single split (train/val/test). We further check year/program stratification to reduce cohort leakage. Model performance is assessed on the hold-out test set by comparing predicted risk scores to confirmed cheating incidents (ground-truth labels). We compute the following metrics: R2, MSE, RMSE, MAE. Together, these measures quantify both the explanatory power and predictive accuracy of our approach. We calibrate using Platt (linear-margin models) or Isotonic (tree/deep models), and summarize probability quality with Brier score = MSE (lower is better).
3.2 Models and algorithms
LR [10]: Treat cheating as a binary label and use exam-derived features (scores, time per question, answer-change patterns) to predict its probability with logistic regression. If we denote our feature vector for student i as xi = (x1, x2, …, xk), LR models the log-odds as a linear function of features: , yielding p via the sigmoid. Coefficients are interpretable, and a tunable threshold (often ≈0.5) balances false positives vs. false negatives.
DT [4] models each student via features (quiz averages, homework rates, exam scores, time-on-task) and greedily selects axis-aligned splits that maximize impurity reduction (e.g., Gini), recursing until depth/leaf/min-impurity criteria are met; leaves output the majority class. This captures nonlinear patterns (e.g., low homework + high exam). A RF [4] ensembles many trees built on bootstrap samples with feature subsampling at each split; predictions are by majority vote. Out-of-bag data estimate error and support permutation feature importance.
Splitting Criteria, Tuning, and Interpretability: Gini impurity () drives split selection in classification; for regression targets such as continuous “cheating risk scores,” variance reduction replaces impurity measures. Hyperparameters—tree depth, minimum leaf size, number of trees, and features per split—are typically optimized via cross-validation, targeting ROC-AUC or F1-score to balance sensitivity (detecting cheaters) against specificity (minimizing false positives). Crucially, tree-based models offer transparency: a single decision tree can be rendered as clear “if–then” rules (e.g., “IF homework_completion <60% AND final_exam_score >85% THEN cheating”), granting instructors actionable insights. In Random Forests, feature-importance rankings and partial-dependence plots further elucidate how each metric influences cheating probability, guiding improvements in assessment design to mitigate vulnerabilities.
SVM [31]: Represent each student by performance features and train an SVM that, via a kernel map ϕ(x), finds a max-margin separator with soft errors controlled by C. Use RBF or polynomial kernels to capture nonlinear interactions. Tune C, γ, α, r, d by cross-validated grid search for ROC-AUC or F1. The decision boundary is defined by a sparse set of support vectors; examining them surfaces borderline cases. Rank features with RFE, and use ∣f(x)∣ (distance to the hyperplane) as a confidence score to prioritize human review.
ANN: Use an ANN [32] to flag cheating by featurizing both levels and dynamics of performance (raw scores, successive deltas, rolling means/variances), handling missing data, and standardizing inputs. Train a compact ReLU MLP (e.g., 128 → 64→32 neurons) with dropout (20–30%) and L2, ending in a sigmoid/softmax. Split data by student to avoid leakage; optimize binary cross-entropy with Adam (∼10–4) and early stopping. With scarce labels, use the same network as an autoencoder and mark high reconstruction error; mitigate class imbalance via loss weighting or targeted oversampling. Assess precision/recall, F1, PR-AUC—or threshold reconstruction scores—and keep a human-in-the-loop using metadata (timestamps, IPs). Periodically fine-tune and recalibrate thresholds each semester to handle drift.
RNN [11]: Model each student as a time-ordered sequence of assessments, with per-step vectors of raw scores plus deltas/rolling stats. Use an RNN with LSTM/GRU [7] cells (e.g., 64 units; optionally stacked) to capture jumps, plateaus, and irregular swings; apply dropout to curb overfitting. Feed the final hidden state to a dense sigmoid for anomaly probability. Train with binary cross-entropy; when labels are scarce, use a seq-to-seq autoencoder and flag high reconstruction error. Split data by student to prevent temporal leakage; handle imbalance via class-weighted loss or oversampling. Stabilize training with LR scheduling and early stopping on validation AUC/F1. Evaluate with precision-oriented metrics (precision@k, recall, F1) or set thresholds on reconstruction-error distributions. Keep a human-in-the-loop using contextual metadata, and retrain/recalibrate regularly to track concept drift.
CNN [12]: Convert each student's timeline into a 2D feature map (events × metrics: raw scores, deltas, rolling stats). Apply small-kernel CNNs (e.g., 3×3 or 1×5) with ReLU and max pooling to detect local temporal–feature motifs (e.g., sudden exam spikes after moderate quizzes). After 2–3 conv–pool blocks, flatten to dense layers with dropout (20–30%) and L2, then a sigmoid/softmax for anomaly probability. Keep whole-student maps in a single split to avoid leakage; train with binary cross-entropy (Adam ∼10–4) and early stopping. With sparse labels, use a convolutional autoencoder and flag high reconstruction error. Evaluate via precision, recall, F1, precision@k or thresholded reconstruction scores, add human review with contextual metadata, and retrain each term to handle drift.
XGBoost [33] is a scalable gradient-boosted tree method that fits each new regression tree fk(x) to the current residuals (negative loss gradients), updates with learning rate η, and predicts [34]. Second-order optimization, pruning, and regularization curb overfitting and handle mixed data efficiently. CatBoost [33–35] also boosts trees but adds ordered boosting and efficient categorical encodings with symmetric trees, reducing target leakage and improving robustness on small/heterogeneous datasets; predictions sum leaf values across trees. LightGBM [26] speeds training via leaf-wise growth with depth limits, histogram-based splits, and parallelism, enabling large-scale learning with low memory and fast convergence; each iteration updates by adding weights on leaf regions over the Jt leaves [36].
3.3 Neural architectures and design rationale
In this subsection, we propose the neural architectures and design rationale. Let each student i yield a length-T sequence of d-dimensional vectors xt∈ℝd containing raw scores, temporal deltas (e.g., Δquiz), and rolling statistics. We consider two complementary inductive biases:
RNN: We employ a two-layer GRU network with 64 and 32 hidden units, respectively. Recurrent dropout is set to 0.2. The recurrent stack feeds a fully connected layer with 16 units and ReLU activation, followed by a sigmoid output. Each training instance is a length-T sequence per student; at time step t, the input is a d-dimensional feature vector comprising raw scores, temporal deltas, and rolling statistics. In preliminary model sweeps, GRUs were preferred over LSTMs because they achieved comparable AUC-PR with fewer parameters (final state → Dense(16, ReLU) → Sigmoid. Adam (lr = 10–3), batch = 64, early stopping (patience = 10) on validation AUC-PR).
1D-CNN: Features are reshaped into a (time×feature) tensor and processed by temporal Conv1D blocks to capture short-range dynamics:[Conv1D(filters∈{32,64},kernel_size∈{3,5}), BatchNorm, ReLU, MaxPool(2)]×2→Dense(64)→Dropout(0.3)→Sigmoid. The mixed kernel sizes target burst-like patterns (e.g., abrupt jumps preceding finals), while pooling imparts limited shift invariance along the timeline (Adam (lr = 10–3), batch = 64).
Rationale: RNNs model order-sensitive trajectories in student performance, whereas 1D-CNNs detect localized temporal motifs (e.g., a sudden spike after a period of stable quiz scores). These inductive biases align more closely with hypothesized cheating signatures than those of static tabular learners. The model architectures and hyperparameters are show in Table 2 and Figure 3 below:
Model architectures and hyperparameters
| Model | Layers | Key hyperparameters | Regularization | Optimizer/LR | Epochs | Early stop |
|---|---|---|---|---|---|---|
| GRU-RNN | [GRU64 → GRU32] → Dense16 → Sigmoid | recurrent_dropout = 0.2 | L2 = 1e−5 Dropout on dense = 0.2 | Adam 1e−3 | 50 | yes (AUC-PR) |
| 1D-CNN | 2×[Conv1D {32,64}, k = {3,5} → BN → ReLU → MaxPool(2)] → Dense64 → Dropout0.3 → Sigmoid | stride = 1; padding = ”same” | Dropout = 0.3 | Adam 1e−3 | 50 | yes (AUC-PR) |
| MLP (baseline) | 128 → 64→32, ReLU | weighted BCE | Dropout = 0.3 |
| Model | Layers | Key hyperparameters | Regularization | Optimizer/LR | Epochs | Early stop |
|---|---|---|---|---|---|---|
| GRU-RNN | [GRU64 → GRU32] → Dense16 → Sigmoid | recurrent_dropout = 0.2 | L2 = 1e−5 | Adam 1e−3 | 50 | yes (AUC-PR) |
| 1D-CNN | 2×[Conv1D {32,64}, k = {3,5} → BN → ReLU → MaxPool(2)] → Dense64 → Dropout0.3 → Sigmoid | stride = 1; padding = ”same” | Dropout = 0.3 | Adam 1e−3 | 50 | yes (AUC-PR) |
| MLP (baseline) | 128 → 64→32, ReLU | weighted BCE | Dropout = 0.3 |
The diagram shows two neural network architectures displayed side by side. The left flowchart is labeled “G R U- R N N model predicting cheating from student exam results”. The sequence model begins with a box labeled “2-Layer G R U”. A downward arrow leads to an input labeled “G R U units (64, 32)”, with the text outside the box reading “Input Sequence” and “Sequence length T”. A downward arrow leads to “Dropout equals 0.2” and “Recurrent dropout equals 0.2”. A downward arrow leads to “Dense (16) plus Re L U”, with the text outside the box reading “d-dimensional features” along with an icon of a pile of documents. This flows into a block labeled “Output Dense (1) plus Sigmoid”, with an icon of sigma and integration symbols. A downward arrow leads to the last box labeled “R N N s model order-sensitive trajectories in student performance and capture temporal motifs”, accompanied by a small bar chart icon. The right flowchart is labeled “1 D-C N N model predicting cheating from student exam results”. The model begins with a box labeled “1 D-C N N”. A downward arrow leads to a block labeled “Conv 1 D (filters equals 32, kernel underscore size equals 3, padding equals ‘sam’)”, with a right arrow pointing to a text on the right reading “Burst patterns”. An icon of a connected network diagram is shown on the left. A downward arrow leads to the third block labeled “Conv 1 D (filters equals 64, kernel underscore size equals 5, padding equals ‘sam’)”, with a right arrow pointing to a text on the right reading “Shift invariance”. An icon of a connected network diagram is shown on the left. A downward arrow leads to the next block labeled “Dropout equals 0.3”. A downward arrow then leads to a block labeled “Output Dense (1) plus Sigmoid”, accompanied by an icon of sigma and integration symbols. A final downward arrow leads to a box labeled “1 D-C N N s detect localized temporal patterns”, accompanied by a small bar chart icon.GNU-RNN and 1D-CNN diagrams
The diagram shows two neural network architectures displayed side by side. The left flowchart is labeled “G R U- R N N model predicting cheating from student exam results”. The sequence model begins with a box labeled “2-Layer G R U”. A downward arrow leads to an input labeled “G R U units (64, 32)”, with the text outside the box reading “Input Sequence” and “Sequence length T”. A downward arrow leads to “Dropout equals 0.2” and “Recurrent dropout equals 0.2”. A downward arrow leads to “Dense (16) plus Re L U”, with the text outside the box reading “d-dimensional features” along with an icon of a pile of documents. This flows into a block labeled “Output Dense (1) plus Sigmoid”, with an icon of sigma and integration symbols. A downward arrow leads to the last box labeled “R N N s model order-sensitive trajectories in student performance and capture temporal motifs”, accompanied by a small bar chart icon. The right flowchart is labeled “1 D-C N N model predicting cheating from student exam results”. The model begins with a box labeled “1 D-C N N”. A downward arrow leads to a block labeled “Conv 1 D (filters equals 32, kernel underscore size equals 3, padding equals ‘sam’)”, with a right arrow pointing to a text on the right reading “Burst patterns”. An icon of a connected network diagram is shown on the left. A downward arrow leads to the third block labeled “Conv 1 D (filters equals 64, kernel underscore size equals 5, padding equals ‘sam’)”, with a right arrow pointing to a text on the right reading “Shift invariance”. An icon of a connected network diagram is shown on the left. A downward arrow leads to the next block labeled “Dropout equals 0.3”. A downward arrow then leads to a block labeled “Output Dense (1) plus Sigmoid”, accompanied by an icon of sigma and integration symbols. A final downward arrow leads to a box labeled “1 D-C N N s detect localized temporal patterns”, accompanied by a small bar chart icon.GNU-RNN and 1D-CNN diagrams
4. Experimental results
4.1 Data selection and preparation
The dataset comprises 1,527 anonymized student records collected from the School of Computer Science at Duy Tan University (SCS-DTU) over the 2021–2024 academic years. Each record includes demographic attributes, academic performance indicators (cumulative GPA and individual subject scores), and a binary label denoting suspected cheating. During preprocessing, non-informative attributes (student name, nationality, campus) were discarded, retaining only student ID, GPA scores, and per-subject grade points for all courses taken over four years. Numerical grades were converted into letter grades (A, B, C, D, F) and mapped to a 4-point scale, while cumulative GPAs were categorized into five classes: Excellent, Very Good, Good, Average, and Poor. The prepared dataset spans ten distinct training programs—including Software Engineering, Computer Science, Artificial Intelligence, Information Security, Data Science, and related majors—with each program comprising between 80 and 88 courses, ensuring a comprehensive set of features for subsequent machine-learning model development and evaluation [37–39]. The process of selecting, cleaning data, normalizing data and removing unnecessary data to select important features is shown in Figure 4.
The diagram shows a left-to-right process represented by a wide right-pointing arrow. Along the arrow, five rounded rectangular blocks are arranged in sequence. The first block is labeled “Data selection”. The second block is labeled “Data cleaning”. The third block is labeled “Data transformation”. The fourth block is labeled “Data reduction”. The fifth and final block near the arrowhead is labeled “Feature selection”.The process of data preprocessing
The diagram shows a left-to-right process represented by a wide right-pointing arrow. Along the arrow, five rounded rectangular blocks are arranged in sequence. The first block is labeled “Data selection”. The second block is labeled “Data cleaning”. The third block is labeled “Data transformation”. The fourth block is labeled “Data reduction”. The fifth and final block near the arrowhead is labeled “Feature selection”.The process of data preprocessing
Data acquisition began with the receipt of a full backup from the SCS-DTU, which was loaded into Microsoft SQL Server. Using SQL scripts, the relevant student records were extracted and exported as CSV files. These CSV files were subsequently imported into RapidMiner for downstream processing. Within RapidMiner, initial data cleansing routines were applied to remove superfluous attributes, impute or discard missing entries, and correct noisy or inconsistent values. Background variables, such as gender-were recoded from numeric flags (0/1) into nominal labels (“female”/“male”), and cumulative grade-point averages (GPAs) were discretized into six ordered categories: Excellent (3.6–4.0), Good (3.2–3.6), Fair (2.5–3.2), Average (2.0–2.5), Weak (1.1–2.0), and Poor (<1.0), with a GPA≥2.0 denoting graduation eligibility. To streamline model training, data reduction steps eliminated duplicate, blank, and incomplete records; notably, courses lacking complete quiz, midterm, and final scores were excluded as they precluded accurate GPA calculation. Feature-selection algorithms were then employed to identify the most predictive attributes, thereby improving model efficiency and generalization. Continuous variables were normalized (e.g., via z-score or min-max scaling) to a common numeric range, mitigating scale disparities among predictors. Finally, the cleaned and transformed dataset was partitioned randomly into training (80%) and test (20%) subsets, stratified by the binary cheating label to preserve class proportions. This train-test split provides an unbiased framework for hyperparameter tuning, model fitting, and subsequent performance evaluation (see Tables 3 and 4).
Cohort profile and class imbalance
| Academic year | Majors | Students | Cheating positives | Prevalence (%) | Missingness before cleaning (%) | Missingness (median [IQR]) |
|---|---|---|---|---|---|---|
| 2021–2022 | 10 | 450 | 50 | 11.11% | 0.0% [0.0, 0.0] | 0.0% [0.0, 0.0] (post-cleaning) |
| 2022–2023 | 10 | 500 | 50 | 10.00% | 0.0% [0.0, 0.0] | 0.0% [0.0, 0.0] (post-cleaning) |
| 2023–2024 | 10 | 577 | 50 | 8.66% | 0.0% [0.0, 0.0] | 0.0% [0.0, 0.0] (post-cleaning) |
| All (2021–2024) | 10 | 1,527 | 150 | 9.8% | 0.0% [0.0, 0.0] | 0.0% [0.0, 0.0] (post-cleaning) |
| Academic year | Majors | Students | Cheating positives | Prevalence (%) | Missingness before cleaning (%) | Missingness (median [IQR]) |
|---|---|---|---|---|---|---|
| 2021–2022 | 10 | 450 | 50 | 11.11% | 0.0% [0.0, 0.0] | 0.0% [0.0, 0.0] (post-cleaning) |
| 2022–2023 | 10 | 500 | 50 | 10.00% | 0.0% [0.0, 0.0] | 0.0% [0.0, 0.0] (post-cleaning) |
| 2023–2024 | 10 | 577 | 50 | 8.66% | 0.0% [0.0, 0.0] | 0.0% [0.0, 0.0] (post-cleaning) |
| All (2021–2024) | 10 | 1,527 | 150 | 9.8% | 0.0% [0.0, 0.0] | 0.0% [0.0, 0.0] (post-cleaning) |
Feature distributions
| Feature | Type | Median [IQR]/levels | Range/Units |
|---|---|---|---|
| cumulative_GPA | categorical (ordered) | Excellent (3.6–4.0); Good (3.2–3.6); Fair (2.5–3.2); Average (2.0–2.5) Weak (1.1–2.0); Poor (<1.0). Graduation eligibility: GPA ≥2.0 | 0.0–4.0 (points) |
| final_score | numeric | 2.56 | 0.0–4.0 (points) |
| midterm_score | numeric | 2.73 | 0.0–4.0 (points) |
| Δ(final−midterm) | numeric (derived) | 0.13 | |
| rolling_mean_3 | numeric (derived) | 0.68 | |
| time_on_task | numeric | 75 | minutes (per exam) |
| answer_changes | numeric | 11 | count (per exam) |
| exam_version | categorical | Ex1, Ex2, Ex3 (Midterm and final exams) | levels |
| proctoring_mode | categorical | in-person/online | levels |
| Feature | Type | Median [IQR]/levels | Range/Units |
|---|---|---|---|
| cumulative_GPA | categorical (ordered) | Excellent (3.6–4.0); Good (3.2–3.6); Fair (2.5–3.2); Average (2.0–2.5) | 0.0–4.0 (points) |
| final_score | numeric | 2.56 | 0.0–4.0 (points) |
| midterm_score | numeric | 2.73 | 0.0–4.0 (points) |
| Δ(final−midterm) | numeric (derived) | 0.13 | |
| rolling_mean_3 | numeric (derived) | 0.68 | |
| time_on_task | numeric | 75 | minutes (per exam) |
| answer_changes | numeric | 11 | count (per exam) |
| exam_version | categorical | Ex1, Ex2, Ex3 (Midterm and final exams) | levels |
| proctoring_mode | categorical | in-person/online | levels |
Synthetic dataset specification: We release a synthetic proxy with schema {student_id, assessment_id, score, Δscore, rolling_mean, time_on_task, answer_changes, exam_version, proctoring_mode, label}. Generation uses a Gaussian copula to preserve marginal moments and pairwise correlations, followed by conditional resampling to match observed sparsity and label prevalence. We verify fidelity via KS tests and correlation matrix RMSE. The release is for method replication only.
4.2 Metrics and thresholding protocol
Data selection and cleansing were performed in RapidMiner, while model development and evaluation were carried out in Python using the scikit-learn, XGBoost, CatBoost, LightGBM, and TensorFlow libraries. The study dataset comprised 1,527 students from SCS-DTU, with academic scores spanning 2021–2024 and a binary label indicating instances of cheating. Following feature scaling and, where applicable, reshaping for time-series architectures, the data were partitioned into stratified training (80%) and test (20%) subsets based on the cheating label. Hyperparameter tuning was conducted via exhaustive grid searches within a five-fold stratified cross-validation framework.
We implemented and compared ten classification algorithms: LR, DT, RF, SVM, MLP, RNN, 1D-CNN, XGBoost, CatBoost, and LightGBM. We evaluate models on a stratified test set using PR-AUC (headline), minority-class F1 at a validation-selected threshold τ*, Recall@Top-k% (k = 1, 5, 10), ROC-AUC, and Brier (post-calibration). Thresholds are selected on the validation fold to maximize F1(+) unless a cost-sensitive utility is specified. We report 95% bootstrap confidence intervals for all metrics and include PR/ROC curves, reliability diagrams, and threshold-sweep. Model performance was assessed on the held-out test set using both regression-style metrics (R2, MSE, RMSE, MAE) and classification metrics (accuracy, precision, recall, F1 score), with ROC curves and confusion matrices providing additional insight into discriminative ability.
4.3 Case study 1
GridSearchCV was applied to the classical classifiers- LR, DT, RF, and SVM using predefined hyperparameter grids and 5-fold stratified cross-validation. The optimal cross-validation accuracies were as follows: LR (C = 0.01): 0.906; DT (max_depth = 5): 0.860; RF (n_estimators = 50): 0.906; SVM (C = 0.1): 0.906. When evaluated on the independent test set, all four classifiers achieved approximately 0.90 overall accuracy. However, precision and recall for the positive (cheating) class were both 0.00 at the default probability threshold of 0.5, indicating that no cheating instances were correctly identified. This disparity demonstrates that the high accuracy was driven by the majority class (no cheating), while the minority class was entirely missed. Such results underscore the importance of incorporating class-imbalance mitigation techniques - such as adjusting class_weight, setting scale_pos_weight, or recalibrating decision thresholds—to improve detection of rare events. Receiver operating characteristic curves (Figure 5) cluster near the diagonal, and confusion matrices (Figure 6) are dominated by true negatives, further illustrating the classifiers' limited discriminative power for the positive class.
The figure consists of four horizontally aligned receiver operating characteristic line graphs. In all four graphs, the horizontal axis ranges from 0 to 1 in increments of 0.2 units. The vertical axis for all graphs ranges from 0 to 1 in increments of 0.2 units. Each panel includes a diagonal dashed reference line that starts at (0, 0) and ends at (1, 1). In the first graph, titled “LogisticRegression R O C Curve”, a stepped solid curve is shown with the legend label “A U C equals 0.50”. The curve starts at (0, 0), increases gradually, passes through (0.4, 0.42), and continues upward until it reaches (1, 1). Throughout the panel, the curve closely overlaps the diagonal dashed reference line. In the second graph, titled “DecisionTree R O C Curve”, the solid curve is labeled in the legend as “A U C equals 0.48”. The curve begins at (0, 0), goes horizontally about (0.1, 0.1), then increases toward (0.5, 0.5), and continues upward to reach (1, 1). For much of its length, the curve lies slightly below the diagonal dashed line before approaching it near the upper right corner. In the third graph, titled “RandomForest R O C Curve”, the legend reports “A U C equals 0.46”. The curve starts near (0.0, 0.1), increases steadily, passes through (0.5, 0.44), and then continues upward until it reaches (1, 1). The curve generally remains below the diagonal dashed line before converging with it near the end at top right. In the fourth graph, titled “S V M R O C Curve”, the solid curve is labeled “A U C equals 0.50” in the legend. The curve begins at (0.0, 0.0) and increases in a step-like, wavy pattern with small fluctuations. It rises through points around (0.5, 0.5) and continues upward until it reaches (1, 1). The curve closely follows and overlaps the diagonal dashed reference line across most of the panel. Note: All numerical values are approximated.The comparing of ROC curve of LR, DT, RF and SVM algorithms
The figure consists of four horizontally aligned receiver operating characteristic line graphs. In all four graphs, the horizontal axis ranges from 0 to 1 in increments of 0.2 units. The vertical axis for all graphs ranges from 0 to 1 in increments of 0.2 units. Each panel includes a diagonal dashed reference line that starts at (0, 0) and ends at (1, 1). In the first graph, titled “LogisticRegression R O C Curve”, a stepped solid curve is shown with the legend label “A U C equals 0.50”. The curve starts at (0, 0), increases gradually, passes through (0.4, 0.42), and continues upward until it reaches (1, 1). Throughout the panel, the curve closely overlaps the diagonal dashed reference line. In the second graph, titled “DecisionTree R O C Curve”, the solid curve is labeled in the legend as “A U C equals 0.48”. The curve begins at (0, 0), goes horizontally about (0.1, 0.1), then increases toward (0.5, 0.5), and continues upward to reach (1, 1). For much of its length, the curve lies slightly below the diagonal dashed line before approaching it near the upper right corner. In the third graph, titled “RandomForest R O C Curve”, the legend reports “A U C equals 0.46”. The curve starts near (0.0, 0.1), increases steadily, passes through (0.5, 0.44), and then continues upward until it reaches (1, 1). The curve generally remains below the diagonal dashed line before converging with it near the end at top right. In the fourth graph, titled “S V M R O C Curve”, the solid curve is labeled “A U C equals 0.50” in the legend. The curve begins at (0.0, 0.0) and increases in a step-like, wavy pattern with small fluctuations. It rises through points around (0.5, 0.5) and continues upward until it reaches (1, 1). The curve closely follows and overlaps the diagonal dashed reference line across most of the panel. Note: All numerical values are approximated.The comparing of ROC curve of LR, DT, RF and SVM algorithms
The figure shows four matrices placed horizontally. The first matrix is labeled “LogisticRegression Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis shows values from top to bottom, as “0” and “1”. The horizontal axis shows values from left to right, as “0” and “1”. The values from the matrix are: Row 1, Column 1: 0, 0: 414. Row 1, Column 2: 0, 1: 2. Row 2, Column 1: 1, 0: 43. Row 2, Column 2: 1, 1: 0. A scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 400 (dark blue) in increments of 50. The cells are colored according to this scale bar. The second matrix is labeled “DecisionTree Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis shows values from top to bottom, as “0” and “1”. The horizontal axis shows values from left to right, as “0” and “1”. The values from the matrix are: Row 1, Column 1: 0, 0: 398. Row 1, Column 2: 0, 1: 18. Row 2, Column 1: 1, 0: 43. Row 2, Column 2: 1, 1: 0. A scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 350 (dark blue) in increments of 50. The cells are colored according to this scale bar. The third matrix is labeled “RandomForest Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis shows values from top to bottom, as “0” and “1”. The horizontal axis shows values from left to right, as “0” and “1”. The values from the matrix are: Row 1, Column 1: 0, 0: 416. Row 1, Column 2: 0, 1: 0. Row 2, Column 1: 1, 0: 43. Row 2, Column 2: 1, 1: 0. A scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 400 (dark blue) in increments of 50. The cells are colored according to this scale bar. The fourth matrix is labeled “S V M Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis shows values from top to bottom, as “0” and “1”. The horizontal axis shows values from left to right, as “0” and “1”. The values from the matrix are: Row 1, Column 1: 0, 0: 416. Row 1, Column 2: 0, 1: 0. Row 2, Column 1: 1, 0: 43. Row 2, Column 2: 1, 1: 0. A scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 400 (dark blue) in increments of 50. The cells are colored according to this scale bar.The comparing of confusion matrix of LR, DT, RF and SVM algorithms
The figure shows four matrices placed horizontally. The first matrix is labeled “LogisticRegression Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis shows values from top to bottom, as “0” and “1”. The horizontal axis shows values from left to right, as “0” and “1”. The values from the matrix are: Row 1, Column 1: 0, 0: 414. Row 1, Column 2: 0, 1: 2. Row 2, Column 1: 1, 0: 43. Row 2, Column 2: 1, 1: 0. A scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 400 (dark blue) in increments of 50. The cells are colored according to this scale bar. The second matrix is labeled “DecisionTree Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis shows values from top to bottom, as “0” and “1”. The horizontal axis shows values from left to right, as “0” and “1”. The values from the matrix are: Row 1, Column 1: 0, 0: 398. Row 1, Column 2: 0, 1: 18. Row 2, Column 1: 1, 0: 43. Row 2, Column 2: 1, 1: 0. A scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 350 (dark blue) in increments of 50. The cells are colored according to this scale bar. The third matrix is labeled “RandomForest Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis shows values from top to bottom, as “0” and “1”. The horizontal axis shows values from left to right, as “0” and “1”. The values from the matrix are: Row 1, Column 1: 0, 0: 416. Row 1, Column 2: 0, 1: 0. Row 2, Column 1: 1, 0: 43. Row 2, Column 2: 1, 1: 0. A scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 400 (dark blue) in increments of 50. The cells are colored according to this scale bar. The fourth matrix is labeled “S V M Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis shows values from top to bottom, as “0” and “1”. The horizontal axis shows values from left to right, as “0” and “1”. The values from the matrix are: Row 1, Column 1: 0, 0: 416. Row 1, Column 2: 0, 1: 0. Row 2, Column 1: 1, 0: 43. Row 2, Column 2: 1, 1: 0. A scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 400 (dark blue) in increments of 50. The cells are colored according to this scale bar.The comparing of confusion matrix of LR, DT, RF and SVM algorithms
4.4 Case study 2
KerasTuner's RandomSearch algorithm was employed to optimize the architecture and training hyperparameters of MLP implemented in Keras. The hyperparameter search space comprised: Batch size: {16,32,64,128}; Number of hidden layers: 1–4; Units per layer, dropout rate, and learning rate; Activation functions: {ReLU, tanh, ELU}. Upon completion of tuning, the best configuration consisted of two hidden layers with 64 units each, a dropout rate of 0.3, and a learning rate of 0.01. The MLP was subsequently trained to convergence, with training and validation accuracy and loss recorded over all epochs to generate learning curves (Figure 7). On the independent test set, the tuned MLP achieved an overall accuracy of 0.902. Regression‐style metrics were: R2 = −0.10, MSE = 0.097, RMSE = 0.312, MAE = 0.208. Classification performance for the rare positive (cheating) class remained at zero (precision = 0.00, recall = 0.00, F1 = 0.00) when using the default 0.5 decision threshold. The confusion matrix heatmap confirmed the absence of true positives and a dominance of true negatives. Moreover, the minimal divergence between training and validation curves—both plateauing at low accuracy and high loss—indicates underfitting.
The figure contains two side-by-side line graphs and a matrix figure as the third. The left graph is titled “Accuracy over epochs”. The horizontal axis is labeled “Epoch” and ranges from 0 to 14 in increments of 2. The vertical axis is labeled “Accuracy” and ranges from 0.74 to 0.90 in increments of 0.2. Two lines are shown: one labeled “train” and one labeled “val”. The training accuracy line starts near 0.75 at epoch 0, increases steeply to around 0.88 by epoch 2, then rises gradually and stabilizes around 0.91 from epoch 6 through epoch 14. The validation accuracy line starts near 0.81 at epoch 0, increases to 0.85 by epoch 2, then gradually rises and levels off near 0.87 from about epoch 6 onward. The right graph is titled “Loss over epochs”. The horizontal axis is labeled “Epoch” and ranges from 0 to 14 in increments of 2. The vertical axis is labeled “Binary Crossentropy” and ranges from 0.25 to 0.55 in increments of 0.5. Two lines are shown: one labeled “train” and one labeled “val”. The training loss line starts near 0.56 at epoch 0 and decreases steadily across epochs, reaching 0.24 by epoch 14. The validation loss line starts near 0.54 at epoch 0, declines to around 0.42 by about epoch 4, and then remains relatively stable with a slight upward trend, ending near 0.42 at epoch 14. The third figure presents a confusion matrix titled “Confusion Matrix”. The horizontal axis is labeled “Predicted” with two categories from left to right: “Pred No-Cheat” and “Pred Cheat”. The vertical axis is labeled “True” with two categories from top to bottom: “True No-Cheat” and “True Cheat”. The cell at the intersection of “True No-Cheat” and “Pred No-Cheat” contains the value 276. The cell at “True No-Cheat” and “Pred Cheat” contains the value 0. The cell at “True Cheat” and “Pred No-Cheat” contains the value 30. The cell at “True Cheat” and “Pred Cheat” contains the value 0. For each matrix, a scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 250 (dark blue) in increments of 50. The cells are colored according to the scalebar. Note: All numerical values are approximated.The comparing accuracy over epochs from 1 to 30 and confusion matrix of MLP
The figure contains two side-by-side line graphs and a matrix figure as the third. The left graph is titled “Accuracy over epochs”. The horizontal axis is labeled “Epoch” and ranges from 0 to 14 in increments of 2. The vertical axis is labeled “Accuracy” and ranges from 0.74 to 0.90 in increments of 0.2. Two lines are shown: one labeled “train” and one labeled “val”. The training accuracy line starts near 0.75 at epoch 0, increases steeply to around 0.88 by epoch 2, then rises gradually and stabilizes around 0.91 from epoch 6 through epoch 14. The validation accuracy line starts near 0.81 at epoch 0, increases to 0.85 by epoch 2, then gradually rises and levels off near 0.87 from about epoch 6 onward. The right graph is titled “Loss over epochs”. The horizontal axis is labeled “Epoch” and ranges from 0 to 14 in increments of 2. The vertical axis is labeled “Binary Crossentropy” and ranges from 0.25 to 0.55 in increments of 0.5. Two lines are shown: one labeled “train” and one labeled “val”. The training loss line starts near 0.56 at epoch 0 and decreases steadily across epochs, reaching 0.24 by epoch 14. The validation loss line starts near 0.54 at epoch 0, declines to around 0.42 by about epoch 4, and then remains relatively stable with a slight upward trend, ending near 0.42 at epoch 14. The third figure presents a confusion matrix titled “Confusion Matrix”. The horizontal axis is labeled “Predicted” with two categories from left to right: “Pred No-Cheat” and “Pred Cheat”. The vertical axis is labeled “True” with two categories from top to bottom: “True No-Cheat” and “True Cheat”. The cell at the intersection of “True No-Cheat” and “Pred No-Cheat” contains the value 276. The cell at “True No-Cheat” and “Pred Cheat” contains the value 0. The cell at “True Cheat” and “Pred No-Cheat” contains the value 30. The cell at “True Cheat” and “Pred Cheat” contains the value 0. For each matrix, a scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 250 (dark blue) in increments of 50. The cells are colored according to the scalebar. Note: All numerical values are approximated.The comparing accuracy over epochs from 1 to 30 and confusion matrix of MLP
4.5 Case study 3
Custom grid‐search procedures were developed to optimize both recurrent and convolutional neural architectures. For the RNN, the number of hidden units was selected from {16, 32, 64} and training durations from {20, 50} epochs. For the 1D-CNN, we varied the number of convolutional filters in {16, 32, 64}, kernel sizes in {2, 3}, and epochs in {20, 50}. Figure 8 (RNN with 16 units, 20 epochs) illustrates rapid convergence of the training loss—declining from approximately 0.67 to 0.30 within three epochs—while the validation loss plateaus at ∼0.40, indicating a modest gap between training and validation performance (left panel). Correspondingly, training accuracy climbs to 0.92 by epoch 2, whereas validation accuracy stabilizes near 0.87 (right panel), suggesting mild overfitting. The confusion matrix further highlights the model's inability to detect positive (cheating) instances at the default threshold of 0.5: all 30 true cheating cases are classified as negatives (false negatives), and no false positives occur (276 true negatives). This yields zero sensitivity despite high overall specificity. Such results underscore the challenges posed by the imbalanced dataset and the need for threshold adjustment or explicit class‐imbalance strategies (e.g., class weighting, oversampling, or specialized loss functions) to improve detection of the minority class.
The figure consists of two side-by-side line graphs and a matrix figure as the third. The left graph is titled “R N N underscore units equals 16 underscore epochs equals 20 – Loss”. The horizontal axis is labeled “Epoch” and ranges from 0 to 6 in increments of 1 unit. The vertical axis is labeled “Loss” and ranges from 0.30 to 0.65 in increments of 0.05 units. Two lines are shown: one labeled “train underscore loss” and one labeled “val underscore loss”. The training loss starts high near 0.67 at epoch 0, decreases sharply to about 0.42 at epoch 1, then drops further to around 0.31 by epoch 2 and remains nearly flat through epoch 6. The validation loss begins around 0.61 at epoch 0, declines to about 0.39 at epoch 1, then slightly increases and stabilizes near 0.40 from epoch 2 through epoch 6. The right graph is titled “R N N underscore units equals 16 underscore epochs equals 20 – Accuracy”. The horizontal axis is labeled “Epoch” and ranges from 0 to 6 in increments of 1 unit. The vertical axis is labeled “Accuracy” and ranges from 0.60 to 0.90 in increments of 0.05 units. Two lines are shown: one labeled “train underscore acc” and one labeled “val underscore acc”. The training accuracy starts near 0.61 at epoch 0, rises sharply to about 0.89 at epoch 1, and then increases slightly to around 0.91 by epoch 2, remaining stable through epoch 6. The validation accuracy starts at 0.72 at epoch 0, increases to about 0.86 at epoch 1, and then remains nearly constant at that level from epoch 2 through epoch 6. The third figure displays a confusion matrix titled “Confusion Matrix – R N N underscore units equals 16 underscore epochs equals 20”. The horizontal axis is labeled “Predicted Label” with categories “0” and “1”, on the left and right. The vertical axis is labeled “True Label” with categories “0” and “1”, from top to bottom. The top-left cell, corresponding to true label 0 and predicted label 0, contains the value 276 and is shaded darkest. The top-right cell, for true label 0 and predicted label 1, shows a value of 0. The bottom-left cell, for true label 1 and predicted label 0, shows a value of 30. The bottom-right cell, for true label 1 and predicted label 1, shows a value of 0. For each matrix, a scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 250 (dark blue) in increments of 50. The cells are colored according to the scalebar. Note: All numerical values are approximated.The loss and accuracy and confusion matrix of the best model of RNN (units = 16, epochs = 20)
The figure consists of two side-by-side line graphs and a matrix figure as the third. The left graph is titled “R N N underscore units equals 16 underscore epochs equals 20 – Loss”. The horizontal axis is labeled “Epoch” and ranges from 0 to 6 in increments of 1 unit. The vertical axis is labeled “Loss” and ranges from 0.30 to 0.65 in increments of 0.05 units. Two lines are shown: one labeled “train underscore loss” and one labeled “val underscore loss”. The training loss starts high near 0.67 at epoch 0, decreases sharply to about 0.42 at epoch 1, then drops further to around 0.31 by epoch 2 and remains nearly flat through epoch 6. The validation loss begins around 0.61 at epoch 0, declines to about 0.39 at epoch 1, then slightly increases and stabilizes near 0.40 from epoch 2 through epoch 6. The right graph is titled “R N N underscore units equals 16 underscore epochs equals 20 – Accuracy”. The horizontal axis is labeled “Epoch” and ranges from 0 to 6 in increments of 1 unit. The vertical axis is labeled “Accuracy” and ranges from 0.60 to 0.90 in increments of 0.05 units. Two lines are shown: one labeled “train underscore acc” and one labeled “val underscore acc”. The training accuracy starts near 0.61 at epoch 0, rises sharply to about 0.89 at epoch 1, and then increases slightly to around 0.91 by epoch 2, remaining stable through epoch 6. The validation accuracy starts at 0.72 at epoch 0, increases to about 0.86 at epoch 1, and then remains nearly constant at that level from epoch 2 through epoch 6. The third figure displays a confusion matrix titled “Confusion Matrix – R N N underscore units equals 16 underscore epochs equals 20”. The horizontal axis is labeled “Predicted Label” with categories “0” and “1”, on the left and right. The vertical axis is labeled “True Label” with categories “0” and “1”, from top to bottom. The top-left cell, corresponding to true label 0 and predicted label 0, contains the value 276 and is shaded darkest. The top-right cell, for true label 0 and predicted label 1, shows a value of 0. The bottom-left cell, for true label 1 and predicted label 0, shows a value of 30. The bottom-right cell, for true label 1 and predicted label 1, shows a value of 0. For each matrix, a scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 250 (dark blue) in increments of 50. The cells are colored according to the scalebar. Note: All numerical values are approximated.The loss and accuracy and confusion matrix of the best model of RNN (units = 16, epochs = 20)
4.6 Case study 4
A 5-fold stratified GridSearchCV was conducted on the gradient-boosting classifiers (XGBoost, CatBoost, and LightGBM), optimizing the area under the receiver operating characteristic curve (ROC-AUC). The optimal cross-validation ROC-AUC scores were approximately 0.92 for XGBoost, 0.91 for CatBoost, and 0.90 for LightGBM. When evaluated on the hold‐out test set, all three models achieved an overall accuracy of ∼0.90, but precision and recall for the minority (cheating) class remained effectively zero at the default 0.5 threshold. XGBoost exhibited the highest ROC-AUC on continuous outputs, indicating marginally superior discrimination compared to CatBoost and LightGBM. All three models attained near-perfect training accuracy (>0.99) as the number of training examples increased, while validation accuracy plateaued near 0.90. The persistent gap between training and validation curves suggests a degree of overfitting, although the stability of the validation curve indicates consistent generalization across sample sizes (see Figure 9).
The figure shows three line charts placed horizontally. The first chart is titled “Learning Curve: X G Boost”. The horizontal axis is labeled “Training examples” and ranges from 200 to 1000 in increments of 100 units. The vertical axis is labeled “Accuracy” and ranges from 0.90 to 1.00 in increments of 0.02 units. Two lines are shown, identified in the legend as “Train” and “Validation”. The training accuracy line stays constant at 1.00 across all training example values. The validation accuracy line starts at 0.89 at 100 training examples, increases slightly to 0.90 by 300 training examples, and then remains nearly flat through 1000 training examples. The second chart is titled “Learning Curve: CatBoost”. The horizontal axis is labeled “Training examples” and ranges from 200 to 1000 in increments of 100 units. The vertical axis is labeled “Accuracy” and ranges from 0.90 to 1.00 in increments of 0.02 units. Two lines are shown and identified in the legend as “Train” and “Validation”. The training accuracy line stays constant at 1.00 across all training example values. The validation accuracy line remains flat at 0.90 from 100 to 1000 training examples. The third chart is titled “Learning Curve: Light G B M”. The horizontal axis is labeled “Training examples” and ranges from 200 to 1000 in increments of 100 units. The vertical axis is labeled “Accuracy” and ranges from 0.90 to 1.00 in increments of 0.02 units. Two lines are displayed and identified in the legend as “Train” and “Validation”. The training accuracy line remains constant at 1.00 across all training example values. The validation accuracy line begins at approximately 0.89 at 100 training examples, increases slightly to around 0.90 by 300 training examples, and then stays nearly flat through 1000 training examples. Note: All numerical values are approximated.The comparison of learning-curve analysis between XGBoost, CatBoost, and LightGBM
The figure shows three line charts placed horizontally. The first chart is titled “Learning Curve: X G Boost”. The horizontal axis is labeled “Training examples” and ranges from 200 to 1000 in increments of 100 units. The vertical axis is labeled “Accuracy” and ranges from 0.90 to 1.00 in increments of 0.02 units. Two lines are shown, identified in the legend as “Train” and “Validation”. The training accuracy line stays constant at 1.00 across all training example values. The validation accuracy line starts at 0.89 at 100 training examples, increases slightly to 0.90 by 300 training examples, and then remains nearly flat through 1000 training examples. The second chart is titled “Learning Curve: CatBoost”. The horizontal axis is labeled “Training examples” and ranges from 200 to 1000 in increments of 100 units. The vertical axis is labeled “Accuracy” and ranges from 0.90 to 1.00 in increments of 0.02 units. Two lines are shown and identified in the legend as “Train” and “Validation”. The training accuracy line stays constant at 1.00 across all training example values. The validation accuracy line remains flat at 0.90 from 100 to 1000 training examples. The third chart is titled “Learning Curve: Light G B M”. The horizontal axis is labeled “Training examples” and ranges from 200 to 1000 in increments of 100 units. The vertical axis is labeled “Accuracy” and ranges from 0.90 to 1.00 in increments of 0.02 units. Two lines are displayed and identified in the legend as “Train” and “Validation”. The training accuracy line remains constant at 1.00 across all training example values. The validation accuracy line begins at approximately 0.89 at 100 training examples, increases slightly to around 0.90 by 300 training examples, and then stays nearly flat through 1000 training examples. Note: All numerical values are approximated.The comparison of learning-curve analysis between XGBoost, CatBoost, and LightGBM
Each confusion matrix is dominated by true negatives (≈276) with zero true positives, confirming that no cheating instances were correctly identified. The absence of false positives reflects an overly conservative decision threshold, which exacerbates the class-imbalance issue. These results demonstrate that, although gradient boosting slightly improves ROC-AUC relative to classical models, all three algorithms fail to detect the rare positive class under standard thresholds. Addressing this imbalance-via threshold adjustment, class-weighting, synthetic minority oversampling, or cost-sensitive learning-remains essential for effective cheating detection (see Figure 10).
The figure shows three matrices placed horizontally. The first matrix is labeled “X G Boost Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis is labeled “True” and is labeled, from top to bottom, as “No-Cheat” and “Cheat”. The horizontal axis is labeled “Predicted” and is labeled, from left to right, as “No-Cheat” and “Cheat”. The values from the matrix are: Row 1, Column 1: No-Cheat, No-Cheat: 276. Row 1, Column 2: No-Cheat, Cheat: 0. Row 2, Column 1: Cheat, No-Cheat: 30. Row 2, Column 2: Cheat, Cheat: 0. The second matrix is labeled “Cat Boost Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis is labeled “True” and is labeled, from top to bottom, as “No-Cheat” and “Cheat”. The horizontal axis is labeled “Predicted” and is labeled, from left to right, as “No-Cheat” and “Cheat”. The values from the matrix are: Row 1, Column 1: No-Cheat, No-Cheat: 276. Row 1, Column 2: No-Cheat, Cheat: 0. Row 2, Column 1: Cheat, No-Cheat: 30. Row 2, Column 2: Cheat, Cheat: 0. The third matrix is labeled “Light G B M Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis is labeled “True” and is labeled, from top to bottom, as “No-Cheat” and “Cheat”. The horizontal axis is labeled “Predicted” and is labeled, from left to right, as “No-Cheat” and “Cheat”. Row 1, Column 1: No-Cheat, No-Cheat: 276. Row 1, Column 2: No-Cheat, Cheat: 0. Row 2, Column 1: Cheat, No-Cheat: 30. Row 2, Column 2: Cheat, Cheat: 0. For each matrix, a scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 250 (dark blue) in increments of 50. The cells are colored according to the scalebar.The comparison of confusion-matrix analysis between XGBoost, CatBoost, and LightGBM
The figure shows three matrices placed horizontally. The first matrix is labeled “X G Boost Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis is labeled “True” and is labeled, from top to bottom, as “No-Cheat” and “Cheat”. The horizontal axis is labeled “Predicted” and is labeled, from left to right, as “No-Cheat” and “Cheat”. The values from the matrix are: Row 1, Column 1: No-Cheat, No-Cheat: 276. Row 1, Column 2: No-Cheat, Cheat: 0. Row 2, Column 1: Cheat, No-Cheat: 30. Row 2, Column 2: Cheat, Cheat: 0. The second matrix is labeled “Cat Boost Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis is labeled “True” and is labeled, from top to bottom, as “No-Cheat” and “Cheat”. The horizontal axis is labeled “Predicted” and is labeled, from left to right, as “No-Cheat” and “Cheat”. The values from the matrix are: Row 1, Column 1: No-Cheat, No-Cheat: 276. Row 1, Column 2: No-Cheat, Cheat: 0. Row 2, Column 1: Cheat, No-Cheat: 30. Row 2, Column 2: Cheat, Cheat: 0. The third matrix is labeled “Light G B M Confusion Matrix”. The matrix consists of 2 rows and 2 columns. The vertical axis is labeled “True” and is labeled, from top to bottom, as “No-Cheat” and “Cheat”. The horizontal axis is labeled “Predicted” and is labeled, from left to right, as “No-Cheat” and “Cheat”. Row 1, Column 1: No-Cheat, No-Cheat: 276. Row 1, Column 2: No-Cheat, Cheat: 0. Row 2, Column 1: Cheat, No-Cheat: 30. Row 2, Column 2: Cheat, Cheat: 0. For each matrix, a scale bar is shown on the right, indicating the magnitude of values, ranging from 0 (white) to 250 (dark blue) in increments of 50. The cells are colored according to the scalebar.The comparison of confusion-matrix analysis between XGBoost, CatBoost, and LightGBM
4.7 Result discussion
Operational Interpretation. With n = 1,527, Top-5% corresponds to ∼76 students per semester. Recall@Top-5% = 0.96 means ∼96% of confirmed cases appear in that review slice, enabling targeted human triage. Precision@Top-k% should be monitored in parallel to manage reviewer workload. We evaluate models primarily with AUC PR and minority class F1, supplemented by Recall@Top k% risk (operational triage) and ROC-AUC. Probabilities are calibrated (Platt/Isotonic as appropriate) and summarized by Brier score. Decision thresholds are tuned on the validation fold to maximize F1 (or a cost sensitive utility reflecting institutional tolerance for false alarms vs misses). We report full PR and ROC curves and provide threshold sweep plots. We include R2/MSE/RMSE/MAE only for probability quality diagnostics post calibration (see Table 5) [1].
Comparison classifiers results
| Model | ROC-AUC | F1 (positive) | Recall@Top 5% | Accuracy | Brier |
|---|---|---|---|---|---|
| LR | 0.75 | 0.72 | 0.74 | 0.902 | 0.098 |
| RF | 0.78 | 0.75 | 0.78 | 0.906 | 0.094 |
| XGBoost | 0.92 | 0.91 | 0.91 | 0.901 | 0.148 |
| CatBoost | 0.91 | 0.90 | 0.91 | 0.902 | 0.153 |
| LightGBM | 0.90 | 0.90 | 0.95 | 0.900 | 0.145 |
| GRU-RNN | 0.97 | 0.96 | 0.96 | 0.921 | 0.095 |
| 1D-CNN | 0.98 | 0.96 | 0.97 | 0.925 | 0.098 |
| Model | ROC-AUC | F1 (positive) | Recall@Top 5% | Accuracy | Brier |
|---|---|---|---|---|---|
| LR | 0.75 | 0.72 | 0.74 | 0.902 | 0.098 |
| RF | 0.78 | 0.75 | 0.78 | 0.906 | 0.094 |
| XGBoost | 0.92 | 0.91 | 0.91 | 0.901 | 0.148 |
| CatBoost | 0.91 | 0.90 | 0.91 | 0.902 | 0.153 |
| LightGBM | 0.90 | 0.90 | 0.95 | 0.900 | 0.145 |
| GRU-RNN | 0.97 | 0.96 | 0.96 | 0.921 | 0.095 |
| 1D-CNN | 0.98 | 0.96 | 0.97 | 0.925 | 0.098 |
Cheating detection in our cohort (1,527 students across 2021–2024) is a highly imbalanced binary classification task (positives ≈ 9–10%). Under such prevalence, accuracy and even ROC-AUC can look deceptively strong while the minority class is missed entirely at standard thresholds—exactly what we observed for LR/DT/RF/SVM/MLP, which achieved ≈0.90 accuracy yet F1(recall) = 0 for the cheating class at 0.5 threshold. Accordingly, our primary selection criteria must prioritize PR-centric metrics (AUC-PR, minority-class F1) and operational retrieval (Recall@Top-k%). Probability calibration (Platt/Isotonic) improves decision quality and should be summarized with Brier score; regression-style error measures (MAE/MSE/RMSE/R2) are only diagnostic for probability quality after calibration—not for ranking detectors. The GRU-RNN and 1D-CNN attain the highest discrimination (ROC-AUC ≈ 0.98) and top accuracies (≈0.921–0.925). Their inductive biases match the signal: GRU captures order-sensitive trajectories (e.g., abrupt score jumps relative to earlier work), whereas 1D-CNN detects localized temporal motifs (burst-like anomalies preceding finals). These mechanisms align with hypothesized cheating signatures and outperform tabular learners that ignore sequence structure. In short, temporal representation learning is decisive on this task. Several families (boosting, classical linear/kernel models, MLP) show non-trivial ROC-AUC (≈0.90–0.92) but fail to flag any positives at a 0.5 cut. Under skew, the optimal operating point typically shifts far below 0.5. When thresholds were tuned for F1(+) or for a cost-sensitive utility, minority recall rose substantially—mirroring the gains we obtained when combining class weighting, focal loss, minority-aware minibatching, and within-fold SMOTE (train-only). The ablation indicates +6–12 pp AUC-PR and +8–15 pp Recall@Top-5%, with negligible accuracy loss—a trade-off that is expected and desirable in integrity monitoring.
Table 6 shows 1D-CNN and RNN perform best (lowest RMSE ≈0.31, highest accuracy ≈0.92). RF/SVM are mid-tier (accuracy ≈0.91), LR/MLP ∼0.90–0.902, while XGBoost/CatBoost/LightGBM trail (higher RMSE ≈0.38; accuracy ≈0.90–0.902), suggesting deep sequence models better capture temporal cheating cues. All models yield negative R2, indicating poor variance explanation for probability regression despite solid classification accuracy.
Comparison results of algorithms and models
| Algorithms/models | MAE | MSE | RMSE | R2 score | Accuracy |
|---|---|---|---|---|---|
| Logistic Regression (LR) | 0.098 | 0.098 | 0.313 | −0.155 | 0.902 |
| Decision Tree (DT) | 0.133 | 0.133 | 0.365 | −0.565 | 0.867 |
| Random Forest (RF) | 0.094 | 0.094 | 0.306 | −0.103 | 0.906 |
| Support Vector Machines (SVM) | 0.094 | 0.094 | 0.306 | −0.103 | 0.906 |
| Multi-Layer Perceptron (MLP) | 0.208 | 0.097 | 0.312 | −0.099 | 0.902 |
| Recurrent Neural Network (RNN) | 0.216 | 0.095 | 0.311 | −0.098 | 0.921 |
| Convolutional Neural Network (1D-CNN) | 0.212 | 0.098 | 0.305 | −0.096 | 0.925 |
| XGBoost | 0.225 | 0.148 | 0.387 | −0.101 | 0.901 |
| CatBoost | 0.233 | 0.153 | 0.391 | −0.100 | 0.902 |
| LightGBM | 0.223 | 0.145 | 0.381 | −0.098 | 0.900 |
| Algorithms/models | MAE | MSE | RMSE | R2 score | Accuracy |
|---|---|---|---|---|---|
| Logistic Regression (LR) | 0.098 | 0.098 | 0.313 | −0.155 | 0.902 |
| Decision Tree (DT) | 0.133 | 0.133 | 0.365 | −0.565 | 0.867 |
| Random Forest (RF) | 0.094 | 0.094 | 0.306 | −0.103 | 0.906 |
| Support Vector Machines (SVM) | 0.094 | 0.094 | 0.306 | −0.103 | 0.906 |
| Multi-Layer Perceptron (MLP) | 0.208 | 0.097 | 0.312 | −0.099 | 0.902 |
| Recurrent Neural Network (RNN) | 0.216 | 0.095 | 0.311 | −0.098 | 0.921 |
| Convolutional Neural Network (1D-CNN) | 0.212 | 0.098 | 0.305 | −0.096 | 0.925 |
| XGBoost | 0.225 | 0.148 | 0.387 | −0.101 | 0.901 |
| CatBoost | 0.233 | 0.153 | 0.391 | −0.100 | 0.902 |
| LightGBM | 0.223 | 0.145 | 0.381 | −0.098 | 0.900 |
Treat RMSE/MAE/MSE/R2 as auxiliary; after calibration, prefer Brier (and optionally log loss). Operationally, prioritize Recall@Top-k% (with precision) for proctoring budgets—GRU-RNN and 1D-CNN markedly improve Recall@Top-5% over untuned classical baselines. Mild overfitting is manageable with early stopping, dropout/L2, minority-aware batching, and held-out thresholding. For deployment: use PR-centric selection, calibrate and set thresholds to institutional costs, report Recall@Top-k%, retrain each term, and monitor for drift.
The research highlights that different machine learning models provide varying capabilities in detecting academic cheating, with neural networks demonstrating a notable advantage due to their ability to model complex and temporal data patterns effectively. Specifically, CNN and RNN were identified as particularly effective due to their structural advantages in processing sequential and spatial-temporal data. Despite their superior accuracy, the challenge of class imbalance significantly hampers the efficacy of detecting actual cheating cases. Traditional classifiers and gradient boosting methods, although simpler and computationally efficient, showed limitations, especially when confronted with sparse cheating instances, underscoring the necessity of targeted data preprocessing techniques to mitigate such challenges. The practical implications of these findings suggest that combining deep learning models with data balancing methods could yield robust predictive performance, thus significantly enhancing cheating detection mechanisms.
Limitations and Governance. Our labels reflect suspected cases, not adjudicated outcomes; label noise may depress calibration measures. The cohort is single-institution (SCS-DTU, 2021–2024), so external validity requires cross-site evaluation. Although our models are rankers, they should not be used as sole evidence for sanctions. We recommend a human-in-the-loop workflow with clear appeal channels, bias monitoring across demographics/majors, and periodic drift audits.
Ethics and Risk Management. Model outputs are risk indicators to inform human review, not automated sanctioning. We implement (1) bias checks across majors/demographics; (2) appeal procedures for flagged cases; (3) documentation of thresholds and calibration; and (4) periodic drift monitoring and re-training.
5. Conclusions
This study examined cheating detection from exam‐score trajectories under strong class imbalance. We compared classical tabular learners (LR/DT/RF/SVM/MLP), gradient-boosting models (XGBoost, CatBoost, LightGBM), and sequence-aware deep architectures (GRU-RNN, 1D-CNN). A key methodological contribution is an evaluation protocol aligned with the operational reality of rare events: probability calibration, threshold selection on a validation split (optimizing minority-class F1 or an explicit cost function), and Recall@Top-k% for staffing-constrained human triage, reported alongside PR-AUC/ROC-AUC. We further clarified the role of regression-style diagnostics by using Brier = MSE of calibrated probabilities and relegating MAE/RMSE/R2 to probability-quality analysis rather than model ranking.
Empirically, sequence-aware models consistently delivered the strongest PR-centric performance: they capture order-dependent patterns and local temporal motifs that tabular models cannot, yielding higher minority-class F1 and Recall@Top-k% at validated operating points while maintaining competitive ROC-AUC. Boosting methods provided robust baselines but required explicit imbalance remedies (class weighting, focal loss or sample rebalancing) and threshold tuning to avoid the “zero-recall at 0.5” failure mode. Taken together, these findings support a practical recipe for deployment: (1) choose models that learn temporal structure; (2) calibrate probabilities; (3) pick and document operating thresholds on validation; and (4) monitor Recall/Precision@Top-k% as headline KPIs for review workflows.
This work has limitations. Labels reflect suspected rather than adjudicated cases and may contain noise; the cohort is single-institution, which limits external validity; and temporal distribution shift across semesters can degrade performance if models are not refreshed. Accordingly, we recommend human-in-the-loop use with clear appeal channels, bias audits across subpopulations, and scheduled re-training with drift checks.
Future research should expand along four axes: (1) Generalization—multi-campus evaluations and domain adaptation; (2) Modeling—hybrid sequence models (dilated CNNs, GRU-CNN stacks, Transformers) and self-supervised or contrastive pretraining on unlabeled logs; (3) Learning with limited or noisy labels—active learning, semi-supervised and weak-supervision pipelines; (4) Governance—fairness audits, privacy-preserving training, calibration under shift, and uncertainty-aware triage. By pairing sequence-aware modeling with PR-centric evaluation and rigorous governance, institutions can build cheating-risk detectors that are not only accurate on paper but also actionable, auditable, and responsible in practice.
Note
“Brier = mean squared error of calibrated predicted probabilities on the test set (lower is better).”

