The science of systems that learn from data and experience. Its scope runs from the statistical and mathematical foundations of learning, including generalization, expressivity, and computational limits, through the design of learning algorithms, deep learning, reinforcement learning, and probabilistic methods, to the empirical study of large models and the trustworthiness, interpretability, and societal impact of learned systems.
Every time a streaming service guesses what you would rate a film you have never seen, it is solving a matrix completion problem: fill in the missing entries of a vast user-by-item table from the few that are observed. The question became famous during the Netflix Prize (2006-2009), and it looks hopeless - infinitely many matrices fit the observed entries - until one assumes the structure that makes recommendation possible: the table is essentially low rank, because tastes are governed by a few latent factors. In their landmark 2009 paper 'Exact Matrix Completion via Convex Optimization' (Foundations of Computational Mathematics), Emmanuel Candes and Benjamin Recht proved that an n-by-n matrix of rank r can be recovered exactly, with high probability, from only about n^1.2 * r * log n randomly observed entries - not by the NP-hard route of minimizing rank, but by minimizing the nuclear norm, a convex surrogate (the sum of the singular values) solvable efficiently. The proof, in the lineage of Candes-Romberg-Tao compressed sensing, turns on two ideas: an incoherence condition ensuring the singular vectors are spread out rather than spiky, and a dual certificate witnessing optimality, whose existence rests on delicate random-matrix concentration. It transformed a practical engineering puzzle into rigorous theory and seeded a decade of work across machine learning, signal processing, computer vision, and sensor localization. This mission formalizes the Candes-Recht exact-recovery theorem in Lean, decomposed into its dual-certificate construction and the probabilistic concentration reductions at its core.
Bandit Algorithms III: Asymptotic and Minimax Optimality of UCBTextbook
The basic UCB regret bound of Mission II is logarithmic but not tight: its leading constant 16/Δi is eight times the information-theoretic limit, and its worst-case rate carries a spurious logn. Chapters 8–9 of Lattimore–Szepesvári close both gaps. A refined confidence schedule f(t)=1+tlog2t yields the asymptotically optimal limsupn→∞Rn/logn≤∑i:Δi>02/Δi — exactly matching the instance-dependent lower bound of Mission VII for Gaussian noise. The MOSS index μ^i+Ti4log+(kTin) achieves minimax regret Rn≤39kn+∑iΔi, matching the Ω(kn) lower bound up to a constant. These two theorems are the gold standard for finite-armed stochastic bandits.
Bandit Algorithms II: Stochastic Bandits and the UCB AlgorithmTextbook
A learner repeatedly chooses one of k slot machines, observes only the reward of the chosen arm, and wants to earn almost as much as the best arm in hindsight. This is the stochastic multi-armed bandit, the canonical model of the exploration–exploitation dilemma. This mission formalizes the model (environments, policies, regret, and the regret decomposition Rn=∑iΔiE[Ti(n)]) and the two classical algorithms of Chapters 6–7 of Lattimore–Szepesvári: Explore-Then-Commit and the Upper Confidence Bound algorithm built on the optimism principle. The goal theorem is the instance-dependent UCB regret bound Rn≤3∑iΔi+∑i:Δi>016log(n)/Δi — logarithmic regret with explicit constants, the single most cited result of bandit theory — together with its distribution-free companion Rn≤8nklogn+3∑iΔi.
Bandit Algorithms XIV: Bayesian Bandits, the Gittins Index and Thompson SamplingTextbook
The oldest bandit algorithm (Thompson, 1933) is also the most modern: sample a parameter from the posterior and act greedily. Chapters 34–36 of Lattimore–Szepesvári develop the Bayesian view in two crowning results. The Gittins index theorem: for infinite-horizon discounted Markov bandits, the seemingly intractable dynamic program is solved exactly by an index policy — each arm gets a retirement-value index computable arm-by-arm, and playing the largest index is Bayesian optimal. And the frequentist analysis of Thompson sampling — the goal theorem: with Gaussian posteriors, Thompson sampling on 1-subgaussian bandits achieves limn→∞Rn/logn=∑i:Δi>02/Δi, exactly asymptotically optimal, alongside the minimax-grade Rn≤Cknlogn. Together they explain why posterior sampling is both principled and practically dominant.
Bandit Algorithms XI: Lower Bounds for Stochastic Linear BanditsTextbook
Is the dn regret of LinUCB (Mission X) an artifact of the algorithm or a law of nature? Chapters 24–25 of Lattimore–Szepesvári prove it is essentially unimprovable. On the unit ball there is a parameter θ with ∥θ∥22=d2/(48n) forcing Rn≥163dn — the goal theorem — and the hypercube gives the same Ω(dn) rate. The asymptotic chapter is more striking still: for fixed finite action sets, the instance-optimal constant c(A,θ) is characterized by an allocation program, and optimism itself is provably suboptimal — LinUCB and Thompson sampling cannot achieve it, because exploration must sometimes deliberately play actions optimism would never touch. These lower bounds define the targets for the entire linear-bandit literature.
Bandit Algorithms V: Adversarial Bandits and Exp3Textbook
What if the rewards are not random at all, but chosen by an adversary who knows your algorithm? Remarkably, a randomized learner can still compete with the best fixed arm in hindsight. Chapters 11–12 of Lattimore–Szepesvári develop the adversarial k-armed bandit: rewards xti∈[0,1] are an arbitrary fixed matrix, the learner samples At∼Pt, and regret is measured against maxi∑txti. The exponential-weights algorithm Exp3, fed by importance-weighted loss estimates X^ti=1−1{At=i}(1−Xt)/Pti, achieves Rn≤2nklogk — the goal theorem. The companion Exp3-IX, which deliberately biases its estimator, upgrades this to a bound holding with high probability rather than only in expectation. These results are the foundation of all adversarial online learning with partial feedback.
Learning and Optimization with Seasonal Patterns 2: Regret Lower Bound for Periodic Multi-Armed Bandits with Known PeriodsResearch Paper
Motivation
Demand, prices and click-through rates often follow seasonal cycles: weekly, monthly, yearly. A decision maker who repeatedly chooses among options whose expected payoffs repeat in time faces a non-stationary multi-armed bandit. Chen, Wang and Wang (arXiv:2005.08088v4, Operations Research, 2025) study the model in which each arm's mean reward is periodic with its own unknown period. They give a two-stage policy whose regret is O~(T∑kTk), and a lower bound showing that no policy, even one told the periods, does better than Ω(TmaxkTk). This mission formalizes that lower bound, Theorem 3 of the paper. A companion mission of the same series covers the upper bound.
Lower bounds of this kind fix what is achievable. For the stationary K-armed bandit, the minimax regret is of order TK (Auer, Cesa-Bianchi, Freund and Schapire, 2002; Lattimore and Szepesvári, Bandit Algorithms, 2020, Ch. 15). In the periodic model a learner effectively faces ∑kTk "effective arms" (arm k at each phase of its cycle), which suggests the rate T∑kTk. The paper proves the weaker but general bound TmaxkTk/(32e), which matches the upper bound up to a factor K and logarithms.
Setting
There are K≥2 arms and epochs t=1,…,T. An instanceν is a family of mean rewards μk,t. If arm k is pulled at epoch t, the reward is drawn from the Gaussian law N(μk,t,1), independently of the past given the arm.
A sequence f has minimum periodp≥1 if f(t+p)=f(t) for all t≥1 and no smaller positive integer has this property.
A policyπ chooses the arm πt at epoch t at random, from a distribution that depends on the previous arms and rewards. Its expected regret on ν is display (1) of the paper:
Eν[RTπ]=t=1∑T(kmaxμk,t−E[μπt,t]).
Three classes of instances appear, all with unit-variance Gaussian rewards and means in [0,1] (Assumption 1):
E1: every arm is stationary;
E2(T1): the first arm's mean has minimum period T1, the other arms are stationary;
E3(T1,…,TK): arm k's mean has minimum period Tk.
The proof compares instances through the metric (B-48),
The periods satisfy Assumption 3 with the parameters (7): with n=⌊T/K⌋ and g=⌈n⌉, T>4K, g≥2 and Tk<n/(2g) for every k.
Formalization targets
Goal: Theorem 3 (p. 24)
For every policy π there is ν∈E3(T1,…,TK) with
Eν[RTπ]≥32e1TkmaxTk.
Milestones
Lemma 12 (p. 23). Stationary Gaussian bandits whose arms other than the first have known mean 1/2: every policy has regret at least T/(82e) on some choice of μ1∈[0,1].
Lemma 13 (p. 23). One periodic arm: every policy has regret at least TT1/(16e) on some instance of E2(T1) whose other arms have mean 1/2.
Fact 1 (proof of Theorem 3, p. 48). E3 is dense in E2∪E3 under M: an instance of E2(T1) can be perturbed on the arms other than the first into an instance of E3(T1,…,TK) at M-distance less than any δ>0.
Fact 2 (proof of Theorem 3, p. 48). The regret fπ(ν)=Eν[RTπ] is continuous in ν on E2∪E3 for M, uniformly in π.
Significance
Theorem 3 shows that knowing the periods does not remove the cost of learning them phase by phase: the regret must grow like TmaxkTk, not like TK. Together with the paper's upper bound O~(T∑kTk) and the inequality TmaxkTk≥T∑kTk/K, it shows the two-stage policy is near-optimal when K is fixed. The intermediate results are reusable: Lemma 12 is a clean two-point Gaussian lower bound, Lemma 13 is a phase-decomposition argument that applies to any bandit model with known cyclic structure, and facts 1 and 2 give a continuity framework for the regret as a function of the instance.
The result has a published proof. As far as the platform's records show, none of it has been machine-checked. The stationary minimax lower bound of Lattimore and Szepesvári (a different family and constant) exists on the platform, together with the Bretagnolle–Huber inequality, the Gaussian KL divergence and the divergence decomposition, which this mission lists as tools.
Difficulty
The two-point argument of Lemma 12 is standard, but its constant needs Δ=1/(2T)<1/2, so small horizons must be handled separately. Lemma 13 decomposes the problem into T1 phase problems; making this rigorous means showing that a policy on the periodic problem induces, on each phase, a policy on a stationary problem of horizon ⌊T/T1⌋, with a history law that factorizes accordingly. The phase-wise instances produced by the decomposition need not have minimum period exactly T1 (all phases may coincide), so a perturbation is needed to land in E2(T1); the constant has slack for it. Fact 2 asks for continuity of an integral against the law of a T-round adaptive process, uniformly over all policies, which needs a domination argument that does not depend on the policy's kernels.
Formalization scope
Arms are Fin K, 0-based: the paper's first arm is index 0. Epochs are 1-based natural numbers; history coordinate i is epoch i+1. Values of the means at t=0 are ignored.
A policy is the platform's BanditAlgorithm.BanditPolicy K: one Markov kernel per round from the history to the next arm. The history law periodicMeasure reuses the platform's step kernel with the Gaussian bandit of the current epoch's means; with constant means it is BanditAlgorithm.banditMeasure of BanditAlgorithm.gaussianBandit. Lemma 12 is stated with the platform's banditRegret, Tμ∗−E[∑tXt], which equals (1) for Gaussian arms; the other statements use (1) directly.
Assumption 1 is stated on epochs 1,…,T. Minimum periods are exact.
Every lower bound has the order "for every policy there is an instance", with the instance in the class, so the policy may depend on K, T and the periods; this is the paper's "policy that knows Tk". The supremum over the class is replaced by an attained instance, which is the stronger form.
Added hypotheses, all disclosed: K≥2 (with one arm the regret is 0); Tk≥1 (the paper's Tk∈N+); in Lemmas 12 and 13, the means of the arms other than the first are 1/2, which is how "a policy that knows μk for k≥2" is read (with arbitrary known values the printed lemmas fail). Assumption 3 with (7) is the standing assumption of Appendix B.4.
A statement with the quantifiers in the order "there is an instance for every policy", without the [0,1] clause, or without exact minimum periods, would be false or trivial; none of these is the target.
Welcome contributions: the bridge between periodicMeasure with constant means and banditMeasure, the phase decomposition of the history law, and the uniform continuity of the regret in the means.
Selected references
N. Chen, C. Wang, L. Wang, Learning and Optimization with Seasonal Patterns, Operations Research, 2025; arXiv:2005.08088v4, 2021. https://arxiv.org/abs/2005.08088v4
P. Auer, N. Cesa-Bianchi, Y. Freund, R. E. Schapire, The nonstochastic multiarmed bandit problem, SIAM Journal on Computing 32(1), 2002. https://doi.org/10.1137/S0097539701398375
Improved Algorithms for Linear Stochastic Bandits I: High-Probability Regret Bound for the OFUL AlgorithmResearch Paper
Motivation
In a linear stochastic bandit, a learner repeatedly chooses an action from a set of vectors and receives a noisy reward whose mean is linear in the action. The model underlies contextual recommendation, adaptive routing and dynamic pricing, where each option is described by features and the payoff of a feature vector must be learned while it is exploited. The quality of a strategy is measured by its regret: the reward lost, relative to always playing the best action, over the first n rounds.
The optimism-in-the-face-of-uncertainty principle (play as if the most favourable parameter consistent with the data were true) was introduced for linear bandits by Auer (2002), and developed by Dani, Hayes and Kakade (2008) (ConfidenceBall, regret O(dnlog3/2n) with confidence sets from a union bound over time) and Rusmevichientong and Tsitsiklis (2010). Abbasi-Yadkori, Pál and Szepesvári (NIPS 2011) replaced the union bound by a self-normalized martingale inequality that holds uniformly in time. It gives smaller confidence ellipsoids and, through them, a high-probability regret bound for the resulting algorithm, OFUL, that improves the earlier ones by logarithmic factors. The inequality became the standard tool for linear and kernelized bandits and for linear reinforcement learning.
Setting
Fix a dimension d≥1 and an unknown parameter θ∗∈Rd. In round t=1,2,… the learner is given a nonempty decision setDt⊆Rd, chooses Xt∈Dt, and observes the reward
Yt=⟨Xt,θ∗⟩+ηt.
There is a filtration {Ft}t≥0 such that Xt is Ft−1-measurable and ηt is Ft-measurable and conditionally R-sub-Gaussian: E[eληt∣Ft−1]≤exp(λ2R2/2) for all λ∈R, with R≥0 fixed.
For a regularization parameter λ>0 let Vt=λI+∑s=1tXsXs⊤ and let θt=Vt−1∑s=1tYsXs be the regularized least-squares estimate. With ∥v∥A=v⊤Av and a known bound ∥θ∗∥2≤S, the confidence ellipsoid is
The OFUL algorithm chooses, in round t, a pair (Xt,θt) maximizing ⟨x,θ⟩ over Dt×Ct−1. The pseudo-regret is Rn=∑t=1n⟨xt∗−Xt,θ∗⟩, where ⟨xt∗,θ∗⟩=maxx∈Dt⟨x,θ∗⟩.
Formalization targets
Goal: Theorem 3, the regret of OFUL
If ∥Xt∥2≤L, ⟨x,θ∗⟩∈[−1,1] for all x∈Dt, and λ≥max(1,L2), then for every δ>0, with probability at least 1−δ,
For any positive definite V, Vt=V+∑s≤tXsXs⊤ and St=∑s≤tηsXs: with probability at least 1−δ, for all t≥0,
∥St∥Vt−12≤2R2log(det(Vt)1/2det(V)−1/2/δ).
Milestones: Theorem 2, the confidence ellipsoids
With probability at least 1−δ, θ∗∈Ct for all t≥0 (first claim). If ∥Xt∥2≤L, then with probability at least 1−δ, for all t, ∥θt−θ∗∥Vt≤Rdlog((1+tL2/λ)/δ)+λ1/2S (second claim, stated here for d≥2).
Significance
Theorem 3 bounds the regret of OFUL by O(dnlogn) with high probability, uniformly over the horizon, so it holds for an unknown horizon without restarting. The bound applies to arbitrary, even adversarially changing, decision sets. Theorem 1 is the ingredient that makes this possible: a deviation bound for a vector-valued martingale, normalized by its own random covariance, that holds for all times simultaneously and whose logarithmic term is a determinant rather than a union-bound count. The same inequality underlies regret analyses of generalized linear bandits, kernelized bandits, linear Markov decision processes and many confidence-sequence constructions.
All three results are proved in the paper's appendices (not included in the source file used here). None of them is formalized in the stated generality. Prove2Me holds the special cases V=λI, R=1, δ<1 of Theorems 1 and 2 (from the Bandit Algorithms textbook series), a pathwise LinUCB regret lemma that assumes the confidence event, and the elliptical potential lemma. A formal proof of Theorem 3 would be the first machine-checked high-probability regret bound for OFUL with the paper's confidence sets.
Difficulty
The actions are chosen adaptively, by an argmax over a data-dependent set, so the sequence Xt has no independence structure and the least-squares estimate is not a sum of independent terms. A fixed-design concentration bound followed by a union bound over time and over a covering of the sphere loses logarithmic factors and does not produce the determinant in the radius; that is the route of the earlier work that Theorem 1 improves. Theorem 1 must hold for all times at once for a quantity normalized by the random matrix Vt, which is itself built from the adaptively chosen actions; a bound for each fixed t does not give it.
Formalization scope
Vectors are Fin d → ℝ, matrices Matrix (Fin d) (Fin d) ℝ, and ∥x∥A is Real.sqrt (x ⬝ᵥ A *ᵥ x). Rounds are indexed t+1 for t:N, so sums over s≤t are sums over Finset.range t at index s + 1, and the time-0 objects are empty sums. The probability space is standard Borel, as Mathlib's conditional sub-Gaussianity (HasCondSubgaussianMGF, variance proxy R2) requires; this is an added hypothesis. Every "with probability at least 1−δ, for all t" is stated as an outer-measure bound ≤δ on the failure event, with the time quantifier inside the event. det(⋅)1/2 is the real square root of the determinant; the matrices inverted are positive definite, so Lean's junk inverse never occurs.
OFUL is a predicate on the whole process: in every round the chosen pair maximizes ⟨x,θ⟩ over Dt×Ct−1, with any tie-breaking. Runs exist whenever the decision sets are nonempty and compact. The measurability of the actions is assumed, as in Theorem 1. The optimal reward ⟨xt∗,θ∗⟩ is the supremum over Dt, finite because of the reward bound.
Two corrections to the printed Theorem 3 are made and disclosed. The printed nL/d is replaced by nL2/d, which is what the determinant–trace bound detVn≤(λ+nL2/d)d gives; for L≤1 the corrected bound implies the printed one. The hypothesis λ≥max(1,L2) is added: for λ<1 the printed logarithm can be negative, and the printed bound would then assert Rn≤0. In the second claim of Theorem 2, d≥2 is added, because at d=1 the claim does not follow from the first claim and the paper's proof is not available.
The goal is a probability bound over the noise, not the pathwise statement "if θ∗∈Ct−1 for all t then Rn≤…". The pathwise statement assumes the confidence event instead of proving it, and is already on the platform. The confidence sets inside the OFUL predicate use the same δ as the conclusion.
Needed infrastructure: maximal inequalities for nonnegative supermartingales, Gaussian integrals of quadratic forms on Rd, log-determinant bounds for sums of rank-one updates, and the determinant–trace inequality. The platform rows BanditAlgorithm.self_normalized_martingale_bound, BanditAlgorithm.least_squares_confidence_ellipsoid and BanditAlgorithm.elliptical_potential_lemma are referenced as tools. Proofs of the milestones, generalizations of the existing special cases to general V and R, and reusable determinant lemmas are all welcome.
Selected references
Y. Abbasi-Yadkori, D. Pál, Cs. Szepesvári, Improved Algorithms for Linear Stochastic Bandits, Advances in Neural Information Processing Systems 24 (NIPS), 2011. https://proceedings.neurips.cc/paper/2011
P. Rusmevichientong, J. N. Tsitsiklis, Linearly Parameterized Bandits, Mathematics of Operations Research 35(2), 2010. https://doi.org/10.1287/moor.1100.0446
High-Dimensional Statistics I: Gaussian Concentration of Lipschitz FunctionsTextbook
Motivation
A recurring question in high-dimensional statistics is how tightly a scalar quantity built from
many random inputs concentrates around its mean, even as the number of inputs grows without
bound. Two classical answers organize the whole toolkit: martingale methods, which control a
sum of dependent increments one conditional step at a time, and Gaussian-specific isoperimetry,
which shows that essentially any regular (Lipschitz) function of a high-dimensional Gaussian
vector concentrates as tightly as a single Gaussian coordinate, regardless of dimension. This
mission formalizes one representative theorem from each line: the general martingale
Bernstein bound (Wainwright, High-Dimensional Statistics, 2019, Theorem 2.19) and the
Gaussian concentration of Lipschitz functions (Theorem 2.26), following Chapter 2 of the same
book.
Setting
A random variable X with mean μ=E[X] is sub-Gaussian with parameter σ
(Definition 2.2) if E[eλ(X−μ)]≤eσ2λ2/2 for all
λ∈R; it is sub-exponential with parameters (ν,α) (Definition 2.7,
a strictly milder condition) if the same bound holds only for ∣λ∣<1/α, with the
convention 1/0=+∞ so that α=0 recovers the sub-Gaussian case exactly.
A sequence {Dk}k≥1, adapted to a filtration {Fk}, is a martingale
difference sequence if each Dk is Fk-measurable and
E[Dk∣Fk−1]=0. Such sequences arise throughout statistics via the
Doob martingale construction: given a function f of independent variables
X1,…,Xn, setting Dk:=E[f(X)∣X1,…,Xk]−E[f(X)∣X1,…,Xk−1] telescopes to f(X)−E[f(X)]=∑kDk, converting a deviation
question about f(X) into a martingale concentration question.
A function f:Rn→R is L-Lipschitz with respect to the Euclidean norm
if ∣f(x)−f(y)∣≤L∥x−y∥2 for all x,y (Eq. (2.38)).
Formalization targets
Goal — Theorem 2.26 (Gaussian concentration of Lipschitz functions)
Let (X1,…,Xn) be i.i.d. standard Gaussian and f be L-Lipschitz with respect to the
Euclidean norm. Then f(X)−E[f(X)] is sub-Gaussian with parameter at most L, and
hence
P[∣f(X)−E[f(X)]∣≥t]≤2e−t2/2L2for all t≥0.
The bound is dimension-free: it depends on n only through f's Lipschitz constant, not the
ambient dimension itself.
For any differentiable f and convex φ,
E[φ(f(X)−E[f(X)])]≤E[φ(2π⟨∇f(X),Y⟩)] for X,Y∼N(0,In) independent — the interpolation identity Theorem 2.26's
proof is built on.
Given a martingale difference sequence with a per-index sub-exponential conditional
moment-generating-function bound E[eλDk∣Fk−1]≤eλ2νk2/2 for ∣λ∣<1/αk, the sum ∑kDk is itself sub-exponential
with parameters (∑kνk2,maxkαk), and satisfies the two-regime
concentration inequality of Eq. (2.28): sub-Gaussian for small deviations, sub-exponential for
large ones. This is the chapter's central general-purpose martingale concentration tool.
Significance
Theorem 2.19 is the source of two of the most-cited concentration inequalities in the field —
the Azuma–Hoeffding inequality (Corollary 2.20) and the bounded-differences/McDiarmid
inequality (Corollary 2.21), both already faithfully covered elsewhere on the platform
(azuma_hoeffding_two_sided, bounded_diff_martingale_two_sided) and included here as
kind: reference milestones rather than redrafted. Theorem 2.26's Gaussian Lipschitz
concentration is separately significant: it is the tool behind dimension-free operator-norm
bounds for random matrices, concentration of the empirical spectral distribution, and much of
the machinery of Chapters 5 and 6 of the same book.
Formalizing it. No faithful prior art exists on the platform for either the martingale
Bernstein bound or Lipschitz-Gaussian concentration itself (a fresh search for "martingale
Bernstein," "sub-exponential martingale," "Gaussian interpolation," and "Lipschitz
concentration" returned no hits; the existing Vershynin-book item
HighDimProb.Isoperimetry.lipschitz_concentration_sphere concentrates a Lipschitz function on
the sphere, a different underlying space and a different proof from Theorem 2.26's Gaussian
vector). Both goal-adjacent theorems and the Gaussian interpolation lemma are drafted here as
open goals (:= by sorry); the two Azuma–Hoeffding/bounded-differences corollaries are reused
from the platform's existing, already-proved formalizations.
Difficulty
The naive approach to Theorem 2.26 — try to bound f(X)−E[f(X)] directly via a
Lipschitz-type argument in Rn — has no obvious route to a dimension-free bound,
since a union bound over coordinates (or over an ε-net of the domain) picks up a
factor that grows with n. The resolution, Lemma 2.27's interpolation identity, instead
exploits a special structural fact about the Gaussian distribution — its rotation invariance —
to replace the nonlinear quantity f(X)−E[f(X)] with the linear, and hence exactly
computable, Gaussian quantity ⟨∇f(X),Y⟩, at the mild cost of a
non-optimal constant. Theorem 2.19's difficulty is bookkeeping rather than a conceptual
obstruction: the recursive conditioning step (Eq. (2.29)) must be iterated exactly n times
while keeping track of the interplay between the two parameters νk,αk per
difference, and Proposition 2.9's two-regime tail bound (small-deviation sub-Gaussian behavior,
large-deviation sub-exponential behavior) must be carried through unchanged into the final
statement — dropping either regime understates what the theorem proves.
Formalization scope
Expectations are Bochner integrals against an explicit probability measure, with integrability
required as an explicit hypothesis in IsSubGaussian and IsSubExponential (Mathlib's Bochner
integral silently returns 0 for a non-integrable function, which this mission's definitions
rule out as a trivializing formalization). The sub-exponential condition's domain restriction
|λ| < 1/α is realized as the disjunction α = 0 ∨ |λ| < 1/α, since Lean's real division
convention 1/0 = 0 is exactly backwards from the book's own stated 1/0 = +\infty convention
for the degenerate sub-Gaussian case.
"X,Y∼N(0,In) independent" (Lemma 2.27, Theorem 2.26) is formalized via Mathlib's
HasGaussianLaw predicate together with explicit coordinatewise mean-zero and
identity-covariance hypotheses, which together pin down the standard multivariate normal law,
plus IndepFun. The inner product ⟨∇f(X),Y⟩ is realized as
fderiv ℝ f (X ω) (Y ω), the Fréchet derivative applied to Y(ω) — equal to
⟨∇f(X(ω)),Y(ω)⟩ by the Riesz representation of the gradient on a
Hilbert space, avoiding the need to separately construct a gradient vector field.
Theorem 2.19's printed parameter pair for part (a), "(∑kνk2,α∗)," is formalized
as (∑kνk2,α∗): Definition 2.7 parametrizes the sub-exponential MGF bound
by ν (with ν2 appearing in the exponent), so a literal transcription of the printed
pair's first entry would silently square the effective parameter and make part (a), read
literally, inconsistent with part (b)'s own tail-bound formula (which the book derives from
part (a) via the general sub-exponential tail bound, Proposition 2.9). The corrected pairing is
the one the book's own proof actually establishes; see MODERATION_NOTES.md for the full
derivation.
Out of scope for this mission: Proposition 2.5 (plain Hoeffding for a sum of independent
sub-Gaussians), used in the book only as background for Theorem 2.26's proof and not redrafted,
since the goal theorem's own statement does not depend on it once Lemma 2.27 is in hand.
Selected references
M. J. Wainwright, High-Dimensional Statistics: A Non-Asymptotic Viewpoint, Cambridge
University Press, 2019. DOI: 10.1017/9781108627771.
Chapter 2.
K. Azuma, "Weighted sums of certain dependent random variables," Tôhoku Mathematical
Journal, 19:357–367, 1967.
W. Hoeffding, "Probability inequalities for sums of bounded random variables," Journal of the
American Statistical Association, 58:13–30, 1963.
Foundations of Reinforcement Learning III: Structured Bandits and the Decision-Estimation CoefficientTextbook
Motivation
Every algorithm in the first three chapters of Foster and Rakhlin's Foundations of
Reinforcement Learning and Interactive Decision Making — ε-Greedy and UCB for the
multi-armed bandit, Inverse Gap Weighting and SquareCB for contextual bandits — is a
special case of the same two-step recipe: estimate a model of the world with an online
regression oracle, then convert the estimate into a decision that trades exploration
against exploitation. Chapter 4 asks whether this recipe can be made generic: given
any structured decision-making problem, specified only by a function class F and a
decision space Π, is there a single quantity that governs the best achievable
regret, the way A/γ governs the multi-armed bandit and d/γ
governs the linear bandit? The chapter's answer is the Decision-Estimation
Coefficient (DEC), introduced by Foster, Kakade, Qian, and Rakhlin [40] as a
complexity measure that both upper- and lower-bounds achievable regret for a general
decision-making protocol, unifying results that were previously proved from scratch,
case by case, for each structured setting. This mission formalizes the chapter's
central upper bound (Proposition 13) together with the machinery that makes it
computable in two concrete cases — the multi-armed bandit (Proposition 14) and the
linear bandit (Propositions 16–17).
Setting
Fix a finite decision space Π and a class F⊆RΠ of candidate
mean-reward functions, with a ground-truth f⋆∈F (realizability). Over T
rounds, at each round t the learner observes an estimate f^t produced by an
online regression oracle, plays a decision distribution pt∈Δ(Π)
(possibly depending on f^t and the history), and the regret is
Reg:=t=1∑Tf⋆(π⋆)−t=1∑TEπ∼pt[f⋆(π)],
where π⋆=argmaxπf⋆(π). The oracle's cumulative estimation
error is assumed bounded: ∑t=1TEπ∼pt[(f^t(π)−f⋆(π))2]≤EstSq(F,T,δ) with probability at least 1−δ
(Definition 7). Writing πf:=argmaxπf(π), the DEC game value at a
reference model f^ and scale γ>0 is the min-max quantity
and the DEC of F itself is decγ(F):=supf^∈co(F)decγ(F,f^). The Estimation-to-Decisions (E2D) algorithm plays,
at each round, a pt certifying (i.e. attaining or beating) the value of this min-max
game at f^t.
Formalization targets
Goal — Proposition 13 (E2D regret bound)
Reg≤decγ(F)⋅T+γ⋅EstSq(F,T,δ)
with probability at least 1−δ, for any exploration parameter γ>0. This
is the weakest stable statement the chapter proves about E2D: it holds for an
arbitrary function class and an arbitrary regression oracle, with no structural
assumption on F beyond realizability, and the chapter's later sections instantiate
it rather than strengthen it.
Milestones
Lemma 9 (Decoupling), general form: for any distribution ν over a finite
model class and any fˉ, Ef∼ν[f(πf)−fˉ(πf)]≤A⋅Ef∼νEπ∼p[(f(π)−fˉ(π))2] —
the estimation-to-decisions bridge the whole chapter's approach rests on, decoupling
the model index from the played decision.
Proposition 14 (IGW minimizes the DEC): for the multi-armed bandit (Π=[A],
F=RA), Inverse Gap Weighting is the exact minimizer of the DEC game,
giving decγ(F)=(A−1)/(4γ) — the first concrete computation of
an abstract quantity, recovering Chapter 3's rate from Proposition 13 alone.
Proposition 16 (G-optimal design): existence, for any compact
full-dimensional-span set Z⊆Rd, of a distribution p with
supz∈Z⟨Σp−1z,z⟩≤d — the classical convex-analysis
primitive Proposition 17 needs.
Proposition 17 (DEC for linear bandits): combining the G-optimal design with
inverse gap weighting gives decγ(F)≲d/γ for the linear
bandit function class, leading via Proposition 13 to a dT regret bound.
Significance
The Decision-Estimation Coefficient is, in the book's own words, "the main result" of
this line of work: Foster, Kakade, Qian, and Rakhlin [40] show it is not merely an
upper bound but (in a suitable localized form, developed further in Chapter 6) a
tight characterization of the minimax regret for structured bandits and, more
generally, for the interactive decision-making protocol the rest of the book studies.
Proposition 13 is the mechanism that makes this useful in practice: it reduces regret
analysis for a new structured problem to a single, purely convex-analytic computation
of decγ(F), in place of a bespoke exploration argument. Propositions
14–17 are the demonstration that this reduction is not vacuous — they recompute, via
the DEC alone, the two rates (multi-armed and linear bandit) that earlier chapters of
the book derived by direct, setting-specific arguments, and the match is exact.
Formalizing this chapter therefore captures the book's unifying abstraction itself,
not just one more instance of it. No formalization of the Decision-Estimation
Coefficient, in any form, currently exists on the platform (see Formalization scope).
Difficulty
The obvious formalization mistake is to state Proposition 13's conclusion with
decγ(F) left as an unconstrained free real-number parameter satisfying
only the inequality the theorem asserts — a formalization under which the "theorem"
would be a triviality about an arbitrary real number, since nothing about the actual
min-max game would ever be checked. The chapter's content is precisely the opposite:
that this specific minimax quantity can be computed (Proposition 14) or bounded
via a concrete strategy (Proposition 17), and — as Chapter 6 shows for a lower bound
outside this chunk's scope — that no smaller quantity would do. A second difficulty is
proof-theoretic rather than notational: the book's own proof of Proposition 13 bounds
regret by an unconstrained supremum over all reference functions f^:Π→R, and only identifies this with the official, co(F)-restricted
decγ(F) of Eq. (4.16) via Proposition 24 — a fact stated on p. 80,
outside this chapter's numbered range, whose own proof the book defers to an exercise.
A formalization that quietly imports Proposition 24 to close this gap would rest the
goal theorem on an unverified fact; this mission instead states the hypothesis the
book's own text uses to motivate restricting to co(F) in the first place
(online estimation algorithms produce f^t∈co(F)), so the goal is
faithful to what is actually established within the chapter's own pages.
Formalization scope
Every item fixes a finite decision space (Fin A, Fin n, or a generic Fintype S)
and states the DEC as the literal sInf-of-sSup transcription of the min-max game
(Eqs. (4.15)–(4.16)), never as an opaque bound — this is the trivializing
formalization the chunk's own reading of the chapter rules out (see Difficulty).
piStar : (S → ℝ) → S is a hypothesized global maximizer selector throughout,
constrained to be a genuine argmax only on the function class in scope (F or
Set.univ), matching how the book treats πf as a fixed but arbitrary
tie-breaking choice. The goal theorem (Proposition 13) adds the explicit hypothesis
hfhat : ∀ t, fhat t ∈ convexHull ℝ F, replacing an appeal to the out-of-range
Proposition 24 (see Difficulty); this is the one place this mission's statement is not
a line-by-line transcription of the book's own displayed proof steps, and it is
recorded here and in MODERATION_NOTES.md. Proposition 14's and Proposition 17's
≲ are replaced by the explicit constants the book's own proofs establish
((A−1)/(4γ) exactly, and (4d+1)/(2γ) respectively — the latter obtained
by summing the three terms the proof of Proposition 17 isolates). Proposition 14's
Lean statement splits the book's single equality decγ(F,f^)=(A−1)/(4γ) into an upper bound on the literal decGf, a lower bound restricted
to full-support distributions, and IGW's own exact game value, because the book's
min over the whole simplex is not provable as a literal Lean equality: a
distribution with a zero-weight arm makes the inner supremum genuinely unbounded, and
Lean's total Real.sSup returns a junk value smaller than (A−1)/(4γ) there
(caught in moderation, MODERATION_NOTES.md); the three-conjunct statement recovers
exactly the book's real content without asserting that false literal equality. Lemma 9
is restated
inside FoundationsRL.Structured rather than imported from the Chapter 2 mission,
since draft items across chunks cannot import one another; its source citation still
points to its original location (p. 32). Proposition 16 is not drafted: the
platform's existing BanditAlgorithm.kiefer_wolfowitz_equivalence (Lattimore &
Szepesvári, Theorem 21.1) states the identical existence claim — compact set with
full-dimensional span, a design with G-value at most d — as one clause of a larger
equivalence, and is reused as a reference item rather than redrafted. Proposition 22
(primal/dual DEC equivalence, §4.4) is deliberately excluded: the book states it "under
mild regularity conditions" it does not pin down in the statement itself, which is
exactly the kind of unquantified hypothesis this series' faithfulness standard
excludes from a goal or milestone. Contributions extending this mission with Chapter
6's lower bound (matching decγ(F) from below, establishing tightness)
or with a formalization of Proposition 24 itself (removing this mission's hfhat
hypothesis) are welcome.
Selected references
D. Foster, S. Kakade, J. Qian, and A. Rakhlin, The Statistical Complexity of Interactive Decision Making, arXiv:2112.13487, 2021. https://arxiv.org/abs/2112.13487
D. Foster and A. Rakhlin, Foundations of Reinforcement Learning and Interactive Decision Making, arXiv:2312.16730, 2023. https://arxiv.org/abs/2312.16730
T. Lattimore and C. Szepesvári, Bandit Algorithms, Cambridge University Press, 2020.
J. Kiefer and J. Wolfowitz, The Equivalence of Two Extremum Problems, Canadian Journal of Mathematics, 1960.
Bandit Algorithms XVI: Markov Decision Processes and UCRL2Textbook
The final step from bandits to reinforcement learning: actions now change the state of the world. Chapter 38 of Lattimore–Szepesvári studies online learning in an unknown Markov decision process with S states, A actions and rewards in [0,1]. The optimism principle of Mission II scales up: UCRL2 maintains confidence sets over transition kernels, solves an extended MDP by extended value iteration, and recomputes only when a state-action count doubles. The goal theorem: with probability 1−δ, R^n<CD(M)SAnlog(nSA/δ), where D(M) is the diameter of the MDP — sublinear regret with no prior knowledge of the dynamics. The matching lower bound E[R^n]≥C′DSAn brackets the true complexity of tabular reinforcement learning up to DS.
Bandit Algorithms VII: Lower Bounds for Finite-Armed BanditsTextbook
How well can any algorithm possibly do? Chapters 13–17 of Lattimore–Szepesvári answer with three matching impossibility results. The divergence decomposition identifies the information a policy collects: D(Pνπ,Pν′π)=∑iE[Ti(n)]D(Pi,Pi′). Feeding it into the Bretagnolle–Huber inequality of Mission VI yields the goal theorem — the minimax lower bound Rn≥271(k−1)n over Gaussian bandits, showing MOSS (Mission III) is optimal up to a constant. The same machinery gives the instance-dependent bound of Lai–Robbins type: every consistent policy suffers liminfnRn/logn≥∑i:Δi>0Δi/dinf(Pi,μ∗,Mi), certifying the asymptotic optimality of the UCB of Mission III and KL-UCB of Mission IV, and a high-probability lower bound showing the Exp3-IX guarantees of Mission V cannot be improved.
Contextual Inverse Optimization: Offline and Online Learning 3: EllipsoidalCones Has Worst-Case Regret O(d² ln(T tan α(C₀))) from a Pointed Initial Knowledge SetResearch Paper
Motivation
In inverse optimization one observes the decisions of an expert who solves an optimization problem and tries to recover the expert's objective. The contextual version of Besbes, Fonseca and Lobel (arXiv:2106.14015; Operations Research, 2025) asks a sharper question: not how well the cost can be estimated, but how well a decision-maker who imitates the expert performs, measured by regret against the expert. Settings of this kind arise whenever a firm wants to reproduce the routine choices of an experienced planner, for example routing, assortment or scheduling decisions, without being told the planner's cost trade-offs.
The online setting is adversarial. Nature fixes an unknown cost vector and then presents a sequence of problems; after each decision it reveals the expert's optimal action, which cuts the set of costs still consistent with the observations. A companion mission of this series shows that the natural greedy policy, which plays the center of this consistent set, can incur regret linear in the horizon. This mission formalizes the paper's first positive answer: an algorithm, EllipsoidalCones, with regret logarithmic in the horizon when the initial set of costs lies inside a pointed cone.
Setting
Costs lie on the unit sphere Sd of Rd, d=m+1≥2. In period t nature presents a feasible set Xt⊆Rn, compact, nonempty and of diameter at most 1, and a 1-Lipschitz contextft:Rn→Rd. The forward problem with cost c has solution set ψ(c,Xt,ft)=argminx∈Xtft(x)′c. The decision-maker plays some xtπ, nature reveals an expert action xt⋆∈ψ(c⋆,Xt,ft) for the true cost c⋆, and the regret over T periods is
RTπ=t=1∑T(ft(xtπ)−ft(xt⋆))′c⋆.
The knowledge setC(It) is the set of costs in the initial set C0⊆Sd under which every expert action seen so far was optimal.
The angle between nonzero vectors is θ(c,c^)=arccos(c′c^/∥c∥∥c^∥), and the uncertainty angle of a set C is α(C)=infc^∈Sdsupc∈Cθ(c,c^); a minimizing c^ is a circumcenter. The initial set is pointed when α(C0)<π/2.
For a positive diagonal W=diag(w1,…,wm) and an orthonormal U, the ellipsoidal cone is
a cone with axis Ue1 and elliptical cross-section. The effective difference of period t is δtπ=(ft(xtπ)−ft(xt⋆))/∥ft(xtπ)−ft(xt⋆)∥ (and 0 if the difference vanishes).
EllipsoidalCones with parameter ϵ starts from the circular cone W1=tan2α(C0)I, U1e1 a circumcenter of C0, and plays the axis Ute1 as proxy cost. With δˉ=Ut−1δtπ it keeps the cone if
δˉ[2:d]′Wtδˉ[2:d]≤ϵ2(12)
and otherwise replaces it by the output of ConeUpdate, which performs one central-cut step of the ellipsoid method on the cone's cross-section and re-centres the result by a spectral decomposition.
Formalization targets
Goal: Theorem 3 (p. 20)
With ϵ=d/T there is an absolute constant K such that for every d≥2, T≥1, closed pointed C0 with α(C0)>0, c⋆∈C0, admissible instances and run,
RTπ≤Kd2ln(Ttanα(C0)) if α(C0)>arctan(d/T),RTπ≤Kd otherwise.
The constant is left unfixed: the shape d2ln(Ttanα(C0)) is the claim, and the proof's explicit constant rests on a step that fails at d=2.
Milestones
Proposition 1:α(E(W,U))=arctanλmax(W).
Lemma 4: if C(It)⊆E(Wt,Ut) and (12) holds, the period's regret is at most ϵ.
Lemma B-1 (η=0): ConeUpdate's output contains E(W,U)∩{δ′c≥0} and ∏λi shrinks by e−1/(2(d−1)).
Lemma 3:C(It+1)⊆E(Wt,Ut)∩{δtπ′c≥0}⊆E(Wt+1,Ut+1), and for d≥3∏λi shrinks by e−1/(d−1).
Cohen et al. (2020), Lemmas 3 and 4 (already posed on the platform, referenced): an ellipsoid-method step shrinks the smallest eigenvalue by at most (d−1)2/d2, and does not shrink it at all below the threshold (ϵ/(10(d−1)))2 when the cut is long.
Lemma 5: for d≥3, λmax(Wt)≤(10d/ϵ)2(d−2)λmax(W1)d−1e−Itπ/(d−1), or no update ever happens.
Significance
The theorem shows that in contextual inverse optimization the decision-maker can force nature to reveal information: a cone update happens only after a large mistake, and each such mistake shrinks the cone geometrically, so large mistakes are logarithmically rare. It is the building block of the paper's general-case algorithm, which restarts the same construction in growing subspaces, and it isolates why the greedy policy fails: it never trades instantaneous regret for information.
The paper proves the result; nothing in it is machine-checked. The mission produces a checked version of the pointed-case analysis, an ellipsoidal-cone layer on top of the platform's ellipsoid-method definitions, and corrected statements where the printed ones are imprecise: the dimension restriction d≥3 for the rate of Lemma 3 and the first part of Lemma 5, the column order of the spectral decomposition in ConeUpdate, and the η=0 case of Lemma B-1.
Difficulty
Two steps resist the obvious argument. First, ConeUpdate is not the ellipsoid method: the minimum-volume ellipsoid of the cut cross-section has the right volume but its cone is no longer centred, and the new cone's axis is an eigenvector of an indefinite bordered matrix. Showing that this re-centred cone still contains the cut cone requires matching the eigenvalues of the bordered matrix with those of the ellipsoid matrix in the right order. Second, a volume decrease alone does not bound the regret: by Proposition 1 the uncertainty angle is governed by the largest eigenvalue, and the ellipsoid method can make that eigenvalue grow. The lazy test (12) and the smallest-eigenvalue lemmas of Cohen et al. are what keep the cones from becoming long and thin; without them the product of eigenvalues falls while the angle does not.
Formalization scope
Costs live in EuclideanSpace ℝ (Fin (m+1)) with d=m+1, actions in EuclideanSpace ℝ (Fin n), so all norms are Euclidean. W is stored as its diagonal w, U as a real orthogonal matrix, and U−1 is written U⊤. Lean periods are 0-based: Lean period t is the paper's period t+1. This text is the arXiv v3 preprint (1 July 2023), not the published 2025 article; all numbering is the preprint's.
Standing assumptions carried by the statements: d≥2, n≥1, feasible sets nonempty, compact and of diameter at most 1, contexts 1-Lipschitz, C0 a closed subset of Sd, ϵ>0. Pinned and disclosed choices:
ConeUpdate is a relation, not a function. The printed U−1Bδ[2:p] is read as (U−1δ)[2:p]. The first column of V is the eigenvector of the positive eigenvalue of M with positive first entry, and larger ∣μ∣ is paired with larger eigenvalues of N. The page leaves this order open, and Algorithm 1 and Lemma B-1 disagree on it.
The printed "U1←Id−1" is read as: U1 orthonormal with a circumcenter of C0 as its first column.
α(C0)>0 on the goal. At α(C0)=0 the algorithm's W1=0 is not positive definite.
The rate of Lemma 3 and the first part of Lemma 5 are stated for d≥3 only; the rest of both lemmas for every d≥2. At d=2 ConeUpdate halves W, and 1/2>e−1.
Lemma B-1 is stated only for η=0. Its η>0 inclusion is false as printed.
The worst-case regret is stated as a bound for every true cost, every adaptive instance sequence and every run, not as a real supremum. The polynomial running-time claim is not formalized.
The run predicate is satisfiable: a circumcenter of a closed nonempty C0 exists, the bordered matrix has exactly one positive eigenvalue, and a cone update is required only when the tail of δˉ is nonzero. A formalization in which no run exists, or in which K may depend on d or T, does not state Theorem 3.
A complete development needs the angle and uncertainty-angle layer, ellipsoidal cones and their cross-sections, the published ellipsoid-method update (LinearOptimization_EllipsoidMethod) with its volume and eigenvalue facts, Cauchy–Schwarz over ellipsoids, and eigenvalue interlacing for bordered matrices. The ellipsoidal-cone layer and the cone-update lemma can be reused in the paper's general case and in other cone-based learning algorithms. Proofs of any milestone, and of the referenced Cohen et al. lemmas, are welcome.
Selected references
O. Besbes, Y. Fonseca, I. Lobel, Contextual Inverse Optimization: Offline and Online Learning, arXiv:2106.14015v3, 2023; Operations Research, 2025. https://arxiv.org/abs/2106.14015
Fast Convergence of Regularized Learning in Games II: If Every Player Runs Optimistic FTRL with Step Size (n−1)^{−1/2}T^{−1/4}, Each Player's Regret Is at Most (R+4)√(n−1)·T^{1/4}Research Paper
Motivation
When every participant in a repeated game runs a no-regret learning algorithm, the time-averaged play converges to a coarse correlated equilibrium, and in smooth games the average welfare approaches the price-of-anarchy guarantee. These facts are a standard justification for equilibrium predictions in auctions, routing and markets, where agents learn rather than solve for equilibria. Against an adversary, the regret of any learning algorithm can grow like T, which gives average regret of order T−1/2. Syrgkanis, Agarwal, Luo and Schapire (arXiv:1507.00407, NIPS 2015) asked whether learning is faster when the opponents are themselves learning, and so their play changes slowly.
Timeline. Daskalakis, Deckelbaum and Kim (2011) gave a decentralized algorithm with O(logT) regret in two-player zero-sum games. Rakhlin and Sridharan (2013) obtained the same rate in zero-sum games with optimistic mirror descent, which uses the last observed utility as a prediction of the next one. Syrgkanis et al. (2015) extended recency-biased learning to general multi-player normal-form games. They isolated the RVU property, proved that it gives constant total regret (Theorem 4, Corollary 8), and proved that each individual player's regret is O(T1/4) (Theorem 11, Corollary 12). This mission formalizes the second result. Later work by Chen and Peng (2020) and by Daskalakis, Fishelson and Golowich (2021) improved the individual rate to polylogarithmic in T.
Setting
A static game has n players. Each player i has a finite strategy set Si with ∣Si∣=d, and a utility function ui:S1×⋯×Sn→[0,1]. A mixed strategy is a point wi of the probability simplex Δ(Si)⊂Rd, and the players randomize independently. The game is played at times t=1,…,T. At time t each player picks wit∈Δ(Si) and then observes its utility vector
uit=(Es−i∼w−it[ui(x,s−i)])x∈Si,
so its expected utility is ⟨wit,uit⟩. The regret of player i after T steps is
ri(T)=wi∗∈Δ(Si)supt=1∑T⟨wi∗−wit,uit⟩.
Strategies are measured in the norm ∥⋅∥=∥⋅∥1. Utilities are measured in its dual, ∥v∥∗=∥v∥∞. A regularizerR is 1-strongly convex if R(2a+b)≤2R(a)+R(b)−8∥a−b∥12. The rangeR is maxi(supΔRi−infΔRi).
Optimistic follow the regularized leader (OFTRL) with step size η>0 and predictor sequence Mit plays wi0=argminΔR and
wiT=w∈Δ(Si)argmax⟨w,t=1∑T−1uit+MiT⟩−ηR(w).
With one-step recency bias, MiT=uiT−1. An algorithm has the RVU property (Regret bounded by Variation in Utilities) with parameters α,β,γ if its regret against every w∗ is at most α+β∑t∥ut−ut−1∥∗2−γ∑t∥wt−wt−1∥2. It has the stability property with constant κ if ∥wit−wit+1∥≤κ.
Formalization targets
Goal: Corollary 12 (p. 5)
If n≥2, every player runs OFTRL with Mit=uit−1 and η=(n−1)−1/2T−1/4, then for every player i and every wi∗∈Δ(Si),
t=1∑T⟨wi∗−wit,uit⟩≤(R+4)n−1⋅T1/4.
Milestones
Lemma 20 (supp. p. 4): ∥wt−gt∥≤η∥Mt−ut∥∗ and ∥gt−wt+1∥≤η∥Mt+1∥∗, where gT=argmaxΔ⟨g,∑t≤Tut⟩−R(g)/η is the leader sequence.
Stability with κ=2η (p. 5, before Corollary 12): OFTRL with Mt=ut−1 and utilities in [0,1] moves by at most 2η per step.
Proposition 7 (p. 5): OFTRL with Mt=ut−1 has the RVU property with α=R/η, β=η, γ=1/(4η).
The display in the proof of Theorem 11 (p. 5): ∥uit−uit−1∥∗2≤(n−1)∑j=i∥wjt−wjt−1∥2≤(n−1)2κ2.
Theorem 11 (p. 5): RVU with α>0,β>0,γ≥0 plus stability κ give each player regret at most α+βκ2(n−1)2T.
Significance
The result. Corollary 12 shows that recency-biased regularized learning in a game guarantees each player average regret O(T−3/4). The adversarial optimum is O(T−1/2). The bound covers every player separately and needs no coordination beyond a common step size. Corollary 8, by contrast, bounds only the sum of the regrets. Theorem 11 isolates the two algorithmic properties responsible, RVU and stability. Any algorithm with both, including optimistic mirror descent and the other recency biases of the paper, inherits the bound.
Formalizing it. The result is proved in the paper, with the proofs of Proposition 7 and Lemma 20 in its supplementary material. No machine-checked version of optimistic FTRL, of the RVU property or of fast rates in games appears to exist. A formalization makes the following explicit: the strong-convexity form the proofs actually use, the convention for the first utility vector u0, and the exact range of times on which stability is needed. A later mission of this series formalizes the companion bound on total regret (Corollary 8) on the same objects.
Difficulty
The game-theoretic part is short. Once Proposition 7 and the stability bound are available, Theorem 11 is a summation and Corollary 12 is arithmetic with real powers. The weight lies in the single-player analysis of optimistic FTRL. The first idea, to compare consecutive argmax points through first-order optimality conditions, fails because the regularizer need not be differentiable and the maximizers can lie on the boundary of the simplex. The arguments must work with the optimality of argmax points over Δ and with a quadratic-growth consequence of strong convexity. That consequence is not immediate from the midpoint condition in footnote 3. It needs continuity of the regularizer to pass from midpoints to the whole segment. Bounding the variation of utility vectors by the movement of the other players needs a total-variation estimate for product distributions.
Formalization scope
Players are Fin n and strategies Fin d, both indexed from 0. Time keeps the paper's indices: w0 is the initial play and sums run over t∈{1,…,T}. The vector ui0 is the utility vector against w0; the paper never defines it. Mixed profiles and expectations use the published agt_games definitions. ∥⋅∥ is ℓ1 and ∥⋅∥∗ is Mathlib's sup norm. Optimistic FTRL and the leader sequence are predicates on trajectories, so every maximizer is a valid run. The goal quantifies over every run of the dynamics.
Standing assumptions and disclosed additions:
utilities lie in [0,1] (§2);
each player may have its own regularizer (footnote 2);
strong convexity is footnote 3's midpoint condition on Δ plus continuity on Δ;
R is any common upper bound on the ranges, which includes the exact range;
n≥2 and T≥1, since otherwise the step size is undefined;
the regret's supremum is expressed as a bound for every comparator w∗∈Δ;
Theorem 11 assumes RVU only on the realized play, which strengthens the theorem, and stability only for the steps t=0,…,T−1.
Trivializations are ruled out. The regularizer class is nonempty: (d/2)∑xwx2 satisfies it. Removing the [0,1] bound on utilities would make the stability claim false. The goal fixes the paper's step size and constant, so it is not a statement about a chosen run, an unbounded supremum or a vacuous parameter range.
A complete development needs the following:
quadratic growth of strongly convex functions at constrained maximizers, for ℓ1 on the simplex;
the be-the-leader induction of Theorem 19;
a total-variation bound for product measures on finite sets.
These are reusable well beyond this mission, for online learning and for learning in games, and contributions of each are welcome.
Selected references
V. Syrgkanis, A. Agarwal, H. Luo, R. E. Schapire, Fast Convergence of Regularized Learning in Games, NIPS 2015; arXiv:1507.00407v5. https://arxiv.org/abs/1507.00407
A. Rakhlin, K. Sridharan, Optimization, Learning, and Games with Predictable Sequences, NIPS 2013. https://arxiv.org/abs/1311.1869
C. Daskalakis, A. Deckelbaum, A. Kim, Near-Optimal No-Regret Algorithms for Zero-Sum Games, SODA 2011. https://arxiv.org/abs/1011.1489
X. Chen, B. Peng, Hedging in Games: Faster Convergence of External and Swap Regrets, NeurIPS 2020. https://arxiv.org/abs/2006.04953
C. Daskalakis, M. Fishelson, N. Golowich, Near-Optimal No-Regret Learning in General Games, NeurIPS 2021. https://arxiv.org/abs/2108.06924
Fast Convergence of Regularized Learning in Games III: In a (λ,µ)-Smooth Cost Game, First-Order Regret Bounds Make the Average Cost at Most λ(1+µ)/(µ(1−µ))·OPT′ + A·n log d/TResearch Paper
Motivation
When independent players repeatedly choose actions, a low-regret guarantee controls each player's performance against a fixed alternative. It does not directly control the social cost, the sum of all players' costs. For a class of games satisfying smoothness, however, the individual comparisons can imply a bound on the cost of the joint play. This connection makes learning guarantees useful for studying decentralized systems in which nobody selects the entire outcome. The paper uses the smooth-game framework for this purpose.
Syrgkanis, Agarwal, Luo, and Schapire study regularized learning in games and, in Appendix H, show a separate consequence for first-order regret in cost-minimization games. A first-order bound scales with the comparator's cumulative cost, rather than only with the number of rounds. Their Theorem 23 turns that guarantee into an average-cost estimate whose additive term falls as 1/T. The result is proved in the paper's supplementary material. This mission targets a machine-checked statement and proof of that known result.
Setting
There are n players. Each player has d pure strategies; a pure profiles=(si)i chooses one for each player. Player i incurs a cost ci(s)∈[0,1]. The total cost is C(s)=∑ici(s), and OPT′=minsC(s) is the minimum over pure profiles. A mixed profilew=(wi)i gives each player a probability distribution over its d strategies. Players randomize independently, and C(w)=Es∼w[C(s)] is the expected total cost.
In round t, player i uses wit and sees a vector cit. Its coordinate ci,xt is the expected cost of playing pure strategy x against the other players' round-t mixed strategies. Thus ⟨wit,cit⟩ is that player's actual expected cost. Equation (21) assumes that, for every fixed pure strategy x,
The constants A1,A2 are the constants of the learning guarantee. The logarithm is natural.
A game is (λ,μ)-smooth when one pure profile s∗ works against every pure profile s:
i∑ci(si∗,s−i)≤λOPT′+μC(s).
Here (si∗,s−i) replaces only player i's strategy in s. This is the paper's cost-game smoothness condition (20). It concerns a single witness s∗, rather than a different witness chosen for each s.
Formalization targets
The goal is Appendix H's Theorem 23. For 0<μ<1 and T≥1, if the game is (λ,μ)-smooth and every player's realized play obeys (21), then
The constants are exactly those printed in the paper. The milestone list follows the displayed statements in its proof: the identity expressing social cost through cost-vector inner products, the aggregate first-order inequality (22), the expected smoothness inequality, the quadratic inequality for the auxiliary quantity x, and the final root bound. Together they describe the intermediate claims available to solvers without changing the theorem's target.
Significance
The theorem says that the joint play of agents with individual first-order regret guarantees approaches an explicit benchmark in average social cost. The 1/T term separates the learning transient from the benchmark determined by λ, μ, and OPT′. It is a statement about arbitrary realized trajectories that meet the regret condition; the players need not coordinate on a common learning rule.
The paper already proves Theorem 23. The remaining work here is to formalize its cost model, the product-distribution expectations, the first-order bound, the smoothness transfer, and the real inequality behind the exact constants. A proof of this mission would give reusable finite-game definitions and a checked route from player-level cost comparisons to a social-cost statement. Its intermediate claims may also be useful when an alternative algorithm supplies the same first-order guarantee.
Difficulty
Summing the individual regret bounds leaves a sum of square roots of comparator costs. The cost benchmark in the conclusion is a single quantity, so the player-level terms cannot simply be replaced one at a time by OPT′. The cost-vector identity and smoothness condition must also agree on the same independent product distribution: changing only one player's pure action while keeping the opponents mixed is easy to encode with a different, unintended distribution. After those choices, the exact coefficient A depends on the signs and domains in a quadratic bound. In particular, the intermediate Cauchy–Schwarz inequality uses A1≥0, while the final theorem has no sign restriction on A1.
Formalization scope
Players are Fin n, and each player's strategies are Fin d. The condition d≥1 makes the pure optimum a minimum over a nonempty finite set. Mixed profiles are probability weights on each strategy set, with the players' choices independent. The published agt_games definition supplies product-profile probabilities and finite expectations; its expectedPayoff operation is used as an expectation of costs. The paper-specific layer defines pure and expected social cost, the pure optimum, cost vectors, smoothness, and the realized first-order regret condition.
Time has the paper's indices: every round sum runs from t=1 through T, without shifting the sequence. The goal uses T≥1 and 0<μ<1 so its divisions represent the paper's formula. Every cost is in [0,1]; every played profile is mixed; (21) is required for every player and every fixed pure comparator. The formulation does not restrict the game to a chosen example, assume the goal inequality, replace OPT′ with a default-valued infimum, or allow an empty strategy set. The first-order condition is stated on the realized sequence, which follows from any algorithm that guarantees (21) for all sequences.
The intermediate inequality (22) and the quadratic milestone explicitly require A1≥0, the sign used in their Cauchy–Schwarz step. Theorem 23 itself allows any real A1,A2,λ, as the paper's stated formula does. The source prints ∑T=1T in (22); the Lean statement uses the intended round index t. The source also switches between logd and lnd; both denote the natural logarithm. Contributions on finite product expectations, the mixed smoothness consequence, and the final real inequalities are all within scope.
Selected references
Vasilis Syrgkanis, Alekh Agarwal, Haipeng Luo, and Robert E. Schapire, Fast Convergence of Regularized Learning in Games, Advances in Neural Information Processing Systems 28, 2015. arXiv:1507.00407v5, Appendix H, supplementary pp. 10–11.
Regularization via Mass Transportation II: Wasserstein-Robust Linear Regression Equals Empirical Loss plus a Dual-Norm PenaltyResearch Paper
Motivation
Ordinary least squares and its robust relatives (least absolute deviations, Huber, ϵ-insensitive and quantile regression) fit a linear predictor by minimizing an average loss over the training sample. That average estimates the out-of-sample error poorly when the sample is small, and the usual fix is to add a regularization penalty such as a norm of the weights. Choosing that penalty is often a heuristic. Distributionally robust optimization gives one principled derivation of it: hedge against every data-generating distribution that is close to the empirical one in a Wasserstein distance, and minimize the worst-case expected loss.
S. Shafieezadeh-Abadeh, D. Kuhn and P. Mohajerin Esfahani, Regularization via Mass Transportation (arXiv:1710.10016v3, 2019; JMLR 20, 2019), show that for linear regression with a convex Lipschitz loss this worst-case problem is exactly a regularized empirical risk problem, with an explicit dual-norm penalty. The same identity appears for the logistic loss in Shafieezadeh-Abadeh, Mohajerin Esfahani and Kuhn, Distributionally Robust Logistic Regression (NeurIPS 2015, arXiv:1509.09259). For a general convex loss it is Theorem 10 of the tutorial by Kuhn, Mohajerin Esfahani, Nguyen and Shafieezadeh-Abadeh (arXiv:1908.08729). Blanchet, Kang and Murthy (arXiv:1610.05627) obtain the square-root Lasso from a type-2 transport cost. This mission formalizes the regression theorem of the first paper, Theorem 3.1(ii) on p. 8. The source PDF is the arXiv preprint v3, and every page and label cited here refers to it.
Setting
Data. Inputs are x∈Rn and outputs are y∈R. A pair is written ξ=(x,y)∈Rn+1. The training samples are (x^i,y^i) for i=1,…,N, and P^N=N1∑i=1Nδ(x^i,y^i) is their empirical distribution (§2.1, p. 5).
Hypotheses and loss. A weight vector w∈Rn defines the linear hypothesis h(x)=⟨w,x⟩. A univariate convex loss L:R→R defines the regression loss ℓ(⟨w,x⟩,y)=L(⟨w,x⟩−y) (§3.1, p. 8). The residual ⟨w,x⟩−y is the pairing of ξ with the vector (w,−1).
Norms.Rn+1 carries an arbitrary norm ∥⋅∥ with dual norm∥z∥∗=sup{⟨z,ξ⟩:∥ξ∥≤1} (§1.1, p. 4). The Lipschitz modulus of L is lip(L)=supz=z′∣L(z)−L(z′)∣/∣z−z′∣, and the convex conjugate is L∗(θ)=supzθz−L(z) (§1.1).
Wasserstein ball. With transport cost d(ξ,ξ′)=∥ξ−ξ′∥, the Wasserstein distance (Definition 2.2, p. 7) is
W(Q,Q′)=Πinf∫d(ξ,ξ′)Π(dξ,dξ′),
where the infimum runs over couplings Π of Q and Q′. The ball of radius ρ≥0 is Bρ(P^N)={Q:Q(Ξ)=1,W(Q,P^N)≤ρ} (3). The distributionally robust regression problem is
winfQ∈Bρ(P^N)supEQ[L(⟨w,x⟩−y)].(4)
Formalization targets
Goal: Theorem 3.1(ii), p. 8
Let Ξ=Rn+1, ρ≥0, N≥1, and let L be nonnegative, convex and Lipschitz continuous. Then for every w
so (4) and the regularized problem (6) have the same optimal value.
Milestones
Lemma A.1 (p. 28). The worst case over the ball equals infλ≥0λρ+N1∑isupξ[I(ξ)−λd(ξ,ξ^i)]. This milestone is linked to the published, proved WassersteinDRO.Duality.strong_duality_worst_case_risk_v2, which covers it for ρ>0.
Lemma A.3 (p. 29). For γ>0, supζL(⟨β,ζ⟩)−γ∥ζ−ζ^∥ equals L(⟨β,ζ^⟩) if lip(L)∥β∥∗≤γ, and +∞ otherwise.
End of the proof of Lemma A.3 (p. 30). lip(L)=sup{∣θ∣:L∗(θ)<∞}.
The inf-over-λ display in the proof of Theorem 3.1(ii) (pp. 30–31). The worst case equals inf{λρ+N1∑iL(⟨w,x^i⟩−y^i):λ≥0,lip(L)∥(w,−1)∥∗≤λ}.
Significance
The result. The theorem turns an infinite-dimensional min–max problem over distributions into a finite convex program over w: an empirical average plus a norm penalty. The penalty is not ∥w∥∗ but ∥(w,−1)∥∗, because the transport cost moves outputs as well as inputs. It scales with the radius ρ and with the Lipschitz modulus of the univariate loss. With the Huber, ϵ-insensitive and pinball losses it gives Corollaries 3.2–3.4 of the paper. It also gives a statistical reading of the regularization weight as a Wasserstein radius, which the paper's §4 turns into generalization bounds.
Formalizing it. The theorem is proved in the paper; none of its items is formalized yet except Lemma A.1. Lemma A.1 is available in a general form (any p≥1, any closed support set, ρ>0) as a proved platform theorem from the mission Wasserstein Distributionally Robust Optimization. The per-w identity also follows from the proved WassersteinDRO.Duality.convex_loss_p1_v2 (Theorem 10 of the tutorial, also present as WassersteinDRO.Regularization.convex_loss_p1_v2; the public mission Wasserstein Distributionally Robust Optimization IV has it as goal) once one shows that lip(L∘⟨(w,−1),⋅⟩)=lip(L)∥(w,−1)∥∗. No platform item states that factorization. What this mission adds is the paper's explicit form of the penalty, Lemma A.3 (the inner supremum for an arbitrary convex Lipschitz loss and an arbitrary norm), and the conjugate description of the Lipschitz modulus. The last two are reusable in classification (mission I of this series states them again in its own namespace) and in every Wasserstein-regularization result with linear hypotheses.
Difficulty
The obvious argument is to write the worst case as a supremum over perturbations of the samples and bound each perturbed loss by its Lipschitz constant. That gives only the inequality "≤", and only for one Lipschitz constant. The identity needs two further steps that this bound does not supply. First, the bound must be attained: the supremum over distributions has to be shown to reach ρlip(L)∥(w,−1)∥∗, a matching lower bound that a Lipschitz estimate cannot give. Second, lip(L) must be the modulus, not any constant: an upper Lipschitz constant K>lip(L) makes the identity false. In Lemma A.3 the case lip(L)∥β∥∗=γ sits exactly on the boundary between a finite and an infinite supremum. A further point is that the duality of Lemma A.1 needs ρ>0, while the theorem includes ρ=0.
Formalization scope
Space.Rn+1 with an arbitrary norm is an abstract finite-dimensional real normed space E with its Borel σ-algebra. A continuous linear isomorphism e : E ≃L[ℝ] (Fin n → ℝ) × ℝ reads off inputs and outputs. The theorem therefore covers every norm on Rn+1; fixing a concrete norm would prove it for that norm only. The functional (w,−1) is pairing e w, and ∥(w,−1)∥∗ is its operator norm.
Published objects. The ball, the worst case, the empirical distribution and the Lipschitz modulus are the published definitions WassersteinDRO_Duality_ambiguitySet, _worstCaseRisk_v2, _empiricalDistribution and _lipschitzModulus. The worst case is taken in the extended reals, and expectations follow the convention of nominalRisk_v2 (no junk value for non-integrable losses). lip(L) takes values in [0,∞].
Standing assumptions.L is nonnegative, convex on R and Lipschitz continuous; N≥1; ρ≥0; Ξ=Rn+1. The paper's L≥0 (§2.1) is imposed as the standing loss-function assumption.
Reading of "equivalent". The paper does not define the word. The goal states the per-w identity that the proof establishes plus equality of optimal values. Attainment of the infimum is not claimed.
Corrections. The printed display on p. 30 omits a factor N1 and carries a stray ∣θ∣ in the constraint; milestone 4 states the corrected identity. The printed Lemma A.1 omits ρ>0, which its proof uses; the linked platform theorem has it.
Ruled out. The goal is not a restatement of Theorem 10 with Lip(ℓw) of the composite loss: its right-hand side is the explicit average of (6) with lip of the univariate L and the dual norm of (w,−1).
Welcome contributions. The factorization lip(L∘β)=lip(L)∥β∥∗ for a linear functional β on a finite-dimensional normed space, conjugate calculus for univariate convex functions, and the evaluation of the worst case at ρ=0.
Selected references
S. Shafieezadeh-Abadeh, D. Kuhn, P. Mohajerin Esfahani, Regularization via Mass Transportation, arXiv:1710.10016v3, 2019; J. Mach. Learn. Res. 20(103), 2019. https://arxiv.org/abs/1710.10016
D. Kuhn, P. Mohajerin Esfahani, V. A. Nguyen, S. Shafieezadeh-Abadeh, Wasserstein Distributionally Robust Optimization: Theory and Applications in Machine Learning, INFORMS TutORials, 2019. https://arxiv.org/abs/1908.08729
P. Mohajerin Esfahani, D. Kuhn, Data-driven distributionally robust optimization using the Wasserstein metric, Math. Program. 171, 2018. https://arxiv.org/abs/1505.05116
S. Shafieezadeh-Abadeh, P. Mohajerin Esfahani, D. Kuhn, Distributionally Robust Logistic Regression, NeurIPS, 2015. https://arxiv.org/abs/1509.09259
J. Blanchet, Y. Kang, K. Murthy, Robust Wasserstein Profile Inference and Applications to Machine Learning, J. Appl. Probab. 56(3), 2019. https://arxiv.org/abs/1610.05627
D. P. Bertsekas, Convex Optimization Theory, Athena Scientific, 2009 (cited as [9] in the paper for Propositions 5.3.1 and 5.5.4).
Online Learning and Online Convex Optimization 4: Finite Littlestone Dimension Gives Agnostic Online Regret 2√(Ldim·ln(eT/Ldim)·T)Research Paper
Motivation
Online classification is the simplest model of learning from a stream: on every round a learner sees an instance, guesses its binary label, and is then told the truth. Cover's impossibility argument shows that against an adversary no deterministic learner can compete with even two fixed hypotheses, so two relaxations are studied. In the realizable case the labels are generated by some unknown hypothesis of a known class H, and the question is how many mistakes a learner must make. In the agnostic (unrealizable) case the labels are arbitrary, the learner may randomize, and it is judged by its regret: its expected number of mistakes minus that of the best hypothesis in H in hindsight.
For the realizable case Littlestone (1988) introduced the dimension now named after him and showed that it is exactly the optimal mistake bound. Whether the same quantity also governs the agnostic case was settled by Ben-David, Pál and Shalev-Shwartz (COLT 2009, link): a class is agnostically online learnable if and only if its Littlestone dimension is finite. This mission formalizes the upper-bound half of that result, in the form given by Shalev-Shwartz's survey Online Learning and Online Convex Optimization (Found. Trends Mach. Learn. 4(2), 2011, doi:10.1561/2200000018), Theorem 3.6.
Timeline:
1988 — Littlestone defines the dimension and the Standard Optimal Algorithm (SOA) for the realizable case (doi:10.1007/BF00116827).
1994 — Littlestone and Warmuth publish the Weighted Majority algorithm for prediction with expert advice (doi:10.1006/inco.1994.1009).
2009 — Ben-David, Pál, Shalev-Shwartz: agnostic online learnability is characterized by finite Littlestone dimension, regret O(Ldim(H)TlnT), with a lower bound Ω(Ldim(H)T).
2011 — the survey presents the expert construction and the upper bound in the form formalized here.
Setting
Let X be a set of instances and let a hypothesis classH be a set of functions h:X→{0,1}. On round t=1,…,T the learner receives xt, outputs pt∈[0,1] (the probability of predicting 1), receives yt∈{0,1} and pays ∣pt−yt∣, the probability of a wrong randomized guess. The prediction pt may depend on (x1,y1),…,(xt−1,yt−1) and xt, never on yt. The regret against h∈H is
t=1∑T∣pt−yt∣−t=1∑T∣h(xt)−yt∣.
A shattered tree of depth d assigns an instance to every node of a complete binary tree of depth d, so that for every labelling (y1,…,yd) some h∈H agrees with yt at the t-th node of the root-to-leaf path that the labels select. The Littlestone dimensionLdim(H) is the largest depth of a shattered tree.
The Standard Optimal Algorithm keeps the version space Vt of hypotheses consistent so far, splits it as Vt(r)={h∈Vt:h(xt)=r}, and predicts the r with the larger Ldim(Vt(r)). For a set I={i1<⋯<iL} of rounds, Expert(i1,…,iL) runs the same rule but flips its prediction on the rounds of I and updates Vt+1=Vt(y^t) using its own prediction instead of the label. Weighted Majority over d experts with costs zt∈[0,1]d keeps weights wt[i]∝exp(−η∑s<tzs[i]) and pays ⟨wt,zt⟩.
Formalization targets
Goal: Theorem 3.6, upper bound
For every class H with L=Ldim(H)<∞ and every horizon T≥L there is an online learner with pt∈[0,1] such that for every h∈H and every sequence of T examples
t=1∑T∣pt−yt∣−t=1∑T∣h(xt)−yt∣≤2Lln(eT/L)T.
The paper states O(Ldim(H)ln(T)T); the constant 2 and the factor ln(eT/L) are what its proof gives.
Milestones
Lemma 3.3: SOA makes at most Ldim(H) mistakes in the realizable case (already proved on the platform).
Lemma 3.7: for every h∈H and instance sequence some Expert(i1,…,iL) with L≤Ldim(H) predicts h(xt) on every round.
Corollary 3.8: some such expert makes at most minh∈H∑t∣h(xt)−yt∣ mistakes.
Theorem 3.1: Weighted Majority has regret at most log(d)/η+ηT, and 2log(d)T for η=log(d)/T.
Significance
The theorem shows that finiteness of one combinatorial parameter suffices for sublinear agnostic regret, for classes that may be infinite and need not have any parametric structure. With the lower bound it characterizes agnostic online learnability, parallel to the role of the VC dimension in PAC learning; the gap between the two dimensions (thresholds on [0,1] have VC dimension 1 and infinite Littlestone dimension) separates online from batch learnability. The expert construction, which simulates a realizable-case learner on guessed mistake rounds, is a reusable reduction from agnostic to realizable online learning.
On the formal side, the realizable theory (Littlestone dimension, SOA and its mistake bound) is already on the platform as published definitions and theorems from the textbook Understanding Machine Learning; the agnostic step, the expert construction, the bound on the number of experts and the Weighted Majority bound with bounded costs are not. To our knowledge no proof assistant library contains a machine-checked proof of the agnostic upper bound.
Difficulty
The obvious approach, running Weighted Majority with one expert per hypothesis, fails for infinite H and gives log∣H∣ instead of Ldim(H) for finite H. The class must be replaced by a finite family that reproduces every h∈H on the given instance sequence while having only poly(T)Ldim(H) members. The experts in the family never see labels, so the argument that each h is reproduced must track two runs, the expert's and SOA's on labels generated by h, and show their version spaces coincide round by round. This depends on the tie-breaking rule and on how an empty side Vt(r) is ranked, details the informal argument leaves implicit. The bound on the number of experts is a Sauer-type binomial estimate that holds only in the range 1≤Ldim(H)≤T.
Formalization scope
The development reuses the published definitions UnderstandingML_Online (online algorithms as maps from the history and the current instance to a prediction, history, mistakes, cumLoss, cumLossHyp, ldim, ldimBot, soa, wmWeights) and the proved theorem UnderstandingML.soa_mistake_bound (Lemma 3.3). Sequences are indexed by Fin T and rounds from 0; the paper's indices 1≤i1<⋯<iL≤T become a set I⊆{0,…,T−1}. A learner of type OnlineAlgR X cannot read the current label, so the existential goal cannot be satisfied by a learner that peeks at yt; the goal also requires pt∈[0,1]. Logarithms are natural.
Committed conventions:
The O(·) of Theorem 3.6 is instantiated as 2Lln(eT/L)T under T≥L, where L=Ldim(H) is a natural number. At L=0 the bound reads regret ≤0.
Expert(i1,…,iL) breaks ties toward 0, as its box on p. 166 says; the SOA box (p. 163) breaks ties toward 1. Lemma 3.7 holds with the expert's rule. An empty Vt(r) ranks below every nonempty class (ldimBot).
Theorem 3.1's bound 1 is stated for every η>0, bound 2 for η∈(0,1), the tuned bound for every d≥1, T≥0; minima over experts and over H are stated as "for every expert", "for every h".
Equation (3.1) is stated with 1≤L≤T, together with the count of index sets.
Excluded: the lower bound Ω(Ldim(H)T) (the survey cites it to Ben-David, Pál and Shalev-Shwartz), Theorem 3.5 (VC dimension at most Littlestone dimension), Corollary 3.4. Lemma 3.2 and the Halving bound are already on the platform.
Welcome contributions: proofs of the milestones, a proof of SOA's mistake bound with ties toward 0, and a reusable version of the Weighted Majority bound for general finite expert types.
Selected references
S. Shalev-Shwartz, Online Learning and Online Convex Optimization, Foundations and Trends in Machine Learning 4(2), 2011, 107–194. https://doi.org/10.1561/2200000018
N. Littlestone, Learning Quickly When Irrelevant Attributes Abound: A New Linear-Threshold Algorithm, Machine Learning 2, 1988, 285–318. https://doi.org/10.1007/BF00116827
N. Littlestone, M. K. Warmuth, The Weighted Majority Algorithm, Information and Computation 108(2), 1994, 212–261. https://doi.org/10.1006/inco.1994.1009
What Can We Learn Privately? I: The Exponential Mechanism Privately Agnostically Learns Any Finite Concept Class from O((ln|C| + ln(1/β))·max{1/(εα), 1/α²}) ExamplesResearch Paper
Motivation
Many learning tasks run on data about people: medical records, purchase histories, survey answers. A learner whose output reveals too much about any one training example is unusable on such data, however accurate it is. Differential privacy (Dwork, McSherry, Nissim and Smith, 2006) makes "reveals too much" precise: changing one individual's record may change the probability of any outcome of the algorithm by at most a factor eε.
Kasiviswanathan, Lee, Nissim, Raskhodnikova and Smith asked which concept classes can be learned under this constraint, and at what cost in examples. This mission formalizes their first answer (§3 of the paper): every finite concept class can be learned privately in the agnostic model, with a sample size that grows only logarithmically in the size of the class. The result is the private counterpart of the cardinality version of Occam's razor (Kearns and Vazirani, An Introduction to Computational Learning Theory, 1994, p. 35). The paper appeared in SIAM J. Comput. 40(3) (2011) 793–826, DOI 10.1137/090756090; all indices below are those of arXiv:0803.0924v3.
The learner is an instance of the exponential mechanism of McSherry and Talwar (FOCS 2007), which selects an output with probability exponentially weighted by a quality score.
Setting
A database of size n over a domain D is a vector z=(z1,…,zn)∈Dn. Databases z,z′ are neighbors if they differ in exactly one entry. A randomized algorithm A is described by the law A(z) of its output on each input. It is ε-differentially private (Definition 2.1) if
Pr[A(z)∈S]≤eεPr[A(z′)∈S]
for all neighbors z,z′ and all measurable sets S of outputs.
Learning uses a finite set X of examples and labels {0,1}. A hypothesis is a map h:X→{0,1}, and H is a finite nonempty class of hypotheses; in the theorem H=Cd, the concept class itself. Data are drawn i.i.d. from an arbitrary probability distribution P on X×{0,1} (the agnostic setting, Definition 2.5). The error of h is err(h)=Pr(x,y)∼P[h(x)=y], and OPT=minf∈Herr(f). On a database z=((xi,yi))i≤n the score is q(z,h)=−∣{i:yi=h(xi)}∣ and the training error is errT(h)=−q(z,h)/n.
The exponential mechanismAqε outputs h∈H with probability proportional to exp(εq(z,h)/2). In Lean these objects are IsDP, err, OPT, score, errT and expMech in the namespace PrivLearn.Generic.
Formalization targets
Goal: Theorem 3.4 (generic private learner)
There is an absolute constant C>0 such that for every finite X, every finite nonempty H, every ε>0 and all α,β∈(0,1/2):
Aqε is ε-differentially private on databases of every size;
for every distribution P on X×{0,1} and every
n≥C(ln∣H∣+lnβ1)max{εα1,α21},
the mechanism run on n i.i.d. samples from P satisfies
Pr[err(Aqε(z))>OPT+α]≤β.
The constant is left unspecified, as the paper's O(⋅) leaves it. The qualitative first sentence of Theorem 3.4 (classes of size exp(poly(d)) are privately agnostically learnable) follows at once, because then ln∣H∣ is polynomial in d.
Milestones
Lemma 3.3: Aqε is ε-differentially private.
Theorem A.2: for i.i.d. Xi∈[a,b] with mean μ, Pr[∣n1∑iXi−μ∣≥δ]≤2exp(−2δ2n/(b−a)2).
Uniform convergence (proof of Theorem 3.4, p. 11): Pr[∃h∈H:∣err(h)−errT(h)∣≥ρ]≤2∣H∣e−2nρ2.
Exponential weights (proof of Theorem 3.4, p. 12): if ∣err(h)−errT(h)∣<ρ for all h∈H, then Aqε(z) outputs some h with errT(h)≥OPT+2ρ with probability at most ∣H∣e−εnρ/2.
Corollary 3.6: the same guarantee with ln∣H∣ replaced by VCDIM(Cd)⋅ln∣X∣.
Significance
Theorem 3.4 shows that privacy costs little in sample complexity when computation is unconstrained. Any class that some learner with polynomially long output can learn is learnable privately from polynomially many examples. The non-private agnostic bound is O((log∣Cd∣+log(1/β))/α2); the private one is larger by a factor O(α/ε) when α>ε, and is the same otherwise (Remark, p. 12). The theorem is the baseline the rest of the paper refines: an efficient private learner for parity (§4) and the equivalence of local privacy with statistical query learning (§5). It also underlies later work on the sample complexity of private learning, such as Beimel, Kasiviswanathan and Nissim (TCC 2010).
The result is proved in the paper; to our knowledge it has not been machine-checked. The mission produces a checked statement of the exponential mechanism's privacy and of its learning guarantee with an explicit dependence on ∣H∣, ε, α and β. It also provides a reusable bounded i.i.d. Hoeffding inequality (Theorem A.2) in the form learning theory uses. Formalizing the proof with any explicit constant is welcome. The constant 6 printed at the end of the proof does not, by itself, make the last inequality close (the two error terms add to more than β), so a correct constant is a genuine contribution.
Difficulty
Privacy and accuracy pull in opposite directions. Empirical risk minimization, the obvious learner, is not private: replacing one example can change which hypothesis minimizes the training error, and the output then changes with probability one. Smoothing the choice must not destroy accuracy, and the accuracy argument must hold simultaneously for every hypothesis in the class, with the mechanism's randomness and the sample's randomness combined into one failure probability. In Lean the measure-theoretic bookkeeping is non-trivial: the output is a measure on a function space, the failure probability integrates the mechanism's law against the product measure Pn, and the privacy bound must hold for every measurable event, not only for single hypotheses.
Formalization scope
Databases are vectors Fin n → X × Bool, indexed from 0; neighbors differ in exactly one coordinate.
Algorithms are output laws: Aqε(z) is the measure ∑h∈Hpz(h)δh on X → Bool with the product σ-algebra. It is a probability measure because H is nonempty and all weights are positive; the theorems assume H nonempty, and the definition does not build normalization in.
Privacy is required on every measurable set of outputs and for every database size n, with no assumption on how the data arise.
The example domain is a finite type with the discrete σ-algebra (the paper's Xd has representations of size at most d). The goal quantifies it over Type so that the absolute constant can precede it.
Labels are Bool, standing for {0,1}. Logarithms are natural (Real.log), as in the paper's ln.
Parameters are ε>0 and α,β∈(0,1/2) as in Definitions 2.4 and 3.1. The failure event is the strict err(h)>OPT+α, as on p. 11, and OPT is the minimum over H=Cd.
O(⋅) is formalized as one absolute constant C>0 quantified before the domain, the class, the distribution and all parameters. A constant chosen after them would make the statement trivially true.
Theorem A.2 assumes measurable, mutually independent, pairwise identically distributed variables with values in [a,b] at every point. It makes no assumption n≥1 or a<b; in those cases the right side is 2 by Lean's division conventions.
Corollary 3.6 uses the platform's VC dimension ComputationalLearning.vcDim (Kearns–Vazirani's Definition 9, the same notion as the paper's Definition 3.5) and references the platform's proved Sauer lemma.
Contributions welcome: proofs of any milestone, a measurable-structure library for finite output laws, and a general exponential-mechanism privacy theorem from which Lemma 3.3 follows.
Selected references
S. P. Kasiviswanathan, H. K. Lee, K. Nissim, S. Raskhodnikova, A. Smith, What Can We Learn Privately?, SIAM J. Comput. 40(3) (2011) 793–826. arXiv:0803.0924v3, DOI 10.1137/090756090
C. Dwork, F. McSherry, K. Nissim, A. Smith, Calibrating Noise to Sensitivity in Private Data Analysis, TCC 2006. DOI 10.1007/11681878_14
F. McSherry, K. Talwar, Mechanism Design via Differential Privacy, FOCS 2007. DOI 10.1109/FOCS.2007.66
W. Hoeffding, Probability Inequalities for Sums of Bounded Random Variables, J. Amer. Statist. Assoc. 58 (1963) 13–30. DOI 10.1080/01621459.1963.10500830
A. Beimel, S. P. Kasiviswanathan, K. Nissim, Bounds on the Sample Complexity for Private Learning and Private Data Release, TCC 2010. DOI 10.1007/978-3-642-11799-2_26
Additive Logistic Regression: A Statistical View of Boosting 2: Population Discrete AdaBoost Is a Newton-Like Step with Exact Line Search c = ½ log((1−err)/err)Research Paper
Motivation
AdaBoost (Freund and Schapire, 1997) combines many weak classifiers fm(x)∈{−1,1} into a weighted vote sign(∑mcmfm(x)). At each round it fits a classifier to a reweighted training sample, sets cm=log((1−errm)/errm) from the weighted error errm, and multiplies the weights of the misclassified points by ecm. The algorithm was introduced as a procedure for boosting the accuracy of weak learners in the PAC sense, and its good test performance was at first not explained by any statistical model.
Friedman, Hastie and Tibshirani (Ann. Statist. 28 (2000) 337–407) recast boosting as stagewise fitting of an additive model by minimizing the exponential criterionJ(F)=E(e−yF(x)). Their Lemma 1 shows that J is minimized at half the log-odds, so an additive F fitted this way estimates a logistic regression model. Result 1 of that paper, the subject of this mission, shows that Discrete AdaBoost is exactly such a fitting procedure: in its population version, each round is a Newton-like step on J followed by an exact line search, and the resulting reweighting is AdaBoost's. This view underlies gradient boosting, LogitBoost and the loss-function analysis of boosting in later work.
Setting
Let X be a measurable space of features and let ν be a probability law of the pair (x,y) on X×{−1,1}. The conditional law of y given x is the disintegration ν(⋅∣x), and E(⋅∣x) is the integral against it. A current estimate is a measurable F:X→R, and the exponential criterion is
J(F)=E(e−yF(x))∈[0,∞].
With weight w(x,y)=e−yF(x), the weighted conditional expectation and the weighted expectation are
so the population reweighting coincides with Algorithm 1's after renormalization.
Milestones
the second-order expansion of Ew(e−cys∣x) in c (proof of Result 1, p. 346);
(17)–(18): for c>0, (16) is maximization of Ew[yf(x)], solved by (18) (p. 347);
(19): −Ew[yf(x)]=Ew[(y−f(x))2]/2−1 (p. 347);
(20): c=21log((1−err)/err) is the unique minimizer of Ewe−cyf(x) (p. 347);
−yf(x)=2×1[y=f(x)]−1 and the weight identity (p. 347);
Corollary 2 (p. 348), referenced from the platform in its sample form;
(21)–(22): the first-order condition of the line search (p. 348).
Significance
The result. Result 1 turns AdaBoost from a combinatorial reweighting scheme into an optimization method with a known objective. It identifies the population quantity the algorithm estimates (via Lemma 1, half the log-odds), explains the specific constant log((1−err)/err) as an exact line search, explains why the weights of misclassified points are inflated by (1−err)/err, and shows that a weak learner worse than chance is used with reversed sign. It is the template for the other population results of the paper (Real AdaBoost, LogitBoost, Gentle AdaBoost) and for later functional-gradient views of boosting.
Formalizing it. The result is proved on paper; it has not been machine-checked. The paper's argument is stated informally ("≈", "arg min", "Newton-like"); this mission fixes each claim as a precise statement about measures on X×{−1,1}. The sample version of Corollary 2 is already proved on the platform (UnderstandingML.adaboost_reweighted_error) and is reused; everything else is new.
Difficulty
Each step is elementary on its own; the work is in carrying the population statements through measure theory. The weighted expectations are ratios of integrals that are meaningful only under integrability of e−yF(x), and the line-search claim compares J(F+c′f) for all real c′, all of which must be shown finite before they can be compared as real numbers. The pointwise claim (16)–(18) uses the conditional law of y given x, which exists only as a disintegration kernel; passing from it to the unconditional (17) requires the disintegration identity and measurability of (18) in x. The tempting shortcut of treating J(F+cf) as its quadratic approximation does not prove the line-search claim: (20) is the exact minimizer of J, not of the approximation.
Formalization scope
Labels are Booleans, with y=1 for true and y=−1 for false. J is a lower Lebesgue integral in [0,∞] and is compared there; Ew and Ew(⋅∣x) are ratios of Bochner integrals, the conditional one against ν.condKernel x. Disclosed additions, each needed for the page's formulas to be defined: F and f measurable, J(F)<∞, and 0<err<1 wherever log((1−err)/err) is used. No condition err<1/2 is imposed. Ties Ew(y∣x)=0 go to −1, as printed. The "≈" of the expansion is read as "second-order Taylor polynomial in c", and derivatives are stated with HasDerivAt.
Trivializing formalizations are excluded: J is not defined by a conditional formula or as a Bochner integral (which would return 0 on non-integrable weights and make every step optimal), the logarithm is never applied outside its domain, and the class of f stays {−1,1}-valued (over real-valued f the Newton step is Gentle AdaBoost's, a different result).
Everything is in the namespace AddLogReg.Discrete, with one definitions file (Setting) holding J, Ew, the classifier (18), err and the coefficient (20). Reusable pieces are lemmas on disintegrations of laws on X×Bool and on weighted expectations. Proofs of any milestone are welcome independently.
Selected references
J. Friedman, T. Hastie and R. Tibshirani, Additive logistic regression: a statistical view of boosting, Ann. Statist. 28(2):337–407, 2000. https://doi.org/10.1214/aos/1016218223
Y. Freund and R. E. Schapire, A decision-theoretic generalization of on-line learning and an application to boosting, J. Comput. System Sci. 55(1):119–139, 1997. https://doi.org/10.1006/jcss.1997.1504
R. E. Schapire and Y. Singer, Improved boosting algorithms using confidence-rated predictions, Machine Learning 37:297–336, 1999. https://doi.org/10.1023/A:1007614523901
The Weighted Majority Algorithm 8: With η Anomalies Every Deterministic Learner Can Be Forced to opt(F, 0) + 2η MistakesResearch Paper
Learning with anomalies
In the on-line prediction model studied by Littlestone and Warmuth (Inform. and Comput. 108, 1994), a learner receives instances one at a time, predicts a binary label for each, and is then told the correct label. No probabilistic assumption is made: performance is measured by the number of mistakes on a worst-case sequence. When every label is given by some function of a known target classF, the best attainable worst-case number of mistakes is a combinatorial invariant of F, characterized by Littlestone (Mach. Learn. 2, 1988) through mistake trees, now called the Littlestone dimension.
Real data are rarely labelled exactly by a member of the class. Littlestone and Warmuth measure the departure by the number of anomalies: the number of trials inconsistent with the best function in the class. Their Weighted Majority algorithm makes a number of mistakes that grows linearly in the number of anomalies, and the rate of growth can be pushed towards 2 mistakes per anomaly by the choice of the update factor β. Theorem 8.1 of the paper shows that no deterministic learner does better: rate 2 is forced for every class with at least two members. This mission formalizes that lower bound.
Setting
Let X be an arbitrary set (the domain) and F a set of functions X→{0,1}, possibly infinite. A trial is a pair (x,y)∈X×{0,1}, and a sequence of trialsS=((x1,y1),…,(xT,yT)) is finite; an instance may occur several times with different labels.
The number of anomalies of S with respect to f:X→{0,1} is anomf(S)=#{t≤T:f(xt)=yt}.
Sη is the collection of sequences with at most η anomalies with respect to F: those S with minf∈Fanomf(S)≤η.
A deterministic on-line algorithmA maps the list of earlier trials and the current instance to a predicted label. MA(S) is the number of trials t at which A's prediction, computed from (x1,y1),…,(xt−1,yt−1) and xt, differs from yt.
opt(F,η)=minAsupS∈SηMA(S)∈{0,1,…}∪{∞}.
A complete k-mistake tree for F is a complete binary tree of height k whose internal nodes carry instances of X such that every labelling of a root-to-leaf path is realized by some f∈F; in Lean this is the predicate ShattersTree F k v.
Formalization targets
Goal: Theorem 8.1 (p. 250)
For every class F with ∣F∣>1 and every η≥0,
opt(F,η)≥opt(F,0)+2η.
The constant 2 is the paper's, and the statement holds for infinite classes, including those with opt(F,0)=∞.
Milestones
The [Lit88, Lit89b] mistake tree (§8, p. 251): for any pool F there is a complete opt(F,0)-mistake tree. This is published and proved on the platform as UnderstandingML.soa_mistake_bound (MSOA(F)≤Ldim(F)), referenced together with UnderstandingML.ldim_le_mistake_bound; together they identify opt(F,0) with the Littlestone dimension.
The first stage of the adversary (§8, p. 251): if F has a complete k-mistake tree with k≥1, then against any deterministic algorithm there are k−1 trials on which it errs every time, an instance x, and two functions f1,f2∈F consistent with all these trials with f1(x)=f2(x).
Significance
The theorem settles the optimal rate at which mistakes must grow with the number of anomalies for deterministic learners: 2 per anomaly, uniformly over all non-trivial classes. Together with the upper bounds of §2 of the paper, it shows that the deterministic Weighted Majority algorithm is within an arbitrarily small factor of optimal in its dependence on the anomalies, and it explains the factor of 2 that separates deterministic from randomized prediction (Theorem 8.2, a separate mission of this series, gives rate 1 for randomized learners).
The result is proved in the paper; to the best of a search of the platform it has not been formalized. The mission produces a machine-checked lower bound in a general setting: arbitrary domain, arbitrary (possibly infinite) class, extended-natural-valued optimal bounds. The noise-free side, the Littlestone dimension and the Standard Optimal Algorithm, is already published on the platform from Shalev-Shwartz and Ben-David's textbook, and this mission connects it to the noisy setting.
Difficulty
The obvious argument forces one mistake per trial and stops: an adversary that always contradicts the learner leaves no function of F with few anomalies once the trials outnumber the mistake tree. The difficulty is to spend the anomaly budget so that each anomaly buys two mistakes while the sequence stays in Sη, and to do so against an arbitrary algorithm, which may behave adversarially itself. The definition of Sη through a minimum over the class (not a fixed target) and the possibility that opt(F,0) is infinite both have to be handled.
Formalization scope
The domain is any type X; the class is F : Set (X → Bool); the hypothesis ∣F∣>1 is F.Nontrivial, which is meaningful for infinite classes (a cardinality ncard would be 0 there).
Algorithms, histories and mistake counts are those of the published module UnderstandingML_Online (OnlineAlg X Bool, history, mistakes, ShattersTree); an algorithm is an arbitrary function of the history and the current instance, so every deterministic learner is covered.
Sequences are Fin T → X × Bool for every finite length T, with repeated instances and conflicting labels allowed.
opt F η : ℕ∞ is an infimum over algorithms of a supremum over sequences in Sη; η is a natural number. Membership in Sη is "some f∈F has at most η anomalies", which is the paper's minimum over F.
Swapping the quantifiers (one sequence hard for all algorithms) or replacing opt by a finite toNat truncation would give a different statement; neither is used.
Needed infrastructure: the identification opt(F,0)=minAMA(F) (sequences with no anomaly are exactly the realizable ones), the existence of a shattered tree of depth opt(F,0) from the published Littlestone-dimension results, and the two-stage adversary. Contributions of general lemmas about ShattersTree (subtrees, truncation to smaller depth) are reusable for the randomized mission and beyond.
Selected references
N. Littlestone, M. K. Warmuth, The Weighted Majority Algorithm, Information and Computation 108(2), 212–261, 1994. https://doi.org/10.1006/inco.1994.1009
N. Littlestone, Learning quickly when irrelevant attributes abound: a new linear-threshold algorithm, Machine Learning 2, 285–318, 1988. https://doi.org/10.1007/BF00116827
S. Shalev-Shwartz, S. Ben-David, Understanding Machine Learning: From Theory to Algorithms, Cambridge University Press, 2014, Chapter 21. https://doi.org/10.1017/CBO9781107298019
Stable Function Approximation in Dynamic Programming II: In a Nondiscounted MDP with All Policies Proper, Value Iteration Through a Self-Weighted Averager ConvergesResearch Paper
Motivation
Value iteration solves a Markov decision process by repeatedly applying the Bellman backup to a table of state values. When the state space is too large for a table, the values are stored in a function approximator instead: each backup is followed by a fit, and the fitted function is backed up again. This approximate value iteration is the basis of fitted dynamic programming and of much of reinforcement learning, and it was known in the early 1990s to diverge for common approximators such as linear regression and neural networks (Boyan and Moore, 1995, cited in the report as [BM95]).
Gordon's technical report (Stable Function Approximation in Dynamic Programming, CMU-CS-95-103, 1995; conference version at ICML 1995) identified a class of approximators, the averagers (k-nearest-neighbour, kernel averaging, linear and bilinear interpolation on a grid), for which the combination is stable. For discounted problems stability follows from a norm argument. For nondiscounted problems, where the backup contracts only in a weighted norm, that argument fails, and the report proves convergence by a different route: approximate value iteration is exact value iteration on a derived MDP. This mission formalizes that route and its conclusion, Corollary 4.1 (p. 11).
Setting
A finite nondiscounted Markov decision process M has non-goal states x=1,…,n, a finite nonempty action set U(x) at each state, an expected one-step cost cxa, and transition probabilities paxy≥0 with ∑ypaxy≤1. The missing mass 1−∑ypaxy is the probability of moving to a cost-free absorbing goal state, which the paper calls state 1. A value function is a vector V∈Rn; the goal's value is 0. The parallel value backup is
(TMV)(x)=a∈U(x)min[cxa+y∑paxyV(y)],
and with a discount factor γ the sum is multiplied by γ. A policyμ chooses an admissible action at each state; it is proper if, from every starting state, the probability of having reached the goal after t steps tends to 1.
An averagerA consists of constants ky, nonnegative weights βy and nonnegative weights βyz (z ranging over all states, the goal included) with βy+∑zβyz=1. Its mapping is MA(V)(y)=βyky+∑zβyzV(z). Two operators MF and T are compatible if (MF∘T)k(V0) converges from every initial guess V0.
The derived MDPM′ has the same states and actions, transition probabilities paxz′=∑ypaxyβyz to non-goal states z, the rest of the mass going to the goal, and costs cxa′=cxa+γ∑ypaxyβyky.
The states are layered by distance from the goal: S1={1}, Uk=⋃j<kSj, and Sk is the set of states outside Uk from which every action reaches Uk with positive probability. Write k(x) for the index of the layer containing x. The averager is self-weighted for M if for every state y either βy>0 or βyx>0 for some state x with k(x)≤k(y).
Formalization targets
Goal: Corollary 4.1 (p. 11)
If every policy of M is proper and A is self-weighted for M, then for every initial guess V0 both
(MA∘TM)k(V0)and(TM∘MA)k(V0)
converge as k→∞. The first is the page's "TM and MA are compatible", the second its "approximate value iteration … will converge".
Milestones
Theorem 4.1 (Derived MDP, p. 7): TM′=TM∘MA, for every discount factor and for the nondiscounted backup.
Exhaustiveness (p. 10, after [BT89]): if all policies of M are proper, every state lies in some layer Sk.
M′ terminates (p. 10): if all policies of M are proper and A is self-weighted, there is m>0 such that under every policy of M′ and from every state the goal is reached within m steps with positive probability; in particular every policy of M′ is proper.
Theorem 2.3 (p. 4), in its all-policies-terminate case: value iteration converges from any initial guess. This is the published, proved item BertsekasDP.ssp_main_theorem.
A further item, Theorem 6.1 (p. 12), says the limit is V∗ when the optimal value function is a fixed point of MA, under self-weighting or under positive costs with k≥0.
Significance
The corollary gives a checkable, purely combinatorial condition on the approximator and the transition graph under which fitted value iteration on a stochastic shortest path problem cannot diverge. Weight decay (βy>0 everywhere) and any averager that only ever averages a state's value with values no farther from the goal satisfy it. The derived-MDP identity, Theorem 4.1, is the tool the rest of the report builds on (§5, §6.1, Appendix A): it reduces questions about approximate dynamic programming with averagers (convergence, the limit's error, the online Q-learning variant) to questions about an ordinary MDP.
The results are proved in the report, in prose. They have no machine-checked version. The formal development needed here is a finite stochastic shortest path theory (proper policies, the layer partition, termination bounds) on top of the published model BertsekasSSPModel and the proved BertsekasDP.ssp_main_theorem.
Difficulty
The obvious argument fails. If all policies of M are proper, TM contracts in some weighted max norm, and MA is a nonexpansion in the unweighted max norm, but MA can expand the weighted norm, so the composition need not contract. Figure 4 of the report gives a three-state chain and an averager, replacing V(2) by V(3), for which approximate value iteration diverges. The proof therefore has to go through M′ and show that self-weighting makes M′ terminate under every policy. That is a probabilistic statement about all, possibly nonstationary, policies of a new process, derived from a layer structure of the old one. Exhaustiveness of the layers itself needs a policy-construction argument.
Formalization scope
The MDP is the published BertsekasSSPModel n C: non-goal states Fin n, a finite action type C with admissible sets M.U i, costs M.g, and transition probabilities M.p whose deficit is the probability of termination. Gordon's state 1 is this implicit termination state: its value is 0 by construction, and where the partition and the averager's goal weight need to name it, it is none : Option (Fin n). This realises the page's two "without loss of generality" steps (state 1 cost-free and absorbing; β1=1, k1=0). The nondiscounted backup is BertsekasSSPBellmanOp M, the discounted one BertsekasDiscountedBellmanOp M γ.
The averager acts on whole value functions over the state set. An approximator that fits only a sample X0 of states is the special case βyx=0 for x∈/X0. "All policies proper" quantifies over stationary policies μ : Fin n → C with admissible actions, as on p. 3. Self-weightedness is stated through layer membership, which agrees with the page's k(x) whenever the partition is exhaustive (milestone 2). Theorem 6.1 takes V∗ as a real fixed point of the backup; under either of its cases such a fixed point is unique and optimal. Convergence is in Rn, which is the same as convergence in max norm.
The goal's hypotheses concern M and A only. A formalization that assumed properness of the derived MDP, dropped the "every action" quantifier in the layers (making the partition and self-weighting weaker notions), or dropped k(x)≤k(y) (making every averager with a positive weight self-weighted) would be a different theorem.
Reusable beyond this mission: the layer partition and its exhaustiveness, and the termination bound for the derived model, are general finite stochastic shortest path facts. Contributions to any milestone are welcome, and so are proofs of the definitions-level lemmas they rely on.
Selected references
G. J. Gordon, Stable Function Approximation in Dynamic Programming, Technical Report CMU-CS-95-103, Carnegie Mellon University, 1995. Conference version in Proceedings of the 12th International Conference on Machine Learning (ICML), 1995. https://doi.org/10.1016/B978-1-55860-377-6.50040-2
D. P. Bertsekas and J. N. Tsitsiklis, Parallel and Distributed Computation: Numerical Methods, Prentice Hall, 1989. ([BT89] in the report.)
D. P. Bertsekas and J. N. Tsitsiklis, An Analysis of Stochastic Shortest Path Problems, Mathematics of Operations Research 16(3), 1991. https://doi.org/10.1287/moor.16.3.580
J. A. Boyan and A. W. Moore, Generalization in Reinforcement Learning: Safely Approximating the Value Function, Advances in Neural Information Processing Systems 7, Morgan Kaufmann, 1995. ([BM95] in the report.)
Rademacher and Gaussian Complexities: Risk Bounds and Structural Results 3: Monotone, Absolutely Homogeneous, Convex-Hull Invariant and Subadditive R_n, and the 2L Bound for L-Lipschitz Maps Fixing 0Research Paper
Motivation
Generalization bounds in statistical learning theory control the gap between the empirical and the true performance of a predictor chosen from a class F. Bartlett and Mendelson (JMLR 3, 2002) showed that this gap is governed by the Rademacher complexity of a loss class built from F, and that, unlike the VC dimension, this complexity can be estimated from data. Their bounds are useful only if the Rademacher complexity of a large class can be bounded in terms of simpler classes: voting methods take convex combinations of base classifiers, neural networks compose fixed squashing functions with linear combinations, and loss classes compose a class with a Lipschitz loss. Section 3.1 of the paper collects the elementary rules for such constructions in one statement, Theorem 12. These rules are used throughout the later literature on margin bounds, boosting, and neural-network generalization (e.g. Koltchinskii and Panchenko, Ann. Statist. 2002).
This mission formalizes those rules, as stated in the published journal article (J. Mach. Learn. Res. 3 (2002), pp. 463–482), Theorem 12 on p. 469.
Setting
Let X be a measurable space, μ a probability measure on it, and n≥0 an integer. A class is a set F of functions f:X→R; no boundedness is assumed.
For a sample x=(x1,…,xn)∈Xn, the empirical Rademacher complexity of F (Definition 2, p. 464) is
the expectation over independent uniform random signs σ1,…,σn. The Rademacher complexity is its expectation over an i.i.d. sample from μ:
Rn(F)=∫XnR^n(F)(x)dμ⊗n(x).
In the Lean development these are empiricalRademacher n F x and rademacherComplexity μ n F, both valued in [0,∞].
The constructions on classes are those of §2, p. 467: convF is the class of convex combinations of functions from F; −F={−f:f∈F}; absconvF is the class of convex combinations of functions from F∪−F; cF={cf:f∈F}; for ϕ:R→R, ϕ∘F={ϕ∘f:f∈F}; and ∑i=1kFi={f1+⋯+fk:fi∈Fi}.
The goal is the conjunction, each part quantified over its own classes.
Milestones
Parts 1, 3, 2 and 4 are milestones, in the order of the paper's proof. Part 7 remains a theorem item in the mission and a conjunct of the goal.
Significance
The result. Theorem 12 reduces the complexity of a constructed class to the complexities of its building blocks. Part 2 says that convex combinations come for free, which is why margin bounds for boosting depend only on the base class. Part 4 converts a complexity bound for a class into one for a Lipschitz loss composed with it, which is how the risk bounds of the paper's §2 are applied to concrete losses. Parts 1, 3 and 7 are the bookkeeping rules that combine the others; the paper uses parts 1 and 3 to show part 7 is tight.
Formalizing it. The results are proved, and part 4 is due to Ledoux and Talagrand (Probability in Banach Spaces, Springer 1991, Corollary 3.17). The platform already has a one-sided contraction lemma in a different normalization (UnderstandingML.contraction_lemma, Shalev-Shwartz and Ben-David's Lemma 26.9: per-coordinate Lipschitz maps, 1/m normalization, no absolute value, bounded nonempty sets of vectors), referenced here as supporting material, and a convex-hull identity for Mohri et al.'s empirical complexity without absolute value. None of the five parts is formalized for Definition 2's complexity (factor 2/n, absolute value inside the supremum, expectation over the sample, arbitrary classes). The mission provides machine-checked versions in exactly that normalization, reusable by the other missions of this paper and by any later development of Rademacher-complexity bounds.
Difficulty
Parts 1, 2, 3 and 7 are pointwise statements about finite averages of suprema, but they have to be carried through the expectation over the sample and through suprema that may be infinite. Part 2 needs the supremum of the absolute value of a linear functional over a convex hull to equal its supremum over the set, including the symmetric hull conv(F∪−F).
Part 4 is the substantial one. The natural first idea, comparing ∣∑iσiϕ(f(xi))∣ with Lϕ∣∑iσif(xi)∣ term by term for each sign vector, fails: the inequality holds only after averaging over the signs, and only with the factor 2 caused by the absolute value. Without the hypothesis ϕ(0)=0 the statement fails (a constant ϕ≡1 has Lϕ=0 but Rn(ϕ∘F)>0 for nonempty F and n≥1), so any argument must use it.
Formalization scope
Representation. Classes are Set (X → ℝ). Signs are averaged over all 2n vectors Fin n → Bool (encoded ±1). Complexities take values in ℝ≥0∞: a class unbounded on the sample has complexity +∞, the empty class has complexity 0, and Rn is the lower Lebesgue integral against Measure.pi (fun _ : Fin n => μ) with μ a probability measure. This is the paper's meaning for arbitrary classes. A real-valued formalization would not be: a real supremum over an unbounded set and a Bochner integral of a non-integrable function both evaluate to 0. That would make the upper bounds false for unbounded classes and every comparison trivially satisfied for others. No boundedness hypothesis is added anywhere, since Theorem 12 is about arbitrary classes. convF is convexHull ℝ F (finite convex combinations, no closure); absconvF is convexHull ℝ (F ∪ -F); cF is c • F; ∣c∣ enters as ENNReal.ofReal |c|, and 0⋅∞=0 makes c=0 consistent. The Lipschitz map in part 4 is LipschitzWith Lφ φ on all of R, as printed. The classes of part 7 are indexed by Fin k and summed pointwise. At n=0 all complexities are 0 and every part holds trivially, matching the page.
Added hypothesis. Part 7 (and its conjunct in the goal) assumes that each sample function x↦R^n(Fi)(x) is almost-everywhere measurable for μ⊗n. The paper does not discuss measurability. The lower integral is monotone and positively homogeneous without it, so parts 1–4 need no such hypothesis, but it is not additive, and part 7 needs additivity. The hypothesis holds whenever each Fi is a countable class of measurable functions.
Omitted parts. Parts 5 and 6 of Theorem 12 are not posed. Part 5, Rn(F+h)≤Rn(F)+∥h∥∞/n, is false as printed: its proof (p. 470) bounds Esup∣∑σi(f+h)(xi)∣ correctly but drops the factor 2/n of Definition 2, and the correct conclusion is Rn(F)+2∥h∥∞/n. For F={0}, h≡1, n=1 the left side is 2 and the printed right side is 1. Part 6 is derived from part 5. A corrected part 5 would not be the paper's statement, so neither is included. The remark after the theorem that parts 1–3 hold for the Gaussian complexity, and the others with an extra lnn factor, is also not included.
Contributions welcome. Proofs of the milestones; a proof of part 4 from the platform's one-sided contraction lemma, converted to Definition 2's normalization; and general lemmas about empiricalRademacher (finiteness on finite classes, measurability for countable classes) that the other missions of this paper can reuse.
Selected references
P. L. Bartlett and S. Mendelson, Rademacher and Gaussian Complexities: Risk Bounds and Structural Results, Journal of Machine Learning Research 3 (2002), 463–482. https://www.jmlr.org/papers/v3/bartlett02a.html
V. Koltchinskii and D. Panchenko, Empirical margin distributions and bounding the generalization error of combined classifiers, Annals of Statistics 30 (2002), 1–50. https://doi.org/10.1214/aos/1015362183
S. Shalev-Shwartz and S. Ben-David, Understanding Machine Learning: From Theory to Algorithms, Cambridge University Press, 2014, Lemma 26.9. https://doi.org/10.1017/CBO9781107298019
Rademacher and Gaussian Complexities: Risk Bounds and Structural Results 2: With Probability 1 − δ, Every ±1 Classifier in F Has Error ≤ Training Error + R_n(F)/2 + √(ln(1/δ)/(2n))Research Paper
Motivation
A binary classifier is judged by its misclassification probability, the chance that it mislabels a fresh example. That probability is unknown; what a learner sees is the training error, the fraction of mistakes on the sample it was trained on. A generalization bound controls the gap between the two simultaneously for every classifier in a class, so that a classifier chosen by looking at the data still has a guaranteed error. The classical bounds of Vapnik and Chervonenkis measure the class by its VC dimension, a fixed combinatorial quantity that does not depend on the data distribution, and the paper recalls evidence that such fixed complexity penalties cannot be universally effective for model selection (p. 464).
Bartlett and Mendelson's paper in the Journal of Machine Learning Research (2002) studies two data-dependent alternatives, the Rademacher and Gaussian complexities, and proves risk bounds and structural rules for them. This mission formalizes the paper's classification bound, Theorem 5(b), which replaces the VC penalty by half the Rademacher complexity of the class. The source is the published JMLR article (jmlr.org/papers/v3/bartlett02a), pp. 463–482; all page numbers below are the journal's.
Timeline. Uniform convergence of empirical frequencies with a VC-dimension rate is due to Vapnik and Chervonenkis (1971). Rademacher penalties for model selection were introduced by Koltchinskii (2001) and by Bartlett, Boucheron and Lugosi (2002), who also proved Theorem 5(a) with the maximum discrepancy. The 2002 paper proves part (b) in Appendix B (p. 480) by adapting its proof of the general risk bound, Theorem 8 (p. 467).
Setting
Let X be a measurable space and P a probability distribution on X×{±1}; a pair (X,Y)∼P is an example X with label Y. Let F be a set of {±1}-valued functions on X (the classifiers), and let (Xi,Yi)i=1n be an i.i.d. sample drawn from Pn.
The misclassification probability of f∈F is P(Y=f(X)).
The training error of f is P^n(Y=f(X))=n1#{i:Yi=f(Xi)}, where P^n is the empirical measure of the sample (p. 463).
Let σ1,…,σn be independent uniform {±1}-valued random variables, independent of the sample. The empirical Rademacher complexity of a class G of real functions on X at the points x1,…,xn is
R^n(G)(x)=Eσg∈Gsupn2i=1∑nσig(xi),
and the Rademacher complexity is Rn(G)=ER^n(G)(X1,…,Xn) with X1,…,Xn i.i.d. (Definition 2, p. 464). Note the factor 2/n and the absolute value.
In Theorem 5, Rn(F) is the Rademacher complexity of F itself, viewed as a class of real functions with values ±1, under the marginal of P on X.
With the 0–1 lossL(Y,f(X))=1(Y=f(X)), the largest gap on a sample is Φ=suph∈L∘F(Eh−E^nh)=supf∈F(P(Y=f(X))−P^n(Y=f(X))).
Formalization targets
Goal: Theorem 5(b), p. 465
For every 0<δ<1, with probability at least 1−δ over the sample, every f∈F satisfies
P(Y=f(X))≤P^n(Y=f(X))+2Rn(F)+2nln(1/δ).
The bound is uniform: the probability that some f violates it is at most δ.
Milestones (Appendix B, p. 480)
Bounded differences. Replacing one example changes Φ by at most 1/n.
McDiarmid step. With probability at least 1−δ, every f∈F satisfies
P(Y=f(X))≤P^n(Y=f(X))+EΦ+2nln(1/δ).
Symmetrization.EΦ≤Rn(F)/2.
Supporting platform items, referenced and not restated: McDiarmid's inequality for i.i.d. samples (StabGen.Uniform.mcdiarmid_inequality, open) and the symmetrization lemma in the normalization of Shalev-Shwartz and Ben-David (UnderstandingML.representativeness_le_rademacher, proved).
Significance
The result. Theorem 5(b) is a distribution-dependent risk bound for classification whose complexity term can be estimated from a single sample. The paper shows (Theorem 6, p. 465) that it is never much worse than the VC bound, since the empirical Rademacher complexity of a {±1} class is O(d/n) in terms of its empirical VC dimension d, and it can be much better. It also serves as the template for margin bounds for large-margin classifiers, kernel machines and voting methods in Section 4 of the paper.
Formalizing it. The result is proved in the paper; no machine-checked proof is known to exist. A formal proof needs McDiarmid's inequality, a symmetrization argument with a ghost sample, and the reduction of the 0–1 loss class to the classifier class through the identity 1(Y=f(X))=(1−Yf(X))/2. The milestones split these, so they can be proved independently. The symmetrization milestone also corrects a printed slip (see Formalization scope).
Difficulty
The bounded-difference step and the final union of the two halves are short. The difficulty is in the symmetrization. The supremum over an uncountable class is not automatically measurable, so the expectations in the paper's chain need not exist as written; a proof must work with exactly the random variables it integrates and use only the measurability it is given. The exchange of a sample point with its ghost copy, the conditioning on the sample, and the replacement of σiYi by σi must each be justified on finite sign averages and product measures. The obvious shortcut, bounding the gap by the Rademacher complexity of the loss class L∘F, does not give the stated term: with the absolute value of Definition 2, the loss class's complexity is not Rn(F)/2, because (1−Yif(Xi))/2 contributes a sign sum that does not depend on f.
Formalization scope
Lean conventions:
Labels {±1} are the units Z×={1,−1}, coerced to R; the real class of F is {x↦(f(x):R)}. Sample indices are 0,…,n−1; the sample law is the product measure Pn.
The Rademacher complexities take values in [0,∞]: the sign expectation is the exact average over the 2n sign vectors, and the sample expectation is a lower Lebesgue integral. An unbounded class therefore has complexity +∞, not a junk 0. The goal compares values in [0,∞] with ENNReal.ofReal on the real terms.
"With probability at least 1−δ, every f in F" is encoded as: the outer Pn-measure of {S:∃f∈F,the bound fails} is at most δ.
Hypotheses added relative to the page, each necessary or a reading of the page:
n≥1: at n=0 Lean's 2/0=0 makes R0(F)=0 and the bound false.
0<δ<1, as in Theorem 8 (p. 467).
Every f∈F is measurable, and three random variables are measurable: the gap supremum Φ, the double-sample supremum (S,S′)↦supf(P^n′−P^n)(Y=f(X)), and x↦R^n(F)(x). The paper does not discuss measurability; these are exactly the variables the proof integrates, as in the proved platform item UnderstandingML.representativeness_le_rademacher. They hold, for instance, for countable classes.
The two intermediate milestones assume F nonempty, so that the real suprema in them are genuine.
Corrected slip: the chain on p. 480 ends "=Esupfn1∑iσif(Xi)=Rn(F)/2". The last equality is only "≤": Rn carries an absolute value, and for F={f} the left side of that step is 0 while Rn(F)/2>0. The symmetrization milestone states the chain's conclusion with ≤, which is all Theorem 5(b) needs.
Not a valid formalization: a version with a real-valued supremum or Bochner integral for Rn (which would read 0 on an unbounded or non-integrable class and make the bound false or free), one that drops the absolute value or uses the 1/n normalization, or one that measures Rn of the loss class instead of F. Theorem 5(a) (maximum discrepancy, due to Bartlett, Boucheron and Lugosi) is not part of this mission.
Contributions welcome: proofs of the three milestones and the goal; the general McDiarmid inequality (the referenced open item); and reusable lemmas on measurable suprema of classifier families.
Selected references
P. L. Bartlett, S. Mendelson, Rademacher and Gaussian Complexities: Risk Bounds and Structural Results, Journal of Machine Learning Research 3 (2002) 463–482. https://www.jmlr.org/papers/v3/bartlett02a.html
V. Koltchinskii, Rademacher penalties and structural risk minimization, IEEE Transactions on Information Theory 47(5) (2001) 1902–1914. https://doi.org/10.1109/18.930926
V. N. Vapnik, A. Ya. Chervonenkis, On the uniform convergence of relative frequencies of events to their probabilities, Theory of Probability and its Applications 16 (1971) 264–280. https://doi.org/10.1137/1116025
Rademacher and Gaussian Complexities: Risk Bounds and Structural Results 1: With Probability 1 − δ, Expected Loss ≤ Empirical Dominating Cost + R_n(φ̃∘F) + √(8 ln(2/δ)/n) for Every f in FResearch Paper
Motivation
A learning algorithm picks a predictor f from a class F after looking at n training examples, and the quantity of interest is its expected loss on a fresh example. Since f depends on the data, its training error is a biased estimate of that loss, and a risk bound quantifies the bias uniformly over the class: with high probability, every f∈F has expected loss at most an observable sample average plus a complexity penalty. Classical penalties (VC dimension, covering numbers, fat-shattering dimension) are fixed in advance; data-dependent penalties such as the Rademacher complexity can be estimated from the sample itself and adapt to the distribution, which matters for model selection by complexity regularization.
Bartlett and Mendelson, Rademacher and Gaussian Complexities: Risk Bounds and Structural Results (JMLR 3, 2002, jmlr.org), state such a bound in a general decision-theoretic setting (Theorem 8, p. 467): the loss is any [0,1]-valued function of an outcome and an action, and the sample average is taken of a dominating costϕ≥L, which may be Lipschitz and hence amenable to the structural results of the paper's §3 even when L is a discontinuous 0–1 loss. Multiclass classification with error-correcting output codes is the paper's example. The ingredients — McDiarmid's inequality (McDiarmid, 1989) and symmetrization — are standard; the paper combines them for this setting with a centred cost class. The source is the published JMLR article (pp. 463–482); all page numbers below are its printed pages.
Setting
Let X be an input space, Y an output space and A an action space with a distinguished action 0∈A, each a measurable space. A probability measure P on X×Y generates independent examples (X1,Y1),…,(Xn,Yn), and (X,Y) is a fresh draw from P.
A loss is L:Y×A→[0,1]; a costϕ:Y×A→[0,1]dominates it if ϕ(y,a)≥L(y,a) for all y,a.
F is a class of maps X→A; the expected loss of f is EL(Y,f(X))=∫L(y,f(x))dP(x,y).
The empirical mean of h:X×Y→R is E^nh=n1∑i=1nh(Xi,Yi).
The centred cost class is ϕ~∘F={(x,y)↦ϕ(y,f(x))−ϕ(y,0):f∈F}.
For a class G of real functions on a space Z and a sample z1,…,zn, the empirical Rademacher complexity is
R^n(G)=Eσg∈Gsupn2i=1∑nσig(zi),
with σ1,…,σn independent uniform signs, and the Rademacher complexity is Rn(G)=ER^n(G), the expectation over an i.i.d. sample (Definition 2, p. 464). The factor is 2/n and the absolute value is inside the supremum.
In Lean these are RadGauss.RiskBound.empiricalRademacher, rademacherComplexity, empMean, supDev (the uniform deviation suph∈G(Eh−E^nh)), doubleSupDev and phiTildeComp.
Formalization targets
Goal: Theorem 8 (p. 467)
For every integer n≥1 and every 0<δ<1, with probability at least 1−δ over the sample, every f∈F satisfies
EL(Y,f(X))≤E^nϕ(Y,f(X))+Rn(ϕ~∘F)+n8ln(2/δ).
The event is uniform over F. The constant 8 is the printed one.
Milestones (the steps of the paper's proof)
Theorem 9 (McDiarmid's inequality), p. 467: for independent, not necessarily identically distributed Xi and f with bounded differences ci,
P{f(X1,…,Xn)−Ef(X1,…,Xn)≥t}≤e−2t2/∑ici2.
Bounded differences, p. 467: replacing one example changes suph∈ϕ~∘F(Eh−E^nh) by at most 2/n.
Concentration, p. 467: with probability at least 1−δ/2,
hsup(Eh−E^nh)≤Ehsup(Eh−E^nh)+2ln(2/δ)/n.
Combined bound, p. 468: with probability at least 1−δ, for all f∈F,
The platform's Proved UnderstandingML.representativeness_le_rademacher (Shalev-Shwartz and Ben-David, Lemma 26.2) is included as a supporting reference: it is the middle inequality of milestone 5 in the 1/n, absolute-value-free normalization.
Significance
Theorem 8 bounds the expected loss of any predictor in the class, including one selected by the data, by quantities the learner can compute or estimate: the empirical cost and the Rademacher complexity of the centred cost class, which concentrates around its empirical version. Combined with the paper's structural results (§3: contraction by Lipschitz maps, convex hulls, sums), it yields margin bounds for voting methods, neural networks and kernel machines (§4) from complexities of simple base classes. Centring by ϕ(y,0) matters: with the absolute value inside the supremum, Rn(ϕ∘F) can exceed Rn(ϕ~∘F) by order 1/n even for a single function.
The result is proved in the paper. As far as the platform's catalog shows, it has no machine-checked proof in this normalization; McDiarmid's inequality in the non-identically distributed form, the 2/n bounded-difference step and the absolute-value Rademacher symmetrization with factor 2/n are not on the platform. A Lean proof would give a reusable risk-bound template for any loss dominated by a bounded cost.
Difficulty
The obvious argument — bound Eϕ(Y,f(X))−E^nϕ(Y,f(X)) for a fixed f by Hoeffding's inequality — does not survive the choice of f after seeing the data: the deviation must be controlled uniformly over F, and F is typically infinite. The supremum over F is a single random variable, but its expectation is not obviously small, and relating it to the Rademacher complexity requires a ghost sample and a sign-swap symmetry of the product measure. In a formal development, the supremum of an uncountable family is not automatically measurable; the expectations the argument manipulates must be genuine, not lower integrals, which is why the statements carry explicit measurability hypotheses. McDiarmid's inequality itself, for independent but not identically distributed coordinates, is a martingale-difference concentration argument that is not in Mathlib in this form.
Formalization scope
Sample and law. A sample is S : Fin n → X × Y, its law the product Measure.pi (fun _ => P). "With probability at least 1−δ, every f∈F satisfies …" is encoded by bounding the (outer) measure of {S∣∃f∈F,¬bound} by δ, with 0<δ<1.
Complexities in [0,∞].R^n is an average over all 2n sign vectors (Fin n → Bool, true=+1) of [0,∞]-valued suprema; Rn is its Lebesgue integral. The goal is stated in [0,∞], with the real terms embedded by ENNReal.ofReal (all are nonnegative). A real supremum or a Bochner integral would silently return 0 on an unbounded class or a non-integrable integrand, making the bound trivially false or vacuous; this is ruled out by the choice of codomain.
Added hypotheses, all implicit in the paper.n≥1 (at n=0 Lean's 1/0=0 makes the bound read EL≤0, which is false); measurability of L, ϕ and every f∈F; a distinguished action 0 ([Zero A]).
Measurability guard. The goal and milestones 3–5 assume measurability of exactly the random variables the proof integrates: the uniform deviation S↦suph∈ϕ~∘F(Eh−E^nh), the empirical Rademacher complexity S↦R^n(ϕ~∘F), and the double-sample deviation (S,S′)↦suph(n1∑ih(Si′)−E^nh). These hold, for example, for countable F; the same guard is used by UnderstandingML.representativeness_le_rademacher.
No printed slip was found in the statements formalized here. In Theorem 9, if every ci=0 the printed exponent has a zero denominator; Lean reads it as 0 and the bound as 1, which is true.
Not formalized: the paper's Theorem 10 and the applications; this mission covers Definition 2 and §2 up to the end of the proof of Theorem 8.
Contributions welcome: a proof of McDiarmid's inequality for Measure.pi (reusable well beyond this mission), the ghost-sample symmetrization with the sign-swap invariance of Pn⊗Pn, and the final assembly.
Selected references
P. L. Bartlett and S. Mendelson, Rademacher and Gaussian Complexities: Risk Bounds and Structural Results, Journal of Machine Learning Research 3 (2002), 463–482. https://www.jmlr.org/papers/v3/bartlett02a.html
C. McDiarmid, On the method of bounded differences, Surveys in Combinatorics 1989, London Math. Soc. Lecture Note Series 141, Cambridge University Press, 148–188. https://doi.org/10.1017/CBO9781107359949.008
S. Shalev-Shwartz and S. Ben-David, Understanding Machine Learning: From Theory to Algorithms, Cambridge University Press, 2014, Chapter 26. https://doi.org/10.1017/CBO9781107298019
V. Koltchinskii and D. Panchenko, Empirical margin distributions and bounding the generalization error of combined classifiers, Annals of Statistics 30 (2002), 1–50. https://doi.org/10.1214/aos/1015362183
Optimal Best Arm Identification with Fixed Confidence IV: Asymptotic Optimality of the Track-and-Stop StrategyResearch Paper
Motivation
In best arm identification with fixed confidence, a learner samples K unknown distributions (arms) sequentially and must, as early as possible, name the arm with the largest mean, while being wrong with probability at most a prescribed risk δ. The problem models adaptive A/B/n testing, clinical and simulation-based selection among alternatives, and the "ranking and selection" problem of operations research and simulation optimization. The quantity of interest is the sample complexityEμ[τδ], the expected number of samples a strategy takes before stopping.
Timeline of the question this mission formalizes:
Chernoff (1959) introduced sequential tests based on generalized likelihood ratios for adaptive design of experiments, with a finite set of hypotheses (doi:10.1214/aoms/1177706205).
Kaufmann, Cappé and Garivier (2016, JMLR) proved a change-of-measure lower bound on the sample complexity of every δ-PAC strategy (arXiv:1407.4443).
Garivier and Kaufmann (COLT 2016) identified the exact constant T∗(μ) in that lower bound and gave the first strategy, Track-and-Stop, whose sample complexity matches it asymptotically as δ→0 (arXiv:1602.04589). This mission covers the upper-bound half of that paper.
Setting
A canonical one-parameter exponential family is a family of laws νθ, θ∈Θ, on R with density exp(θx−b(θ)) with respect to a reference measure ξ; b is twice differentiable and strictly convex, and νθ has mean b˙(θ). Bernoulli, Poisson and Gaussian laws with known variance are examples. The divergence d(μ,μ′) is the Kullback–Leibler divergence between the members with means μ and μ′.
A bandit model μ=(μ1,…,μK) assigns a member of the family to each arm. The class S consists of models with a unique optimal arm a∗(μ). At each round t=1,2,… the learner picks an arm At as a function of past observations, observes a reward drawn from that arm's law, and at a stopping time τδ recommends an arm. Na(t) is the number of draws of arm a in the first t rounds and μ^a(t) its empirical mean.
With Alt(μ)={λ∈S:a∗(λ)=a∗(μ)} and ΣK the probability simplex, the characteristic time is
T∗(μ)−1=w∈ΣKsupλ∈Alt(μ)infa=1∑Kwad(μa,λa),
and the maximizer w∗(μ) are the optimal proportions of arm draws.
Track-and-Stop combines two ingredients:
a sampling rule that tracks the plug-in proportions w∗(μ^(t)) while forcing each arm to be drawn about t times: C-Tracking tracks the cumulated sum of projections of w∗(μ^(s)) onto ΣKϵs={w∈ΣK:wa≥ϵs}, ϵs=(K2+s)−1/2/2; D-Tracking draws an under-sampled arm when some Na(t)<t−K/2, and otherwise the arm maximizing twa∗(μ^(t))−Na(t);
Chernoff's stopping rule, which stops at the first t at which some arm a beats every other arm b in a generalized likelihood ratio test, Za,b(t)>β(t,δ), here with β(t,δ)=log(r(t)/δ).
Formalization targets
Goal: Theorem 14 (p. 13)
For α∈[1,e/2] and r(t)=O(tα), Chernoff's stopping rule with β(t,δ)=log(r(t)/δ) combined with C-Tracking or D-Tracking satisfies
Lemma 7 (p. 7): C-Tracking ensures Na(t)≥t+K2−2K and maxa∣Na(t)−∑s<twa∗(μ^(s))∣≤K(1+t).
Lemma 8 (p. 7): D-Tracking ensures Na(t)≥(t−K/2)+−1, and proportions within 3(K−1)ϵ of w∗(μ) after a time tϵ that does not depend on the trajectory, once the plug-in targets are within ϵ.
Proposition 9 (p. 8): under either rule, Na(t)/t→wa∗(μ) almost surely.
Lemma 18 (p. 27): an explicit x with c1x≥log(c2xα) for α∈[1,e/2].
Proposition 13 (p. 11): with any sampling rule whose proportions converge almost surely to w∗, τδ<∞ almost surely and limsupδ→0τδ/log(1/δ)≤αT∗(μ) almost surely.
Significance
Theorem 1 of the same paper shows Eμ[τδ]≥T∗(μ)kl(δ,1−δ) for every δ-PAC strategy, and kl(δ,1−δ)∼log(1/δ). Theorem 14 with α=1 therefore shows that the lower bound is attained: T∗(μ) is the exact asymptotic sample complexity of best arm identification in exponential family models, and Track-and-Stop is asymptotically optimal.
The result is proved in the paper. What is not available is a machine-checked proof for exponential families. The platform already holds a Lean development of the Gaussian case following Lattimore and Szepesvári, Bandit Algorithms, Ch. 33, stated for one existentially chosen policy with a different threshold. This mission asks for the universal statement: every run of either tracking rule, for every exponential family, with the paper's thresholds. The tracking lemmas (Lemmas 15, 7, 8) are deterministic combinatorics and reusable by any tracking-based algorithm.
Difficulty
The obvious argument plugs the almost-sure behaviour of Proposition 13 into an expectation. That step fails: almost-sure convergence of τδ/log(1/δ) does not control E[τδ], because on the rare events where the empirical means are far from μ the stopping time may be very large. Theorem 14 needs a quantitative concentration of μ^(t) on events whose complements have summable probability, which in turn relies on the forced exploration guaranteed by the t lower bounds on Na(t) (the concentration step of App. D, Lemmas 19–20).
A second obstacle is the regularity of w∗: the tracking lemmas only transfer convergence of μ^(t) to convergence of Na(t)/t through the continuity of μ↦w∗(μ) on S, proved from the characterization of w∗ in §2.2 (Proposition 6). The GLR statistic also needs its closed form (7) near μ, which requires the empirical means to lie in the interior of the mean space.
Formalization scope
Model. The exponential family is a structure (ξ,Θ,b) with Θ a nonempty open interval, each νθ normalized, b twice continuously differentiable and b¨>0 on Θ. Openness and b¨>0 are added to the paper's "convex, twice differentiable"; strict convexity is what makes νμ unique. Bandit models are parameter vectors θ∈ΘK with K≥2; arms are indexed 0,…,K−1. S is the set of parameter vectors with a unique arm of largest mean b˙(θa).
Protocol. Policies, the trajectory law Pμ, pull counts, empirical means and T∗(μ) are the platform's published definitions (BanditPolicy, BanditTrajectory, TrackAndStop). T∗ uses Kullback–Leibler divergences of the arm laws over the class S and takes values in [0,∞]. Trajectory coordinate t is round t+1. An arm never drawn has empirical mean 0.
Target map.w∗(μ^(t)) is undefined in the paper when μ^(t)∈/S (an unsampled arm, ties, a mean outside b˙(Θ)). Every tracking statement quantifies over every target map with values in ΣK that returns optimal proportions on S, over every choice of L∞ projections, and over every tie-breaking, including randomized ones.
Stopping rule. The two maxima in Za,b(t) are suprema over Θ in the extended reals. Za,b(t)>β is written without subtracting infinities. The stopping time is the first t≥1 at which the test succeeds, +∞ if none. "r(t)=O(tα)" is r(t)≤Dtα for t≥1; r>0 is added so that log(r(t)/δ) is defined.
Values in [0,∞]. Expectations of τδ, the ratios and T∗ live in [0,∞]. No statement converts them to reals, so an infinite expected stopping time is never read as 0.
Corrections, disclosed. Proposition 9's printed Pw is Pμ. Lemma 18 adds c2/c1α>1 and x>0, without which its expressions are undefined.
Ruled out. Specializing to Gaussian arms, or asserting that some sampling policy achieves the bound, would restate existing platform results and is not this theorem: the goal is about every C-Tracking or D-Tracking run in every exponential family.
Welcome contributions. Exponential-family facts (b˙ is the mean, the KL formula, concentration of empirical means); the continuity of w∗ (Proposition 6, App. A.3); Lemma 17 (App. B.2), from which Lemma 8 follows; the closed form (7) of the GLR statistic.
Selected references
A. Garivier, E. Kaufmann, Optimal Best Arm Identification with Fixed Confidence, COLT 2016, JMLR W&CP 49. arXiv:1602.04589v2
E. Kaufmann, O. Cappé, A. Garivier, On the Complexity of Best-Arm Identification in Multi-Armed Bandit Models, JMLR 17, 2016. arXiv:1407.4443