The Bayesian Learning Rule Beyond KL Geometry

cs.LG Sophia Sklaviadis, Wu Lin, Andre Martins, Mario Figueiredo, Thomas Moellenhoff, Emtiyaz Khan · Oct 22, 2026
Local to this browser
What it does
The paper attacks a structural restriction in the Bayesian Learning Rule (BLR): its reliance on KL geometry ties it to fixed-support exponential families, ruling out posteriors with adaptive tails or bounded support. The authors show that...
Why it matters
The authors show that mixed-space Bregman discrepancies—specifically Fenchel-Young discrepancies—generalize the BLR while preserving its efficient additive canonical-coordinate updates. This matters because it yields a principled framework...
Main concern
The paper presents a theoretically sound extension of the BLR beyond KL geometry. The identification of Fenchel-Young discrepancies as the structure that preserves additive updates is novel, and the recovery of the classical BLR as a...
Community signal
0
0 up · 0 down
Sign in to vote with arrows
AI Review AI reviewed
Plain-language introduction

The paper attacks a structural restriction in the Bayesian Learning Rule (BLR): its reliance on KL geometry ties it to fixed-support exponential families, ruling out posteriors with adaptive tails or bounded support. The authors show that mixed-space Bregman discrepancies—specifically Fenchel-Young discrepancies—generalize the BLR while preserving its efficient additive canonical-coordinate updates. This matters because it yields a principled framework for scalable variational optimizers, such as $q$-VON, that use heavy-tailed or bounded-support $q$-Gaussian posteriors and can be implemented by modifying only a few constants in existing code.

Critical review
Verdict
Bottom line

The paper presents a theoretically sound extension of the BLR beyond KL geometry. The identification of Fenchel-Young discrepancies as the structure that preserves additive updates is novel, and the recovery of the classical BLR as a special case is clean. However, the empirical evaluation is explicitly positioned as preliminary; the authors defer large-scale benchmarks and state-of-the-art comparisons to future work, which limits the practical impact of the contribution.

“drawing these connections is beyond the scope of this paper; as is demonstrating state-of-the-art or scale. We show experiments up to CIFAR-100 (App. G), and defer large experiments, for example on quantization-aware training, pruning, and language model training with bounded variance to future work.”
Sklaviadis et al., Sec. 1 · Section 1
“While the CIFAR-100 experiment we show in App. G is not a large-scale benchmark, there is no structural obstacle to scaling further, and applying q-VON to transformer- or LLM-scale training is an empirical engineering question, which we leave to future work.”
Sklaviadis et al., Sec. 4 · Section 4
What holds up

The convex-analytic derivation is rigorous. Theorem 1 shows that the canonical gradient under a mixed-space Bregman divergence retains the same additive form as the BLR, and Propositions 1 and 2 correctly establish the induced maxent family and its dual coordinates. The instantiation for Tsallis negentropy yielding bounded-support and heavy-tailed $q$-Gaussians is a compelling concrete payoff.

“When \Omega=-\mathcal{H} we recover the classical BLR.”
Sklaviadis et al., Sec. 3.4 · Section 3.4
“the generalized canonical gradient satisfies \mathbf{F}(\boldsymbol{\lambda})^{-1}\nabla_{\boldsymbol{\lambda}}[\mathbb{E}_{\pi_{\boldsymbol{\lambda}}}[\ell(\boldsymbol{\theta})]+\mathbb{D}_{\Omega}(\eta_{\boldsymbol{\lambda}_0},\pi_{\boldsymbol{\lambda}})] = \nabla_{\boldsymbol{\mu}}\mathbb{E}_{\pi_{\boldsymbol{\lambda}(\boldsymbol{\mu})}}[\bar{\ell}(\boldsymbol{\theta})]|_{\boldsymbol{\mu}=\boldsymbol{\mu}(\boldsymbol{\lambda})} + \boldsymbol{\lambda}”
Sklaviadis et al., Sec. 3.4 · Theorem 1, Eq. 16
Main concerns

