cs.LGOct 9, 2026

The Sample Complexity of Multicalibration with an ERM Oracle

Authors: Princewill Okoroafor

Organizations: Harvard University

Abstract

A predictor is multicalibrated with respect to a class of subpopulations (formally, of functions on the inputs, called auditors) if, for every subpopulation in the class and every value the predictor outputs, the average outcome among the members of that subpopulation who receive that value is close to the value itself. For finite classes, about ε−3\varepsilon^{-3} samples are necessary and sufficient to learn a predictor with multicalibration error ε\varepsilon, but the algorithms that attain this rate enumerate the class, which is impossible for the large classes of interest. We ask how many samples suffice when the learner accesses the class only through an empirical risk minimization (ERM) oracle. We give an algorithm that uses O~(d/ε3)\tilde{O}(d/\varepsilon^{3}) samples and polynomially many ERM calls, where dd is the VC dimension of the class. Its dependence on ε\varepsilon is optimal for all learners. When dd is at least about 1/ε1/\varepsilon, every learner that can only access the class through polynomially many ERM calls needs Ω~(d/ε3)\tildeΩ(d/\varepsilon^{3}) samples. Learners with unlimited computation need only O~(ε−3+dε−2)\tilde{O}(\varepsilon^{-3}+d\varepsilon^{-2}) samples, so the extra factor 1/ε1/\varepsilon on the dimension term is the price of oracle efficiency. For swap multicalibration, which allows a different function from the class for each predicted value, Θ~(d/ε3)\tildeΘ(d/\varepsilon^{3}) samples are necessary and sufficient. For estimating the error of a given predictor, the ordinary notion is the harder one: the swap error can be estimated with an ERM oracle, whereas the ordinary error is NP-hard to approximate even for classes with linear-time ERM. Consequently, there are classes for which it is computationally easier to learn a multicalibrated predictor than to compute the multicalibration error of a given predictor.

Explore similar work

Apr 23, 2026

The Sample Complexity of Multicalibration

We study the minimax sample complexity of multicalibration in the batch setting. A learner observes nn i.i.d. samples from an unknown distribution and must output a (possibly randomized) predictor whose population multicalibration error, measured by Expected Calibration Error (ECE), is at most ε\varepsilon with respect to a given family of groups. For every fixed κ>0κ> 0, in the regime ∣G∣≤ε−κ|G|\le \varepsilon^{-κ}, we prove that Θ~(ε−3)\widetildeΘ(\varepsilon^{-3}) samples are necessary and sufficient, up to polylogarithmic factors. The lower bound holds even for randomized predictors, and the upper bound is realized by a randomized predictor obtained via an online-to-batch reduction. This separates the sample complexity of multicalibration from that of marginal calibration, which scales as Θ~(ε−2)\widetildeΘ(\varepsilon^{-2}), and shows that mean-ECE multicalibration is as difficult in the batch setting as it is in the online setting, in contrast to marginal calibration which is strictly more difficult in the online setting. In contrast we observe that for κ=0κ= 0, the sample complexity of multicalibration remains Θ~(ε−2)\widetildeΘ(\varepsilon^{-2}) exhibiting a sharp threshold phenomenon. More generally, we establish matching upper and lower bounds, up to polylogarithmic factors, for a weighted LpL_p multicalibration metric for all 1≤p≤21 \le p \le 2, with optimal exponent 3/p3/p. We also extend the lower-bound template to a regular class of elicitable properties, and combine it with the online upper bounds of Hu et al. (2025) to obtain matching bounds for calibrating properties including expectiles and bounded-density quantiles.
Aug 4, 2026

Sample Complexity of Multicalibration for Multilevel Properties

Calibration requires a predictor to be unbiased after conditioning on its own predictions. Multicalibration asks for this guarantee simultaneously across a collection of groups. Many prediction tasks ask for several related features of the same conditional outcome distribution: variance is defined relative to the mean, skewness relative to both mean and variance, and conditional value at risk relative to a quantile. We study multicalibration for a sequence of kk properties in which each property is identifiable once the preceding properties are fixed. This framework includes Bayes pairs but does not require the properties to arise from a single loss. For every fixed k≥2k\ge2, we establish matching upper and lower sample-complexity bounds up to logarithmic factors under regularity conditions. Even with only polylogarithmically many binary groups, achieving multicalibration error ε\varepsilon requires Ω~(ε−(k+2))\widetildeΩ(\varepsilon^{-(k+2)}) samples. Conversely, for any finite group family G\mathcal G, we give a randomized learner using O(ε−(k+2)+ε−2log⁡∣G∣)O(\varepsilon^{-(k+2)}+\varepsilon^{-2}\log|\mathcal G|) samples. Thus the sample complexity is Θ~(ε−(k+2))\widetildeΘ(\varepsilon^{-(k+2)}) for polynomial-size group families. We instantiate the theory for three canonical examples.
Oct 11, 2026

Sample-Optimal Multi-Group Learning with an ERM Oracle

In multi-group learning, a single predictor must perform nearly as well as the best hypothesis in a class on every group in a rich, possibly infinite and overlapping, collection of subpopulations. The known algorithms with the optimal sample complexity keep track of every pair of a group and a hypothesis, which is impossible for the large classes of interest, and the known oracle-based algorithms have a worse rate. We ask whether the optimal rate can be attained by a learner that accesses the classes only through an empirical risk minimization (ERM) oracle. We give an algorithm that uses O~((dG+dH)/ε2)\tilde{O}((d_G+d_H)/\varepsilon^2) samples and polynomially many calls to a weighted ERM oracle over pairs of a group and a hypothesis, where dGd_G and dHd_H are the VC dimensions of the group and hypothesis classes and ε\varepsilon bounds the excess risk on every group, weighted by its probability. We also show a matching lower bound of Ω((dG+dH)/ε2)Ω((d_G+d_H)/\varepsilon^2) samples for every learner even with unlimited computation, so oracle efficiency costs no additional samples. We also consider a natural special case in which every input belongs to at most kk groups and these groups can be listed. In this case, an ERM oracle for the hypothesis class alone suffices, with O~((dH+log⁡k)/ε2)\tilde{O}((d_H+\log k)/\varepsilon^2) samples in the agnostic setting and O~((dH+log⁡k)/ε)\tilde{O}((d_H+\log k)/\varepsilon) samples when every group has a hypothesis that is correct on it.