Tutorials › Statistical Learning (Graduate) › Peeling Localized Excess-Risk Classes

Localized Complexity and Stability · Tutorial 175 of 989

Peeling Localized Excess-Risk Classes

Partition excess risks into dyadic shells and combine localized deviation bounds to control the whole class with an explicit probability guarantee.

Advanced 12 min read

What You'll Learn

  • Define dyadic shells for a class of functions with nonnegative excess risk.
  • Combine shellwise deviation guarantees using a summable allocation of failure probabilities.
  • Prove a peeling bound for the supremum over the full class.
  • Apply Bernstein’s inequality to obtain local bounds for finite classes.
  • Distinguish absolute control at small risk from relative control away from zero.

Why Peel a Class by Risk?

A deviation bound for functions with small risk is often sharper than a bound that treats every function in a class as equally complex. Localization makes this possible: one first controls functions whose risk is at most a specified radius. The difficulty is that a class of interest usually contains functions at many different risk levels. A single localized bound need not control all of them at once.

Peeling resolves this by splitting the class into geometrically increasing risk ranges, or dyadic shells. Each shell is contained in a localized class, so a deviation bound at the shell’s upper radius applies to every function in that shell. Allocating the total failure probability across the shells then gives a simultaneous guarantee over the full class.

Definitions and the Peeling Bound

Let \(Z_1,\ldots,Z_n\) be independent observations with common distribution \(P\), and let \(P_n g=n^{-1}\sum_{i=1}^n g(Z_i)\) be the empirical average of a measurable function \(g\). Write \(Pg=E[g(Z)]\). Consider a class \(\mathcal G\) of measurable functions whose risks satisfy \(0\le Pg\le R\). In excess-risk applications, \(Pg\) is the risk in excess of a reference risk and is nonnegative; the functions themselves need not be pointwise nonnegative.

Definition: For a radius \(r>0\), the localized class is \(\mathcal G(r)=\{g\in\mathcal G: Pg\le r\}\). Given a base radius \(r_0>0\), the dyadic radii are \(r_j=2^j r_0\), for integers \(j\ge0\). The first shell is \(\{g:Pg\le r_0\}\); for \(j\ge1\), shell \(j\) is \(\{g:r_{j-1}<Pg\le r_j\}\).

Choose a finite integer \(J\) with \(R\le r_J\). These shells cover the class, including functions of zero risk. Assume the suprema below are measurable; alternatively, the probability statements can be interpreted using outer probability when the usual measurability conditions for empirical-process suprema are not available.

Suppose that for each \(j=0,\ldots,J\) and each \(t>0\), a localized deviation bound is available:

$$ P\left\{\sup_{g\in\mathcal G(r_j)} |P_n g-Pg|>b(r_j,t)\right\}\le e^{-t}. $$

Here \(b(r,t)\) is a nonnegative bound, and the statement is an assumption about the localized classes. The peeling argument does not require a particular formula for \(b\).

Theorem (Dyadic peeling bound): Let \(0<\delta<1\), and set \(\delta_j=\delta\,2^{-(j+1)}\) and \(t_j=\log(1/\delta_j)=\log(1/\delta)+(j+1)\log 2\). Under the localized deviation assumption, with probability at least \(1-\delta\), every \(g\in\mathcal G\) satisfies $$ |P_n g-Pg|\le b(r_j,t_j) $$ for the unique shell index \(j\) containing \(g\). In particular, on this event, $$ \sup_{g\in\mathcal G}|P_n g-Pg|\le \max_{0\le j\le J} b(r_j,t_j). $$

Proof. For each \(j\), apply the assumed localized bound with \(t=t_j\). The probability that the bound fails at radius \(r_j\) is at most \(e^{-t_j}=\delta_j\). By the union bound, the probability of any failure among the \(J+1\) radii is at most

$$ \sum_{j=0}^{J}\delta_j =\delta\sum_{j=0}^{J}2^{-(j+1)} \le \delta. $$

Therefore, with probability at least \(1-\delta\), every localized bound holds simultaneously. On this event, take any \(g\in\mathcal G\). If \(Pg\le r_0\), then \(g\in\mathcal G(r_0)\). Otherwise, there is a shell index \(j\ge1\) with \(r_{j-1}<Pg\le r_j\), so \(g\in\mathcal G(r_j)\). In either case, the deviation bound at the corresponding radius applies to \(g\). Taking the supremum gives the final inequality. \(\square\)

