Tutorials › Statistical Learning (Graduate) › VC Entropy Integrals and Rademacher Bounds

VC and Entropy Complexity · Tutorial 166 of 989

VC Entropy Integrals and Rademacher Bounds

Learn how multiscale covering numbers control Rademacher complexity and how VC-subgraph entropy yields a uniform complexity bound.

Intermediate 10 min read

What You'll Learn

  • Define empirical Rademacher complexity and the sample-dependent \(L_2\) metric.
  • Derive a Dudley entropy-integral bound by chaining finite covers.
  • Convert VC-subgraph entropy bounds into \(L_2\) covering bounds for bounded functions.
  • Integrate the VC entropy envelope to obtain a dimension-dependent rate.
  • Check numerical entropy calculations and avoid common scale and normalization errors.

Why Covering Numbers Enter Rademacher Bounds

A class can contain many functions and still have modest statistical complexity if, at the observed sample, its functions are close to one another. Covering numbers quantify how many distinct approximations are needed at each resolution. A Dudley entropy integral adds the contributions across resolutions: coarse scales cost little, while fine scales matter only to the extent that many distinguishable functions remain.

Here the metric is computed on the observed sample. This makes the resulting bound conditional on the data and connects the entropy calculation directly to empirical Rademacher complexity. For a uniformly bounded VC-subgraph class, a uniform covering-number envelope makes the integral finite and gives the familiar square-root dependence on dimension and sample size.

Definition: For a fixed sample \(S=(X_1,\ldots,X_n)\) and a class \(\mathcal F\) of real-valued functions, define the empirical \(L_2\) pseudometric $$ d_S(f,g)=\left(\frac{1}{n}\sum_{i=1}^n(f(X_i)-g(X_i))^2\right)^{1/2}. $$ A covering number \(N(\epsilon,\mathcal F,d_S)\) is the smallest number of \(d_S\)-balls of radius \(\epsilon\), with centers in \(\mathcal F\), whose union contains \(\mathcal F\). It is a pseudometric because distinct functions can agree on all observed sample points. The empirical Rademacher complexity is $$ \widehat{\mathfrak R}_S(\mathcal F) = E_\sigma\left[\sup_{f\in\mathcal F}\frac{1}{n}\sum_{i=1}^n \sigma_i f(X_i)\right], $$ where \(\sigma_1,\ldots,\sigma_n\) are independent random signs, each equal to \(1\) or \(-1\) with probability \(1/2\).

The Dudley Entropy Integral

For each fixed \(f\), the signed sum in the definition has mean zero. Its increments across functions are controlled by the empirical distance: for every real \(\lambda\), $$ E_\sigma\exp\left\{\lambda\left[\frac{1}{n}\sum_{i=1}^n\sigma_i(f(X_i)-g(X_i))\right]\right\} \leq \exp\left(\frac{\lambda^2 d_S(f,g)^2}{2n}\right). $$ This follows by independence and \(\cosh(u)\leq \exp(u^2/2)\). Thus close functions have small random increments, and a sequence of increasingly fine covers can control the supremum.

Theorem (Dudley entropy-integral bound): Suppose the vectors \((f(X_1),\ldots,f(X_n))\), \(f\in\mathcal F\), form a bounded, separable subset of \(\mathbb R^n\). Let \(D\) be its \(d_S\)-diameter. If the entropy integral is finite, then $$ \widehat{\mathfrak R}_S(\mathcal F) \leq \frac{12}{\sqrt n}\int_0^D \sqrt{\log N(\epsilon,\mathcal F,d_S)}\,d\epsilon. $$

Here is the chaining argument for this constant. If \(D=0\), all functions have the same sample vector, so the expected supremum is zero. Otherwise, let \(\epsilon_k=D2^{-k}\) for \(k=0,1,\ldots\). Choose a finite \(\epsilon_k\)-cover \(T_k\) of \(\mathcal F\), with centers in the class. At \(k=0\), one center suffices because the diameter is \(D\). For each \(f\), choose an approximant \(\pi_k(f)\in T_k\) within \(\epsilon_k\) of \(f\). Then $$ d_S(\pi_k(f),\pi_{k-1}(f))\leq \epsilon_k+\epsilon_{k-1}=3\epsilon_k. $$ The signed process converges along these approximants to its value at \(f\), since its sample-vector linear functional is continuous in \(d_S\).

