The Illusion of AI Expertise Under Uncertainty: A Probabilistic Paradigm to Interpret Results Under Ground Truth Variation
Authors: Aparna Elangovan, Lei Xu, Mahsa Elyasi, Ismail Akdulum, Mehmet Aksakal, Enes Gurun, Brian Hur, Saab Mansour, Ravid Shwartz Ziv, Karin Verspoor, Dan Roth
Organizations: Independent Researcher · Amazon · Gazi University · University of Colorado · Samsun University · The University Of Melbourne · New York University · RMIT University · University of Pennsylvania
Benchmarking the capabilities of AI systems, including Large Language Models (LLMs) and Vision Models, typically ignores the impact of uncertainty in the underlying ground truth answers from experts. This ambiguity is not just limited to human preferences, but is also consequential even in safety critical domains such as medicine where uncertainty is pervasive. In this paper, we introduce a probabilistic paradigm to theoretically explain how high certainty in ground truth answers is almost always necessary for even an expert to achieve high scores. In contrast, in datasets with high variation in ground truth answers, the expert performance drops substantially resulting in almost no difference between a random labeller and an expert. This characteristic also manifests when comparing models, where uncertainty obfuscates differences between poor and high performing models. Therefore, ignoring uncertainty in ground truth evaluation data can result in the misleading conclusion that a non-expert has similar performance to that of an expert. Given ground truth answer variability and commonly used F1 and accuracy metrics, we bring forth a simple and explainable probabilistic coin toss model to estimate expected accuracyand expected F1 and compare the estimated scores with those that experts (human or model) achieve across multiple domains -- including radiology and models such as GPT-5 and Gemini-3. The results lead to the recommendation that stratification by the probability of the ground truth answer becomes critical when expert performance is relatively low. Under stratified evaluation, performance comparison becomes more reliable in high certainty bins, mitigating the effect of the key confounding factor --- uncertainty.
Figures & tables
N=100
N=10
pd
rc
P(Score>=rc)
rc
P(Accurcy>=rc)
pd=pr=0.25
Random labeller baseline for 4 classes
0.25
50
6.6e-8
5
0.08
0.25
80
1.2e-30
8
0.00
pd=pr=0.5
Random labeller baseline for 2 classes
0.5
50
0.54
5
0.62
Table 1: Examples of how given probability pd of following the majority label, the corresponding probability of scores over 50% or 80% accuracy changes out of N 100 vs. 10 examples.
pd (Expected Recall)
m
Expected Precision
Expected F1
pd=pr=0.5
Random labeller baseline for 2 classes
0.5
0.01
0.01
0.02
0.5
0.10
0.10
0.17
0.5
0.50
0.50
0.50
0.6
0.01
0.01
0.03
0.6
0.10
0.14
0.23
Table 2: Examples showing how given probability pd following the majority label, the positive class ratio m , the expected precision and F1 increases with increasing m .
Figure 1: Simulated impact of agreement for the majority label (X axis). Y axis is the performance of a simulated human (H) who follows the majority label distribution ( pd ) vs. the random labeler (R) performance across different positive label ratios m∈0.1,0.3,0.5 . At low certainty, the difference between human labeler and random labeler is relatively low compared to high certainty. Recall and Accuracy are unaffected by the positive label ratio. The vertical bars indicate standard deviation.
Figure 2: Comparison of theoretical vs. empirical performance of Models ( M ) and human ( H ) radiologists, compared against ground truth (majority label across 5 human radiologists). H is the average performance of additional 3 radiologists. Dashed lines are theoretical expected performance at given pd and m according to equations 2 and 14 . Theoretical performance closely matches human expert performance ( pd→1.0 is approximated as pd=0.97 ). Humans achieve a relatively high F1-score ( >0.8 ) as pd→1.0 , whereas pd→0.6 the peak performance drops close to a random labeller baseline. At low positive class ratio ( m<0.01 ), even at pd→1.0 (see Pneumonia and Consolidation), F1 is <0.35 . Δ=H−M is the vertical distance illustrated by the length of the vertical line connecting the human and the model performance. See Appendix Tables 8 and 11 for detailed results with sample size. Similar results for Gemini-2.5 & Gemini-3.1 are in Appendix Figure 10 and Tables 10 9 .
∣GT∣=5,∣L∣ =2
F1
Accuracy
Pathology
pd
S
m
AI
H
E
ΔH−AI
ΔE−H
ΔE−M
AI
H
E
ΔH−AI
ΔE−H
ΔE−M
Gemini-3.0 Preview
Cardio- megaly
0.60
142
0.507
0.67
0.64 ± 0.10
0.60
-0.03
-0.04
-0.07
0.63
0.66 ± 0.03
0.60
0.03
-0.06
-0.03
0.80
196
0.372
0.69
0.71 ± 0.04
0.75
0.02
0.04
0.06
0.72
0.79 ± 0.03
0.80
0.07
0.01
0.08
1.00
162
0.037
0.35
0.70 ± 0.12
0.71
0.35
0.01
0.36
0.91
0.98 ± 0.02
0.97
0.07
-0.01
0.06
All
500
0.302
0.66
0.68 ± 0.07
0.68
0.02
0.00
0.02
0.75
0.81 ± 0.03
0.80
0.06
-0.01
0.05
Table 3: Examples of models ( M ) Gemini-3.0 and GPT-5.1 performance variation when stratified by pd compared to average human ( H ) performance of 3 radiologists and the corresponding standard deviation, sample size S . For Lung Opacity, the overall ( pd=All ) performance is over 80%, and even at stratified pd , the models achieves reasonable F1 and accuracy scores compared to humans. On the other hand, for Cardiomegaly and Support Devices the overall ( pd=All ) performance falls below 80% and and Gemini even achieves a Δ=H−M of just 2 points for Cardiomegaly, however as pd→1 , the Δ is over 20 points. Expected ( E ) performance given ⟨pd,m⟩ per equations 2 and 14 (note: pd→1.0 is approximated as pd=0.97 ). Scores (Model or Human) over 0.80 are highlighted and they all occur at pd∈0.8,1.0 but never at 0.60 . Full results in Appendix Tables 8 , and 11 .
Figure 3: (Left and Mid) Distribution of difference in Human vs. Model Δ=H−M : Each point Δi,o,pd=Ho,pd−Mi,o,pd where i∈{Gemini-3.1,Gemini-3.0,Gemini-2.5,GPT-5.1} , pd∈{0.6,0.8,1.0,All} , o∈Pathologies . Ho,pd is the human score (F1 or accuracy) against the pathology o at pd . Similarly, Mi,o,pd is the performance of the model i against the pathology o at pd . Average Δ is generally higher as pd→1.0 compared to pd→0.6 where Δ F1 is close to zero. (Right) Distribution of absolute model accuracy: All 4 model accuracies (agnostic of m ) across different pd . Empirical vs. (expected) mean annotated in text. Note: Significance Welch t-test for means ( H0 Null hypothesis: there is no difference between the means) and Levene’s test (trimmed, proportion to cut=0.1) for variance ( H0 : there is no difference in variance): **** ( p<=0.0001 ), *** ( p<=0.001 ), ** ( p<=0.01 ), * ( p<=0.05 ), ns (not significant). Blue : Indicates significance of variance.
∣GT∣∈⟨2,3⟩ , ∣L∣ =2
F1
Accuracy (A)
pd
S
m
AI
H′
E
ΔH′−AI
ΔE−H′
ΔE−M
AI
H′
E
ΔH′−AI
ΔE−H′
ΔE−M
0.67
15834
0.229
0.39
0.48 ± 0.0045
0.48
0.09
0.00
0.09
0.42
0.67 ± 0.0045
0.67
0.25
0.00
0.25
1.00
185245
0.029
0.11
1.00 ± 0.0000
1.00
0.89
-0.00
0.89
0.58
1.00 ± 0.0000
1.00
0.42
0.00
0.42
ALL
201079
0.045
0.15
0.75 ± 0.0000
0.75
0.60
-0.00
0.60
0.56
0.97 ± 0.0000
0.97
0.41
0.00
0.41
Table 4: Stratified performance on existing public results [ 5 ] on Mammogram. S is the sample size, H′ is simulated human performance based on GT distribution and hence matches expected performance. ΔE−M narrows as certainty lowers as shown when pd→0.67 compared to pd→1.0 .
Figure 4: Model-2-Model comparison: Left and Mid: Strong (SM) vs. weak model (WM) : The mean of ΔSM−WM is less than zero at pd→0.6 , indicating that the weaker model is better than the stronger model, demonstrating the unreliability of comparison when ground truth is not certain. At pd→1.0 , ΔSM−WM of stronger model is higher than zero. Right: Comparison of 2 SMs: pd→1.0 shows that there is only marginal difference between two models but pd→0.6 is not a clear indicator as it is not any different to weak vs. strong model comparison. Paired Student’s T-test for statistical significance ( H0 : No difference in performance between a pair of models). Note: Improvements may be marginal between models, hence not statistically significant.
Appendix figures & tables13 assets
Supplementary material from the paper’s appendix.
Appendix
Figure 5: Proportion of samples (CheXpert test set with 500 samples) by agreement level as annotated by 5 radiologists independently. The number of samples with 100% agreement ( pd=1.0 ) on positive findings drops dramatically to less than 0.01 for 4 of 5 the pathologies. The exception is lung opacity, where the proportion of samples approximately remains the same regardless of agreement levels. Agreement level for all pathologies are in Appendix 6 .
Figure 6: Chest x ray agreement for binary classification across 13 labels. The agreement is across 5 radiologists. The x-axis shows the percentage agreement for the majority label ( yes/ no) and the corresponding y-axis indicate the total number of records with that agreement. For instance, for Atelectasis, 100% agreement ( 5 out of 5 radiologists) is observed for 32 records for a positive finding (yes), while 100% agreement for negative finding for 200 records.
Figure 7: MNLI Matched, MNLI Mismatched and SNLI performance (Accuracy) distribution by pd on 4 models - Mistral, Mistral Large, Sonnet and Llama 3.1. These are existing results from Elangovan et. al. [ 8 ] . The variation in scores is low when pd→0.6 compared to pd→1.0 , as demonstrated by the statistical significance using Levene’s test for variance. This shows the lack of differentiating capability between strong vs. weak models when pd→0.6 . The mean performance at pd→1.0 is also higher (statistically significant using Welch’s t-test) compared to pd→0.6 Detailed tabular results are available in Table 5 .
Dataset
pd
S
E(A)
H′
ΔEH′
Model
A
ΔEM
MNLI M ⟨ GT = 5, L=3 ⟩
0.60
1599
0.60
0.61 ± 0.02
-0.01
Llama3.1
0.55
0.05
Mistral
0.50
0.10
Mistral-L
0.59
0.01
Sonnet
0.61
-0.01
0.80
2457
0.80
0.81 ± 0.01
-0.01
Llama3.1
0.68
0.12
Mistral
0.63
0.17
Appendix
Table 5: Expected accuracy ( A ) vs empirical accuracy. Here M is the existing result on Models by Elangovan et al [ 8 ] . GT is the number of ground truth annotations and L is the number of labels. H ’ is the simulated performance of another Human. ΔEM=E(A)−M , ΔEH′=E(A)−H′ . † Only 7 samples (in pd→0.75 ) have 4 annotations in the SNLI dataset, the rest have 5.
Figure 8: Difference in expected vs. actual model performance of G-Eval original results [ 22 ] for 4 criteria - Coherence, relevance, consistency and fluency. Ground truth is the majority label. Model performance is the most frequent label from the model. ΔEM is close to zero pd→0.33 , the gap between expected and the model’s substantially lower performance becomes apparent as pd→1.0 . Detailed tabular performance in Appendix Table 6 . Examples of differences between G-Eval and human experts are shown in Appendix Table 7 . Paired (expected vs. actual) T-test for significance.
Criteria
pd
S
E(A)
Score
ΔEM
Coherence ⟨GT =3, L=5 ⟩
0.33
341
0.33
0.185
0.145
0.67
1059
0.67
0.291
0.379
1.00
200
1.00
0.165
0.835
Consistency ⟨GT =3, L=5 ⟩
0.33
19
0.33
0.368
-0.038
0.67
241
0.67
0.241
0.429
1.00
1340
1.00
0.305
0.695
Appendix
Table 6: Accuracy of GEval: Model performance is the most frequent label from the model. Human ground truth is also the majority label. ΔEM narrows down the most as pd→0.33
Table 7: Examples of differences at GEval vs Human experts for coherence.
Figure 9: Model-2-Model accuracy comparison: Paired Student’s T-test to measure the statistical significance where the null hypothesis ( Ho )is there is no difference between the pair of models. Significance: * ( p<=0.05 ), ns (not significant). Left and Mid : Comparison of strong and weak model: Right : Comparison of 2 strong models.
F1
Accuracy
Pathology
pd
S
m
AI
H
E
ΔH−AI
ΔE−H
ΔE−M
AI
H
E
ΔH−AI
ΔE−H
ΔE−M
Atelectasis
0.60
130
0.500
0.57
0.66 ± 0.03
0.60
0.09
-0.06
0.03
0.53
0.60 ± 0.02
0.60
0.07
-0.00
0.07
0.80
140
0.400
0.63
0.70 ± 0.03
0.76
0.07
0.06
0.13
0.67
0.73 ± 0.03
0.80
0.06
0.07
0.13
1.00
230
0.139
0.67
0.74 ± 0.03
0.90
0.07
0.16
0.23
0.88
0.91 ± 0.02
0.97
0.03
0.06
0.09
All
500
0.306
0.62
0.69 ± 0.02
0.73
0.07
0.04
0.11
0.73
0.78 ± 0.02
0.83
0.05
0.05
0.10
Cardio- megaly
0.60
142
0.507
0.67
0.64 ± 0.10
0.60
-0.03
-0.04
-0.07
0.63
0.66 ± 0.03
0.60
0.03
-0.06
-0.03
Appendix
Table 8: Performance of Gemini Pro 3.0 Preview ( GP ) model vs. human ( H ) radiologists, compared against ground truth (majority label across 5 human radiologists). The data is stratified by pd probability of observed agreement among the 5 radiologists for the majority label. Size (S) is the number of samples corresponding to that pd , positive class ratio is m . H is the average performance of additional 3 radiologists and the corresponding standard deviation. Humans achieve a peak F1-score of over 0.98 @ pd=1.0 , while the best performance at pd=0.6 is 0.72 across all pathologies. Similar trend is seen for recall, precision and accuracy. At low positive class ratio ( m<0.01 ), even at pd=1.0 (see Pneumonia and Consolidation), F1 is <0.35 . Δ=H−GP is the difference between the average of human performance and the model. All the scores (Model or Human) over 0.80 are highlighted and they all occur at pd∈0.8,1.0 but never at 0.60 given the sample size ( S ).
F1
Accuracy
Pathology
pd
S
m
AI
H
E
ΔH−AI
ΔE−H
ΔE−M
AI
H
E
ΔH−AI
ΔE−H
ΔE−M
Atelectasis
0.60
130
0.500
0.61
0.66 ± 0.03
0.60
0.05
-0.06
-0.01
0.52
0.60 ± 0.02
0.60
0.08
-0.00
0.08
0.80
140
0.400
0.63
0.70 ± 0.03
0.76
0.07
0.06
0.13
0.63
0.73 ± 0.03
0.80
0.10
0.07
0.17
1.00
230
0.139
0.58
0.74 ± 0.03
0.90
0.16
0.16
0.32
0.82
0.91 ± 0.02
0.97
0.09
0.06
0.15
All
500
0.306
0.61
0.69 ± 0.02
0.73
0.08
0.04
0.12
0.69
0.78 ± 0.02
0.83
0.09
0.05
0.14
Cardio- megaly
0.60
142
0.507
0.43
0.64 ± 0.10
0.60
0.21
-0.04
0.17
0.62
0.66 ± 0.03
0.60
0.04
-0.06
-0.02
Appendix
Table 9: Performance on CheXpert using Gemini-3.1 pro
Figure 10: Performance of Gemini-2.5, Gemini-3.1 ( M ) vs. human ( H ) radiologists, compared against ground truth (majority label across 5 human radiologists). H is the average performance of additional 3 radiologists who are different from the 5 radiologists used to determine the ground truth. The data is stratified by pd probability of observed agreement among the 5 radiologists for the majority label, corresponding positive class ratio is m . Dashed lines are theoretical expected performance at given pd and m according to equations 2 and 14 , ( pd→1.0 is approximated as pd=0.97 ). Humans are able to achieve a relatively high F1-score ( >0.8 ) as pd→1.0 , whereas pd→0.6 the peak performance drops close to a random labeller baseline. At low positive class ratio ( m<0.01 ), even at pd→1.0 (see Pneumonia and Consolidation), F1 is <0.35 . Δ=H−M is the vertical distance illustrated by the length of the vertical line connecting the human and the model performance.
F1
Accuracy
Pathology
pd
S
m
AI
H′
E
ΔH−AI
ΔE−H
ΔE−M
AI
H
E
ΔH−AI
ΔE−H
ΔE−M
Atelectasis
0.60
130
0.500
0.65
0.66 ± 0.03
0.60
0.01
-0.06
-0.05
0.52
0.60 ± 0.02
0.60
0.08
-0.00
0.08
0.80
140
0.400
0.62
0.70 ± 0.03
0.76
0.08
0.06
0.14
0.53
0.73 ± 0.03
0.80
0.20
0.07
0.27
1.00
230
0.139
0.40
0.74 ± 0.03
0.90
0.34
0.16
0.50
0.60
0.91 ± 0.02
0.97
0.31
0.06
0.37
All
500
0.306
0.56
0.69 ± 0.02
0.73
0.13
0.04
0.17
0.56
0.78 ± 0.02
0.83
0.22
0.05
0.27
Cardio- megaly
0.60
142
0.507
0.69
0.64 ± 0.10
0.60
-0.05
-0.04
-0.09
0.56
0.66 ± 0.03
0.60
0.10
-0.06
0.04
Appendix
Table 10: Performance on CheXpert using Gemini-2.5 pro
An LLM judge evaluates outputs at scale. Experts should label only where it is least sure. Its natural escalation signal conflates two uncertainties: aleatoric, real disagreement in the expert pool, which labels cannot reduce, and epistemic, the judge's ignorance, which labels do reduce. A small Bayesian model separates them: a regression on labels already collected learns how far to trust a black-box judge's prediction. Both components follow as simple formulas, with no sampling or further judge calls. The components isolate on a real LLM judge against exactly known truth, and stated confidence is no guide to its actual error. On real human disagreement (ChaosNLI) the epistemic ranking removes 83% more error than total uncertainty for the same expert labels, though simply escalating the least-labelled items does as well there. We demonstrate we can estimate where a judge is ignorant rather than where experts genuinely disagree, and propose using this to direct expert labelling. Code and data are available at https://github.com/composo-ai/judge-uncertainty-decomposition.
Do expert personas improve language model performance? The Wharton Generative AI Lab reports that they do not, broadcasting to millions via social media the recommendation that practitioners abandon a technique recommended by Anthropic, Google, and OpenAI. We demonstrate that this null finding was structurally predictable. Five core mechanisms precluded detection before data collection began: baseline contamination elevating the starting point to near-ceiling, system prompt hierarchy subordinating experimental manipulation, impossible expert specifications collapsing to generic competence, format constraints suppressing reasoning processes, and provider exclusion limiting generalizability. Controlled trials correcting these limitations reveal what the original design obscured. To test this, we selected the GPQA Diamond hardest questions to prevent baseline pattern matching, forcing reliance on genuine expert reasoning. On items with valid key answers, expert personas achieve ceiling accuracy. They eliminated all baseline errors through confidence amplification. Furthermore, forensic examination of model divergence identified that half of the hardest GPQA items contain chemically or logically indefensible answers. The model's CoT revealed reasoning away from impossible answers, yielding penalization for accurate chemistry. These findings recontextualize the original null results. Methodologically sound persona research faces measurement constraints imposed by benchmark validity limitations. Answering the persona question requires evaluation infrastructure the field does not yet possess.
Drake Mullens, Stella Shen
Tarleton State University · University of Texas at Tyler
Reliable Large Language Models (LLMs) should abstain when confidence is insufficient. However, prior studies often treat refusal as a generic "I don't know'', failing to distinguish input-level ambiguity (data uncertainty) from capability limitations (model uncertainty). This lack of distinction limits downstream action decisions like requesting clarification or invoking external tools. In this work, we introduce UA-Bench, a benchmark of over 3,500 questions drawn from six datasets spanning knowledge-intensive and reasoning-intensive tasks, designed to evaluate explicit uncertainty attribution. An evaluation of 18 frontier LLMs shows that even state-of-the-art models struggle to reliably discriminate between data uncertainty and model uncertainty, and that high answer accuracy does not necessarily imply strong uncertainty attribution ability. To narrow this gap, we propose a lightweight data synthesis and reinforcement learning strategy. Experiments on both Qwen3-4B-Instruct-2507 and Qwen3-8B in thinking mode show that the proposed method improves uncertainty attribution while preserving answer accuracy. Our code and data are publicly available now.
Jingyi Ren, Ante Wang, Yunghwei Lai +5
Dept. of Comp. Sci. & Tech., Institute for AI, Tsinghua University, Beijing, China · Institute for AI Industry Research (AIR), Tsinghua University, Beijing, China