Tutorials › Statistical Learning (Graduate) › Fast Rates from Bernstein Conditions for Kernel Scores

Kernel Classification Guarantees · Tutorial 776 of 989

Fast Rates from Bernstein Conditions for Kernel Scores

Learn how a Bernstein condition localizes stochastic error around near-optimal kernel scores and improves the resulting excess-risk rate.

Advanced 12 min read

What You'll Learn

  • Define excess losses and the Bernstein condition for a kernel-score class.
  • Derive a high-probability bound for empirical risk minimization.
  • Solve the resulting self-bounding inequality for different Bernstein exponents.
  • Compare localized rates with a global uniform-deviation bound.
  • Check when a kernel-score loss is bounded and why the Bernstein assumption needs justification.

Why Variance Near the Optimum Matters

A uniform-deviation argument controls empirical risk error equally across a whole score class. That is safe, but it can be wasteful: a score with large excess risk need not be estimated as precisely as one that is nearly optimal. A Bernstein condition formalizes a useful form of localization. It says that the variance of an excess loss decreases as its mean excess decreases. For empirical risk minimization, this makes the stochastic error smaller precisely where the selected score must be compared with an optimum.

We develop the argument for a finite set of kernel scores, where the complexity contribution can be stated explicitly. The same variance calculation is the core of more general localized empirical-process bounds, in which a local complexity measure replaces the logarithm of the number of candidates. The finite-class setting also makes clear which assumptions produce the fast rate and which do not.

Kernel Scores and Excess Losses

Let \(Z=(X,Y)\) be an observation, and let \(Z_1,\ldots,Z_n\) be independent observations with the same distribution as \(Z\). Let \(\mathcal F\) be a finite set of real-valued scores from a reproducing kernel Hilbert space (RKHS), with \(|\mathcal F|=M\). For a kernel satisfying \(K(x,x)\leq \kappa^2\), every score with \(\|f\|_{\mathcal H}\leq R\) obeys \(|f(x)|\leq \kappa R\). We use a loss \(\ell(y,f(x))\) taking values in \([0,1]\), and write \(P\ell_f=E[\ell(Y,f(X))]\) and \(P_n\ell_f=n^{-1}\sum_{i=1}^n\ell(Y_i,f(X_i))\).

Assume that some \(f_0\in\mathcal F\) minimizes population risk over \(\mathcal F\). For each \(f\), its excess loss relative to \(f_0\) is the random variable \(g_f(Z)=\ell(Y,f(X))-\ell(Y,f_0(X))\). Its mean is the excess risk \(\mu_f=Pg_f=P\ell_f-P\ell_{f_0}\), which is nonnegative. An empirical risk minimizer is any \(\hat f\in\mathcal F\) satisfying \(P_n\ell_{\hat f}=\min_{f\in\mathcal F}P_n\ell_f\).

Definition: The class satisfies a Bernstein condition with constants \(B>0\) and \(0\leq\beta\leq1\) if, for every \(f\in\mathcal F\), $$ E[g_f(Z)^2]\leq B\mu_f^\beta. $$ For \(\beta=0\), interpret \(\mu_f^0=1\), including when \(\mu_f=0\). Since \(\operatorname{Var}(g_f)\leq E[g_f^2]\), this condition also bounds the variance by \(B\mu_f^\beta\). It is an assumption about the distribution, loss, and score class—not a consequence of using a kernel.

Here \(g_f\in[-1,1]\), because each loss is in \([0,1]\). The condition says that when a score has small excess risk, its excess loss also has small second moment. The case \(\beta=1\) is especially consequential: the variance is at most proportional to the excess risk, rather than merely being bounded by a constant.

Worked Example: Verifying a Bernstein Condition

Suppose a candidate score has excess loss \(g_f(Z)=\mathbf 1_A(Z)\) for an event \(A\) with probability \(p\), and the reference score has zero loss. Then \(g_f\) is a valid loss difference in \([0,1]\), and $$ \mu_f=E[g_f]=p,\qquad E[g_f^2]=E[\mathbf 1_A^2]=p. $$ Thus \(E[g_f^2]=\mu_f\), so the Bernstein condition holds with \(B=1\) and \(\beta=1\). This calculation verifies the condition for this candidate; to use it for a class, it must hold for every candidate in that class.