Write \(Z_f=n^{-1}\sum_i\sigma_i f(X_i)\). The expected maximum of \(M\) mean-zero variables whose moment generating functions are bounded by \(\exp(\lambda^2s^2/2)\) is at most \(s\sqrt{2\log M}\). This follows by bounding the exponential of the maximum by the sum of the \(M\) exponentials, taking logarithms and expectations, and minimizing over \(\lambda>0\). For the increments \(Z_{\pi_k(f)}-Z_{\pi_{k-1}(f)}\), the increment scale is at most \(3\epsilon_k/\sqrt n\), and there are at most \(|T_k||T_{k-1}|\) possible pairs. With \(H(\epsilon)=\log N(\epsilon,\mathcal F,d_S)\), monotonicity of covering numbers gives $$ E_\sigma\sup_f\left(Z_{\pi_k(f)}-Z_{\pi_{k-1}(f)}\right) \leq \frac{6\epsilon_k}{\sqrt n}\sqrt{H(\epsilon_k)}. $$ The initial center has expected process value zero. Summing the increments, and using that \(H(u)\geq H(\epsilon_k)\) for \(\epsilon_k/2\leq u\leq\epsilon_k\), yields $$ \sum_{k\geq1}\epsilon_k\sqrt{H(\epsilon_k)} \leq 2\int_0^D\sqrt{H(u)}\,du. $$ Combining these inequalities proves the theorem. Boundedness and separability ensure the approximating process converges to the supremum; finiteness of the integral justifies the sum. \(\square\)

Worked Example: A Two-Function Class

Suppose that on a fixed sample the class consists of \(f_0\) and \(f_1\), with \(d_S(f_0,f_1)=1\). Its diameter is \(D=1\). For \(0<\epsilon<1\), two balls are required, while at radius \(1\) one ball centered at either function covers both. Hence, apart from endpoints that do not affect the integral, $$ N(\epsilon,\mathcal F,d_S)=2 \quad (0<\epsilon<1), \qquad \int_0^1\sqrt{\log N(\epsilon,\mathcal F,d_S)}\,d\epsilon =\sqrt{\log 2}. $$ Dudley’s bound is therefore \(12\sqrt{\log 2}/\sqrt n\), approximately \(9.991/\sqrt n\). The exact complexity can be smaller: the integral bound pays for a chain that works generally, not for a class-specific exact calculation.

From VC-Subgraph Entropy to a Complexity Rate

For a real-valued function \(f\), its subgraph is the set of pairs \((x,t)\) satisfying \(t<f(x)\). The VC-subgraph dimension of \(\mathcal F\) is the VC dimension of the family of these sets. The standard uniform VC entropy bound says that a class of sets with VC dimension \(V\) has covering numbers bounded by a power of the reciprocal scale, uniformly over probability measures. For bounded functions, the subgraph representation transfers this control to function covers.

Theorem (VC-subgraph entropy integral): Let \(\mathcal F\) be a measurable class of functions taking values in \([0,1]\), with VC-subgraph dimension \(V\geq1\). Suppose its uniform \(L_2(Q)\) covering numbers satisfy $$ N(\epsilon,\mathcal F,L_2(Q)) \leq \left(\frac{A}{\epsilon}\right)^{2V}, \qquad 0<\epsilon\leq1, $$ for every probability measure \(Q\), where \(A\geq1\) is a universal constant. Then, for every fixed sample, $$ \widehat{\mathfrak R}_S(\mathcal F) \leq \frac{12\sqrt V}{\sqrt n} \left(\sqrt{2\log A}+\sqrt{\frac{\pi}{2}}\right). $$

The entropy premise is a standard consequence of the VC-subgraph covering bound. To see why its scale has this form, for functions in \([0,1]\) associate to \(f\) the indicator of its subgraph restricted to \(0\leq t\leq1\). Under \(Q\) on \(x\) and an independent uniform \(t\) on \([0,1]\), the \(L_1\) distance between these indicators is \(E_Q|f-g|\). Since \(|f-g|^2\leq|f-g|\), an \(L_1\) cover at radius \(\epsilon^2\) gives an \(L_2\) cover at radius \(\epsilon\). The uniform VC set entropy bound therefore gives an \(L_2\) covering bound with exponent \(2V\), absorbing fixed constants into \(A\).

To derive the displayed rate, \(d_S\) has diameter at most \(1\), because all sample values lie in \([0,1]\). The uniform covering premise applies to the empirical measure on the observed sample, so $$ \log N(\epsilon,\mathcal F,d_S) \leq 2V\log(A/\epsilon). $$ Substitute this envelope in Dudley’s theorem and use \(\sqrt{a+b}\leq\sqrt a+\sqrt b\): $$ \int_0^1\sqrt{2V\log(A/\epsilon)}\,d\epsilon \leq \sqrt{2V}\left(\sqrt{\log A}+\int_0^1\sqrt{\log(1/\epsilon)}\,d\epsilon\right). $$ With \(u=\sqrt{\log(1/\epsilon)}\), the last integral is \(2\int_0^\infty u^2e^{-u^2}\,du=\sqrt{\pi}/2\). Thus the right side is \(\sqrt V(\sqrt{2\log A}+\sqrt{\pi/2})\); equivalently, retaining the sharper separation before multiplying constants gives the stated bound. More directly, \(\sqrt{2V\log(A/\epsilon)}\leq\sqrt{2V\log A}+\sqrt{2V\log(1/\epsilon)}\), which leads to \(12\sqrt V(\sqrt{2\log A}+\sqrt{\pi/2})/\sqrt n\). We use this latter, directly justified expression as the bound: $$ \widehat{\mathfrak R}_S(\mathcal F) \leq \frac{12\sqrt V}{\sqrt n} \left(\sqrt{2\log A}+\sqrt{\frac{\pi}{2}}\right). $$