What the Shell Bound Does—and Does Not—Say

The allocation \(\delta_j=\delta 2^{-(j+1)}\) spends more of the failure-probability budget on the first shells and less on distant ones, while ensuring that the total expenditure is at most \(\delta\). The price of looking across shells appears in the confidence parameter \(t_j\): it increases by \(\log 2\) each time the radius doubles. A different summable allocation is possible, but its corresponding confidence parameters must also be included in the localized bounds.

The theorem gives an absolute deviation bound on each shell. It also gives relative control on shells bounded away from zero, provided the shell’s deviation is small enough. For example, if for some \(0<\eta<1\) we have \(b(r_j,t_j)\le\eta r_{j-1}\) for every \(j\ge1\), then on the same simultaneous event, every function in those shells satisfies

$$ |P_n g-Pg|\le\eta r_{j-1}<\eta Pg, \qquad P_n g\ge(1-\eta)Pg. $$

This condition cannot be applied to the first shell by the same reasoning: its members can have risk arbitrarily close to zero, so there is no positive lower bound \(r_{-1}\) to compare with the absolute deviation. Near zero risk, the useful conclusion is typically an absolute bound, not a relative one.

Worked Example: Allocating Error Across Four Shells

Suppose the class has risks at most \(1\), take \(r_0=0.125\), and use four shells with upper radii \(0.125,0.25,0.5,1\). Let \(\delta=0.10\). Assume the localized bounds at the allocated confidence levels give the following values.

Shell \(j\)Risk range\(r_j\)\(\delta_j\)\(b(r_j,t_j)\)
0\(Pg\le0.125\)0.1250.0500.03
1\(0.125<Pg\le0.25\)0.250.0250.05
2\(0.25<Pg\le0.5\)0.50.01250.08
3\(0.5<Pg\le1\)10.006250.12

The four failure probabilities sum to \(0.09375\), which is no greater than \(0.10\). Thus, with probability at least \(1-0.09375=0.90625\), the deviation is at most \(0.03,0.05,0.08,\) or \(0.12\), according to the function’s shell. In particular, the supremum deviation over the class is at most \(0.12\) on this event.

For instance, a function with risk \(Pg=0.31\) lies in shell 2 because \(0.25<0.31\le0.5\); its deviation is bounded by \(0.08\). A function with risk \(Pg=0.02\) lies in shell 0 and is bounded by \(0.03\). The absolute bound in that shell is larger than this function’s risk, which is not a contradiction: peeling alone does not promise relative accuracy near zero.

Obtaining Shellwise Bounds for a Finite Class

A common way to supply the localized bounds is to combine Bernstein’s inequality with a union bound over functions. Suppose \(\mathcal G\) is finite, \(|g(Z)|\le B\) almost surely for every \(g\in\mathcal G\), and \(\operatorname{Var}(g(Z))\le vPg\) for some \(v\ge0\). For each fixed \(g\), Bernstein’s inequality for independent observations gives

$$ P\left\{|P_n g-Pg|> \sqrt{\frac{2\operatorname{Var}(g(Z))t}{n}}+ \frac{2Bt}{3n}\right\} \le 2e^{-t}. $$

The boundedness assumption is stated for \(g\); consequently \(|g(Z)-Pg|\le2B\), as required by this version of Bernstein’s inequality. If \(Pg\le r\), the variance condition implies \(\operatorname{Var}(g(Z))\le vr\). If there are \(N(r)\) functions in \(\mathcal G(r)\), a union bound therefore gives, for every \(t>0\),

$$ P\left\{ \sup_{g\in\mathcal G(r)}|P_n g-Pg|> \sqrt{\frac{2vr t}{n}}+\frac{2Bt}{3n} \right\} \le 2N(r)e^{-t}. $$

In particular, set \(t=\log(2N(r)/\varepsilon)\), for \(0<\varepsilon<1\). Then the right side is \(\varepsilon\), giving an explicit localized bound. Combining these bounds at the dyadic radii with the peeling theorem yields a full-class guarantee. The calculation illustrates how the sample-size term and the number of functions in each localized class both enter the shellwise control.

Worked Example: A Localized Bernstein Bound

Suppose \(n=500\), \(|g(Z)|\le1\), and \(\operatorname{Var}(g(Z))\le2Pg\). At radius \(r=0.10\), assume there are at most \(N(r)=20\) functions in the localized class. To obtain failure probability \(\varepsilon=0.05\) at this radius, use

