We present a complete theoretical characterization of Latent Posterior Factors (LPF), a principled framework for aggregating multiple heterogeneous evidence items in probabilistic prediction tasks. Multi-evidence reasoning arises pervasively in high-stakes domains including healthcare diagnosis, financial risk assessment, legal case analysis, and regulatory compliance, yet existing approaches either lack formal guarantees or fail to handle multi-evidence scenarios architecturally. LPF encodes each evidence item into a Gaussian latent posterior via a variational autoencoder, converting posteriors to soft factors through Monte Carlo marginalization, and aggregating factors via exact Sum-Product Network inference (LPF-SPN) or a learned neural aggregator (LPF-Learned). We prove seven formal guarantees spanning the key desiderata for trustworthy AI: Calibration Preservation (ECE <= epsilon + C/sqrt(K_eff)); Monte Carlo Error decaying as O(1/sqrt(M)); a non-vacuous PAC-Bayes bound with train-test gap of 0.0085 at N=4200; operation within 1.12x of the information-theoretic lower bound; graceful degradation as O(epsilondeltasqrt(K)) under corruption, maintaining 88% performance with half of evidence adversarially replaced; O(1/sqrt(K)) calibration decay with R^2=0.849; and exact epistemic-aleatoric uncertainty decomposition with error below 0.002%. All theorems are empirically validated on controlled datasets spanning up to 4,200 training examples. Our theoretical framework establishes LPF as a foundation for trustworthy multi-evidence AI in safety-critical applications.
Figures & tables
Figure 1: Dependency graph of LPF theoretical results. Assumptions (top) support lemmas and intermediate results, which enable the seven main theorems. Arrows indicate logical dependence. Note that different theorems use different subsets of assumptions: Theorem 3.3 (Generalization) is data-dependent and does not directly rely on Assumptions A1–A6, while Theorems 3.6 and 3.7 build on the results of Theorems 3.1 , 3.2 , and 3.4 rather than their assumptions alone.
Theorem
Key Implementation Details
Verification Experiment
Dataset
Key Metric
Code Variable
T1: Calibration
Does NOT use σmax ; only A1, A3, A4
10-bin calibration
Synthetic ( N=700 )
ECE
epsilon , delta_theoretical
T2: MC Error
Uses A2 for bounded decoder inputs
M -ablation study
20 posteriors
Max error
M , errors
T3: Generalization
Uses deff , NOT σmax
Train/test split
Dedicated ( N=4200 )
Gap vs bound
vc_dim , empirical_gap
T4: Info-Theoretic
Uniform weighting
MI computation
Synthetic ( N=100 )
ECE vs bound
I_E_Y , noise
T5: Robustness
Uses A1, A6
Corruption injection
Synthetic ( N=100 )
L1 distance
corruption_levels , l1_distances
T6: Sample Compl.
K∈{1,…,20} for scaling
K -ablation
Synthetic ( N=100 )
ECE vs K
evidence_counts , lpf_ece
Table 1: Mapping from theoretical guarantees to implementation and empirical verification. All experiments use K≤5 evidence items for main results (extended to K=20 for Theorem 3.6 scaling studies), except Theorem 3.3 which uses a dedicated dataset with N=4200 training examples to achieve non-vacuous generalization bounds.
Figure 2: Calibration verification (Theorem 1). Left: ECE for individual evidence ( 0.140 ), LPF-SPN ( 0.185 ), and LPF-Learned ( 0.058 ), with Hoeffding ( 0.772 ) and Bernstein ( 0.459 ) tight bounds annotated. Centre and right: reliability diagrams for LPF-SPN and LPF-Learned showing confidence vs. accuracy against the perfect-calibration diagonal.
M
Mean Error
Std Error
95th Percentile
Theoretical Bound
4
0.019±0.044
0.044
0.080
0.774
8
0.016±0.030
0.030
0.069
0.547
16
0.013±0.018
0.018
0.053
0.387
32
0.010±0.012
0.012
0.037
0.274
64
0.008±0.009
0.009
0.025
0.193
Table 2: Monte Carlo error bounds: empirical results vs. theoretical guarantees (Theorem 2).
Figure 3: Monte Carlo error bounds (Theorem 2). Left: log-log plot of mean error, 95th-percentile error, and theoretical bound vs. M∈{4,8,16,32,64} ; all empirical curves remain well below the bound. Right: normalised error scaling confirms the empirical rate closely tracks O(1/M) theory.
N
Train Loss
Test Loss
Gap
Bound
2002
0.0407
0.0496
0.0089
0.278
3003
0.0393
0.0455
0.0062
0.253
4200
0.0379
0.0463
0.0085
0.228
Table 3: Generalization bound verification across training sizes (Theorem 3).
Figure 4: Generalization bound verification (Theorem 3). Top-left: train and test loss learning curves with confidence intervals across N∈{2002,3003,4200} . Top-right: empirical gap (near zero) vs. VC bound (loose) and data-dependent PAC-Bayes bound (tight, 0.228 at N=4200 ). Bottom-left: bound-to-gap ratio on a log scale. Bottom-right: test loss vs. N with effective dimension deff=1335 marked.
Metric
Value
Interpretation
Hˉ(Y∣E) (uniform)
0.158 bits
Reported value
H(Y)
1.399 bits
Maximum possible
Reduction
88.7%
Evidence is highly informative
Evidence noise
0.317 bits
Moderate conflicts exist
Table 4: Theorem 4 approximation quality.
Figure 5: Information-theoretic lower bound (Theorem 4). Top-left: decomposition of total uncertainty H(Y)=1.399 bits into evidence information I(E;Y)=1.399 and residual H(Y∣E)≈0 . Top-right: ECE comparison — theoretical lower bound ( 0.158 ), achievable bound including MC term ( 0.317 ), and LPF-SPN empirical ECE ( 0.178 ). Bottom-left: evidence quality distribution (mean ≈1.0 ). Bottom-right: scatter of calibration error vs. evidence conflict (KL divergence), with trend y=0.248x+0.137 .
Figure 6: Robustness to evidence corruption (Theorem 5). Left: empirical L1 distance ∥pclean−pcorrupted∥ (blue) remains near zero while the theoretical O(ϵK) bound (red dashed) grows linearly; the safe region is shaded. Right: bound-to-empirical ratio (up to 6×107 at ϵ=0.1 ), confirming the bound is highly conservative in practice.
Figure 7: Sample complexity scaling (Theorem 6). Top-left: LPF-Learned ECE (blue) and baseline uniform ECE (green) both lie far below the theoretical O(1/K) bound (red dashed) for K∈{1,…,20} . Top-right: bound-to-empirical ECE ratio. Bottom-left: O(1/K) fit ( 0.25/K+0.12 , R2=0.849 ) with empirical ECE plateauing at K≈7 . Bottom-right: LPF vs. uniform baseline at K∈{1,2,3,5} ; baseline available only for K≥5 .