Beyond Performance Metrics: Uncertainty Mapping of Label Ambiguity in Fazekas Score Prediction
Authors: Susanne Schmid, Johanna Ospel, Richard Frayne, Roberto Souza
Organizations: Schulich School of Engineering, Department of Electrical and Software Engineering, University of Calgary, Calgary, T2N 1N4, AB, Canada · Hotchkiss Brain Institute, University of Calgary, Calgary, T2N 1N4, AB, Canada · Department of Radiology and Clinical Neurosciences, University of Calgary, Calgary, T2N 1N4, AB, Canada · Calgary Image Processing and Analysis Centre, Foothills Medical Centre, Calgary, T2N 2T9, AB, Canada · Seaman Family MR Research Centre, Foothills Medical Centre, Calgary, T2N 2T9, AB, Canada
Reference labels used to train medical image classification models are not always as certain as they may appear, and this uncertainty has implications on performance metrics. In this study, we propose a framework to analyze model performance for periventricular Fazekas score prediction that goes beyond conventional metrics. The Fazekas score is an ordinal visual rating scale used to assess the severity of white matter hyperintensities and is known to be affected by inter-rater variability. While the best Fazekas score prediction model achieved a Matthews correlation coefficient (MCC) of 0.70, performance varied across data splits and loss functions, making interpretation of model capabilities difficult. Rather than interpreting epistemic uncertainty of a model's prediction as an isolated scalar value, our approach of uncertainty mapping relates uncertainty to its position within the learned feature representation. This highlights regions of class-boundary transitions where cases appear more ambiguous and misclassifications are more likely. It also identifies potential label disagreement, including low-uncertainty misclassified cases that expert review found to be inconsistent with the original reference Fazekas score. Therefore, uncertainty mapping allows model behaviour to be examined in relation to class separation and potential model-label disagreement. Loss function choice also influenced the uncertainty profile, with some models showing clearer class separation and more localized uncertainty in ambiguous regions than others. These findings suggest that uncertainty mapping for Fazekas score predictions can support model interpretation and targeted dataset review when reference labels are affected by ambiguity/ inter-rater variability.
Figures & tables
Figure 1: Distribution of the periventricular Fazekas score (PVWMH) by age. With increasing age, the WMHs become more prevalent. By the age of 80 years, all individuals have some extent of WMHs.
Figure 2: Two-channel 3D volume input to the convolutional neural network architecture. Multiple dropout layers were included to enable uncertainty estimation using Monte Carlo dropout.
Figure 3: Uncertainty Mapping: Feature Visualization combined with epistemic uncertainty estimation obtained through Monte Carlo dropout.
Figure 4: Visual comparison of Layerwise-Relevance Propagation (LRP) patterns for four representative subjects with periventricular Fazekas scores of 0, 1, 2, and 3. LRP heatmaps are shown for models trained with different loss functions from the fold with the best overall performance. While performance differed substantially between models, these differences were not clearly reflected in the heatmaps. Red regions indicate areas supporting the model prediction, whereas blue regions indicate areas contributing against the prediction or toward an alternative class. The subject with Fazekas score 1 was incorrectly predicted by the GCE-trained model.
Loss function
Accuracy
MCC
Cohen’s κ
F1: Class 0
F1: Class 1
F1: Class 2
F1: Class 3
CCE
0.73 ± 0.05
0.46 ± 0.09
0.45 ± 0.09
0.54 ± 0.11
0.80 ± 0.05
0.65 ± 0.13
0.33 ± 0.44
GCE
0.68 ± 0.06
0.39 ± 0.07
0.37 ± 0.08
0.44 ± 0.19
0.76 ± 0.07
0.61 ± 0.14
0.24 ± 0.43
OLL
0.73 ± 0.04
0.47 ± 0.10
0.46 ± 0.11
0.53 ± 0.16
0.80 ± 0.04
0.66 ± 0.12
0.31 ± 0.46
OLL clip
0.74 ± 0.05
0.49 ± 0.10
0.48 ± 0.10
0.57 ± 0.10
0.81 ± 0.04
0.66 ± 0.14
0.40 ± 0.46
Table 1: Nested cross-validation results: Average of all 20 runs ± the standard deviation.
Figure 5: Distribution of epistemic uncertainty per Fazekas score class for the baseline model trained with CCE, achieving an MCC of 0.59. Misclassified cases are highlighted in red, while correct predictions are gray.
Figure 6: t -SNE dimension reduction of the penultimate layer for CCE baseline model (MCC=0.59). Each point represents one subject from the test set. The fill colour indicated the reference label, while the edge colour indicated the model predictions. Hence, points with different fill and edge colours correspond to misclassified samples, while homogeneous points are correctly predicted.
Figure 7: Uncertainty results combined with the t -SNE analysis. The shapes represent the reference Fazekas score, while the color scale represents epistemic uncertainty. In a), epistemic uncertainty is shown for all samples, while in b), the same results are visualized only for incorrect predictions; correct predictions are shown in grey. Dashed ellipses indicate Fazekas score region and were added to guide visual interpretation.
Figure 8: Epistemic uncertainty estimation visualized combined with the corresponding model feature space representation for the four different loss functions analyzed (a) CCE, (b) GCE, (c) OLL, and (d) OLL clip. The color scale represents epistemic uncertainty, while the shape corresponds to the reference Fazekas score.
Figure 9: Agreement between the predictions of the CCE and OLL clip model is shown for both models uncertainty mapping. Coloured points indicate the epistemic uncertainty of test samples where the two models disagreed. Grey data points reflect subjects both models agreed on. The shape corresponds to the reference label.
Selection Group
Cases Reviewed
Revised
Observed revisions
Inconsistency
4
4
(A) 0 → 2; (B) 0 → 2; (C) 0 → 1; (D)1 → 3
Ambiguous boundary
4
2
(E) 1-; (H) 1 → 2; (G) 1 → 3-
Confident
4
1
(I) 0 → 1; (J) 1+
Table 2: Target re-evaluation of the selected 12 cases summarized by group. Case letters correspond to the samples shown in Figure 10 . Label revisions are reported as original reference label → reassessed label. The symbols + and − indicate directional uncertainty in the reassessment score, where + indicates that the score is at the higher end and − indicates that the score is at the lower end.
Figure 10: Identification of cases selected for targeted re-evaluation using the OLL clip uncertainty map. Three groups were selected: (1) label-inconsistency cases (A-D); (2) ambiguous boundary cases (E-H), selected from regions of elevated uncertainty near transitions between neighbouring Fazekas scores, and (3) confident agreement cases (I-L), defined as low-uncertainty correctly classified samples located within well-defined Fazekas score regions.
Figure 11: Label-inconsistency cases (misclassified despite low epistemic uncertainty) for the OLL clip model based on Figure 10 .
Human-model alignment is critical for trustworthy AI-assisted decision-making systems. Yet, most work evaluates model predictions against single ground-truth labels, overlooking that humans themselves often disagree on labels, a signal of genuine ambiguity. We investigate whether models struggle on the same instances that humans find difficult. We measure this on two vision datasets (FER+ and CIFAR-10H) where multiple human annotations per image capture human disagreement patterns. We evaluate eight pretrained models across three architectures (ResNet, EfficientNet, MobileNetV3) in two parts: first, whether model uncertainty (softmax confidence, entropy) correlates with human disagreement, and second, whether predictive multiplicity measures (inter-model disagreement, Jensen-Shannon divergence) do. We find that it does not: alignment is weak in both dimensions. At the discrete label level, 50.4% of CIFAR-10H images and 33.5% of FER+ images receive multiple valid classifications from humans, while the models converge on only one. These instances represent a critical failure case where humans perceive ambiguity and would request expert review, yet models decide confidently. At the continuous score level, single-model uncertainty correlates weakly with human disagreement (ρ=0.24−−0.55), and predictive multiplicity provides only modest improvement. Widely-used uncertainty quantification methods do not reliably identify instances humans find ambiguous. Model uncertainty should not be treated as a trustworthy signal by default for decision-making in high-stakes scenarios.
Manya Singh, Arjun Pakrashi
School of Computer Science, University College Dublin, Ireland
Medical segmentation uncertainty is often evaluated with calibration error or pixel-level error-detection scores, but these metrics do not show whether uncertainty improves downstream review decisions. We propose a decision-aware evaluation protocol that separates uncertainty estimation from deferral policy and evaluates their combination under fixed review budgets: which pixels, regions, or images should be accepted, and which should be reviewed. The protocol reports fixed-budget error reduction, risk-coverage behaviour, region-level review proxies, structural failure modes, probability-derived baselines, and sensitivity to review assumptions. We instantiate it across retinal vessel, skin lesion, and polyp segmentation with multiple backbones and uncertainty estimators. Under matched inference budgets and canonical multi-layer MC Dropout, test-time augmentation (TTA) gives a stronger error-detection signal across matched MC--TTA comparisons, including architecture and dropout-placement checks. However, estimator rankings are not universal. For example, probability-derived scores remain competitive, deep ensembles can outperform TTA, learned deferral helps mainly when the base estimator is weak, and rankings can change with review unit, uncertainty target, correction assumptions, and cost model. Structural and second-observer analyses further show that boundary and thin-structure failures are partly tied to annotation ambiguity. These findings support evaluating uncertainty as an estimator-policy-review configuration, not as a standalone score map. Code: https://anonymous.4open.science/r/Uncertainity-in-Medical-Image-Segmentation/README.md
Saket Maganti, Shreyank N Gowda
Independent Researcher · School of Computer Science, University of Nottingham, Nottingham, UK
Many machine learning (ML) applications rely on expert labels, and qualified experts may provide different but plausible interpretations of the same observation. Such variation across expert labels may reflect genuine disagreement or ambiguity rather than annotation error. When individual experts additionally report intervals rather than exact values, the supervision contains two distinct sources of label uncertainty: within-label imprecision and between-expert variation. Existing methods treat these forms separately: multi-expert approaches collapse labels to a consensus, interval-target methods often yield a single prediction, and predictive-uncertainty methods rarely validate their uncertainty estimates against observed expert disagreement. To address this problem, we propose an approach that preserves individual expert intervals, separates within-label imprecision from between-expert variation, and validates the corresponding predictive uncertainty components. First, heterogeneous label vocabularies are harmonized into a common probabilistic label space, separating encoding differences from expert judgement. Second, individual label intervals are retained and modeled with a mixture of Beta distributions trained using a proper Cramér-distance objective, preserving distinct expert-reported labels. Third, we decompose predictive uncertainty into within-component, between-component, and model uncertainty, and evaluate whether these components correspond to within-label uncertainty, between-label uncertainty, and model error, respectively. Because this correspondence is not guaranteed, we introduce decomposition matching, which aligns the predictive components to their intended label-side sources. On sea-ice concentration the model reduces MAE by 31% over hard labels and outperforms aggregation, interval-distribution, and interval-regression baselines.
Samira Alkaee Taleghan, Younghyun Koo, Andrew P. Barrett +1
University of Colorado Denver Denver, Colorado, USA · National Snow and Ice Data Center (NSIDC) CIRES, University of Colorado Boulder Boulder, Colorado, USA