$$ t=\log\left(\frac{2N(r)}{\varepsilon}\right) =\log\left(\frac{40}{0.05}\right) =\log(800)\approx6.685. $$

Here \(v=2\) and \(B=1\). The square-root term in the localized bound is

$$ \sqrt{\frac{2vr t}{n}} =\sqrt{\frac{2(2)(0.10)(6.685)}{500}} =\sqrt{0.005348} \approx0.0731. $$

The linear term is

$$ \frac{2Bt}{3n} =\frac{2(1)(6.685)}{3(500)} \approx0.00891. $$

Thus the localized deviation is at most approximately \(0.0731+0.00891=0.0820\), with probability at least \(0.95\). This calculation concerns one radius only. To claim simultaneous control over multiple shells, their failure probabilities must still be allocated and summed.

Worked Example: Applying Peeling to Three Numerical Shells

Take \(n=1000\), \(B=1\), \(v=1\), \(\delta=0.05\), and radii \(r_0=0.125\), \(r_1=0.25\), \(r_2=0.5\). Suppose each localized class contains at most four functions. Allocate \(\delta_j=\delta 2^{-(j+1)}\), giving \(0.025,0.0125,\) and \(0.00625\). The Bernstein confidence parameters are \(t_j=\log(2N/\delta_j)\), so

$$ t_0=\log(320)\approx5.768,\qquad t_1=\log(640)\approx6.461,\qquad t_2=\log(1280)\approx7.155. $$

At each radius, use \(b(r_j,t_j)=\sqrt{2r_jt_j/n}+2t_j/(3n)\). The three calculations are

$$ b(r_0,t_0)\approx\sqrt{0.001442}+0.003845 \approx0.03797+0.00385=0.04182, $$
$$ b(r_1,t_1)\approx\sqrt{0.003231}+0.004307 \approx0.05684+0.00431=0.06115, $$
$$ b(r_2,t_2)\approx\sqrt{0.007155}+0.004770 \approx0.08459+0.00477=0.08936. $$

The allocated probabilities sum to \(0.04375\), so all three bounds hold simultaneously with probability at least \(0.95625\). The full-class absolute deviation is at most \(0.08936\) on that event. Moreover, a function in shell 1 has \(Pg>0.125\), and its deviation is at most \(0.06115<0.49(0.125)\). A function in shell 2 has \(Pg>0.25\), and its deviation is at most \(0.08936<0.36(0.25)\). These comparisons yield relative control on the latter shells, but the first-shell bound remains an absolute one.

Common Pitfall: Treating One Local Bound as a Uniform Bound

A tempting but incorrect argument is: “The localized deviation bound holds with probability at least \(1-\delta\) at every radius, so it holds at all radii with probability at least \(1-\delta\).” Individual probability guarantees do not imply a simultaneous guarantee. If there are \(J+1\) radii and each bound fails with probability at most \(\delta\), the union bound only gives a failure probability at most \((J+1)\delta\), which may be much larger than \(\delta\).

The correction is to assign each shell its own failure probability \(\delta_j\) with \(\sum_j\delta_j\le\delta\), and to use the confidence parameter required by that smaller probability in that shell’s localized bound. Another mistake is to extend a relative deviation conclusion down to zero risk. Since \(Pg\) can be arbitrarily small in the first shell, an absolute error bound there does not automatically imply \(|P_n g-Pg|\le\eta Pg\). Keep the absolute and relative conclusions distinct.

Key takeaway: Peeling is a probability bookkeeping argument paired with localization: cover the risk range by dyadic shells, apply the bound at each shell’s upper radius, and make the shell failure probabilities summable. The resulting simultaneous guarantee is only as strong as the localized bounds supplied to it.

Check Your Understanding

Work through these questions using the shell definitions and probability allocation above.

  1. For \(r_0=0.04\), which shell contains a function with risk \(Pg=0.18\)? What localized radius controls its deviation?
  2. Verify that \(\delta_j=\delta 2^{-(j+1)}\) allocates total failure probability at most \(\delta\) over any finite number of shells.
  3. Why does a deviation bound on \(\mathcal G(r)\) apply to every function in shell \(j\) when \(r=r_j\)?
  4. State a sufficient condition for the peeling bound to imply \(P_n g\ge(1-\eta)Pg\) on shells \(j\ge1\).
  5. In the finite-class Bernstein calculation, identify where the class size \(N(r)\) enters the confidence parameter.
Clinical Biostats