A Localized Excess-Risk Bound

The key concentration step is Bernstein's inequality. In the form needed here, if \(W_1,\ldots,W_n\) are independent copies of a random variable \(W\) with finite variance and \(|W-EW|\leq 2\), then for \(t>0\), $$ P\left(\frac{1}{n}\sum_{i=1}^n W_i<EW-\sqrt{\frac{2\operatorname{Var}(W)t}{n}}-\frac{4t}{3n}\right)\leq e^{-t}. $$ For an excess loss \(g_f\), the centered variable has absolute value at most \(2\), so the stated version applies. Taking \(t=\log(M/\delta)\) and using a union bound over the \(M\) candidates gives, with probability at least \(1-\delta\), the lower inequality below simultaneously for all \(f\). The variance term can then be bounded using the Bernstein condition.

Theorem: Suppose \(\mathcal F\) is finite, \(f_0\) minimizes population risk over \(\mathcal F\), the loss lies in \([0,1]\), and the Bernstein condition holds with \(B>0\) and \(0\leq\beta\leq1\). For an empirical risk minimizer \(\hat f\), let \(s=\log(M/\delta)/n\), where \(0<\delta<1\). With probability at least \(1-\delta\), $$ P\ell_{\hat f}-P\ell_{f_0} \leq \min\left\{1,\ (8Bs)^{1/(2-\beta)}+\frac{8}{3}s\right\}. $$

Proof. Set \(t=\log(M/\delta)\). For each \(f\), apply the one-sided Bernstein inequality to \(g_f(Z_i)\), whose mean is \(\mu_f\). A union bound shows that, with probability at least \(1-\delta\), every \(f\in\mathcal F\) satisfies $$ P_ng_f\geq \mu_f-\sqrt{\frac{2\operatorname{Var}(g_f)t}{n}}-\frac{4t}{3n} \geq \mu_f-\sqrt{2Bs\,\mu_f^\beta}-\frac{4}{3}s. $$ On this event, empirical minimization gives \(P_ng_{\hat f}\leq0\), since \(f_0\in\mathcal F\). Write \(\mu=\mu_{\hat f}\). Combining the two inequalities yields $$ \mu\leq \sqrt{2Bs}\,\mu^{\beta/2}+\frac{4}{3}s. $$ Set \(a=\sqrt{2Bs}\), \(c=4s/3\), and \(\alpha=\beta/2<1\). If \(\mu\leq2c\), then \(\mu\leq2c+a^{1/(1-\alpha)}2^{1/(1-\alpha)}\) immediately. If \(\mu>2c\), then \(\mu-c>\mu/2\); the displayed inequality implies \(\mu/2<a\mu^\alpha\). Dividing by \(\mu^\alpha>0\) and raising both sides to the power \(1/(1-\alpha)\) gives $$ \mu<(2a)^{1/(1-\alpha)}=(8Bs)^{1/(2-\beta)}. $$ In either case, \(\mu\leq (8Bs)^{1/(2-\beta)}+2c\), and \(2c=8s/3\). Finally \(0\leq\mu\leq1\) because both risks lie in \([0,1]\), so we may also take the minimum with \(1\). \(\square\)

Worked Example: A Fast Rate for a Finite Kernel Dictionary

Suppose \(M=20\) candidate kernel scores are compared using a bounded loss, with \(n=1000\), confidence level \(1-\delta=0.95\), and a verified Bernstein condition with \(B=1\), \(\beta=1\). Then $$ s=\frac{\log(20/0.05)}{1000} =\frac{\log(400)}{1000} \approx\frac{5.9914645}{1000} =0.0059914645. $$ The theorem gives $$ (8s)+\frac{8}{3}s =\left(8+\frac{8}{3}\right)(0.0059914645) \approx0.06391. $$ Thus, with probability at least \(0.95\), the empirical minimizer's excess risk over the best candidate is at most \(0.06391\) (or the bound \(1\), if that were smaller). The calculation uses the unrounded value of \(s\) before rounding the final result.

Why This Improves on Global Uniform Deviation