The largest weakness is the limited empirical evidence. The authors state that demonstrating state-of-the-art performance or scale is beyond the scope and that the CIFAR-100 experiment is not a large-scale benchmark, deferring transformer-scale training to future work. This undercuts the claim that the method readily scales to modern deep learning. A second concern is an unresolved uniqueness question: the authors note that whether mixed-space Bregman divergences are the only family that recovers BLR-style additive updates while preserving linear expectations remains open, leaving the theoretical picture incomplete.

“drawing these connections is beyond the scope of this paper; as is demonstrating state-of-the-art or scale.”
Sklaviadis et al., Sec. 1 · Section 1
“An interesting question that we hope to address in the future is whether mixed-space Bregman divergences uniquely recover BLR-style updates when we also wish to preserve the linear expectation of the loss term in Eq. 3.”
Sklaviadis et al., Sec. 2.3 · Section 2.3
Evidence and comparison

The theoretical evidence supports the formal claims: the derivations of the generalized BLR update and the $q$-Gaussian Newton-like optimizer are complete. Comparisons to related work are fair but narrow. The paper correctly distinguishes its contribution from $f$-divergence approaches such as Alquier (2021), noting that only the KLD permits the necessary separation into linear expectations within that class, whereas the proposed Fenchel-Young framework avoids this restriction. Empirically, the method is shown to be competitive with IVON on CIFAR-10 and CIFAR-100, though the improvements are modest and the evaluation metrics standard.

“This separation in terms of a linear expectation of \eta_{\boldsymbol{\lambda}_0} and a generalized negentropy of \pi_{\boldsymbol{\lambda}} is not possible with ff-divergences (except the KLD).”
Sklaviadis et al., Sec. 3.3 · Section 3.3
“On CIFAR, q-VON-MF is competitive with (highly tuned) IVON without additional tuning.”
Sklaviadis et al., Sec. 5 · Section 5
Reproducibility

Reproduction is plausible for an expert reader but not turnkey. The algorithms are specified in detail, and the authors note that only the sampling law, the constants relating $\mathbf{S}$ to the scale, and the corresponding diagonal Stein estimator change relative to IVON. However, no code repository is provided in the text, and key implementation details such as exact hyperparameter schedules for the image-classification experiments are relegated to appendices. The heavy reliance on appendices for derivations and experimental details means independent reproduction would require substantial manual transcription.

“Our new Newton-like optimizers require changing a few lines of code in state-of-the-art Improved Variational Online Newton (IVON) ... to obtain heavy-tailed or bounded-support noise injection, maintain non-Gaussian posteriors, and target a generalized Bayesian objective.”
Sklaviadis et al., Sec. 1 · Section 1
“The difference between q<1 and q>1 lies in the relation between (\mathbf{S},\mathbf{m}) and the location-scale parametrizations, and in the corresponding escort reweighting appearing in the Price identity.”
Sklaviadis et al., Sec. 4 · Section 4
Abstract

The Bayesian learning rule (BLR) exploits the relation between Bayesian updating and the Kullback-Leibler divergence (KLD) to derive efficient additive updates for exponential family variational posteriors. The KL geometry is powerful but restrictive: combined with posterior families with parameter-dependent support it can lead to undefined objectives. We ask which non-KL geometries preserve the additive BLR computation. Previous works that use $f$-divergences suggest that this is only possible with the KLD. Here, we show that this is not true. We show that additive information-geometric optimization algorithms can be derived for the class of mixed-space Bregman or Fenchel-Young discrepancies. These discrepancies induce generalized maximum-entropy families, dual coordinates, and generalized-BLR canonical gradient updates, recovering the classical BLR. The generalized BLR leads to new algorithms that \textit{cannot} be derived from the KL-based BLR. For the Tsallis negentropy with quadratic statistics, the generalized BLR yields bounded-support and heavy-tailed $q$-Gaussian posteriors and an easy Newton-like implementation. Experiments up to CIFAR-100 illustrate controllable support and tail behavior and competitive performance.

Challenge the Review

Pick a starting point or write your own. Challenges run in the background, so you can keep reading while the AI investigates.

No challenges yet. Disagree with the review? Ask the AI to revisit a specific claim.