The VC dimension enters as \(\sqrt V\), and the sample size as \(1/\sqrt n\). These are rate statements; the universal constants can make the numerical bound loose. If functions instead take values in \([-B,B]\), an affine rescaling to \([0,1]\) gives the corresponding bound multiplied by \(2B\). The additive constant introduced by the rescaling has zero expected Rademacher contribution.

Worked Example: Applying a VC Entropy Envelope

Suppose \(\mathcal F\) takes values in \([0,1]\), has VC-subgraph dimension \(V=3\), and its uniform covers obey the stated premise with \(A=4\). For \(n=10{,}000\), the bound just derived gives $$ \widehat{\mathfrak R}_S(\mathcal F) \leq \frac{12\sqrt3}{100} \left(\sqrt{2\log4}+\sqrt{\frac{\pi}{2}}\right). $$ Using \(\log4\approx1.3863\), the two square roots are approximately \(1.6651\) and \(1.2533\). The prefactor is \(12\sqrt3/100\approx0.20785\), so the bound is approximately \(0.607\). Every scale from zero to the sample diameter contributes to this calculation; replacing the entropy by its value at one arbitrary scale would not give this integral.

Worked Example: Using an Integrated Entropy Calculation

Suppose a scale-dependent covering calculation on a sample of size \(n=100\) gives a diameter of at most \(1\) and $$ \int_0^1\log N(\epsilon,\mathcal F,d_S)\,d\epsilon=2.0986. $$ Cauchy–Schwarz on the interval of length one implies $$ \int_0^1\sqrt{\log N(\epsilon,\mathcal F,d_S)}\,d\epsilon \leq \sqrt{\int_0^1\log N(\epsilon,\mathcal F,d_S)\,d\epsilon} =\sqrt{2.0986}. $$ Dudley’s bound is therefore at most $$ \frac{12}{\sqrt{100}}\sqrt{2.0986} =1.2\sqrt{2.0986} \approx1.73839. $$ To three decimal places this is \(1.738\), or \(1.74\) to two decimal places. This is a valid but loose upper bound; Rademacher complexities can be bounded more sharply by other arguments or by using a tighter entropy calculation.

Common Pitfall: Confusing Entropy and Its Square Root

The integrand in Dudley’s bound is \(\sqrt{\log N(\epsilon)}\), not \(\log N(\epsilon)\). Integrating the entropy itself and then treating that value as the Dudley integral changes the bound. Cauchy–Schwarz can convert an integrated entropy into an upper bound on the Dudley integral, as in the numerical example, but the square root and the interval length must be handled explicitly. Another frequent error is to use a covering number from a population metric without checking that it controls the empirical metric; the uniform-in-\(Q\) VC entropy premise permits the empirical measure, whereas a bound for only one fixed distribution may not.

The entropy integral is useful because it records how complexity varies with resolution rather than treating the class as simply finite or infinite. A class may have a large fine-scale cover but still have a manageable integral if those scales occupy a short interval. Conversely, a coarse covering estimate alone can miss the accumulation of complexity across finer scales. Uniform VC-subgraph entropy supplies a convenient envelope that makes this accumulation explicit and yields a distribution-free complexity rate.

Check Your Understanding

Answer each question by referring to the relevant metric, scale, or bound.

  1. Why is \(d_S\) a pseudometric rather than necessarily a metric on functions?
  2. In Dudley’s bound, what quantity is integrated over scale: the covering number, its logarithm, or the square root of its logarithm?
  3. Why can the uniform \(L_2(Q)\) covering bound for a function class be applied with \(Q\) equal to the empirical distribution on the sample?
  4. For a class in \([0,1]\), what upper bound can be used for the empirical \(L_2\) diameter, and why?
  5. If \(\int_0^1\log N(\epsilon)\,d\epsilon=H\), state the Cauchy–Schwarz bound for \(\int_0^1\sqrt{\log N(\epsilon)}\,d\epsilon\).
Clinical Biostats