Explore-Then-Commit and the Upper Confidence Bound Algorithms

The core finite-armed stochastic bandit algorithms in order of sophistication. Explore-Then-Commit (ETC) and its regret bound; tuning the commitment time; worst-case versus gap-dependent bounds; anytime algorithms, the doubling trick, ε-greedy and phased elimination. The optimism principle and UCB(δ), the random-sample-size subtlety, the good-event proof of the O(Σ log n/Δᵢ) and O(√(nk log n)) bounds. Asymptotically optimal UCB with confidence level f(t) = 1 + t log²t and the optimal constant 2/Δᵢ. MOSS, Doob's maximal inequality, the peeling device and the O(√(kn)) minimax bound, and MOSS's two weaknesses (suboptimality relative to UCB, instability of the regret distribution). Bernoulli concentration through the relative entropy d(p,q), Pinsker's inequality, the Chernoff bound and KL confidence intervals; KL-UCB and its asymptotic regret Σ Δᵢ/d(μᵢ,μ*). Extensions to known per-arm scales, empirical Bernstein (UCB-V), median-of-means and bounded rewards.

· 35 min read · 6,932 words

Corresponds to Lattimore & Szepesvári, Chapters 6–10.

Motivation: one problem, five increasingly sharp answers

Every algorithm here attacks the same problem: kk arms with unknown means, nn rounds, minimize regret. By the regret decomposition, Rn=∑iΔiE[Ti(n)]R_n = \sum_i\Delta_i\mathbb E[T_i(n)], so the whole game is to pull each suboptimal arm as few times as possible while still being sure it is suboptimal. The algorithms differ in how they decide when they’ve seen enough.

AlgorithmExploration ruleNeeds to knowGap-dependent regretWorst case
ETC (Ch. 6)fixed budget mm per arm, then commitnn and Δ\Delta, for tuningΔ+4Δ(1+log⁡nΔ24)\Delta + \frac4\Delta\big(1 + \log\frac{n\Delta^2}4\big) (k=2k=2)O(n)O(\sqrt n) tuned; O(n2/3)O(n^{2/3}) untuned
UCB(δ\delta) (Ch. 7)optimism with δ=1/n2\delta = 1/n^2nn∑i(3Δi+16log⁡nΔi)\sum_i\big(3\Delta_i + \frac{16\log n}{\Delta_i}\big)8nklog⁡n8\sqrt{nk\log n}
Asymptotically optimal UCB (Ch. 8)optimism with log⁡f(t)\log f(t)nothing (anytime)lim sup⁡Rnlog⁡n≤∑i2Δi\limsup\frac{R_n}{\log n}\le\sum_i\frac2{\Delta_i}, optimalO(nklog⁡n)O(\sqrt{nk\log n})
MOSS (Ch. 9)per-arm confidence log⁡+nkTi\log^+\frac n{kT_i}n,kn, kO(kΔmin⁡log⁡nΔmin⁡2k)O\big(\frac k{\Delta_{\min}}\log\frac{n\Delta_{\min}^2}k\big)39kn39\sqrt{kn}, optimal
KL-UCB (Ch. 10)KL confidence intervalsnothinglim sup⁡Rnlog⁡n≤∑iΔid(μi,μ∗)\limsup\frac{R_n}{\log n}\le\sum_i\frac{\Delta_i}{d(\mu_i,\mu^*)}(see notes)

ConventionStanding assumption for Chapters 6–9

Every bandit is in ESGk(1)\mathcal E^k_{SG}(1): each arm’s reward noise is 1-subgaussian (definition). Other known constants σ\sigma follow by scaling rewards (Lemma 5.4). Two points this hides:

  1. every algorithm here relies on knowing σ\sigma;
  2. arms may have different subgaussian constants: easy if they’re known (Exercise 7.2), harder if not (Exercise 7.7).

Chapter 10 replaces subgaussianity with Bernoulli or [0,1][0,1] rewards.

ConventionNotation

Ti(t)=∑s≤tI{As=i}T_i(t) = \sum_{s\le t}\mathbb I\{A_s = i\} is the number of pulls of arm ii after round tt, and μ^i(t)=1Ti(t)∑s≤tI{As=i}Xs\hat\mu_i(t) = \frac1{T_i(t)}\sum_{s\le t}\mathbb I\{A_s = i\}X_s is its empirical mean. μ^is\hat\mu_{is} is the mean of arm ii‘s first ss samples, in the reward-stack model where each arm has a pre-drawn stack. Arm 1 is optimal without loss of generality.


Explore-then-commit (Chapter 6)

AlgorithmExplore-then-commit, ETC(mm)

Input: the number of exploration pulls per arm, mm.