For comparison, Hoeffding's inequality and a union bound give a global deviation statement: with probability at least \(1-\delta\), every candidate satisfies $$ |P_ng_f-Pg_f|\leq\sqrt{\frac{2\log(2M/\delta)}{n}}. $$ On that event, \(P_ng_{\hat f}\leq0\) implies \(\mu_{\hat f}\leq2\sqrt{2\log(2M/\delta)/n}\). This rate is proportional to \(n^{-1/2}\), regardless of how little variance the excess losses have near the optimum. Under the Bernstein condition, the theorem instead has leading order \(s^{1/(2-\beta)}\). For \(\beta=1\), this is order \(n^{-1}\), up to the logarithm of the class size; for \(\beta=0\), it is order \(n^{-1/2}\). The improvement comes from controlling fluctuations at the scale of the excess risk, not from changing the empirical minimization rule.

Worked Example: Comparing the Two Bounds

Continue with \(M=20\), \(n=1000\), and \(\delta=0.05\). The global bound is $$ 2\sqrt{\frac{2\log(2M/\delta)}{n}} =2\sqrt{\frac{2\log(800)}{1000}} \approx2\sqrt{\frac{13.369223}{1000}} \approx0.23125. $$ Under the additional condition \(B=1,\ \beta=1\), the localized bound from the preceding example is approximately \(0.06391\). The latter is substantially smaller for these inputs. This comparison is between valid upper bounds, not a claim that the realized risk must equal either number.

Worked Example: How the Bernstein Exponent Changes the Rate

Take \(M=10\), \(n=2000\), \(\delta=0.1\), and \(B=2\). Here $$ s=\frac{\log(10/0.1)}{2000} =\frac{\log(100)}{2000} \approx0.0023025851. $$ If \(\beta=1/2\), the bound is $$ (8Bs)^{1/(2-\beta)}+\frac{8}{3}s =(16s)^{2/3}+\frac{8}{3}s \approx(0.0368414)^{2/3}+0.0061402 \approx0.11072+0.00614 \approx0.11686. $$ If instead \(\beta=0\), the corresponding bound is $$ (16s)^{1/2}+\frac{8}{3}s \approx\sqrt{0.0368414}+0.0061402 \approx0.19194+0.00614 \approx0.19808. $$ Both calculations use the same \(B,s\); the smaller excess-risk bound for \(\beta=1/2\) reflects stronger variance decay as the excess risk shrinks.

Reading the Bound Correctly

The theorem is finite-class and compares the empirical minimizer with the best population-risk score in that class. It does not by itself control the gap between that best score and the Bayes classifier. Nor does it establish the Bernstein condition from the RKHS norm constraint. For an infinite score class, a corresponding result needs a uniform localized deviation bound; the class-size term must be replaced by an appropriate local complexity or entropy control. A finite candidate dictionary avoids that additional step while retaining the variance-based rate mechanism.

The earlier material on margin conditions gives one route for understanding when low-noise assumptions can lead to variance control for classification losses. But the implication depends on the chosen surrogate and on how its excess loss relates to classification risk. In an application, verify the relevant condition for the actual loss and comparator rather than inferring it solely from the presence of a margin assumption.

Common pitfall: A small global uniform-deviation bound does not automatically imply a fast rate. The fast-rate argument needs a variance bound that shrinks with excess risk. Conversely, the Bernstein condition alone does not establish that an infinite RKHS class has manageable complexity. Keep the distributional variance assumption and the class-complexity control as distinct requirements.
Key takeaway: Bernstein localization replaces a fixed fluctuation scale by one proportional to \(\mu_f^{\beta/2}\). For an empirical minimizer this yields excess risk of order \(s^{1/(2-\beta)}+s\), with \(s=\log(M/\delta)/n\); the strongest case \(\beta=1\) gives a fast, near-\(n^{-1}\) rate for a finite score class.

Check Your Understanding

Use the definitions and assumptions stated above to answer these questions.

  1. Why does \(g_f\in[-1,1]\) when the loss takes values in \([0,1]\), and what bound does this give on \(|g_f-Eg_f|\)?
  2. For \(\beta=1\), what is the leading dependence of the theorem's bound on \(n\), with \(M,\delta,B\) fixed?
  3. Why is the empirical minimizer's excess loss empirically nonpositive relative to \(f_0\)?
  4. What additional fact, beyond an RKHS norm bound, is needed to apply the finite-class theorem to a proposed score class?
  5. How does the rate exponent change when \(\beta\) decreases from \(1\) to \(0\)?
Clinical Biostats