stat.MLOct 6, 2026

High-dimensional online calibration from harmonic weights

Authors: Maxwell Fishelson, Mehryar Mohri

Organizations: Institute for Advanced Study. · Google Research and Courant Institute of Mathematical Sciences, New York.

Abstract

We study the online calibration of multidimensional forecasts over an arbitrary convex set Y⊆RdY\subseteq\mathbb{R}^d relative to an arbitrary error norm ∥⋅∥L\|\cdot\|_{L}. For forecasting dd binary outcomes simultaneously (Y=[0,1]dY=[0,1]^d), we give the first algorithm that achieves ε\varepsilon-calibration in a number of rounds that is polynomial in dd for every fixed accuracy. It requires dO(1/ε)d^{O(1/\varepsilon)} rounds, exponentially improving the dimension dependence of previous bounds. For multi-class forecasting (Y=ΔdY=Δ_d), we obtain the same dO(1/ε)d^{O(1/\varepsilon)} rate, improving the dO~(1/ε2)d^{\widetilde{O}(1/\varepsilon^2)} bounds of Peng and Fishelson et al. Our algorithm is simple: on each round, it outputs a harmonically weighted distribution over harmonically smoothed past outcomes. The same algorithm works for every forecast set and norm. More generally, it achieves ε\varepsilon-calibration after exp⁡(O(γ(Y,L)/ε))\exp(O(γ(Y,L)/\varepsilon)) rounds, where γ(Y,L)γ(Y,L) is a geometric parameter defined by a matrix discrepancy problem. The harmonic weights are motivated by the fact that the discrete Hilbert transform matrix achieves the optimal discrepancy up to a universal constant, simultaneously for every LL. This optimality result may be of independent interest.

Figures & tables

Explore similar work

Jul 22, 2026stat.ML

Optimal Recalibration of an Online Predictor

We study the problem of recalibrating an online predictor [KE17, OKS24]: given an arbitrary "hint" sequence of forecasts, the learner must output new predictions that are calibrated while incurring small excess error relative to the original forecasts, under a proper loss. We give an online algorithm that achieves (ε,ε2)(\varepsilon, \varepsilon^2)-recalibration for Lipschitz proper losses in T≈ε−3T \approx \varepsilon^{-3} rounds, using an imbalanced extension of the recent simultaneous Blackwell approachability reduction framework of [HTY26]. We show that this tradeoff is optimal by proving a matching lower bound for recalibrating against the squared loss. We also prove a companion K2\mathcal{K}_2-recalibration theorem that obtains the same tradeoffs up to a logarithmic factor. As our main application, we show how our recalibration algorithms can be combined with the online refinement method of [FH23] to obtain simultaneous ε\varepsilon-calibration and ε2\varepsilon^2-calibeating for smooth proper losses at the same asymptotic rate, improving upon prior works that achieved these properties separately or with a worse ε\varepsilon dependence. In particular, the K2\mathcal{K}_2 variant answers a question of [CHJL26] on simultaneously achieving near-optimal calibeating and calibration rates. We also derive extensions to settings with multiple hint sequences. Finally, we empirically evaluate our algorithms on a classification dataset undergoing distribution shift.
Jun 19, 2024cs.LG

Breaking the T2/3T^{2/3} Barrier for Sequential Calibration

A set of probabilistic forecasts is calibrated if each prediction of the forecaster closely approximates the empirical distribution of outcomes on the subset of timesteps where that prediction was made. We study the fundamental problem of online calibrated forecasting of binary sequences under the standard ℓ1\ell_1 calibration error metric, which was initially studied by Foster & Vohra (1998). They derived an algorithm with O(T2/3)O(T^{2/3}) calibration error after TT time steps, and showed a lower bound of Ω(T1/2)Ω(T^{1/2}). These bounds remained stagnant for two decades, until Qiao & Valiant (2021) improved the lower bound to Ω(T0.528)Ω(T^{0.528}) by introducing a combinatorial game called sign preservation and showing that lower bounds for this game imply lower bounds for calibration. In this paper, we give the first improvement to the O(T2/3)O(T^{2/3}) upper bound on calibration error of Foster & Vohra. We do this by introducing a variant of Qiao & Valiant's game that we call sign preservation with reuse (SPR). We prove that the relationship between SPR and calibrated forecasting is bidirectional: not only do lower bounds for SPR translate into lower bounds for calibration, but algorithms for SPR also translate into new algorithms for calibrated forecasting. We then give an improved upper bound for the SPR game, which implies, via our equivalence, a forecasting algorithm with calibration error O(T2/3−ε)O(T^{2/3 - \varepsilon}) for some ε>0\varepsilon > 0, improving Foster & Vohra's upper bound for the first time. Using similar ideas, we then prove a slightly stronger lower bound than that of Qiao & Valiant, namely Ω(T0.54389)Ω(T^{0.54389}). Our lower bound is obtained by an oblivious adversary, marking the first ω(T1/2)ω(T^{1/2}) calibration lower bound for oblivious adversaries.
Jul 14, 2026cs.LG

Efficient Sequential Calibration with O(T2/3−ε)O(T^{2/3-ε}) Error Bound

We study the online binary sequential calibration problem. A recent breakthrough by \citet{dagan2024breaking} overcomes the classical T2/3T^{2/3} barrier for calibration error. Building on this result, we present an efficient randomized forecaster that achieves an expected calibration error O(T2/3−ε)O(T^{2/3-\varepsilon}) for some constant ε>0\varepsilon>0. Our forecaster combines the \textsc{SPR-Calibration} procedure \citep{dagan2024breaking} with an outer Blackwell-style correction layer. The \textsc{SPR-Calibration} procedure controls calibration with respect to a surrogate sequence of conditional-mean estimates, while the correction layer controls the additional error incurred when these surrogates are used to approximate the true outcomes. The analysis decomposes the total calibration error into the surrogate calibration error and the residual discrepancy between the surrogate sequence and the true outcomes. The former is bounded by the \textsc{SPR-Calibration} guarantee in \citet{dagan2024breaking}, and the latter is controlled using a quadratic potential argument together with the sparsity of the \textsc{SPR-Calibration} forecaster.