At={(t mod k)+1t≤mk(round-robin exploration)arg⁡max⁡iμ^i(mk)t>mk(commit to the empirical best)A_t = \begin{cases}(t \bmod k) + 1 & t\le mk \quad\text{(round-robin exploration)}\\ \arg\max_i\hat\mu_i(mk) & t > mk\quad\text{(commit to the empirical best)}\end{cases}

Ties in the argmax are broken arbitrarily.

TheoremTheorem 6.1 — Regret of ETC

On any 1-subgaussian bandit with 1≤m≤n/k1\le m\le n/k,

Rn≤m∑i=1kΔi+(n−mk)∑i=1kΔiexp⁡(−mΔi24)R_n\le m\sum_{i=1}^k\Delta_i + (n - mk)\sum_{i=1}^k\Delta_i\exp\Big(-\frac{m\Delta_i^2}4\Big)

ProofDecompose; exploration is deterministic; commitment errors are subgaussian tails

Step 1: Decomposition. Rn=∑iΔiE[Ti(n)]R_n = \sum_i\Delta_i\mathbb E[T_i(n)].

Step 2: Count pulls. Exploration pulls each arm exactly mm times. After that, arm ii is played in every remaining round iff it was the empirical winner:

E[Ti(n)]=m+(n−mk) P(Amk+1=i)≤m+(n−mk) P(μ^i(mk)≥max⁡j≠iμ^j(mk))\mathbb E[T_i(n)] = m + (n - mk)\,\mathbb P(A_{mk+1} = i)\le m + (n - mk)\,\mathbb P\big(\hat\mu_i(mk)\ge\max_{j\ne i}\hat\mu_j(mk)\big)

Step 3: Being beaten by arm ii means a gap-sized deviation. The winner must at least beat arm 1:

P(μ^i(mk)≥μ^1(mk))=P((μ^i−μi)−(μ^1−μ1)≥Δi)\mathbb P\big(\hat\mu_i(mk)\ge\hat\mu_1(mk)\big) = \mathbb P\big((\hat\mu_i - \mu_i) - (\hat\mu_1 - \mu_1)\ge\Delta_i\big)

Step 4: Subgaussian calculus. Each of μ^i−μi\hat\mu_i - \mu_i and μ^1−μ1\hat\mu_1 - \mu_1 is an average of mm independent 1-subgaussian terms, so each is 1/m1/\sqrt m-subgaussian, and their difference is 2/m\sqrt{2/m}-subgaussian (rigorous version: Exercise 6.1). By Theorem 5.3 the probability is at most exp⁡(−Δi22⋅2/m)=e−mΔi2/4\exp\big(-\frac{\Delta_i^2}{2\cdot2/m}\big) = e^{-m\Delta_i^2/4}.

∎

IntuitionThe exploration–exploitation trade-off, in one formula

The first term is the price of exploring (every arm pulled mm times, needed or not). The second is the price of committing wrongly (with probability about e−mΔ2/4e^{-m\Delta^2/4}, losing Δ\Delta per round for the rest of the game). A larger mm shrinks the second term and grows the first.

Tuning the commitment time (k=2k = 2)

With Δ=Δ2\Delta = \Delta_2, the bound becomes

Rn≤mΔ+(n−2m)Δe−mΔ2/4≤mΔ+nΔe−mΔ2/4(6.4)R_n\le m\Delta + (n - 2m)\Delta e^{-m\Delta^2/4}\le m\Delta + n\Delta e^{-m\Delta^2/4} \tag{6.4}

FormulaOptimal mm and the resulting bound

m=max⁡{1, ⌈4Δ2log⁡nΔ24⌉}(6.5)m = \max\Big\{1,\ \Big\lceil\frac4{\Delta^2}\log\frac{n\Delta^2}4\Big\rceil\Big\} \tag{6.5}

Rn≤min⁡{nΔ, Δ+4Δ(1+max⁡{0, log⁡nΔ24})}(6.6)R_n\le\min\Big\{n\Delta,\ \Delta + \frac4\Delta\Big(1 + \max\Big\{0,\ \log\frac{n\Delta^2}4\Big\}\Big)\Big\} \tag{6.6}

TheoremWorst-case bound (Exercise 6.2)

(6.6) implies Rn≤Δ+CnR_n\le\Delta + C\sqrt n for a universal constant CC, so Rn≤1+CnR_n\le1 + C\sqrt n when Δ≤1\Delta\le1.

DerivationMaximize the gap-dependent term over Δ\Delta

If Δ≤2/n\Delta\le2/\sqrt n, then nΔ≤2nn\Delta\le2\sqrt n. Otherwise set x=Δn/2>1x = \Delta\sqrt n/2 > 1, so 4Δ=2nx\frac4\Delta = \frac{2\sqrt n}x and log⁡nΔ24=2log⁡x\log\frac{n\Delta^2}4 = 2\log x:

4Δ(1+log⁡nΔ24)=2n⋅1+2log⁡xx\frac4\Delta\Big(1 + \log\frac{n\Delta^2}4\Big) = 2\sqrt n\cdot\frac{1 + 2\log x}x

The derivative vanishes where 2−(1+2log⁡x)=02 - (1 + 2\log x) = 0, i.e. x=ex = \sqrt e, giving maximum value 2/e≈1.212/\sqrt e\approx1.21. So the term is at most 2.5n2.5\sqrt n, and C=2.5C = 2.5 works. (Constant worked here.)

DefinitionWorst-case versus gap-dependent bounds

  • Worst-case (problem-free, problem-independent, gap-free): depends only on nn and the environment class, as in (6.7).
  • Gap- (problem-, distribution-, instance-) dependent: depends on the specific instance through the gaps, as in (6.6).

WarningWorst-case bounds need bounded rewards

Without a bound such as Δ≤1\Delta\le1, the worst-case regret of ETC, and of every reasonable algorithm that tries each arm at least once, is infinite: an adversary makes the bad arm arbitrarily bad. “Worst-case” always means “over a class with bounded reward range”.

WarningThe tuning cheats

(6.5) uses Δ\Delta, which is never known. The horizon is sometimes known; the gap almost never is.

  • Horizon-only tuning: m≈n2/3m\approx n^{2/3} gives Rn≤(Δ+C)n2/3R_n\le(\Delta + C)n^{2/3} for all instances, but no constant gives Δ+Cn2/3\Delta + Cn^{2/3} (Exercise 6.5a).
  • Data-dependent commitment: alternate arms until a confidence-interval test separates them, which achieves Δ+Clog⁡nΔ\Delta + \frac{C\log n}\Delta and hence Δ+Cnlog⁡n\Delta + C\sqrt{n\log n}. A sharper stopping test achieves Δ+CΔlog⁡max⁡{e,nΔ2}\Delta + \frac C\Delta\log\max\{e, n\Delta^2\} and hence Δ+Cn\Delta + C\sqrt n (Exercise 6.5b–e).

ExampleExperiment 6.1 — Theory versus practice

Gaussian arms with means 00 and −Δ-\Delta, n=1000n = 1000, ETC tuned by (6.5), averaged over 10510^5 runs. The bound (6.6) tracks the true expected regret closely as Δ\Delta ranges over [0,1][0,1]. The regret peaks at intermediate Δ\Delta: small gaps cost little per mistake, and large gaps are found quickly. ETC’s immediate expected regret is also monotonically non-increasing in tt, though not smoothly.

WarningExpected regret hides a bimodal distribution (Exercise 6.9)

With Δ=110\Delta = \frac1{10} and n=2000n = 2000, both the expected regret and the standard deviation of the pseudo-regret depend sharply on mm. ETC’s regret is nearly two-valued: about mΔm\Delta if it commits correctly, and about nΔn\Delta if it commits wrongly. A standard deviation is a poor summary of such a distribution.

Variants and notes

DefinitionAnytime algorithm and the doubling trick

An algorithm is anytime if it needs no advance knowledge of nn. ETC isn’t. The doubling trick restarts a horizon-nℓn_\ell algorithm on epochs of lengths n1<n2<⋯n_1 < n_2 < \cdots (without resetting the global clock). If the base algorithm has regret at most fn(ν)f_{n}(\nu) up to horizon nn, the restarted algorithm has Rn≤∑ℓ≤ℓmax⁡fnℓ(ν)R_n\le\sum_{\ell\le\ell_{\max}}f_{n_\ell}(\nu) (Exercise 6.6).

  • fn=nf_n = \sqrt n with nℓ=2ℓ−1n_\ell = 2^{\ell-1}: Rn≤CnR_n\le C\sqrt n. The geometric sum is dominated by its last term.
  • fn=glog⁡nf_n = g\log n with nℓ=2ℓ−1n_\ell = 2^{\ell-1}: Rn=O(glog⁡2n)R_n = O(g\log^2n), because every epoch pays a log. Faster-growing epochs such as nℓ=22ℓn_\ell = 2^{2^\ell} recover O(glog⁡n)O(g\log n).
  • Cons: restarts throw away data, the transient behaviour is poor, and it is wasteful in practice. Natively anytime algorithms are preferable. (Nothing here uses stochasticity, so it also applies to adversarial and contextual settings.)

Algorithmε-greedy (Exercise 6.7)

Play each arm once. Then in round tt, play the empirical best with probability 1−εt1 - \varepsilon_t, and otherwise a uniformly random arm.

  • A constant ε\varepsilon gives linear regret: Rn/n→εk∑iΔiR_n/n\to\frac\varepsilon k\sum_i\Delta_i.
  • εt=min⁡{1,CktΔmin⁡2}\varepsilon_t = \min\{1, \frac{Ck}{t\Delta_{\min}^2}\} gives Rn≤C′∑i(Δi+ΔiΔmin⁡2log⁡max⁡{e,nΔmin⁡2k})R_n\le C'\sum_i\big(\Delta_i + \frac{\Delta_i}{\Delta_{\min}^2}\log\max\{e, \frac{n\Delta_{\min}^2}k\}\big), but that requires knowing Δmin⁡\Delta_{\min}.

AlgorithmPhased elimination (Exercise 6.8; Auer–Ortner)

Keep an active set Aℓ\mathcal A_\ell, starting from all arms. In phase ℓ\ell, pull every active arm mℓm_\ell times, using only that phase’s samples for the means μ^i,ℓ\hat\mu_{i,\ell}, and keep

Aℓ+1={i:μ^i,ℓ+2−ℓ≥max⁡j∈Aℓμ^j,ℓ}\mathcal A_{\ell+1} = \{i : \hat\mu_{i,\ell} + 2^{-\ell}\ge\max_{j\in\mathcal A_\ell}\hat\mu_{j,\ell}\}

Phase ℓ\ell aims to eliminate every arm with Δi≥2−ℓ\Delta_i\ge2^{-\ell}. With suitable mℓ≈22ℓlog⁡(⋅)m_\ell\approx2^{2\ell}\log(\cdot):

  • the optimal arm is wrongly eliminated in a phase with probability at most k e−mℓ2−2ℓ/4k\,e^{-m_\ell2^{-2\ell}/4};
  • Rn≤C∑i:Δi>0(Δi+log⁡max⁡{e,nΔi2}Δi)R_n\le C\sum_{i:\Delta_i>0}\big(\Delta_i + \frac{\log\max\{e, n\Delta_i^2\}}{\Delta_i}\big);
  • the worst case is ∑iΔi+C′nklog⁡k\sum_i\Delta_i + C'\sqrt{nk\log k}.

RemarkHistory

Robbins (1952) used “certainty equivalence with forcing”: play the empirical best, except on prescribed exploration times. Anscombe (1963) largely solved the Gaussian clinical-trial version. Later relatives include epoch-greedy (contextual bandits) and forced exploration / PEGE (linear bandits). Garivier et al. (2016) showed that all explore-then-commit and elimination policies are asymptotically suboptimal (Chapter 16).


The upper confidence bound algorithm (Chapter 7)

The optimism principle

IntuitionOptimism in the face of uncertainty

Act as if the environment is as nice as plausibly possible. In a new country, choosing between the unknown local cuisine and a familiar chain, optimism says: the local food might be amazing, so try it. After a few visits the statistics settle the question. Pessimism never tries it and may miss something wonderful forever.

ImportantWhy optimism gives sublinear regret

Give each arm an upper confidence bound that overestimates its mean with high probability. Then:

  1. the optimal arm’s index is at least μ∗\mu^* (with high probability);
  2. a suboptimal arm is played only if its index beats the optimal arm’s, hence only if its index is at least μ∗\mu^*;
  3. every play of that arm shrinks its confidence width, so its index soon falls below μ∗\mu^*.

Optimism is self-correcting: overestimates attract plays, and plays remove overestimates.

DefinitionUCB index (7.2)

From (5.6), P(μ≥μ^+2log⁡(1/δ)/n)≤δ\mathbb P\big(\mu\ge\hat\mu + \sqrt{2\log(1/\delta)/n}\big)\le\delta. So

UCBi(t−1,δ)={∞Ti(t−1)=0,μ^i(t−1)+2log⁡(1/δ)Ti(t−1)otherwise.\mathrm{UCB}_i(t-1,\delta) = \begin{cases}\infty & T_i(t-1) = 0,\\[2pt] \hat\mu_i(t-1) + \sqrt{\dfrac{2\log(1/\delta)}{T_i(t-1)}} & \text{otherwise.}\end{cases}

The added term is the exploration bonus (confidence width). An algorithm that plays the arm maximizing such a per-arm score is an index algorithm.

ComponentDetail
InputsArm ii‘s empirical mean and pull count; the confidence level δ\delta.
OutputA plausible optimistic value for μi\mu_i.
Large whenThe arm looks good (μ^i\hat\mu_i large) or is under-explored (TiT_i small).

AlgorithmUCB(δ\delta) — Algorithm 3

For t=1,…,nt = 1, \ldots, n: play At=arg⁡max⁡iUCBi(t−1,δ)A_t = \arg\max_i\mathrm{UCB}_i(t-1,\delta), observe XtX_t, and update.

IntuitionWhy this particular form

If arm 1 has been played far more often, then μ^1≈μ1\hat\mu_1\approx\mu_1, and to keep playing it the algorithm should be confident that every other arm is worse:

μ^i(t−1)+2log⁡(1/δ)Ti(t−1)≤μ1≈μ^1(t−1)+2log⁡(1/δ)T1(t−1)(7.3)\hat\mu_i(t-1) + \sqrt{\frac{2\log(1/\delta)}{T_i(t-1)}}\le\mu_1\approx\hat\mu_1(t-1) + \sqrt{\frac{2\log(1/\delta)}{T_1(t-1)}} \tag{7.3}

Two questions remain: whether the width can be made much narrower (the lower bounds say not for subgaussian noise), and how to choose δ\delta. If δ\delta is too large, the optimal arm’s index may drop below μ∗\mu^* and it gets abandoned, at a cost of linear regret. If δ\delta is too small, suboptimal arms get over-explored. A first guess is δ≈1/n\delta\approx1/n; the analysis below needs δ=1/n2\delta = 1/n^2.

WarningThe sample size inside the index is random (Exercise 7.1)

(5.6) holds for a fixed number of samples, but Ti(t−1)T_i(t-1) depends on the data.

  • If TT is independent of the samples, the bound still holds.
  • If TT is a stopping time, it can fail with probability 1: by the law of the iterated logarithm, lim sup⁡∑t≤nXt2nlog⁡log⁡n=1\limsup\frac{\sum_{t\le n}X_t}{\sqrt{2n\log\log n}} = 1, so the running mean eventually exceeds any fixed-level bound. Stop then.
  • Fix: pay a union-bound penalty, P(μ^−μ≥2log⁡(T(T+1)/δ)/T)≤δ\mathbb P\big(\hat\mu - \mu\ge\sqrt{2\log(T(T+1)/\delta)/T}\big)\le\delta. log⁡log⁡T\log\log T almost suffices (Garivier 2013; Chapter 20).

The proofs below avoid this by union-bounding over all possible sample counts s∈[n]s\in[n].

Regret analysis

TheoremTheorem 7.1 — Gap-dependent regret of UCB

With δ=1/n2\delta = 1/n^2, on any kk-armed 1-subgaussian bandit,

Rn≤3∑i=1kΔi+∑i:Δi>016log⁡nΔiR_n\le3\sum_{i=1}^k\Delta_i + \sum_{i:\Delta_i>0}\frac{16\log n}{\Delta_i}

ProofThe good-event argument

Fix a suboptimal arm ii and a constant ui∈[n]u_i\in[n] (chosen later).

Step 1: The good event.

Gi={μ1<min⁡t∈[n]UCB1(t,δ)}∩{μ^iui+2log⁡(1/δ)ui<μ1}G_i = \Big\{\mu_1 < \min_{t\in[n]}\mathrm{UCB}_1(t,\delta)\Big\}\cap\Big\{\hat\mu_{iu_i} + \sqrt{\frac{2\log(1/\delta)}{u_i}} < \mu_1\Big\}

The optimal arm is never under-estimated, and after uiu_i samples, arm ii‘s index is already below μ1\mu_1.

Step 2: On GiG_i, arm ii is pulled at most uiu_i times. Suppose not. Then some round tt has Ti(t−1)=uiT_i(t-1) = u_i and At=iA_t = i. But

UCBi(t−1,δ)=μ^iui+2log⁡(1/δ)ui<μ1<UCB1(t−1,δ)\mathrm{UCB}_i(t-1,\delta) = \hat\mu_{iu_i} + \sqrt{\frac{2\log(1/\delta)}{u_i}} < \mu_1 < \mathrm{UCB}_1(t-1,\delta)

so arm ii isn’t the argmax, which is a contradiction. Since Ti(n)≤nT_i(n)\le n always,

E[Ti(n)]≤ui+n P(Gic)(7.5)\mathbb E[T_i(n)]\le u_i + n\,\mathbb P(G_i^c) \tag{7.5}

Step 3: The optimal arm fails rarely. Union bound over its possible sample counts:

P(μ1≥min⁡tUCB1(t,δ))≤∑s=1nP(μ1≥μ^1s+2log⁡(1/δ)s)≤nδ(7.7)\mathbb P\Big(\mu_1\ge\min_t\mathrm{UCB}_1(t,\delta)\Big)\le\sum_{s=1}^n\mathbb P\Big(\mu_1\ge\hat\mu_{1s} + \sqrt{\frac{2\log(1/\delta)}s}\Big)\le n\delta \tag{7.7}

Step 4: Arm ii‘s index falls by uiu_i samples. Choose uiu_i so that Δi−2log⁡(1/δ)/ui≥cΔi\Delta_i - \sqrt{2\log(1/\delta)/u_i}\ge c\Delta_i for some c∈(0,1)c\in(0,1). Then

P(μ^iui+2log⁡(1/δ)ui≥μ1)=P(μ^iui−μi≥Δi−⋅)≤P(μ^iui−μi≥cΔi)≤e−uic2Δi2/2\mathbb P\Big(\hat\mu_{iu_i} + \sqrt{\frac{2\log(1/\delta)}{u_i}}\ge\mu_1\Big) = \mathbb P\Big(\hat\mu_{iu_i} - \mu_i\ge\Delta_i - \sqrt{\cdot}\Big)\le\mathbb P(\hat\mu_{iu_i} - \mu_i\ge c\Delta_i)\le e^{-u_ic^2\Delta_i^2/2}

Step 5: Combine and tune. With the smallest valid ui=⌈2log⁡(1/δ)(1−c)2Δi2⌉u_i = \big\lceil\frac{2\log(1/\delta)}{(1-c)^2\Delta_i^2}\big\rceil and δ=1/n2\delta = 1/n^2,

E[Ti(n)]≤⌈2log⁡n2(1−c)2Δi2⌉+1+n1−2c2/(1−c)2\mathbb E[T_i(n)]\le\Big\lceil\frac{2\log n^2}{(1-c)^2\Delta_i^2}\Big\rceil + 1 + n^{1 - 2c^2/(1-c)^2}

The last term is polynomial in nn unless 2c2/(1−c)2≥12c^2/(1-c)^2\ge1. Taking c=12c = \frac12 makes it at most 1, and gives E[Ti(n)]≤3+16log⁡nΔi2\mathbb E[T_i(n)]\le3 + \frac{16\log n}{\Delta_i^2}. Multiply by Δi\Delta_i and sum.

∎

IntuitionThe template for every optimistic analysis

Split “arm ii was played” into (a) the optimal arm’s index failed (rare by concentration and a union bound) or (b) arm ii‘s index is still high (can only last about log⁡(1/δ)/Δi2\log(1/\delta)/\Delta_i^2 pulls). Every UCB-style proof in the book is a refinement of this split.

TheoremTheorem 7.2 — Worst-case regret of UCB

With δ=1/n2\delta = 1/n^2, on any ν∈ESGk(1)\nu\in\mathcal E^k_{SG}(1),

Rn≤8nklog⁡n+3∑i=1kΔiR_n\le8\sqrt{nk\log n} + 3\sum_{i=1}^k\Delta_i

ProofSplit arms at a threshold gap Δ\Delta

Rn=∑Δi<ΔΔiE[Ti(n)]+∑Δi≥ΔΔiE[Ti(n)]≤nΔ+∑Δi≥Δ(3Δi+16log⁡nΔi)≤nΔ+16klog⁡nΔ+3∑iΔiR_n = \sum_{\Delta_i<\Delta}\Delta_i\mathbb E[T_i(n)] + \sum_{\Delta_i\ge\Delta}\Delta_i\mathbb E[T_i(n)]\le n\Delta + \sum_{\Delta_i\ge\Delta}\Big(3\Delta_i + \frac{16\log n}{\Delta_i}\Big)\le n\Delta + \frac{16k\log n}\Delta + 3\sum_i\Delta_i

The small-gap arms together get at most nn pulls. Then choose Δ=16klog⁡n/n\Delta = \sqrt{16k\log n/n}.

∎

IntuitionSmall gaps are cheap, large gaps are found fast

Arms with tiny gaps cost little per pull even if pulled often. Arms with big gaps are pulled rarely. The worst case sits in between, at Δ≈klog⁡n/n\Delta\approx\sqrt{k\log n/n}. The same balancing argument converts any ∑ilog⁡nΔi\sum_i\frac{\log n}{\Delta_i} bound into nklog⁡n\sqrt{nk\log n}.

RemarkHow good is this?

  • The additive ∑iΔi\sum_i\Delta_i is unavoidable: every reasonable algorithm must try every arm once.
  • No algorithm beats O(nk)O(\sqrt{nk}) in the worst case (lower-bounds-for-finite-armed-bandits, Chapter 15). MOSS removes the log⁡n\sqrt{\log n}.

ExampleExperiment 7.1 — UCB versus ETC

k=2k = 2, Gaussian arms, n=1000n = 1000. ETC with the gap-optimal mm (which requires knowing Δ\Delta) beats parameter-free UCB, but only slightly. ETC with any fixed mm (25, 50, 75 or 100) loses badly for some range of Δ\Delta: a bump, a dip, then linear growth as committing wrongly becomes costly. A variant of UCB from Chapter 8 even beats optimally tuned ETC.

Notes and extensions

RemarkTwo problems with δ=1/n2\delta = 1/n^2

A slightly smaller δ\delta improves both regret and empirical performance, and depending on nn means the algorithm isn’t anytime. Chapter 8 fixes both by making the confidence level depend on tt.

RemarkIndependence can be relaxed

All results hold if, with Ht−1\mathcal H_{t-1} the history,

E[Xt∣Ht−1,At]=μAt,E[eλ(Xt−μAt)∣Ht−1,At]≤eλ2/2(7.11–7.12)\mathbb E[X_t\mid\mathcal H_{t-1},A_t] = \mu_{A_t}, \qquad \mathbb E\big[e^{\lambda(X_t - \mu_{A_t})}\mid\mathcal H_{t-1},A_t\big]\le e^{\lambda^2/2} \tag{7.11–7.12}

i.e. conditionally centred, conditionally subgaussian noise. The proofs then use martingale concentration (Chapter 20).

WarningOptimism isn’t universal

It works well when (i) playing an action reveals that action’s quality and (ii) actions reveal nothing about other actions. If (i) fails, even sublinear regret can be lost. If (ii) fails, optimism may skip a low-reward, high-information action whose information would pay off. An example comes with linear bandits (Chapter 25).

RemarkOptimism bias

Psychologists document an “optimism bias” in humans and animals, among the most robust biases known. The authors note that its mathematical justification as a driver of exploration seems to have escaped notice in those fields.

TheoremKnown per-arm scales (Exercise 7.2)

If arm ii is σi\sigma_i-subgaussian with known σi\sigma_i, scale its bonus to σi2log⁡(1/δ)/Ti\sigma_i\sqrt{2\log(1/\delta)/T_i}. The regret becomes O(∑i:Δi>0(Δi+σi2log⁡nΔi))O\big(\sum_{i:\Delta_i>0}(\Delta_i + \frac{\sigma_i^2\log n}{\Delta_i})\big), independent of the optimal arms’ σ\sigma: an optimal arm’s noise only matters through its lower deviations, which the union bound handles cheaply.

TheoremAdapting to the variance: UCB-V (Exercise 7.6)

For rewards in [0,b][0,b], the empirical Bernstein inequality gives

P(∣μ^−μ∣≥2σ^2nlog⁡3δ+3bnlog⁡3δ)≤δ\mathbb P\Big(|\hat\mu - \mu|\ge\sqrt{\frac{2\hat\sigma^2}n\log\frac3\delta} + \frac{3b}n\log\frac3\delta\Big)\le\delta

with σ^2\hat\sigma^2 the empirical variance. An index built on it achieves

Rn≤C∑i:Δi>0(Δi+(b+σi2Δi)log⁡n)(7.14)R_n\le C\sum_{i:\Delta_i>0}\Big(\Delta_i + \big(b + \frac{\sigma_i^2}{\Delta_i}\big)\log n\Big) \tag{7.14}

This replaces UCB’s bΔi\frac b{\Delta_i} with b+σi2Δib + \frac{\sigma_i^2}{\Delta_i}, a big gain for low-variance arms, and essentially unimprovable (Exercise 16.7). See Bernstein.

TheoremHeavy tails: median of means (Exercise 7.7)

Split nn samples into m≈8log⁡(e1/8/δ)m\approx8\log(e^{1/8}/\delta) blocks and take the median of the block means. With only a finite variance,

P(μ^M+192σ2nlog⁡e1/8δ≤μ)≤δ\mathbb P\Big(\hat\mu_M + \sqrt{\frac{192\sigma^2}n\log\frac{e^{1/8}}\delta}\le\mu\Big)\le\delta

A UCB built on it gets Rn≤C∑i(Δi+σ2log⁡nΔi)R_n\le C\sum_i\big(\Delta_i + \frac{\sigma^2\log n}{\Delta_i}\big) on EVk(σ2)\mathcal E^k_V(\sigma^2). Subgaussianity is relaxed to a finite variance at the cost of constants. Alternatives include truncation with empirical Bernstein (which needs to know roughly where the mean is), Catoni’s estimator, and Huber loss. Bubeck et al. (2013) handle even infinite variance, with worse regret.

RemarkLazy updates (Exercises 7.4–7.5)

Phased UCB recomputes the argmax only at phase boundaries (lengths 1,2,4,…1, 2, 4, \ldots), or keeps playing the chosen arm until its count grows by a factor α\alpha. The regret degrades only mildly. This matters with delayed or batched feedback, communication costs, or expensive index computations.


Asymptotically optimal UCB (Chapter 8)

AlgorithmAlgorithm 6

Play each arm once, then

At=arg⁡max⁡i μ^i(t−1)+2log⁡f(t)Ti(t−1),f(t)=1+tlog⁡2tA_t = \arg\max_i\ \hat\mu_i(t-1) + \sqrt{\frac{2\log f(t)}{T_i(t-1)}},\qquad f(t) = 1 + t\log^2t

It is anytime: the confidence level depends on tt, not nn.

TheoremTheorem 8.1 — Finite-time and asymptotic regret

On any 1-subgaussian bandit,

Rn≤∑i:Δi>0inf⁡ε∈(0,Δi)Δi(1+5ε2+2(log⁡f(n)+πlog⁡f(n)+1)(Δi−ε)2)(8.1)R_n\le\sum_{i:\Delta_i>0}\inf_{\varepsilon\in(0,\Delta_i)}\Delta_i\Big(1 + \frac5{\varepsilon^2} + \frac{2\big(\log f(n) + \sqrt{\pi\log f(n)} + 1\big)}{(\Delta_i - \varepsilon)^2}\Big) \tag{8.1}

lim sup⁡n→∞Rnlog⁡n≤∑i:Δi>02Δi(8.2)\limsup_{n\to\infty}\frac{R_n}{\log n}\le\sum_{i:\Delta_i>0}\frac2{\Delta_i} \tag{8.2}

With ε=Δi/2\varepsilon = \Delta_i/2, Rn≤∑i(Δi+1Δi(8log⁡f(n)+8πlog⁡f(n)+28))R_n\le\sum_i\big(\Delta_i + \frac1{\Delta_i}(8\log f(n) + 8\sqrt{\pi\log f(n)} + 28)\big), so Rn≤C∑i(Δi+log⁡nΔi)R_n\le C\sum_i(\Delta_i + \frac{\log n}{\Delta_i}), and the usual argument gives the worst case C∑iΔi+2Cnklog⁡nC\sum_i\Delta_i + 2\sqrt{Cnk\log n}.

ImportantWhy the messier bound

(8.3) has the right order but the wrong constant: dividing by log⁡n\log n gives 8/Δi8/\Delta_i. Only the ε→0\varepsilon\to0 version (8.2) achieves 2/Δi2/\Delta_i, which Lai and Robbins showed no reasonable policy can beat on any instance (lower-bounds-for-finite-armed-bandits, Chapters 15–16).

LemmaLemma 8.2 — How long an index stays high

Let XtX_t be independent and 1-subgaussian with mean 0, μ^t=1t∑s≤tXs\hat\mu_t = \frac1t\sum_{s\le t}X_s, ε,a>0\varepsilon, a > 0, and u=2a/ε2u = 2a/\varepsilon^2. Then

E[∑t=1nI{μ^t+2at≥ε}]≤1+2ε2(a+πa+1)\mathbb E\Big[\sum_{t=1}^n\mathbb I\Big\{\hat\mu_t + \sqrt{\frac{2a}t}\ge\varepsilon\Big\}\Big]\le1 + \frac2{\varepsilon^2}\big(a + \sqrt{\pi a} + 1\big)

ProofCount the first uu rounds for free, then integrate the subgaussian tail

E[κ]≤u+∑t≥⌈u⌉P(μ^t≥ε−2at)≤u+∑t≥⌈u⌉exp⁡(−t2(ε−2at)2)≤1+u+∫u∞(⋅) dt\mathbb E[\kappa]\le u + \sum_{t\ge\lceil u\rceil}\mathbb P\Big(\hat\mu_t\ge\varepsilon - \sqrt{\frac{2a}t}\Big)\le u + \sum_{t\ge\lceil u\rceil}\exp\Big(-\frac t2\Big(\varepsilon - \sqrt{\frac{2a}t}\Big)^2\Big)\le1 + u + \int_u^\infty(\cdot)\,dt

Substituting s=εt−2as = \varepsilon\sqrt t - \sqrt{2a} evaluates the integral to 2ε2(πa+1)\frac2{\varepsilon^2}(\sqrt{\pi a} + 1), and u=2aε2u = \frac{2a}{\varepsilon^2} supplies the remaining term.

∎

IntuitionThe right order

The index μ^t+2a/t\hat\mu_t + \sqrt{2a/t} can’t drop below ε\varepsilon before t≈2a/ε2t\approx2a/\varepsilon^2, because the bonus alone exceeds ε\varepsilon until then. The lemma says it doesn’t stay up much longer.

ProofProof of Theorem 8.1: allow the optimal index to dip by ε

Step 1: Split.

Ti(n)≤∑tI{μ^1(t−1)+2log⁡f(t)T1(t−1)≤μ1−ε}+∑tI{μ^i(t−1)+2log⁡f(t)Ti(t−1)≥μ1−ε, At=i}T_i(n)\le\sum_t\mathbb I\Big\{\hat\mu_1(t-1) + \sqrt{\frac{2\log f(t)}{T_1(t-1)}}\le\mu_1 - \varepsilon\Big\} + \sum_t\mathbb I\Big\{\hat\mu_i(t-1) + \sqrt{\frac{2\log f(t)}{T_i(t-1)}}\ge\mu_1 - \varepsilon,\ A_t = i\Big\}

Step 2: The optimal arm dips rarely. Union bound over its sample count ss:

E[first sum]≤∑t∑sexp⁡(−s2(2log⁡f(t)s+ε)2)≤∑t1f(t)∑se−sε2/2≤5ε2\mathbb E[\text{first sum}]\le\sum_t\sum_s\exp\Big(-\frac s2\Big(\sqrt{\frac{2\log f(t)}s} + \varepsilon\Big)^2\Big)\le\sum_t\frac1{f(t)}\sum_se^{-s\varepsilon^2/2}\le\frac5{\varepsilon^2}

using (a+b)2≥a2+b2(a+b)^2\ge a^2 + b^2. The function ff was chosen so that ∑t1/f(t)\sum_t1/f(t) converges (Exercise 8.1: the geometric sum is at most 2/ε22/\varepsilon^2, and ∑t11+tlog⁡2t\sum_t\frac1{1 + t\log^2t} is at most about 2.52.5).

Step 3: Arm ii‘s index. Replace f(t)f(t) by f(n)f(n) (which only makes the event larger), index by arm ii‘s own sample count (each count occurs in at most one round with At=iA_t = i), and apply Lemma 8.2 with a=log⁡f(n)a = \log f(n) and gap Δi−ε\Delta_i - \varepsilon.

Step 4: Asymptotics. Choose ε=log⁡−1/4n\varepsilon = \log^{-1/4}n and let n→∞n\to\infty.

∎

ImportantThe trick

Chapter 7 needed the optimal arm’s index to stay above μ1\mu_1, which forced δ≈1/n2\delta\approx1/n^2 and a doubled log⁡\log. Here it only needs to stay above μ1−ε\mu_1 - \varepsilon, and occasional dips are allowed as long as their total expected number is bounded (5/ε25/\varepsilon^2). That permits a narrower confidence width, log⁡f(t)≈log⁡t\log f(t)\approx\log t, and recovers the optimal constant 2.

RemarkNotes

  • Any f(t)=tlog⁡αtf(t) = t\log^\alpha t with α>0\alpha > 0 works, with messier algebra. For Gaussian rewards, f(t)=tf(t) = t or even slightly less suffices.
  • Asymptotics can mislead: lower-order terms can dominate in every practical regime.
  • History: Lai–Robbins (1985) proved (8.2) is optimal and achievable; Lai (1987), Katehakis–Robbins (1995) and Agrawal (1995) gave UCB-type asymptotic analyses; finite-time versions came later.

ExampleOne-armed bandits (Exercises 8.2–8.3)

Arm 2 pays a known mean 0, and arm 1 is N(μ1,1)\mathcal N(\mu_1,1). Standard UCB wastes optimism on arm 2. Use the index μ^1+2log⁡f(t)/T1\hat\mu_1 + \sqrt{2\log f(t)/T_1} against the fixed value 0. Then lim sup⁡Rn/log⁡n\limsup R_n/\log n is 0 if μ1≥0\mu_1\ge0 and −2/μ1=2/Δ-2/\mu_1 = 2/\Delta if μ1<0\mu_1 < 0. Optimism is for uncertainty; don’t spend it on what you know.


MOSS: minimax optimality (Chapter 9)

The problem. UCB’s worst case is O(nklog⁡n)O(\sqrt{nk\log n}), and phased elimination gets O(nklog⁡k)O(\sqrt{nk\log k}). The lower bound is Ω(nk)\Omega(\sqrt{nk}). Removing the log requires a confidence level that depends on each arm’s own count.

AlgorithmMOSS — Algorithm 7

Input n,kn, k. Play each arm once, then

At=arg⁡max⁡i μ^i(t−1)+4Ti(t−1)log⁡+(nkTi(t−1)),log⁡+(x)=log⁡max⁡{1,x}A_t = \arg\max_i\ \hat\mu_i(t-1) + \sqrt{\frac4{T_i(t-1)}\log^+\Big(\frac n{kT_i(t-1)}\Big)},\qquad\log^+(x) = \log\max\{1,x\}

(“Minimax Optimal Strategy in the Stochastic case”, Audibert–Bubeck 2009.)

IntuitionStop being optimistic about well-explored arms

Once Ti≥n/kT_i\ge n/k, the bonus is zero: an arm that has received its “fair share” of the budget gets no exploration help. Arms pulled rarely get a bonus of order log⁡(n/(kTi))/Ti\sqrt{\log(n/(kT_i))/T_i}, which is smaller than UCB’s log⁡n/Ti\sqrt{\log n/T_i}. MOSS explores less and saves the log⁡n\sqrt{\log n} factor.

TheoremTheorem 9.1 — Minimax regret of MOSS

Rn≤39kn+∑i=1kΔiR_n\le39\sqrt{kn} + \sum_{i=1}^k\Delta_i

TheoremTheorem 9.2 — Maximal Hoeffding inequality

If X1,…,XnX_1, \ldots, X_n are independent and σ\sigma-subgaussian and St=∑s≤tXsS_t = \sum_{s\le t}X_s, then

P(∃t≤n:St≥ε)≤exp⁡(−ε22nσ2)(9.1)\mathbb P(\exists t\le n : S_t\ge\varepsilon)\le\exp\Big(-\frac{\varepsilon^2}{2n\sigma^2}\Big) \tag{9.1}

This is the same bound as for SnS_n alone, but now uniform over all t≤nt\le n.

ProofDoob’s submartingale inequality applied to eλSte^{\lambda S_t}

eλSte^{\lambda S_t} is a nonnegative submartingale (Exercise 9.1, by Jensen). Doob’s maximal inequality for submartingales gives

P(max⁡t≤neλSt≥eλε)≤E[eλSn]eλε≤exp⁡(nσ2λ22−λε)\mathbb P\Big(\max_{t\le n}e^{\lambda S_t}\ge e^{\lambda\varepsilon}\Big)\le\frac{\mathbb E[e^{\lambda S_n}]}{e^{\lambda\varepsilon}}\le\exp\Big(\frac{n\sigma^2\lambda^2}2 - \lambda\varepsilon\Big)

Then take λ=ε/(nσ2)\lambda = \varepsilon/(n\sigma^2).

∎

LemmaLemma 9.3 — The peeling device

For independent 1-subgaussian XsX_s with mean 0, δ∈(0,1)\delta\in(0,1) and Δ>0\Delta > 0,

P(∃s≥1:μ^s+4slog⁡+(1sδ)+Δ≤0)≤15δΔ2\mathbb P\Big(\exists s\ge1 : \hat\mu_s + \sqrt{\frac4s\log^+\Big(\frac1{s\delta}\Big)} + \Delta\le0\Big)\le\frac{15\delta}{\Delta^2}

ProofUnion bound over a geometric grid of times

Step 1: Peel. Split s≥1s\ge1 into blocks [2j,2j+1][2^j, 2^{j+1}]. Within block jj, make the event larger by replacing ss with the block’s endpoints, then apply the maximal inequality (9.1) on [0,2j+1][0, 2^{j+1}]:

P(⋅)≤∑j≥0exp⁡(−(2j+2log⁡+12j+1δ+2jΔ)22j+2)\mathbb P(\cdot)\le\sum_{j\ge0}\exp\Big(-\frac{\big(\sqrt{2^{j+2}\log^+\frac1{2^{j+1}\delta}} + 2^j\Delta\big)^2}{2^{j+2}}\Big)

Step 2: Algebra. Using (a+b)2≥a2+b2(a+b)^2\ge a^2 + b^2, each term is at most δ2j+1exp⁡(−Δ22j−2)\delta2^{j+1}\exp(-\Delta^22^{j-2}). The terms are unimodal in jj, so the sum is at most the largest term (8δ/(eΔ2)8\delta/(e\Delta^2)) plus an integral. Both are O(δ/Δ2)O(\delta/\Delta^2), and the total is at most 15δ/Δ215\delta/\Delta^2.

∎

IntuitionPeeling

A union bound over every time ss would cost a factor nn. A union bound over geometrically growing blocks costs only the number of blocks, about log⁡n\log n, and the maximal inequality handles each block at the price of its endpoint. This device also proves the law of the iterated logarithm.

ProofProof of Theorem 9.1: control the expected optimism deficit

Step 1: The deficit. Define

Δ=(μ1−min⁡s≤n(μ^1s+4slog⁡+nks))+\Delta = \Big(\mu_1 - \min_{s\le n}\Big(\hat\mu_{1s} + \sqrt{\frac4s\log^+\frac n{ks}}\Big)\Big)^+

This is how far the optimal arm’s index ever drops below its mean. Arms with gaps above 2Δ2\Delta are rarely played. Arms with gaps below 2Δ2\Delta may be played often, but Δ\Delta is small on average.

Step 2: Split the regret. Rn≤E[2nΔ+8kn+∑i:Δi>max⁡{2Δ,8k/n}ΔiTi(n)]R_n\le\mathbb E\big[2n\Delta + 8\sqrt{kn} + \sum_{i:\Delta_i>\max\{2\Delta, 8\sqrt{k/n}\}}\Delta_iT_i(n)\big].

Step 3: The deficit is small on average. By the tail-integral formula and Lemma 9.3 (with δ=k/n\delta = k/n), E[2nΔ]≤2n∫0∞min⁡{1,15knx2}dx≤16kn\mathbb E[2n\Delta]\le2n\int_0^\infty\min\{1, \frac{15k}{nx^2}\}dx\le16\sqrt{kn}.

Step 4: Large-gap arms. For Δi>2Δ\Delta_i > 2\Delta, the optimal index stays above μi+Δi/2\mu_i + \Delta_i/2, so Ti(n)≤κiT_i(n)\le\kappa_i, the number of samples at which arm ii‘s index exceeds μi+Δi/2\mu_i + \Delta_i/2. Lemma 8.2 gives ΔiE[κi]≤Δi+15n/k\Delta_i\mathbb E[\kappa_i]\le\Delta_i + 15\sqrt{n/k} when Δi≥8k/n\Delta_i\ge8\sqrt{k/n}. Summing over at most kk arms gives 15nk+∑Δi15\sqrt{nk} + \sum\Delta_i.

Total: (16+8+15)kn+∑Δi=39kn+∑Δi(16 + 8 + 15)\sqrt{kn} + \sum\Delta_i = 39\sqrt{kn} + \sum\Delta_i.

∎

ImportantOptimism in expectation, not with high probability

UCB’s analysis needed the optimal arm to be optimistic with high probability. MOSS’s bonus is too small for that. Instead the proof bounds the expected amount of pessimism, and lets the arms it harms (those with gaps below 2Δ2\Delta) be ones that cost little anyway.

Two problems with MOSS

WarningMOSS can be much worse than UCB

Take Gaussian arms with n=k3n = k^3, μ1=0\mu_1 = 0, μ2=−k/n\mu_2 = -\sqrt{k/n}, and μi=−1\mu_i = -1 for i>2i > 2. Then

RnUCB=O(klog⁡k),RnMOSS=Ω(kn)=Ω(k2)R_n^{\mathrm{UCB}} = O(k\log k), \qquad R_n^{\mathrm{MOSS}} = \Omega(\sqrt{kn}) = \Omega(k^2)

Minimax and asymptotic optimality both hide this. Only fully finite-time, instance-dependent bounds reveal it.

WarningPushing expected regret down makes the regret unstable

With two arms and a well-tuned algorithm, the pseudo-regret is roughly bimodal:

R^n≈{nΔwith probability δ (optimal arm misidentified),1Δlog⁡1δotherwise,Rn=O(nΔδ+1Δlog⁡1δ)\hat R_n\approx\begin{cases}n\Delta & \text{with probability }\delta\ (\text{optimal arm misidentified}),\\ \frac1\Delta\log\frac1\delta & \text{otherwise,}\end{cases}\qquad R_n = O\Big(n\Delta\delta + \frac1\Delta\log\frac1\delta\Big)

  • Minimizing the mean (δ≈1/(nΔ2)\delta\approx1/(n\Delta^2)) gives Rn=O(1Δ(1+log⁡nΔ2))R_n = O(\frac1\Delta(1 + \log n\Delta^2)), but on the failure event FF the second moment E[IFR^n2]=Ω(n)\mathbb E[\mathbb I_F\hat R_n^2] = \Omega(n).
  • A slightly smaller δ\delta (=(nΔ)−2= (n\Delta)^{-2}) barely raises the mean, to O(1Δ(1n+log⁡n2Δ2))O(\frac1\Delta(\frac1n + \log n^2\Delta^2)), while E[IFR^n2]=O(1)\mathbb E[\mathbb I_F\hat R_n^2] = O(1).

The last log factor of expected regret is bought with heavy tails.

Notes

RemarkFixes and refinements

  • Asymptotics: MOSS has lim sup⁡Rnlog⁡n≤∑i4Δi\limsup\frac{R_n}{\log n}\le\sum_i\frac4{\Delta_i}. Replacing the 4 with 2 and peeling on finer grids [ξj,ξj+1][\xi^j,\xi^{j+1}] (Garivier 2013) recovers the optimal constant. For Gaussian noise, Brownian-motion boundary-crossing bounds avoid peeling altogether (Exercise 9.4: the Bachelier–Lévy formula via the method of images, and tangent bounds for concave boundaries).
  • UCB*: the index μ^i+4Tilog⁡+nTi\hat\mu_i + \sqrt{\frac4{T_i}\log^+\frac n{T_i}} (no kk) is never worse than UCB, and has worst case O(nklog⁡k)O(\sqrt{nk\log k}) (Exercise 9.3). “Improved UCB” (phased elimination) is the same idea in disguise.
  • Best of both: optimally confident UCB, μ^i+2(1+ε)Tilog⁡nt\hat\mu_i + \sqrt{\frac{2(1+\varepsilon)}{T_i}\log\frac nt}, and AdaUCB, with confidence level log⁡n∑jmin⁡{Ti,TiTj}\log\frac n{\sum_j\min\{T_i, \sqrt{T_iT_j}\}}, are minimax optimal and never worse than UCB. AdaUCB is also asymptotically optimal.
  • Anytime MOSS: Degenne–Perchet (2016).

Bernoulli rewards and KL-UCB (Chapter 10)

The problem. Bernoulli rewards (clicks, conversions) are 12\frac12-subgaussian for every mean, so UCB already applies. But a Bernoulli arm’s variance μ(1−μ)\mu(1-\mu) depends on its mean, and near 0 or 1 the empirical mean concentrates much faster than the worst case. Confidence intervals should shrink accordingly.

Concentration for sums of Bernoulli random variables

DefinitionRelative entropy between Bernoullis (Def. 10.1)

d(p,q)=plog⁡pq+(1−p)log⁡1−p1−qd(p,q) = p\log\frac pq + (1-p)\log\frac{1-p}{1-q}

with limits at the boundary: d(0,q)=log⁡11−qd(0,q) = \log\frac1{1-q}, d(1,q)=log⁡1qd(1,q) = \log\frac1q, and d(p,0)d(p,0), d(p,1)d(p,1) equal 0 if p=qp = q and ∞\infty otherwise. This is the relative entropy D(B(p)∥B(q))D(\mathcal B(p)\Vert\mathcal B(q)), and it is the large-deviation rate of a Bernoulli average.

ComponentDetail
Signatured:[0,1]×[0,1]→[0,∞]d : [0,1]\times[0,1]\to[0,\infty]
First argumentThe observed (or hypothesized) frequency.
Second argumentThe true (or reference) mean.
Asymmetricd(p,q)≠d(q,p)d(p,q)\ne d(q,p) in general.

LemmaLemma 10.2 — Properties of dd

For p,q,ε∈[0,1]p, q, \varepsilon\in[0,1]:

  1. d(⋅,q)d(\cdot,q) and d(p,⋅)d(p,\cdot) are convex, with unique minimizers at qq and pp respectively.
  2. Pinsker: d(p,q)≥2(p−q)2d(p,q)\ge2(p-q)^2.
  3. If p≤q−ε≤qp\le q - \varepsilon\le q, then d(p,q−ε)≤d(p,q)−d(q−ε,q)≤d(p,q)−2ε2d(p, q-\varepsilon)\le d(p,q) - d(q-\varepsilon, q)\le d(p,q) - 2\varepsilon^2.
ProofConvexity from the binary entropy; Pinsker by calculus; (3) by linearity in pp

(1) d(⋅,q)d(\cdot,q) is the negative binary entropy h(p)=plog⁡p+(1−p)log⁡(1−p)h(p) = p\log p + (1-p)\log(1-p), with h′′=1p+11−p>0h'' = \frac1p + \frac1{1-p} > 0, plus a linear function of pp. d(p,⋅)d(p,\cdot) is h(p)h(p) plus the convex functions plog⁡1qp\log\frac1q and (1−p)log⁡11−q(1-p)\log\frac1{1-q}. And d(p,q)>0d(p,q) > 0 unless p=qp = q.

(2) g(x)=d(p,p+x)−2x2g(x) = d(p, p+x) - 2x^2 has g(0)=0g(0) = 0; show by differentiation that it is minimized there on [−p,1−p][-p, 1-p] (Exercise 10.1).

(3) p↦d(p,q−ε)−d(p,q)p\mapsto d(p, q-\varepsilon) - d(p, q) is linear and increasing in pp, so for p≤q−εp\le q-\varepsilon it is at most its value at q−εq - \varepsilon, namely −d(q−ε,q)-d(q-\varepsilon, q). Then apply Pinsker.

∎

LemmaLemma 10.3 — Chernoff’s bound

For i.i.d. Bernoulli(μ\mu) with sample mean μ^\hat\mu:

P(μ^≥μ+ε)≤e−n d(μ+ε,μ)  (ε∈[0,1−μ]),P(μ^≤μ−ε)≤e−n d(μ−ε,μ)  (ε∈[0,μ])(10.1–10.2)\mathbb P(\hat\mu\ge\mu + \varepsilon)\le e^{-n\,d(\mu+\varepsilon,\mu)}\ \ (\varepsilon\in[0,1-\mu]), \qquad \mathbb P(\hat\mu\le\mu - \varepsilon)\le e^{-n\,d(\mu-\varepsilon,\mu)}\ \ (\varepsilon\in[0,\mu]) \tag{10.1–10.2}

ProofCramér–Chernoff with the exact Bernoulli MGF

For λ>0\lambda > 0,

P(μ^≥μ+ε)≤Eeλ∑t(Xt−μ)eλnε=(μeλ(1−μ−ε)+(1−μ)e−λ(μ+ε))n\mathbb P(\hat\mu\ge\mu + \varepsilon)\le\frac{\mathbb E e^{\lambda\sum_t(X_t - \mu)}}{e^{\lambda n\varepsilon}} = \Big(\mu e^{\lambda(1-\mu-\varepsilon)} + (1-\mu)e^{-\lambda(\mu+\varepsilon)}\Big)^n

The minimizing λ=log⁡(μ+ε)(1−μ)μ(1−μ−ε)\lambda = \log\frac{(\mu+\varepsilon)(1-\mu)}{\mu(1-\mu-\varepsilon)} makes the bracket equal e−d(μ+ε,μ)e^{-d(\mu+\varepsilon,\mu)}. The lower tail is symmetric.

∎

RemarkCredit

Chernoff attributed the bound to his friend Herman Rubin, but the name stuck.

IntuitionChernoff versus Hoeffding

By Pinsker, e−nd≤e−2nε2e^{-nd}\le e^{-2n\varepsilon^2}, which recovers Hoeffding. The two agree near μ=12\mu = \frac12. Near 0 or 1, Chernoff is far stronger, because the true deviation rate d(μ±ε,μ)d(\mu\pm\varepsilon,\mu) is much larger than 2ε22\varepsilon^2 there.

CorollaryCorollary 10.4 — KL confidence intervals

For a≥0a\ge0, P(d(μ^,μ)≥a, μ^≤μ)≤e−na\mathbb P(d(\hat\mu,\mu)\ge a,\ \hat\mu\le\mu)\le e^{-na}, and likewise for μ^≥μ\hat\mu\ge\mu. Define

U(a)=max⁡{u∈[0,1]:d(μ^,u)≤a},L(a)=min⁡{u∈[0,1]:d(μ^,u)≤a}U(a) = \max\{u\in[0,1] : d(\hat\mu,u)\le a\}, \qquad L(a) = \min\{u\in[0,1] : d(\hat\mu,u)\le a\}

Then P(μ≥U(a))≤e−na\mathbb P(\mu\ge U(a))\le e^{-na} and P(μ≤L(a))≤e−na\mathbb P(\mu\le L(a))\le e^{-na}. With a=log⁡(1/δ)/na = \log(1/\delta)/n, these are 1−δ1-\delta upper and lower confidence bounds.

ProofMonotonicity of dd on each side of μ^\hat\mu

Lower deviations: d(⋅,μ)d(\cdot,\mu) is decreasing on [0,μ][0,\mu], so {d(μ^,μ)≥a, μ^≤μ}={μ^≤μ−x}\{d(\hat\mu,\mu)\ge a,\ \hat\mu\le\mu\} = \{\hat\mu\le\mu - x\}, where xx solves d(μ−x,μ)=ad(\mu - x,\mu) = a. Apply (10.2).

Upper confidence bound: U≥μ^U\ge\hat\mu and d(μ^,⋅)d(\hat\mu,\cdot) is strictly increasing on [μ^,1][\hat\mu,1] with d(μ^,U)=ad(\hat\mu,U) = a, so {μ≥U}={d(μ^,μ)≥a, μ≥μ^}\{\mu\ge U\} = \{d(\hat\mu,\mu)\ge a,\ \mu\ge\hat\mu\}. Use the first part.

∎

RemarkComputing the bounds

UU and LL have no closed form, but dd is convex in its second argument, so Newton’s method solves d(μ^,u)=ad(\hat\mu,u) = a quickly.

ExampleExample 10.5 — An asymmetric interval

μ^=34\hat\mu = \frac34, n=10n = 10, 90% confidence: a=log⁡(10)/10≈0.23a = \log(10)/10\approx0.23. The interval {x:d(34,x)≤0.23}\{x : d(\frac34, x)\le0.23\} is highly asymmetric: it extends further below 34\frac34 than above. Hoeffding’s half-width, log⁡(1/δ)/2n≈0.34\sqrt{\log(1/\delta)/2n}\approx0.34, gives a similar lower bound but a vacuous upper bound (0.75+0.34>10.75 + 0.34 > 1). KL intervals automatically respect [0,1][0,1] and tighten as μ^→0\hat\mu\to0 or 1.

The KL-UCB algorithm

AlgorithmKL-UCB — Algorithm 8

Play each arm once, then

At=arg⁡max⁡i max⁡{μ~∈[0,1]:d(μ^i(t−1),μ~)≤log⁡f(t)Ti(t−1)},f(t)=1+tlog⁡2tA_t = \arg\max_i\ \max\Big\{\tilde\mu\in[0,1] : d\big(\hat\mu_i(t-1),\tilde\mu\big)\le\frac{\log f(t)}{T_i(t-1)}\Big\},\qquad f(t) = 1 + t\log^2t

The index is the largest mean that is still statistically plausible in KL terms.

TheoremTheorem 10.6 — Regret of KL-UCB

If Xt∼B(μAt)X_t\sim\mathcal B(\mu_{A_t}), then

Rn≤∑i:Δi>0 inf⁡ε1,ε2>0ε1+ε2<ΔiΔi(log⁡f(n)d(μi+ε1, μ∗−ε2)+12ε12+2ε22)R_n\le\sum_{i:\Delta_i>0}\ \inf_{\substack{\varepsilon_1,\varepsilon_2>0\\ \varepsilon_1+\varepsilon_2<\Delta_i}}\Delta_i\Big(\frac{\log f(n)}{d(\mu_i+\varepsilon_1,\ \mu^*-\varepsilon_2)} + \frac1{2\varepsilon_1^2} + \frac2{\varepsilon_2^2}\Big)

lim sup⁡n→∞Rnlog⁡n≤∑i:Δi>0Δid(μi,μ∗)\limsup_{n\to\infty}\frac{R_n}{\log n}\le\sum_{i:\Delta_i>0}\frac{\Delta_i}{d(\mu_i,\mu^*)}

ImportantNever worse than UCB, much better near the edges

UCB tuned for 12\frac12-subgaussian rewards has asymptotic constant 12Δi\frac1{2\Delta_i}. By Pinsker, d(μi,μ∗)≥2Δi2d(\mu_i,\mu^*)\ge2\Delta_i^2, so Δid≤12Δi\frac{\Delta_i}{d}\le\frac1{2\Delta_i}: KL-UCB is never worse asymptotically. When μi≈μ∗\mu_i\approx\mu^*, a Taylor expansion gives

d(μi,μ∗)=Δi22μi(1−μi)+o(Δi2) ⟹ lim sup⁡Rnlog⁡n≈∑i:Δi>02μi(1−μi)Δi(10.5)d(\mu_i,\mu^*) = \frac{\Delta_i^2}{2\mu_i(1-\mu_i)} + o(\Delta_i^2)\ \Longrightarrow\ \limsup\frac{R_n}{\log n}\approx\sum_{i:\Delta_i>0}\frac{2\mu_i(1-\mu_i)}{\Delta_i} \tag{10.5}

The worst-case scale 14\frac14 is replaced by the arm’s actual variance μi(1−μi)\mu_i(1-\mu_i), which helps a lot when means are near 0 or 1 (rare clicks).

WarningThis is not a CLT effect

(10.5) looks like “Bernoulli behaves like a Gaussian with matched variance”, but that approximation’s quality doesn’t improve with nn. As n→∞n\to\infty the confidence level must drive failure probabilities below n−1/2n^{-1/2}, where the CLT says nothing (why the CLT fails for tails). The KL rate is the true large-deviation rate.

LemmaLemma 10.7 — The optimal arm’s index rarely dips

With dˉ(p,q)=d(p,q) I{p≤q}\bar d(p,q) = d(p,q)\,\mathbb I\{p\le q\}, i.i.d. Bernoulli(μ\mu) samples, and

τ=min⁡{t:max⁡1≤s≤n(dˉ(μ^s,μ−ε)−log⁡f(t)s)≤0},E[τ]≤2ε2\tau = \min\Big\{t : \max_{1\le s\le n}\Big(\bar d(\hat\mu_s, \mu - \varepsilon) - \frac{\log f(t)}s\Big)\le0\Big\},\qquad\mathbb E[\tau]\le\frac2{\varepsilon^2}

After time τ\tau, the KL index built from any sample count is at least μ−ε\mu - \varepsilon.

ProofShift the reference by ε, then Chernoff, then integrate the tail

P(τ>t)≤∑s=1nP(dˉ(μ^s,μ−ε)>log⁡f(t)s, μ^s<μ−ε)≤∑sP(d(μ^s,μ)>log⁡f(t)s+2ε2, μ^s<μ)\mathbb P(\tau > t)\le\sum_{s=1}^n\mathbb P\Big(\bar d(\hat\mu_s,\mu-\varepsilon) > \frac{\log f(t)}s,\ \hat\mu_s < \mu - \varepsilon\Big)\le\sum_s\mathbb P\Big(d(\hat\mu_s,\mu) > \frac{\log f(t)}s + 2\varepsilon^2,\ \hat\mu_s < \mu\Big)

using Lemma 10.2(3). By Corollary 10.4 this is at most ∑se−log⁡f(t)−2sε2≤12f(t)ε2\sum_se^{-\log f(t) - 2s\varepsilon^2}\le\frac1{2f(t)\varepsilon^2}. Integrating, E[τ]≤12ε2∫0∞dtf(t)≤2ε2\mathbb E[\tau]\le\frac1{2\varepsilon^2}\int_0^\infty\frac{dt}{f(t)}\le\frac2{\varepsilon^2}.

∎

LemmaLemma 10.8 — A suboptimal arm’s index stays high only briefly

For i.i.d. Bernoulli(μ\mu), Δ>0\Delta > 0, a>0a > 0 and κ=∑s=1nI{d(μ^s,μ+Δ)≤a/s}\kappa = \sum_{s=1}^n\mathbb I\{d(\hat\mu_s,\mu+\Delta)\le a/s\},

E[κ]≤inf⁡ε∈(0,Δ)(ad(μ+ε,μ+Δ)+12ε2)\mathbb E[\kappa]\le\inf_{\varepsilon\in(0,\Delta)}\Big(\frac a{d(\mu+\varepsilon,\mu+\Delta)} + \frac1{2\varepsilon^2}\Big)

ProofEither the mean overshoots or ss is still small

Since d(⋅,μ+Δ)d(\cdot,\mu+\Delta) is decreasing on [0,μ+Δ][0,\mu+\Delta], {d(μ^s,μ+Δ)≤as}⊆{μ^s≥μ+ε}∪{s≤u}\{d(\hat\mu_s,\mu+\Delta)\le\frac as\}\subseteq\{\hat\mu_s\ge\mu+\varepsilon\}\cup\{s\le u\} with u=a/d(μ+ε,μ+Δ)u = a/d(\mu+\varepsilon,\mu+\Delta). So E[κ]≤u+∑se−s d(μ+ε,μ)≤u+1d(μ+ε,μ)≤u+12ε2\mathbb E[\kappa]\le u + \sum_s e^{-s\,d(\mu+\varepsilon,\mu)}\le u + \frac1{d(\mu+\varepsilon,\mu)}\le u + \frac1{2\varepsilon^2}, by Chernoff and then Pinsker.

∎

ProofProof of Theorem 10.6

Fix a suboptimal arm ii and ε1,ε2>0\varepsilon_1,\varepsilon_2 > 0 with ε1+ε2<Δi\varepsilon_1 + \varepsilon_2 < \Delta_i. Before τ\tau (Lemma 10.7 applied to arm 1 with ε2\varepsilon_2), count every round. After τ\tau, arm 1’s index is at least μ∗−ε2\mu^* - \varepsilon_2, so arm ii is played only when its own index reaches μ∗−ε2\mu^* - \varepsilon_2, i.e. when d(μ^i,Ti,μ∗−ε2)≤log⁡f(n)/Tid(\hat\mu_{i,T_i}, \mu^* - \varepsilon_2)\le\log f(n)/T_i. That is the count κ\kappa of Lemma 10.8 with μ+Δ=μ∗−ε2\mu + \Delta = \mu^* - \varepsilon_2 and a=log⁡f(n)a = \log f(n):

E[Ti(n)]≤E[τ]+E[κ]≤2ε22+log⁡f(n)d(μi+ε1,μ∗−ε2)+12ε12\mathbb E[T_i(n)]\le\mathbb E[\tau] + \mathbb E[\kappa]\le\frac2{\varepsilon_2^2} + \frac{\log f(n)}{d(\mu_i+\varepsilon_1,\mu^*-\varepsilon_2)} + \frac1{2\varepsilon_1^2}

Multiply by Δi\Delta_i and sum. For the asymptotics, let ε1,ε2→0\varepsilon_1,\varepsilon_2\to0 slowly with nn (Exercise 10.2).

∎

Notes and extensions

RemarkBernoulli is the worst case on [0,1][0,1]

For any X∈[0,1]X\in[0,1] with mean μ\mu, convexity of x↦eλxx\mapsto e^{\lambda x} gives EeλX≤μeλ+1−μ\mathbb E e^{\lambda X}\le\mu e^\lambda + 1 - \mu, the Bernoulli MGF (Garivier–Cappé). So Chernoff’s bound, and KL-UCB’s guarantee lim sup⁡Rnlog⁡n≤∑Δid(μi,μ∗)\limsup\frac{R_n}{\log n}\le\sum\frac{\Delta_i}{d(\mu_i,\mu^*)}, hold for any [0,1][0,1]-valued rewards. They even hold for adapted sequences with deterministic conditional means (Exercise 10.3), though the bound degrades toward Hoeffding when the means are spread between 0 and 1.

RemarkNot quite tight for general bounded rewards

The optimal asymptotic constant over [0,1][0,1]-supported distributions uses di=inf⁡{D(Pi,P):μ(P)>μ∗, Supp(P)⊆[0,1]}d_i = \inf\{D(P_i,P) : \mu(P) > \mu^*,\ \mathrm{Supp}(P)\subseteq[0,1]\}, a KL projection onto distributions with a better mean. This is never smaller than d(μi,μ∗)d(\mu_i,\mu^*), and is achieved by DMED and IMED (Honda–Takemura).

RemarkExponential families (Exercises 10.4–10.5)

For a one-parameter exponential family arm, replace dd by the family’s KL divergence to get lim sup⁡Rnlog⁡n≤∑iΔidi,inf⁡\limsup\frac{R_n}{\log n}\le\sum_i\frac{\Delta_i}{d_{i,\inf}} (Cappé et al. 2013). For non-canonical families, ties must be broken toward the most-played arm, and the lower-order term 1/di,min⁡1/d_{i,\min} may dominate in practice.

RemarkHistory

Lai (1987) gave the first asymptotic version. Finite-time KL-UCB was proved simultaneously by Garivier–Cappé (2011) and Maillard et al. (2011). KL-UCB-type asymptotic optimality is also known for Thompson sampling (Kaufmann et al.; Korda et al.) and BayesUCB. Ménard–Garivier (2017) combined MOSS and KL-UCB to be minimax and asymptotically optimal for exponential families.

ExampleKL-UCB versus UCB (Exercise 10.6)

With n=10,000n = 10{,}000 and k=2k = 2, compare μ1=12\mu_1 = \frac12, μ2=12+Δ\mu_2 = \frac12 + \Delta (where the two should be similar, since the variance is near its maximum) with μ1=110\mu_1 = \frac1{10} or 910\frac9{10} (where KL-UCB should win clearly, since its intervals shrink near the boundary).


Quick reference

ObjectStatementNote
ETC regretm∑Δi+(n−mk)∑Δie−mΔi2/4m\sum\Delta_i + (n-mk)\sum\Delta_ie^{-m\Delta_i^2/4}difference of means is 2/m\sqrt{2/m}-subgaussian
ETC tuned (k=2k=2)m≈4Δ2log⁡nΔ24m\approx\frac4{\Delta^2}\log\frac{n\Delta^2}4; Rn≤Δ+4Δ(1+log⁡nΔ24)R_n\le\Delta + \frac4\Delta(1 + \log\frac{n\Delta^2}4)needs Δ\Delta
ETC worst caseΔ+Cn\Delta + C\sqrt n tuned; (Δ+C)n2/3(\Delta + C)n^{2/3} untuned
Doubling trickRn≤∑ℓfnℓR_n\le\sum_\ell f_{n_\ell}n\sqrt n preserved; log⁡n\log n needs faster growth
UCB indexμ^i+2log⁡(1/δ)/Ti\hat\mu_i + \sqrt{2\log(1/\delta)/T_i}optimism
UCB, δ=1/n2\delta = 1/n^2ETi(n)≤3+16log⁡nΔi2\mathbb E T_i(n)\le3 + \frac{16\log n}{\Delta_i^2}good-event proof
UCB worst case8nklog⁡n+3∑Δi8\sqrt{nk\log n} + 3\sum\Delta_isplit at Δ≈klog⁡n/n\Delta\approx\sqrt{k\log n/n}
Asymptotically optimal UCBlog⁡f(t)\log f(t), f=1+tlog⁡2tf = 1 + t\log^2t; lim sup⁡Rnlog⁡n≤∑2Δi\limsup\frac{R_n}{\log n}\le\sum\frac2{\Delta_i}Lai–Robbins optimal
Lemma 8.2E#{μ^t+2a/t≥ε}≤1+2ε2(a+πa+1)\mathbb E\#\{\hat\mu_t + \sqrt{2a/t}\ge\varepsilon\}\le1 + \frac2{\varepsilon^2}(a + \sqrt{\pi a} + 1)
Maximal HoeffdingP(∃t≤n:St≥ε)≤e−ε2/2nσ2\mathbb P(\exists t\le n : S_t\ge\varepsilon)\le e^{-\varepsilon^2/2n\sigma^2}Doob
Peelingunion bound over [2j,2j+1][2^j, 2^{j+1}]Lemma 9.3: 15δ/Δ215\delta/\Delta^2
MOSSbonus 4Tilog⁡+nkTi\sqrt{\frac4{T_i}\log^+\frac n{kT_i}}; Rn≤39kn+∑ΔiR_n\le39\sqrt{kn} + \sum\Delta_iminimax optimal
MOSS failuren=k3n = k^3 instance: Ω(k2)\Omega(k^2) versus UCB’s O(klog⁡k)O(k\log k)
Bernoulli KLd(p,q)=plog⁡pq+(1−p)log⁡1−p1−qd(p,q) = p\log\frac pq + (1-p)\log\frac{1-p}{1-q}Pinsker: ≥2(p−q)2\ge2(p-q)^2
ChernoffP(μ^≥μ+ε)≤e−nd(μ+ε,μ)\mathbb P(\hat\mu\ge\mu+\varepsilon)\le e^{-nd(\mu+\varepsilon,\mu)}valid for any [0,1][0,1] rewards
KL intervalU(a)=max⁡{u:d(μ^,u)≤a}U(a) = \max\{u : d(\hat\mu,u)\le a\}asymmetric, respects [0,1][0,1]
KL-UCBlim sup⁡Rnlog⁡n≤∑Δid(μi,μ∗)≈∑2μi(1−μi)Δi\limsup\frac{R_n}{\log n}\le\sum\frac{\Delta_i}{d(\mu_i,\mu^*)}\approx\sum\frac{2\mu_i(1-\mu_i)}{\Delta_i}never worse than UCB
UCB-VC∑(Δi+(b+σi2/Δi)log⁡n)C\sum(\Delta_i + (b + \sigma_i^2/\Delta_i)\log n)empirical Bernstein
Median of meansfinite variance sufficesconstants grow

Facts to keep to hand

  • Every analysis is “pulls of each bad arm”. Bound E[Ti(n)]\mathbb E[T_i(n)] by splitting into “the optimal arm’s index failed” (rare: concentration plus a union bound over sample counts) and “arm ii‘s index is still high” (lasts about log⁡(⋅)/Δi2\log(\cdot)/\Delta_i^2 pulls). Then convert ∑log⁡nΔi\sum\frac{\log n}{\Delta_i} into a worst-case bound by splitting at a threshold gap.
  • ETC’s weakness is its fixed budget. It is nearly optimal when tuned with Δ\Delta and nn, but Δ\Delta is unknown. Adaptive stopping, elimination and optimism all fix this by letting the data decide when exploration is done.
  • The confidence level is the whole design. 1/n21/n^2 (simple proof) leads to log⁡f(t)\log f(t) (anytime, optimal constant 2) and then to per-arm log⁡+nkTi\log^+\frac n{kT_i} (minimax, no log factor). Each refinement trades high-probability optimism for optimism in expectation, and pays in instance-dependent behaviour or regret tails.
  • Optimism must match the noise model. Subgaussian widths are worst-case. KL widths for Bernoulli or bounded rewards, variance-adaptive widths (UCB-V), and robust estimators (median of means) all win by modelling the noise more faithfully.
  • No single criterion tells the whole story. Minimax optimality (MOSS), asymptotic optimality (Lai–Robbins), finite-time instance behaviour (UCB beats MOSS on some instances), and regret tails (instability) can disagree.

References

← All notes