第 8 章精校翻译:链式方法
第 8 章链式方法
本章介绍界定随机过程 $(X_t)_{t\in T}$ 的一些核心方法。第 8.1 节引入链式方法,并用它通过 $T$ 的覆盖数来界定 $\mathbb E\sup_{t\in T}X_t$;这个结果称为 Dudley 不等式。第 8.2 节用 Monte Carlo 积分和经验过程作为应用来说明它。
随后,第 8.3 节引入 VC 理论。它给出的不是度量视角,而是组合视角下对随机过程的理解。第 8.4 节把它应用于统计学习理论。
第 8.5 节把链式方法精炼为泛型链式方法,并用 Talagrand 的 $\gamma_2(T)$ 泛函得到随机过程的最优双侧界(本书这里只证明上界)。一个有用的后果是 Talagrand 比较不等式,它把 Sudakov-Fernique 不等式推广到次高斯过程。
最后,第 8.6 节用 Talagrand 比较不等式推出 Chevet 不等式,这是处理一般集合上随机双线性型的一个方便工具。
学完这些内容后,有些习题会像有趣的小型研究结果:高维中的 Lipschitz 大数定律(Exercise 8.10)、一比特量化(Exercise 8.26)、带有重尾随机矩阵应用的小球方法(Exercise 8.27),以及随机矩阵的 $p\to q$ 范数(Exercise 8.41),这里只列举其中几个。
8.1 Dudley 不等式
对一般高斯过程 $(X_t)_{t\in T}$,Sudakov 不等式(Theorem 7.4.1)给出
$$ \mathbb E\sup_{t\in T}X_t $$
关于 $T$ 的度量熵的下界。现在我们将追求一个上界。而且我们不会只停留在高斯过程上,而是会进一步处理更一般的次高斯过程。
设 $(X_t)_{t\in T}$ 是度量空间 $(T,d)$ 上的随机过程。若存在 $K\ge0$ 使得对所有 $t,s\in T$,
$$\|X_t-X_s\|_{\psi_2}\le K d(t,s), \tag{8.1}$$则称该过程具有次高斯增量。
高斯过程的典范度量为 $d(t,s)=\|X_t-X_s\|_{L^2}$。由于中心化高斯的 $\psi_2$ 范数与标准差等价,高斯过程自动具有次高斯增量。
还有一个平凡例子:任意随机过程都可以通过定义
$$ d(t,s):=\|X_t-X_s\|_{\psi_2} $$而被强制变成具有次高斯增量的过程。
现在给出一个关于一般次高斯随机过程 $(X_t)_{t\in T}$ 的上界,用度量熵
$$ \log\mathcal N(T,d,\varepsilon) $$
来控制它。
设 $(X_t)_{t\in T}$ 是度量空间 $(T,d)$ 上均值为零的随机过程,并且具有 (8.1) 中的次高斯增量。那么
$$ \mathbb E\sup_{t\in T}X_t \le CK\int_0^\infty \sqrt{\log\mathcal N(T,d,\varepsilon)}\,d\varepsilon. $$ 查看学习笔记:Theorem 8.1.3 完整证明在证明 Dudley 不等式之前,先把它与 Sudakov 不等式(Theorem 7.4.1)比较。对高斯过程,Sudakov 不等式说
$$ \mathbb E\sup_{t\in T}X_t \ge c\sup_{\varepsilon\gt 0} \varepsilon\sqrt{\log\mathcal N(T,d,\varepsilon)}. $$
Figure 8.1 同时展示这两个界。它们之间有明显间隙;后面会看到,仅靠度量熵一般无法消除这个间隙。
Dudley 不等式暗示 $\mathbb E\sup_{t\in T}X_t$ 是一个多尺度量:为了控制它,需要在所有尺度 $\varepsilon$ 上观察 $T$。证明也正是这样进行的。我们先用二进尺度 $\varepsilon=2^{-k}$ 证明一个离散版本(类似 Riemann 求和),然后再转到连续版本。
设 $(X_t)_{t\in T}$ 是度量空间 $(T,d)$ 上均值为零的随机过程,并且具有 (8.1) 中的次高斯增量。那么
$$ \mathbb E\sup_{t\in T}X_t \le CK\sum_{k\in\mathbb Z} 2^{-k}\sqrt{\log\mathcal N(T,d,2^{-k})}. \tag{8.2} $$ 查看学习笔记:Theorem 8.1.4 完整证明这个证明使用链式方法技术。它本质上是我们此前用过的 $\varepsilon$-net 论证的多尺度版本,例如 Theorem 4.4.3 和 Theorem 7.6.1 的证明中都出现过类似思想。
在 $\varepsilon$-net 论证中,我们用一个 $\varepsilon$-net $\mathcal N$ 近似 $T$,使每个 $t\in T$ 都接近某个 $\pi(t)\in\mathcal N$,并满足 $d(t,\pi(t))\le\varepsilon$。于是增量条件 (8.1) 给出
$$ \|X_t-X_{\pi(t)}\|_{\psi_2}\le K\varepsilon. \tag{8.3} $$
这导出
$$ \mathbb E\sup_{t\in T}X_t \le \mathbb E\sup_{t\in T}X_{\pi(t)} + \mathbb E\sup_{t\in T}\bigl(X_t-X_{\pi(t)}\bigr). $$
第一项可以在 $|\mathcal N|=\mathcal N(T,d,\varepsilon)$ 个点 $\pi(t)$ 上用并集界处理。但第二项比较麻烦:不清楚如何把 (8.3) 与对所有 $t\in T$ 的并集界结合。为了解决这个问题,我们不在单个网处停止,而是选择越来越小的 $\varepsilon$,得到 $\pi_1(t),\pi_2(t),\ldots$ 这些越来越精细的近似点。这就是链式方法的思想;下面正式写出证明。
步骤 1:链式方法集合。 不失一般性,可以假设 $K=1$(为什么?)并且 $T$ 是有限集(见 Remark 7.2.1)。定义二进尺度
$$ \varepsilon_k=2^{-k}, \qquad k\in\mathbb Z, \tag{8.4} $$并选取 $T$ 的 $\varepsilon_k$-nets $T_k$,使得
$$ |T_k|=\mathcal N(T,d,\varepsilon_k). \tag{8.5} $$只需要二进尺度的一部分。因为 $T$ 有限,存在足够小的整数 $\kappa$(定义最粗网)和足够大的整数 $K$(定义最细网),使得
$$ T_\kappa=\{t_0\}\quad\text{for some }t_0\in T, \qquad T_K=T. \tag{8.6} $$对点 $t\in T$,令 $\pi_k(t)$ 表示 $T_k$ 中离 $t$ 最近的点,于是
$$ d(t,\pi_k(t))\le\varepsilon_k. \tag{8.7} $$由假设 $\mathbb EX_{t_0}=0$,有
$$ \mathbb E\sup_{t\in T}X_t = \mathbb E\sup_{t\in T}(X_t-X_{t_0}). $$把 $X_t-X_{t_0}$ 写成望远镜求和:沿着由 $\pi_k(t)$ 给出的链,从固定点 $t_0$ 逐步走向 $t$:
$$ \begin{aligned} X_t-X_{t_0} &= \bigl(X_{\pi_\kappa(t)}-X_{t_0}\bigr) +\bigl(X_{\pi_{\kappa+1}(t)}-X_{\pi_\kappa(t)}\bigr) +\cdots\\ &\qquad +\bigl(X_t-X_{\pi_K(t)}\bigr). \end{aligned} \tag{8.8} $$见 Figure 8.2。这个和式的第一项与最后一项由 (8.6) 为零,因此
$$ X_t-X_{t_0} = \sum_{k=\kappa+1}^{K} \bigl(X_{\pi_k(t)}-X_{\pi_{k-1}(t)}\bigr). \tag{8.9} $$由于和的上确界不超过上确界之和,得到
$$ \mathbb E\sup_{t\in T}(X_t-X_{t_0}) \le \sum_{k=\kappa+1}^{K} \mathbb E\sup_{t\in T} \bigl(X_{\pi_k(t)}-X_{\pi_{k-1}(t)}\bigr). \tag{8.10} $$
步骤 2:控制增量。 在 (8.10) 中,看起来每个求和项都在整个 $T$ 上取上确界;但实际上,它只在较小的点对集合 $(\pi_k(t),\pi_{k-1}(t))$ 上取上确界。这类点对的数量为
$$ |T_k|\cdot |T_{k-1}|\le |T_k|^2, $$而这个数量可以通过 (8.5) 控制。另一方面,对固定的 $t$,(8.10) 中的增量可按如下方式估计:
$$ \begin{aligned} \|X_{\pi_k(t)}-X_{\pi_{k-1}(t)}\|_{\psi_2} &\le d(\pi_k(t),\pi_{k-1}(t)) &&\text{by (8.1), since }K=1\\ &\le d(\pi_k(t),t)+d(t,\pi_{k-1}(t)) &&\text{by 三角不等式}\\ &\le \varepsilon_k+\varepsilon_{k-1} &&\text{by (8.7)}\\ &\le 2\varepsilon_{k-1}. \end{aligned} $$回忆 (2.22):$N$ 个次高斯随机变量的期望最大值至多为 $CL\sqrt{\log N}$,其中 $L$ 是最大的 $\psi_2$ 范数。把它用于 (8.10) 中每一项,得到
$$ \mathbb E\sup_{t\in T} \bigl(X_{\pi_k(t)}-X_{\pi_{k-1}(t)}\bigr) \le C\varepsilon_{k-1}\sqrt{\log |T_k|}. \tag{8.11} $$步骤 3:累加增量。 我们已经证明
$$ \mathbb E\sup_{t\in T}(X_t-X_{t_0}) \le C\sum_{k=\kappa+1}^{K} \varepsilon_{k-1}\sqrt{\log |T_k|}. \tag{8.12} $$现在代入 (8.4) 中 $\varepsilon_k=2^{-k}$,以及 (8.5) 中对 $|T_k|$ 的界,得到
$$ \mathbb E\sup_{t\in T}(X_t-X_{t_0}) \le C_1\sum_{k=\kappa+1}^{K} 2^{-k}\sqrt{\log\mathcal N(T,d,2^{-k})}. $$Theorem 8.1.4 得证。
现在由离散形式推出 Dudley 不等式的积分形式。
为了把 (8.2) 中的求和转成积分,把 $2^{-k}$ 写成
$$ 2^{-k} = 2\int_{2^{-k-1}}^{2^{-k}}d\varepsilon. $$于是
$$ \begin{aligned} \sum_{k\in\mathbb Z} 2^{-k}\sqrt{\log\mathcal N(T,d,2^{-k})} &= 2\sum_{k\in\mathbb Z} \int_{2^{-k-1}}^{2^{-k}} \sqrt{\log\mathcal N(T,d,2^{-k})}\,d\varepsilon. \end{aligned} $$在积分区间内,$2^{-k}\ge\varepsilon$,所以
$$ \log\mathcal N(T,d,2^{-k}) \le \log\mathcal N(T,d,\varepsilon). $$因此这个求和被下面的量控制:
$$ 2\sum_{k\in\mathbb Z} \int_{2^{-k-1}}^{2^{-k}} \sqrt{\log\mathcal N(T,d,\varepsilon)}\,d\varepsilon = 2\int_0^\infty \sqrt{\log\mathcal N(T,d,\varepsilon)}\,d\varepsilon. $$Theorem 8.1.3 得证。
作为练习,请证明离散和积分形式的 Dudley 不等式实际上在常数因子意义下等价(Exercise 8.3)。
8.1.1 变体与例子
快速回看证明可以看出,链式方法实际上给出增量上确界版本。对任意固定 $t_0\in T$,有
$$ \mathbb E\sup_{t\in T}|X_t-X_{t_0}| \le CK\int_0^\infty \sqrt{\log\mathcal N(T,d,\varepsilon)}\,d\varepsilon. \tag{8.13} $$再由三角不等式可得二元增量版本:
$$ \mathbb E\sup_{t,s\in T}|X_t-X_s| \le CK\int_0^\infty \sqrt{\log\mathcal N(T,d,\varepsilon)}\,d\varepsilon. \tag{8.14} $$Dudley 不等式只给出期望界,但链式方法实际上给出高概率界。假设 $T$ 有限,则对每个 $u\ge0$,都有
$$ \sup_{t,s\in T}|X_t-X_s| \le CK\left[ \int_0^\infty\sqrt{\log\mathcal N(T,d,\varepsilon)}\,d\varepsilon +u\cdot\operatorname{diam}(T) \right] \tag{8.15} $$以至少 $1-2\exp(-u^2)$ 的概率成立(见 Exercise 8.1)。对高斯过程,这也可以直接由高斯集中推出(见 Exercise 8.2)。
注意,(8.14) 和 (8.15) 不需要均值零假设 $\mathbb EX_t=0$;但 Dudley 定理 8.1.3 需要这个假设,否则结论可能失败。你可以试着找出原因。
查看学习笔记:Exercise 8.1 证明虽然 Dudley 积分写成在 $[0,\infty)$ 上积分,但上限可以截到 $\operatorname{diam}(T)$:当 $\varepsilon>\operatorname{diam}(T)$ 时,一个 $\varepsilon$-球就覆盖 $T$,所以 $\log\mathcal N(T,d,\varepsilon)=0$。因此
$$ \mathbb E\sup_{t\in T}X_t \le CK\int_0^{\operatorname{diam}(T)} \sqrt{\log\mathcal N(T,d,\varepsilon)}\,d\varepsilon. \tag{8.16} $$ 查看学习笔记:积分上限截断验证任意有界集合 $T\subset\mathbb R^n$ 的高斯宽度满足
$$ w(T)\le C\int_0^\infty\sqrt{\log\mathcal N(T,\varepsilon)}\,d\varepsilon. \tag{8.17} $$这里 $\mathcal N(T,\varepsilon)$ 表示以 $T$ 中点为中心、半径为 $\varepsilon$ 的欧氏球覆盖 $T$ 所需的最少球数。
查看学习笔记:Theorem 8.1.8 完整证明把 Dudley 不等式测试在欧氏单位球 $T=B_2^n$ 上。由体积法,若 $\varepsilon\in(0,1]$,则
$$ \mathcal N(B_2^n,\varepsilon)\le\left(\frac3\varepsilon\right)^n, $$而若 $\varepsilon\gt 1$,一个球已经覆盖 $B_2^n$。因此
$$ w(B_2^n) \le C\int_0^1\sqrt{n\log(3/\varepsilon)}\,d\varepsilon \le C\sqrt n. $$这个上界与真实量 $w(B_2^n)=\mathbb E\|g\|_2\asymp\sqrt n$ 同阶,所以在欧氏球上 Dudley 不等式是尖锐的。
Dudley 积分在一般情形下可能高估高斯宽度。一个典型例子是
$$ T=\left\{ \frac{e_k}{\sqrt{1+\log k}}:k=1,\ldots,n \right\}. $$Exercise 8.4 会要求你证明:这个集合的高斯宽度有绝对界,但 Dudley 积分会随 $n$ 增大。
这种偏松不是无限坏的:Dudley 不等式和 Sudakov 不等式在最坏情况下只差一个对数因子。Exercise 8.5 会要求你证明这个对数意义下的尖锐性;而第 8.5 节会升级链式方法并去掉这一对数因子。Exercise 8.6 会改进 Dudley 积分的下限;Exercises 8.7 和 8.8 分别给出次指数 Dudley 与局部 Dudley 的版本。
8.2 应用:经验过程
8.2.1 Monte Carlo 方法
给定概率空间 $(\Omega,\mu)$ 和函数 $f:\Omega\to\mathbb R$,目标是计算积分
$$ \int_\Omega f\,d\mu=\mathbb Ef(X). $$
Monte Carlo 方法取独立样本 $X_1,\dots,X_n\sim\mu$。由大数定律(Theorem 1.7.1),
$$ \frac1n\sum_{i=1}^nf(X_i) \to \mathbb Ef(X) \quad\text{almost surely}. \tag{8.18} $$
因此可用
$$ \int_\Omega f\,d\mu \approx \frac1n\sum_{i=1}^nf(X_i) \tag{8.19} $$
来近似该积分。
Monte Carlo 估计 (8.19) 的期望误差是 $O(1/\sqrt n)$。这来自大数定律中的收敛速率 (1.23):
$$ \mathbb E\left| \frac1n\sum_{i=1}^n f(X_i)-\mathbb Ef(X) \right| \le \left[ \operatorname{Var}\left( \frac1n\sum_{i=1}^n f(X_i) \right) \right]^{1/2} = O\left(\frac1{\sqrt n}\right). \tag{8.20} $$Monte Carlo 的基本误差率不直接依赖维度;只要能从 $\mu$ 抽样并计算 $f(X_i)$,就可实施。
8.2.2 Lipschitz 大数定律
一个样本不能同时估计所有函数的积分。若函数可在样本点之间剧烈振荡,则经验平均会完全误判整体积分。
令
$$ \mathcal F=\{f:[0,1]\to\mathbb R:\|f\|_{\mathrm{Lip}}\le L\}. \tag{8.21} $$若 $X,X_1,\dots,X_n$ 是取值于 $[0,1]$ 的 i.i.d. 随机变量,则
$$ \mathbb E\sup_{f\in\mathcal F} \left| \frac1n\sum_{i=1}^nf(X_i)-\mathbb Ef(X) \right| \le \frac{CL}{\sqrt n}. \tag{8.22} $$ 查看学习笔记:Theorem 8.2.3 完整证明在证明之前,先重复关键点:左侧对 $f\in\mathcal F$ 的上确界位于期望内部。
因此,由 Markov 不等式可知,同一组样本 $X_1,\ldots,X_n$ 会以高概率同时对所有 $f\in\mathcal F$ 都工作良好。这里“工作良好”指的是以 $O(1/\sqrt n)$ 的误差近似每一个积分,这与只处理一个固定函数时通常的大数定律速率相同。所以,我们把大数定律变成了一致版本,而且没有损失速率。
为了让 Theorem 8.2.3 的证明更直观,可以把 (8.22) 左侧看成一个由函数 $f\in\mathcal F$ 索引的随机过程的最大幅度。这类随机过程称为经验过程。
设 $\mathcal F$ 是某个集合 $\Omega$ 上的一类实值函数 $f:\Omega\to\mathbb R$。令 $X$ 是按照某个 probability distribution 从 $\Omega$ 中抽取的随机点,并令 $X_1,X_2,\ldots,X_n$ 是 $X$ 的独立副本。由
$$ X_f := \frac1n\sum_{i=1}^nf(X_i)-\mathbb Ef(X) \tag{8.23} $$定义的随机过程 $(X_f)_{f\in\mathcal F}$ 称为由 $\mathcal F$ 索引的经验过程。
不失一般性,只需对函数类
$$ \mathcal F := \{f:[0,1]\to[0,1],\ \|f\|_{\mathrm{Lip}}\le1\} \tag{8.24} $$证明定理。(为什么?)我们要控制 (8.23) 中经验过程 $(X_f)_{f\in\mathcal F}$ 的量
$$ \mathbb E\sup_{f\in\mathcal F}|X_f|. $$步骤 1:检查次高斯增量。 使用 Dudley 不等式(Theorem 8.1.3)。为了应用它,需要检查经验过程关于 $L^\infty$ 度量
$$ d(f,g)=\|f-g\|_{L^\infty} $$具有次高斯增量。固定一对函数 $f,g\in\mathcal F$,写
$$ \|X_f-X_g\|_{\psi_2} = \frac1n\left\|\sum_{i=1}^nZ_i\right\|_{\psi_2}, \qquad Z_i:=(f-g)(X_i)-\mathbb E(f-g)(X). $$因为 $Z_i$ 是独立、均值为零随机变量,Proposition 2.7.1 给出
$$ \|X_f-X_g\|_{\psi_2} \lesssim \frac1n \left( \sum_{i=1}^n\|Z_i\|_{\psi_2}^2 \right)^{1/2}. $$现在用 centering(Lemma 2.7.8)得到
$$ \|Z_i\|_{\psi_2} \lesssim \|(f-g)(X_i)\|_{\psi_2} \lesssim \|f-g\|_{L^\infty}. $$因此
$$ \|X_f-X_g\|_{\psi_2} \lesssim \frac1n\cdot n^{1/2}\|f-g\|_{L^\infty} = \frac1{\sqrt n}\|f-g\|_{L^\infty}. $$步骤 2:应用 Dudley 不等式。 现在以 (8.13) 的形式应用 Dudley 不等式:
$$ \begin{aligned} \mathbb E\sup_{f\in\mathcal F}|X_f| &= \mathbb E\sup_{f\in\mathcal F}|X_f-X_0|\\ &\lesssim \frac1{\sqrt n} \int_0^1 \sqrt{\log\mathcal N(\mathcal F,\|\cdot\|_{L^\infty},\varepsilon)} \,d\varepsilon. \end{aligned} \tag{8.25} $$这里使用了零函数属于 $\mathcal F$,并且由 (8.24) 可知 $\mathcal F$ 在 $L^\infty$ 度量下的直径至多为 $1$。
函数类 $\mathcal F$ 的覆盖数不难估计为
$$ \mathcal N(\mathcal F,\|\cdot\|_{L^\infty},\varepsilon) \le e^{C/\varepsilon}, $$这会在 Exercise 8.9 中验证。把这个界代入积分,得到
$$ \mathbb E\sup_{f\in\mathcal F}|X_f| \lesssim \frac1{\sqrt n} \int_0^1\sqrt{\frac C\varepsilon}\,d\varepsilon \lesssim \frac1{\sqrt n}. $$Theorem 8.2.3 得证。
为了练习,可以尝试 Exercise 8.10,把 Theorem 8.2.3 推广到高维。这题难度较高,但不用担心:主要是计算工作量大,并不是特别取巧。
8.2.3 经验测度
为了获得更宽的视角,再看 Definition 8.2.5 中的经验过程。给定一组 i.i.d. 样本 $X_1,\ldots,X_n$,它们按照某个概率测度 $\mu$ 从 $\Omega$ 中抽取。考虑经验测度 $\mu_n$,它给每个样本点分配相同概率 $1/n$,并计入重数:
$$ \mu_n=\frac1n\sum_{i=1}^n\delta_{X_i}. \tag{8.26} $$
这里 $\delta_x$ 是 $x$ 处的 Dirac 概率测度;也就是说,对任意集合 $A$,若 $x\in A$ 则 $\delta_x(A)=1$,否则为 $0$。因此,$\mu_n(A)$ 是落入集合 $A\subset\Omega$ 的样本点比例。
函数 $f$ 关于原测度 $\mu$ 的积分是 $\mathbb Ef(X)$,也就是 $f$ 的总体平均;而 $f$ 关于经验测度 $\mu_n$ 的积分为
$$ \frac1n\sum_{i=1}^nf(X_i), $$
也就是 $f$ 的样本平均或经验平均。(8.23) 中的经验过程 $X_f$ 记录的是总体期望与经验期望的偏差。
Theorem 8.2.3 中控制的这个偏差,可以看成测度 $\mu$ 与 $\mu_n$ 之间的一种距离,称为 Wasserstein 距离 $W_1(\mu,\mu_n)$。它还有一个等价解释:把一个测度运输成另一个测度的运输成本。这个等价性由 Kantorovich-Rubinstein 对偶定理给出。因此 Theorem 8.2.3 常被称为 Wasserstein 大数定律。
你已经学到的许多工具也适用于经验过程。Exercise 8.11 会要求你证明对称化的经验版本。
8.3 VC 维数
8.3.1 定义与例子
若布尔函数类 $\mathcal F$ 能在有限集 $\Lambda\subset\Omega$ 上实现全部 $2^{|\Lambda|}$ 种二元标记,则称 $\Lambda$ 被 $\mathcal F$ 打散。$\operatorname{vc}(\mathcal F)$ 是可被打散的集合最大基数;若无最大值,则为 $\infty$。
令 $\mathcal F$ 是实线上所有闭区间的指示函数组成的类:
$$ \mathcal F=\{\mathbf 1_I:I\subset\mathbb R\text{ is a closed interval}\}. $$任意两个点都可以由区间实现全部四种二元标记,因此 $\operatorname{vc}(\mathcal F)\ge2$。但对三个有序点 $p\lt q\lt r$,标记 $g(p)=1,g(q)=0,g(r)=1$ 无法由任何区间实现:区间若包含 $p$ 和 $r$,就必然包含中间的 $q$。故不存在被该类打散的三点集,因此 $\operatorname{vc}(\mathcal F)=2$。
令 $\mathcal F$ 是 $\mathbb R^2$ 中所有闭半平面的指示函数组成的类。三个一般位置点可被打散:任意二元标记都可由某条直线分离。另一方面,对任意四个点,总能找到一种二元标记不是线性可分的。因此
$$ \operatorname{vc}(\mathcal F)=3. $$
令 $\Omega=\{1,2,3\}$。可以把 $\Omega$ 上的布尔函数看成长度为 $3$ 的二进制串。例如原书考虑类
$$ \mathcal F=\{001,010,100,111\}. $$集合 $\Lambda=\{1,3\}$ 被 $\mathcal F$ 打散。因为把 $\mathcal F$ 中的函数限制到 $\Lambda$ 上,等于删去第二位,得到 $00,01,10,11$,恰好实现全部二元标记。因此 $\operatorname{vc}(\mathcal F)\ge2$。
另一方面,三点集 $\{1,2,3\}$ 若被打散,则 $\mathcal F$ 中必须出现全部 $8$ 个长度为 $3$ 的二进制串,但事实并非如此,所以 $\operatorname{vc}(\mathcal F)\lt 3$。故 $\operatorname{vc}(\mathcal F)=2$。
这个例子用于展示被打散子集的计数方式,也会在 Pajor 引理的证明后再次出现。
$\mathbb R^n$ 中的半空间是如下形式的集合:
$$ \{x\in\mathbb R^n:\langle a,x\rangle\le b\}. $$半空间的指示函数类具有 VC 维数 $n+1$;证明留给 Exercise 8.17。
函数类的 VC 维数常常与参数个数同阶。例如 $\mathbb R^n$ 中的半空间由 $n+1$ 个参数描述,而 VC 维数也是 $n+1$。不过这只是一个有用启发,不是普遍定理。
8.3.2 Pajor 引理
假设定义域 $\Omega$ 是有限的,并由 $n$ 个点组成。那么 $\Omega$ 上任意布尔函数类 $\mathcal F$ 也是有限的,并且
$$ 2^{\operatorname{vc}(\mathcal F)} \le |\mathcal F| \le 2^n. \tag{8.27} $$
(为什么?)上界通常很松;大多数函数类的大小更接近下界。这一点并不显然。为了准备证明这个结果,先证明:$\Omega$ 中被 $\mathcal F$ 打散的子集至少和 $\mathcal F$ 中的函数一样多。
若 $\mathcal F$ 是有限集合 $\Omega$ 上的布尔函数类,则
$$|\mathcal F|\le |\{\Lambda\subseteq\Omega:\Lambda\text{ 被 }\mathcal F\text{ 打散}\}|.$$ 查看学习笔记:Lemma 8.3.7 完整证明右侧计数中包括空集 $\Lambda=\emptyset$。
<p>证明之前,先用 Example 8.3.4 说明这个结果。在那里 $|\mathcal F|=4$,并且有六个子集 $\Lambda$ 被 $\mathcal F$ 打散,分别是 $\{1\}$、$\{2\}$、$\{3\}$、$\{1,2\}$、$\{1,3\}$ 和 $\{2,3\}$。(检查!)因此 Pajor 引理中的不等式读作 $4\le6$。</p>
对 $\Omega$ 的基数做归纳。情形 $|\Omega|=1$ 是平凡的,因为右侧计数包含空集。假设引理对任意 $n$ 点集合 $\Omega$ 成立,现在证明 $|\Omega|=n+1$ 的情形。
从 $\Omega$ 中取出一个任意点,写作
$$ \Omega=\Omega_0\cup\{x_0\}, \qquad |\Omega_0|=n. $$函数类 $\mathcal F$ 自然分解成两个子类:
$$ \mathcal F_0:=\{f\in\mathcal F:\ f(x_0)=0\}, \qquad \mathcal F_1:=\{f\in\mathcal F:\ f(x_0)=1\}. $$由归纳假设,计数函数
$$ S(\mathcal F) = \left| \{\Lambda\subseteq\Omega:\ \Lambda\text{ 被 }\mathcal F\text{ 打散}\} \right| $$满足
$$ S(\mathcal F_0)\ge|\mathcal F_0|, \qquad S(\mathcal F_1)\ge|\mathcal F_1|. \tag{8.28} $$为了完成证明,只需检查
$$ S(\mathcal F) \ge S(\mathcal F_0)+S(\mathcal F_1). \tag{8.29} $$这样 (8.28) 就会给出
$$ S(\mathcal F) \ge |\mathcal F_0|+|\mathcal F_1| = |\mathcal F|, $$正是所需结论。
不等式 (8.29) 看起来似乎平凡。任何被 $\mathcal F_0$ 或 $\mathcal F_1$ 打散的集合 $\Lambda$,都自动被更大的类 $\mathcal F$ 打散,因此每个被 $S(\mathcal F_0)$ 或 $S(\mathcal F_1)$ 计入的 $\Lambda$ 都自动被 $S(\mathcal F)$ 计入。问题在于重复计数。
假设同一个集合 $\Lambda$ 同时被 $\mathcal F_0$ 和 $\mathcal F_1$ 打散。计数函数 $S(\mathcal F)$ 不会把 $\Lambda$ 计两次。不过,另一个集合会被 $S(\mathcal F)$ 计入,而它既没有被 $S(\mathcal F_0)$ 计入,也没有被 $S(\mathcal F_1)$ 计入;这个集合就是 $\Lambda\cup\{x_0\}$。稍作思考可知,该集合确实被 $\mathcal F$ 打散。(检查!)这证明了 (8.29),从而完成 Pajor 引理的证明。
下面说明 Pajor lemma 证明中的关键点。
回到 Example 8.3.4 中的 $\Omega=\{1,2,3\}$ 和
$$ \mathcal F=\{001,010,100,111\}. $$按照 Pajor 引理的证明,把 $x_0=3$ 从 $\Omega=\{1,2,3\}$ 中取出,于是 $\Omega_0=\{1,2\}$。类 $\mathcal F$ 分解为两个子类:
$$ \mathcal F_0=\{010,100\}, \qquad \mathcal F_1=\{001,111\}. $$被 $\mathcal F_0$ 打散的子集恰好有两个,即 $\{1\}$ 和 $\{2\}$;被 $\mathcal F_1$ 打散的也是同样两个子集。因此 $S(\mathcal F_0)=S(\mathcal F_1)=2$。当然,同样的两个子集也被 $\mathcal F$ 打散,但为了关键不等式 (8.29),还需要另外两个被打散子集,使 $S(\mathcal F)\ge4$。
构造方法是:把 $x_0=3$ 添加到已经计入的子集 $\Lambda$ 中。得到的集合 $\{1,3\}$ 和 $\{2,3\}$ 也被 $\mathcal F$ 打散,而且此前还没有被计数。这样就至少有四个子集被 $\mathcal F$ 打散,使关键不等式 (8.29) 成立。
8.3.3 Sauer-Shelah 引理
现在从 Pajor 引理推出一个重要结论:函数类的基数可由 VC 维数控制。
若 $\mathcal F$ 是 $n$ 点集合 $\Omega$ 上的布尔函数类,则
$$ |\mathcal F| \le \sum_{k=0}^d\binom nk \le \left(\frac{en}{d}\right)^d, \qquad d=\operatorname{vc}(\mathcal F). $$ 查看学习笔记:Lemma 8.3.9 完整证明Pajor 引理说明 $|\mathcal F|$ 被 $\Omega$ 中被 $\mathcal F$ 打散的子集 $\Lambda$ 的数量控制。按 VC 维数的定义,每个这样的 $\Lambda$ 的基数至多为 $d=\operatorname{vc}(\mathcal F)$。因此
$$ |\mathcal F| \le \left| \{\Lambda\subseteq\Omega:\ |\Lambda|\le d\} \right| = \sum_{k=0}^d\binom nk, $$因为右侧的和正是 $n$ 元集合中基数不超过 $d$ 的子集总数。这证明 Sauer-Shelah 引理的第一个不等式。第二个不等式来自 Exercise 0.6 中证明过的二项和界。
Pajor 引理与 Sauer-Shelah 引理通常都是尖锐的;见 Exercise 8.19。
8.3.4 增长函数
Sauer-Shelah lemma 假设定义域 $\Omega$ 有限。那么,当函数类 $\mathcal F$ 定义在无限定义域上,例如 $\Omega=\mathbb R^n$ 时,应如何处理?通常方便的做法是用增长函数衡量 $\mathcal F$ 的复杂度。
设 $\mathcal F$ 是定义域 $\Omega$ 上的布尔函数类。$\mathcal F$ 的增长函数定义为:把 $\mathcal F$ 中所有函数限制到一个 $n$ 点子集上时,最多能得到多少个不同函数:
$$ \Pi_{\mathcal F}(n) = \sup\{ \left|\mathcal F|_{\Lambda}\right|:\ \Lambda\subset\Omega,\ |\Lambda|=n \}. $$从这个角度看,$\mathcal F$ 的 VC 维数就是使 $\Pi_{\mathcal F}(d)=2^d$ 成立的最大 $d$。若 $d=\operatorname{vc}(\mathcal F)\lt \infty$,则增长函数有直接界:
$$ 2^d \le \Pi_{\mathcal F}(n) \le \left(\frac{en}{d}\right)^d \qquad\text{for all }n\ge d. \tag{8.30} $$
下界是 (8.27) 的重新表述,上界来自 Sauer-Shelah Lemma(Lemma 8.3.9)。
为了说明增长函数的用处,下面从 (8.30) 推出 VC 维数对自然运算的稳定性。
令 $\mathcal F$ 和 $\mathcal G$ 是同一定义域上的两个布尔函数类。令
$$ \mathcal F\wedge\mathcal G = \{f\wedge g:\ f\in\mathcal F,\ g\in\mathcal G\}, $$其中 $f\wedge g$ 表示函数 $f$ 与 $g$ 的逐点最小值。那么
$$ \operatorname{vc}(\mathcal F\wedge\mathcal G) \le 10\max\bigl( \operatorname{vc}(\mathcal F), \operatorname{vc}(\mathcal G) \bigr). $$同样的结论也适用于逐点最大值 $\vee$。
查看学习笔记:Proposition 8.3.11 完整证明反证。设
$$ n:=\operatorname{vc}(\mathcal F\wedge\mathcal G)\gt 10d, \qquad d:=\max(\operatorname{vc}(\mathcal F),\operatorname{vc}(\mathcal G)). $$那么
$$ 2^n \le \Pi_{\mathcal F\wedge\mathcal G}(n) \le \Pi_{\mathcal F}(n)\cdot\Pi_{\mathcal G}(n) \le \left(\frac{en}{d}\right)^{2d}. $$第一个和最后一个界来自 (8.30),中间的界由定义直接可见。然而简单计算表明,只要 $n\gt 10d$,就有 $2^n\gt (en/d)^{2d}$。矛盾。
Proposition 8.3.11 可以推广到任何给定的函数类组合方式;见 Exercise 8.21。当我们不想直接计算 VC 维数,而只想给出上界时,它很有用。下面是一个例子。
$\mathbb R^n$ 中的 strip 可写成
$$ \{x\in\mathbb R^n:|\langle a,x\rangle-b|\le c\}. $$见 Figure 8.7。令 $\mathcal F$ 是条带的指示函数类。Example 8.3.5 给出
$$ \operatorname{vc}(\mathcal F)\le20(n+1)\le40n. $$原因如下。每个 strip 都可以表示成两个半空间
$$ \langle a,x\rangle-b\le c, \qquad \langle a,x\rangle-b\ge -c. $$的交。因此,每个带状区域的指示函数都是两个半空间指示函数的逐点最小值。现在应用 VC 稳定性(Proposition 8.3.11)以及 Example 8.3.5 的结果。
为了继续练习增长函数,可以证明一个很惊人的二分现象:增长函数只能多项式增长或指数增长,不存在中间增长(Exercise 8.20);也可以计算并类的 VC 维数(Exercise 8.22)。
8.3.5 通过 VC 维数控制覆盖数
覆盖数通常会随维数指数增长,例如 (4.17)。现在把这个启发再精炼一下:用 VC 维数替代代数维数,这能节省很多。
令 $\mathcal F$ 是某定义域 $\Omega$ 上的布尔函数类,令 $\mu$ 是 $\Omega$ 上任意概率测度。定义函数之间的距离:
$$ d(f,g) = \|f-g\|_{L^2(\mu)} = \bigl(\mathbb E(f-g)(X)^2\bigr)^{1/2}, \tag{8.31} $$
其中 $X$ 的分布为 $\mu$。(如果你没有学过测度论,就把 $X$ 看作任意取值于 $\Omega$ 的随机变量,并把 $\mu$ 理解为 $X$ 的分布名称。)
现在用 VC 维数控制 $\mathcal F$ 关于度量 (8.31) 的覆盖数 $\mathcal N(\mathcal F,L^2(\mu),\varepsilon)$。
设 $\mathcal F$ 是定义域 $\Omega$ 上的布尔函数类,并且 $\Omega$ 上有概率测度 $\mu$。那么,对任意 $\varepsilon\in(0,1)$,
$$ \mathcal N(\mathcal F,L^2(\mu),\varepsilon) \le \left(\frac2\varepsilon\right)^{Cd}, \qquad d=\operatorname{vc}(\mathcal F). $$ 查看学习笔记:Theorem 8.3.13 完整证明作为证明 Theorem 8.3.13 的第一次尝试,先假设 $\Omega$ 是有限集,比如 $|\Omega|=n$。Sauer-Shelah Lemma 8.3.9 给出
$$ \mathcal N(\mathcal F,L^2(\mu),\varepsilon) \le |\mathcal F| \le \left(\frac{en}{d}\right)^d. $$
这还不是 Theorem 8.3.13,但已经接近了。为了强化这个界,需要去掉 $n$;做法是缩小 $\Omega$。下面的引理会帮助我们。
设 $\mathcal F$ 是定义域 $\Omega$ 上有限布尔函数类,且 $\Omega$ 上有概率测度 $\mu$。假设 $\mathcal F$ 中所有函数都是 $\varepsilon$-分离的,即
$$ \|f-g\|_{L^2(\mu)}\gt \varepsilon \qquad \text{for all distinct }f,g\in\mathcal F. $$如果 $n\ge C\varepsilon^{-4}\log|\mathcal F|$,则经验测度 $\mu_n$ 以至少 $0.99$ 的概率满足
$$ \|f-g\|_{L^2(\mu_n)}\gt \varepsilon/2 \qquad \text{for all distinct }f,g\in\mathcal F. $$ 查看学习笔记:Lemma 8.3.14 完整证明由经验测度的定义(见第 8.2.3 节),$\|f-g\|_{L^2(\mu_n)}$ 与 (8.31) 相同,只是把总体平均换成样本平均:
$$ \|f-g\|_{L^2(\mu_n)} = \left( \frac1n\sum_{i=1}^n(f-g)(X_i)^2 \right)^{1/2}, \tag{8.32} $$
其中 $X_i$ 是 $X$ 的 i.i.d. 副本。
注意它与另一个降维结果 Johnson-Lindenstrauss 引理(Theorem 5.3.1)有相似性。证明路线也相同:集中加并集界。
固定一对不同的函数 $f,g\in\mathcal F$,并考虑
$$ \|f-g\|_{L^2(\mu_n)}^2 - \|f-g\|_{L^2(\mu)}^2 = \frac1n\sum_{i=1}^n h(X_i)-\mathbb Eh(X), $$其中 $h=(f-g)^2$。右侧是独立有界随机变量之和,因此也是次高斯;Hoeffding 不等式(Theorem 2.7.3)给出
$$ \mathbb P\left\{ \left| \|f-g\|_{L^2(\mu_n)}^2 - \|f-g\|_{L^2(\mu)}^2 \right| \gt \frac{\varepsilon^2}{4} \right\} \le 2\exp(-cn\varepsilon^4). $$(检查!)因此,以至少 $1-2\exp(-cn\varepsilon^4)$ 的概率,
$$ \|f-g\|_{L^2(\mu_n)}^2 \ge \|f-g\|_{L^2(\mu)}^2-\frac{\varepsilon^2}{4} \gt \varepsilon^2-\frac{\varepsilon^2}{4} \gt \frac{\varepsilon^2}{4}, \tag{8.33} $$其中最后用了引理的假设。现在对所有不同函数对 $f,g\in\mathcal F$ 做并集界。这样的函数对至多有 $|\mathcal F|^2$ 个,因此以至少
$$ 1-|\mathcal F|^2\cdot2\exp(-cn\varepsilon^4) \tag{8.34} $$的概率,界 (8.33) 对所有不同 $f,g\in\mathcal F$ 同时成立。由于假设 $n\ge C\varepsilon^{-4}\log|\mathcal F|$,选择足够大的绝对常数 $C$ 后,可以让 (8.34) 中的数量至少为 $0.99$。
由 packing-covering equivalence(Lemma 4.2.8),可以找到
$$ N=\mathcal N(\mathcal F,L^2(\mu),\varepsilon) $$个 $\mathcal F$ 中的函数,它们在 $L^2(\mu)$ 度量下两两 $\varepsilon$-分离。令
$$ n=\left\lceil C\varepsilon^{-4}\log N\right\rceil $$并把 Lemma 8.3.14 应用于这些函数。以正概率,这些函数在 (8.32) 定义的 $L^2(\mu_n)$ 度量下仍然 $(\varepsilon/2)$-分离,因此它们限制到
$$ \Omega_n=\{X_1,\ldots,X_n\} $$上后仍然全部不同。
固定一个使上述事件成立的 $X_1,\ldots,X_n$ 实现。因此存在子集 $\Omega_n\subset\Omega$,满足
$$ |\Omega_n|\le n\le2C\varepsilon^{-4}\log N, $$并且把所有函数限制到 $\Omega_n$ 后得到的类 $\mathcal F_n=\mathcal F|_{\Omega_n}$ 满足 $|\mathcal F_n|\ge N$。现在对 $\mathcal F_n$ 和 $\Omega_n$ 应用 Sauer-Shelah Lemma 8.3.9,得到
$$ N \le \left(\frac{en}{d_n}\right)^{d_n} \le \left( \frac{2C\varepsilon^{-4}\log N}{d_n} \right)^{d_n}, $$其中 $d_n=\operatorname{vc}(\mathcal F_n)$。化简可得
$$ N\le(2C\varepsilon^{-4})^{2d_n}. $$最后,把 $d_n=\operatorname{vc}(\mathcal F_n)$ 替换为更大的 $d=\operatorname{vc}(\mathcal F)$,并把绝对常数吸收到 $C$ 中,就得到 Theorem 8.3.13。
8.3.6 VC 大数定律
回到第 8.2.2 节,我们对所有 Lipschitz 函数上的一般经验过程做了界。现在把 Lipschitz 替换为任意具有有限 VC 维数的布尔函数类。
令 $\mathcal F$ 是某定义域 $\Omega$ 上具有有限 VC 维数的布尔函数类。令 $X,X_1,X_2,\ldots,X_n$ 是 $\Omega$ 中具有相同分布的独立随机点。那么
$$ \mathbb E\sup_{f\in\mathcal F} \left| \frac1n\sum_{i=1}^nf(X_i)-\mathbb Ef(X) \right| \le C\sqrt{\frac{\operatorname{vc}(\mathcal F)}{n}}. $$ 查看学习笔记:Theorem 8.3.15 完整证明我们把 Dudley 不等式与 Theorem 8.3.13 中的覆盖数界结合起来。首先,用对称化的经验版本(Exercise 8.11)对过程做对称化:
$$ \mathbb E\sup_{f\in\mathcal F} \left| \frac1n\sum_{i=1}^nf(X_i)-\mathbb Ef(X) \right| \le \frac2{\sqrt n} \mathbb E\sup_{f\in\mathcal F} \underbrace{ \left| \frac1{\sqrt n}\sum_{i=1}^n\varepsilon_i f(X_i) \right| }_{Z_f}. $$对 $(X_i)$ 条件化,把所有 randomness 留在 random signs $(\varepsilon_i)$ 中。为了对过程 $(Z_f)_{f\in\mathcal F}$ 使用 Dudley 不等式,需要检查增量是次高斯。三角不等式给出
$$ |Z_f-Z_g| \le \frac1{\sqrt n} \left| \sum_{i=1}^n\varepsilon_i(f-g)(X_i) \right|. $$于是,利用 Proposition 2.7.1 以及显然的事实 $\|\varepsilon_i\|_{\psi_2}\lesssim1$,得到
$$ \begin{aligned} \|Z_f-Z_g\|_{\psi_2} &\le \frac1{\sqrt n} \left\| \sum_{i=1}^n\varepsilon_i(f-g)(X_i) \right\|_{\psi_2}\\ &\lesssim \left( \frac1n\sum_{i=1}^n(f-g)(X_i)^2 \right)^{1/2} = \|f-g\|_{L^2(\mu_n)}, \end{aligned} $$其中 $\mu_n$ 是经验测度;回忆 (8.32)。
现在在给定 $(X_i)$ 的条件下使用 Dudley 不等式(Theorem 8.1.3),再对 $(X_i)$ 取期望解除条件化。得到
$$ \frac2{\sqrt n} \mathbb E\sup_{f\in\mathcal F}Z_f \lesssim \frac1{\sqrt n} \mathbb E \int_0^1 \sqrt{ \log\mathcal N(\mathcal F,L^2(\mu_n),\varepsilon) }\,d\varepsilon. \tag{8.35} $$最后,用 Theorem 8.3.13 控制覆盖数:
$$ \log\mathcal N(\mathcal F,L^2(\mu_n),\varepsilon) \lesssim \operatorname{vc}(\mathcal F)\log(2/\varepsilon). $$把它代入 (8.35),得到 $\sqrt{\log(2/\varepsilon)}$ 的 integral;该 integral 被绝对常数控制。因此
$$ \frac2{\sqrt n} \mathbb E\sup_{f\in\mathcal F}Z_f \lesssim \sqrt{\frac{\operatorname{vc}(\mathcal F)}{n}}, $$正是所需结论。
若 $\mathcal F$ 是有限 VC 维数的布尔类,那么证明 Theorem 8.3.15 时出现的量
$$ \mathbb E_\varepsilon \sup_{f\in\mathcal F} \left| \frac1n\sum_{i=1}^n\varepsilon_i f(x_i) \right| $$称为 $\mathcal F$ 在给定点集 $x_1,x_2,\ldots,x_n\in\Omega$ 上的 Rademacher 复杂度。Rademacher 复杂度反映 $\mathcal F$ 有多丰富。在 Theorem 8.3.15 的证明中,关键步骤就是把它与另一种丰富度度量,即 VC 维数,联系起来:我们证明了对任意 $n$ 点集合,$\mathcal F$ 的 Rademacher 复杂度都被 $C\sqrt{\operatorname{vc}(\mathcal F)/n}$ 控制。
下面把 Theorem 8.3.15 应用到一个经典统计问题:从样本估计随机变量 $X$ 的分布。为了从独立同分布样本 $X_1,\ldots,X_n$ 估计 $X$ 的 CDF
$$ F(x)=\mathbb P\{X\le x\}, $$
一个自然估计量是 empirical CDF,也就是满足 $X_i\le x$ 的样本点比例:
$$ F_n(x):=\frac1n|\{i:\ X_i\le x\}|. $$
令人惊讶的是,$F_n$ 在所有 $x\in\mathbb R$ 上 uniformly 逼近 $F$。
令 $X_1,\ldots,X_n$ 是独立随机变量,并具有共同 cumulative distribution 函数 $F$。那么
$$ \mathbb E\|F_n-F\|_{L^\infty} = \mathbb E\sup_{x\in\mathbb R}|F_n(x)-F(x)| \le \frac C{\sqrt n}. $$ 查看学习笔记:Theorem 8.3.17 完整证明这只是 Theorem 8.3.15 在 $\Omega=\mathbb R$ 以及半无限区间指示函数类
$$ \mathcal F := \{\mathbf 1_{(-\infty,x]}:\ x\in\mathbb R\} $$上的重述。正如 Example 8.3.2 中注意到的,这个函数类的 VC 维数至多为 $2$。
从单位正方形 $[0,1]^2$ 的均匀分布中取 $n$ 个独立同分布样本点,如 Figure 8.8 所示。对该正方形中所有圆的指示函数类 $\mathcal F$ 应用 Theorem 8.3.15。该类的 VC 维数至多为 $3$(见 Exercise 8.13)。于是,以高概率,样本满足
$$ \text{落在 }\mathcal C\text{ 中的点的比例} = \operatorname{Area}(\mathcal C)+O(1/\sqrt n) $$并且这个关系对正方形中所有圆 $\mathcal C$ 同时成立。这是几何偏差理论中的经典结果;它也适用于半平面、矩形、顶点数少的多边形等任何有限 VC 维数的类。
定义在集合 $\Omega$ 上的实值函数类 $\mathcal F$ 称为一致 Glivenko-Cantelli 类,如果对任意 $\varepsilon\gt 0$,都有
$$ \lim_{n\to\infty} \sup_\mu \mathbb P\left\{ \sup_{f\in\mathcal F} \left| \frac1n\sum_{i=1}^nf(X_i)-\mathbb Ef(X) \right| \gt \varepsilon \right\} =0. $$其中上确界对 $\Omega$ 上所有概率测度 $\mu$ 取,并且 $X,X_1,\ldots,X_n$ 是 $\Omega$ 中分布为 $\mu$ 的独立同分布点。Theorem 8.3.15 加 Markov 不等式说明:任意有限 VC 维数的布尔类都是一致 Glivenko-Cantelli 类。反过来也成立(见 Exercise 8.28),因此二者等价。
为了练习,可以做 Exercise 8.24,它给出 Theorem 8.3.15 的一个较弱但更简单版本。也建议继续看几个重要应用:高维分布的一维边缘学习(Exercise 8.25)、一比特量化(Exercise 8.26),以及没有矩假设的随机矩阵(Exercise 8.27;听起来几乎好得不像真的)。
8.4 应用:统计学习理论
统计学习(或机器学习)研究如何从数据中做预测。假设在某个集合 $\Omega$ 上有一个未知函数 $T:\Omega\to\mathbb R$,称为目标函数;我们能在若干样本点 $X_1,\ldots,X_n$ 处观察到它的函数值,其中这些样本点独立地从 $\Omega$ 上某个分布中抽取。因此训练数据是
$$ (X_i,T(X_i)), \qquad i=1,\ldots,n. \tag{8.36} $$
目标是利用这个样本来预测同一分布下新随机点 $X$ 的 $T(X)$;见 Figure 8.9。
一类重要的学习问题是分类,此时函数 $T$ 是布尔的,只取 $0$ 和 $1$,用于把 $\Omega$ 中的点分到两个类别。比如在一项健康研究中,有 $n$ 个病人。对每个病人,我们记录 $d$ 个健康参数,如血压或体温;这就是向量 $X_i\in\mathbb R^d$。假设还知道他们是否患有糖尿病:$T(X_i)=0$ 表示健康,$T(X_i)=1$ 表示患病。
目标是从数据 (8.36) 中学习如何根据健康参数预测糖尿病;也就是学习函数 $T:\mathbb R^d\to\{0,1\}$,以便根据新病人的健康参数做诊断。
8.4.1 风险、拟合与复杂度
给定训练数据 (8.36),我们希望找到函数 $f:\Omega\to\mathbb R$ 来近似 $T$,以便之后通过查看 $f(X)$ 为新病人诊断。目标是最小化误诊新病人的风险,定义为
$$ R(f)=\mathbb E\left(f(X)-T(X)\right)^2. \tag{8.37} $$
在分类问题中,若 $T$ 与 $f$ 都是布尔函数,则风险就是误分类概率(检查!):
$$ R(f)=\mathbb P\{f(X)\ne T(X)\}. \tag{8.38} $$需要多少训练数据?这取决于问题的复杂度。如果认为目标函数 $T(X)$ 的行为很复杂,就需要更多数据。由于通常事先不知道它有多复杂,我们会把候选函数 $f$ 限制在某个函数类 $\mathcal F$ 中,称为假设类。
但如何选择 $\mathcal F$?没有通用规则,不过它需要在拟合和复杂度之间平衡。如果 $\mathcal F$ 太简单,例如只含线性函数,可能会欠拟合(见 Figure 8.10a),错过真实模式。若 $\mathcal F$ 太复杂,可能过拟合,只是记住训练数据,而不是从数据中泛化(Figure 8.10b),并且需要更多数据。理想情况是假设类刚好足够复杂,能捕捉真实模式,又不追逐随机噪声(Figure 8.10c)。
8.4.2 经验风险最小化
一旦选定假设空间 $\mathcal F$,可能会想直接选择其中最好的函数 $f^*$,即最小化风险 (8.37) 的函数:
$$ f^* := \arg\min_{f\in\mathcal F}R(f). $$
问题是我们无法实际计算风险 $R(f)$,因为这需要对整个总体 $\Omega$ 取期望。解决办法是:改为对训练数据取期望。
对函数 $f:\Omega\to\mathbb R$,定义经验风险与经验最小化器为
$$ R_n(f):= \frac1n\sum_{i=1}^n(f(X_i)-T(X_i))^2, \qquad f_n^*:= \arg\min_{f\in\mathcal F}R_n(f). \tag{8.39} $$现在,$R_n(f)$ 和 $f_n^*$ 都能从训练数据中计算出来。学习的输出是 $f_n^*$,它的质量由泛化误差 $R(f_n^*)$ 衡量。
在分类中,$f$ 与 $T$ 都只取 $0$ 或 $1$,所以经验风险 $R_n(f)$ 就是训练集中被 $f$ 分类错误的比例。ERM 会选择在训练数据上犯错最少的 $f\in\mathcal F$。
8.4.3 VC 泛化界
若目标函数 $T$ 是布尔函数,且假设类 $\mathcal F$ 具有有限 VC 维数,则
$$\mathbb E R(f_n^*)\le R(f^*)+C\sqrt{\frac{\operatorname{vc}(\mathcal F)}{n}}.$$ 查看学习笔记:Theorem 8.4.5 完整证明步骤 1:超额风险。 下列逐点界成立:
$$ R(f_n^*)-R(f^*) \le 2\sup_{f\in\mathcal F}|R_n(f)-R(f)|. \tag{8.40} $$为了检查它,记
$$ \varepsilon:= \sup_{f\in\mathcal F}|R_n(f)-R(f)| $$并写出链式不等式:
$$ \begin{aligned} R(f_n^*) &\le R_n(f_n^*)+\varepsilon &&\text{since }f_n^*\in\mathcal F\text{ by construction}\\ &\le R_n(f^*)+\varepsilon &&\text{since }f_n^*\text{ minimizes }R_n\text{ in }\mathcal F\\ &\le R(f^*)+2\varepsilon &&\text{since }f^*\in\mathcal F\text{ by construction}. \end{aligned} $$两边减去 $R(f^*)$,得到 (8.40)。
步骤 2:应用 VC 大数定律。 由 (8.40),只需证明
$$ \mathbb E\sup_{f\in\mathcal F}|R_n(f)-R(f)| \lesssim \sqrt{\frac{\operatorname{vc}(\mathcal F)}{n}}. $$回忆经验风险与真实总体风险的定义 (8.39)、(8.37),可把它改写成
$$ \mathbb E\sup_{\ell\in\mathcal L} \left| \frac1n\sum_{i=1}^n\ell(X_i)-\mathbb E\ell(X) \right| \lesssim \sqrt{\frac{\operatorname{vc}(\mathcal F)}{n}}, \tag{8.41} $$其中
$$ \mathcal L=\{(f-T)^2:\ f\in\mathcal F\}. $$稍作思考可知 $\operatorname{vc}(\mathcal L)=\operatorname{vc}(\mathcal F)$(Exercise 8.29)。因此应用 Theorem 8.3.15 即完成证明。
假设从单位正方形中均匀抽取 $n$ 个训练点,并按照某个固定圆 $C$ 给标签:在圆内标为 $1$,在圆外标为 $0$。目标是从数据中学出这个“疾病圆”。
我们使用经验风险最小化:选一个最匹配标签的圆,也就是使误分类数量最少的圆(圆内健康点和圆外患病点的总数最少)。
效果如何?真实圆 $C$ 给出零误差,即 $R(f^*)=0$;而圆类的 VC 维数至多为 $3$(Exercise 8.13)。因此 Theorem 8.4.5 告诉我们,学到的圆的风险至多为 $O(1/\sqrt n)$。所以新点只需检查是否在学到的圆内即可分类,误诊概率为 $O(1/\sqrt n)$,并且随着数据增多而下降。
VC 泛化界揭示了学习误差的两个来源:
$$ \mathbb E R(f_n^*)\le R(f^*)+C\sqrt{\frac{\operatorname{vc}(\mathcal F)}{n}}. $$偏差项 $R(f^*)$ 来自假设类选择不完美,也就是欠拟合。可以通过在 $\mathcal F$ 中加入更多函数来降低偏差;理想情况下,$\mathcal F$ 足以捕捉真实目标函数 $T$,使偏差等于零。但这样方差项 $O(\sqrt{\operatorname{vc}(\mathcal F)/n})$ 会增大。为了控制它,需要使用更多训练数据,也就是增大 $n$,以避免过拟合。
为了练习,可以把理论扩展到更现实的设置:标签是随机的但与输入相关(Exercise 8.30);离开布尔类,说明如何学习 Lipschitz 函数(Exercise 8.31);并展示如果 VC 维数无限,则学习会失败(Exercise 8.32)。
8.5 泛型链式方法
虽然 Dudley 不等式是一个简单且通用的工具,但它可能比较偏松(见 Exercise 8.4)。覆盖数 $\mathcal N(T,d,\varepsilon)$ 本身并不包含足够信息,无法完全捕捉随机过程 $(X_t)_{t\in T}$ 的大小。
幸运的是,有一种方法可以根据 $T$ 的几何给出 $\mathbb E\sup_{t\in T}X_t$ 的准确的双侧界。这个方法称为泛型链式方法,本质上是对 Dudley 不等式(Theorem 8.1.4)证明中链式方法的强化。
8.5.1 改造 Dudley 不等式
回忆我们通过链式方法得到的界 (8.12):
$$ \mathbb E\sup_{t\in T}X_t \lesssim \sum_{k=\kappa+1}^{\infty} \varepsilon_{k-1}\sqrt{\log |T_k|}, \tag{8.42} $$
其中 $\varepsilon_k=2^{-k}$,$T_k$ 是 $T$ 的 smallest $\varepsilon_k$-nets,所以 $|T_k|=\mathcal N(T,d,\varepsilon_k)$,并且 $\kappa$ 选择为使 $|T_\kappa|=1$。
现在换一种角度:不是固定 $\varepsilon_k$ 并最小化 $|T_k|$,而是固定 $|T_k|$ 并最小化 $\varepsilon_k$。具体地,选取子集 $T_k\subset T$,使得
$$ |T_0|=1, \qquad |T_k|\le2^{2^k}, \qquad k=1,2,\ldots, \tag{8.43} $$
并定义
$$ \varepsilon_k=\sup_{t\in T}d(t,T_k), $$
其中 $d(t,T_k)$ 表示点 $t$ 到集合 $T_k$ 的距离。这样每个 $T_k$ 就是一个 $\varepsilon_k$-net,链式方法界 (8.42) 变成
$$ \mathbb E\sup_{t\in T}X_t \lesssim \sum_{k=1}^{\infty} 2^{k/2}\sup_{t\in T}d(t,T_{k-1}), $$
重新编号后为
$$ \mathbb E\sup_{t\in T}X_t \lesssim \sum_{k=0}^{\infty} 2^{k/2}\sup_{t\in T}d(t,T_k). \tag{8.44} $$
8.5.2 $\gamma_2$ 泛函与泛型链式方法
到目前为止,我们只是把 Dudley 不等式换了一种形式,还没有本质改进。关键步骤现在出现:泛型链式方法允许我们把 (8.44) 中的上确界拉到 sum 的外面。得到的量有一个名字。
设 $(T,d)$ 是度量空间。一列子集 $(T_k)_{k=0}^{\infty}$ 若满足 (8.43),称为可容许序列。$T$ 的 $\gamma_2$ 泛函定义为
$$ \gamma_2(T,d) = \inf_{(T_k)} \sup_{t\in T} \sum_{k=0}^{\infty}2^{k/2}d(t,T_k), \tag{8.45} $$其中 infimum 取遍所有 admissible sequences。
因为 $\gamma_2$ 泛函中的上确界位于 sum 外面,所以它小于 Dudley 求和 (8.44)。这个变化看似很小,但在一些情形下会造成真实差别(见 Exercise 8.34)。
好消息是,可以把 Dudley 不等式(Theorem 8.1.4)中的 Dudley 求和(或积分)替换成更紧的 $\gamma_2$ 泛函。
若 $(X_t)_{t\in T}$ 是均值为零、具有次高斯增量的随机过程,则
$$\mathbb E\sup_{t\in T}X_t\le CK\gamma_2(T,d).$$ 查看学习笔记:Theorem 8.5.2 完整证明我们沿用 Dudley 不等式(Theorem 8.1.4)证明中的链式方法,但这次更细致。
步骤 1:链式方法集合。 和之前一样,可以假设 $K=1$ 且 $T$ 有限,此时 $\gamma_2(T,d)$ 有限(为什么?)。令 $(T_k)$ 是 $T$ 的可容许序列,并且几乎达到 (8.45) 中的上确界:
$$ \sup_{t\in T} \sum_{k=0}^{\infty} 2^{k/2}d(t,T_k) \le 2\gamma_2(T,d) \lt \infty. \tag{8.46} $$记 $T_0=\{t_0\}$。必然存在某个 $K$ 使得 $T_K=T$;否则某个 $t\in T$ 会被无穷多个集合 $T_k$ 遗漏,于是对这些 $k$ 有 $d(t,T_k)\gt \varepsilon$,其中 $\varepsilon\gt 0$ 固定,这会使 (8.46) 中的级数发散。
从 $t_0$ 到一般点 $t\in T$,沿着有限链走:
$$ t_0=\pi_0(t)\to\pi_1(t)\to\pi_2(t)\to\cdots\to t, $$其中 $\pi_k(t)\in T_k$ 是 $T_k$ 中对 $t$ 的最佳近似点,即
$$ d(t,\pi_k(t))=d(t,T_k). $$位移 $X_t-X_{t_0}$ 可写成类似 (8.9) 的望远镜求和:
$$ X_t-X_{t_0} = \sum_{k=1}^{\infty} \bigl(X_{\pi_k(t)}-X_{\pi_{k-1}(t)}\bigr). \tag{8.47} $$步骤 2:控制增量。 这里需要更谨慎。我们希望证明:以高概率,下列事件成立:
$$ |X_{\pi_k(t)}-X_{\pi_{k-1}(t)}| \lesssim 2^{k/2}d(t,T_k) \qquad \forall k\in\mathbb N,\ \forall t\in T. \tag{8.48} $$若能对所有 $k$ 求和,就会得到关于 $\gamma_2(T,d)$ 的所需界。
为了证明 (8.48),先固定 $k$ 和 $t$。次高斯假设给出
$$ \|X_{\pi_k(t)}-X_{\pi_{k-1}(t)}\|_{\psi_2} \le d(\pi_k(t),\pi_{k-1}(t)). $$因此,对每个 $u\ge0$,事件
$$ |X_{\pi_k(t)}-X_{\pi_{k-1}(t)}| \le Cu\,2^{k/2}d(\pi_k(t),\pi_{k-1}(t)) \tag{8.49} $$以至少 $1-2\exp(-8u^2 2^k)$ 的概率成立。(为了得到常数 $8$,把绝对常数 $C$ 选得足够大。)
现在通过并集界去掉对 $t\in T$ 的固定。可出现的点对 $(\pi_k(t),\pi_{k-1}(t))$ 至多有
$$ |T_k|\cdot |T_{k-1}| \le |T_k|^2 \le 2^{2^{k+1}} $$个。再对所有 $k\in\mathbb N$ 做并集界。于是 (8.49) 对所有 $t\in T$ 与 $k\in\mathbb N$ 同时成立的概率至少为
$$ 1- \sum_{k=1}^{\infty} 2^{2^{k+1}}\cdot2\exp(-8u^2 2^k) \ge 1-2\exp(-u^2), $$只要 $u\gt c$。(检查最后一个不等式!)
步骤 3:累加增量。 在 (8.49) 对所有 $t$ 和 $k$ 都成立的事件上,对 $k\in\mathbb N$ 求和,并代入链式方法求和 (8.47),得到
$$ |X_t-X_{t_0}| \lesssim u\sum_{k=1}^{\infty} 2^{k/2}d(\pi_k(t),\pi_{k-1}(t)), \tag{8.50} $$其中 $\lesssim$ 隐藏绝对常数。由三角不等式,
$$ d(\pi_k(t),\pi_{k-1}(t)) \le d(t,\pi_k(t))+d(t,\pi_{k-1}(t)). $$用这个界、重新编号,并代入 (8.46),可知 (8.50) 右侧至多为 $Cu\gamma_2(T,d)$;也就是说
$$ |X_t-X_{t_0}| \lesssim u\gamma_2(T,d). $$(检查!)对 $T$ 取上确界,得到
$$ \sup_{t\in T}|X_t-X_{t_0}| \lesssim u\gamma_2(T,d). $$因为它对任意 $u\gt c$ 都以至少 $1-2\exp(-u^2)$ 的概率成立,所以
$$ \left\| \sup_{t\in T}|X_t-X_{t_0}| \right\|_{\psi_2} \lesssim \gamma_2(T,d). $$这很快推出 Theorem 8.5.2 的结论(检查!)。
和 Dudley 不等式(Remark 8.1.5)类似,泛型链式方法实际上给出
$$ \mathbb E\sup_{t,s\in T}|X_t-X_s| \le CK\gamma_2(T,d). $$这个结论即使没有均值零假设 $\mathbb EX_t=0$ 也成立。
Theorem 8.5.2 只给出期望界,但泛型链式方法实际上也给出高概率界;这个逻辑我们在 Dudley 不等式(Remark 8.1.6)中已经见过。
假设 $T$ 有限。对每个 $u\ge0$,事件
$$ \sup_{t,s\in T}|X_t-X_s| \le CK\bigl(\gamma_2(T,d)+u\,\operatorname{diam}(T)\bigr) $$以至少 $1-2\exp(-u^2)$ 的概率成立。Exercise 8.35 会要求你证明它。对高斯过程,也可以直接从高斯集中推出该结果。
为了练习,可以把泛型链式方法应用于经验过程(Exercise 8.36)。
8.5.3 主控测度与比较定理
$\gamma_2$ 泛函(Definition 8.5.1)通常比 Dudley 不等式(Theorem 8.1.3)中的覆盖数更难计算。但这种努力通常是值得的:泛型链式方法不同于 Dudley,它在常数因子意义下是尖锐的。
设 $(X_t)_{t\in T}$ 是集合 $T$ 上的均值为零高斯过程,并用典范度量
$$ d(t,s)=\|X_t-X_s\|_{L^2} $$装备 $T$,如 (7.1)。那么
$$ c\gamma_2(T,d) \le \mathbb E\sup_{t\in T}X_t \le C\gamma_2(T,d). $$Theorem 8.5.5 的上界直接来自泛型链式方法(Theorem 8.5.2)。下界更难;它的证明可看作 Sudakov 不等式(Theorem 7.4.1)的精细多尺度推广,本书不纳入。
Theorem 8.5.2 告诉我们,上界不仅对高斯过程成立,也对所有次高斯过程成立。因此,把上界和下界结合起来,就可以用高斯过程控制任意次高斯过程。
令 $(X_t)_{t\in T}$ 是集合 $T$ 上的均值为零随机过程,令 $(Y_t)_{t\in T}$ 是均值为零高斯过程。假设
$$\|X_t-X_s\|_{\psi_2}\le K\|Y_t-Y_s\|_{L^2},$$对所有 $t,s\in T$ 成立。那么
$$\mathbb E\sup_{t\in T}X_t\le CK\,\mathbb E\sup_{t\in T}Y_t.$$ 查看学习笔记:Corollary 8.5.6 完整证明考虑 $T$ 上的典范度量
$$ d(t,s)=\|Y_t-Y_s\|_{L^2}. $$现在先用泛型链式方法界(Theorem 8.5.2),再用主控测度定理(Theorem 8.5.5)的下界,得到
$$ \mathbb E\sup_{t\in T}X_t \lesssim K\gamma_2(T,d) \lesssim K\mathbb E\sup_{t\in T}Y_t. $$Corollary 8.5.6 可以看作 Sudakov-Fernique 不等式(Theorem 7.2.8)的次高斯过程版本。代价是右侧多出一个绝对常数因子,而不再局限于高斯过程之间的比较。
下面把 Corollary 8.5.6 应用于集合 $T\subset\mathbb R^n$ 上的典范高斯过程 $Y_x=\langle g,x\rangle$,其中 $g\sim N(0,I_n)$。回忆第 7.5 节中
$$ w(T)=\mathbb E\sup_{x\in T}\langle g,x\rangle $$
是集合 $T$ 的高斯宽度。于是立刻得到下面的几何形式。
令 $(X_x)_{x\in T}$ 是子集 $T\subset\mathbb R^n$ 上的均值为零随机过程。假设
$$\|X_x-X_y\|_{\psi_2}\le K\|x-y\|_2,$$对所有 $x,y\in T$ 成立。那么
$$\mathbb E\sup_{x\in T}X_x\le CKw(T).$$ 查看学习笔记:Corollary 8.5.8 完整证明由典范高斯过程的定义,
$$ \|Y_x-Y_y\|_{L^2} = \|\langle g,x-y\rangle\|_{L^2} = \|x-y\|_2. $$因此 Corollary 8.5.6 的假设正好成立。结论给出
$$ \mathbb E\sup_{x\in T}X_x \le CK\mathbb E\sup_{x\in T}\langle g,x\rangle = CKw(T). $$一个直接后果是:若 $X$ 是 $\mathbb R^n$ 中的次高斯随机向量,则
$$ \mathbb E\sup_{x\in T}\langle X,x\rangle\lesssim \|X\|_{\psi_2}w(T). $$这是把 Corollary 8.5.8 应用于过程 $X_x=\langle X,x\rangle$ 得到的,因为
$$ \|\langle X,x\rangle-\langle X,y\rangle\|_{\psi_2} = \|\langle X,x-y\rangle\|_{\psi_2} \le \|X\|_{\psi_2}\|x-y\|_2. $$Exercise 8.38 会进一步把它用于次高斯 vectors 的 $\ell^p$ 范数。
8.6 Chevet 不等式
Talagrand 比较不等式,更一般地说泛型链式方法,是一个可以用于很多场景的强大工具。现在用它得到一个非常一般的随机双线性型一致界:
$$ \sup_{x\in T,y\in S}\langle Ax,y\rangle\le ? $$
这里 $A$ 是随机矩阵,$T,S$ 是任意有界集。
当 $T$ 和 $S$ 是欧氏球时,一个特殊情形会给出 $A$ 的算子范数;这正是我们已经在 Theorem 4.4.3 和 Theorem 7.3.1 中研究过的对象。现在我们处理一般有界集,目标是只用两个几何量来刻画上界:高斯宽度和半径,其中
$$ \operatorname{rad}(T):=\sup_{x\in T}\|x\|_2. \tag{8.51} $$
令 $A$ 为 $m\times n$ 随机矩阵,其行 $A_i$ 独立、均值为零且次高斯。若 $T\subset\mathbb R^n$、$S\subset\mathbb R^m$ 有界,则
$$\mathbb E\sup_{x\in T,y\in S}\langle Ax,y\rangle \le CK\{w(T)\operatorname{rad}(S)+w(S)\operatorname{rad}(T)\}.$$其中 $K=\max_i\|A_i\|_{\psi_2}$。若把“行”换成“列”,同样结论也成立。
查看学习笔记:Theorem 8.6.1 完整证明证明沿用 Theorem 7.3.1 的思路,只是把 Sudakov-Fernique 不等式换成 Talagrand 比较不等式。
不失一般性,设 $K=1$。我们需要控制如下随机过程:
$$ X_{uv}:=\langle Au,v\rangle,\qquad u\in T,\ v\in S. $$为了验证它的增量是次高斯,固定 $(u,v),(w,z)\in T\times S$,并写成
$$ X_{uv}-X_{wz} =X_{uv}-X_{wv}+X_{wv}-X_{wz} =\langle A(u-w),v\rangle+\langle Aw,v-z\rangle. $$由三角不等式和次高斯假设(见 Exercise 3.34;如果当时跳过了,现在应当补做),得到
$$ \begin{aligned} \|X_{uv}-X_{wz}\|_{\psi_2} &\le \|\langle A(u-w),v\rangle\|_{\psi_2} +\|\langle Aw,v-z\rangle\|_{\psi_2} \\ &\lesssim \|u-w\|_2\|v\|_2+\|v-z\|_2\|w\|_2 \\ &\le \|u-w\|_2\operatorname{rad}(S) +\|v-z\|_2\operatorname{rad}(T). \end{aligned} \tag{8.52} $$接下来为 Talagrand 比较不等式(Corollary 8.5.6)选一个更简单的高斯过程。增量界 (8.52) 指向如下自然选择:
$$ Y_{uv}:=\langle g,u\rangle\operatorname{rad}(S) +\langle h,v\rangle\operatorname{rad}(T), $$其中 $g\sim N(0,I_n)$,$h\sim N(0,I_m)$,且二者独立。这个过程的增量满足
$$ \|Y_{uv}-Y_{wz}\|_{L^2}^2 = \|u-w\|_2^2\operatorname{rad}(S)^2 +\|v-z\|_2^2\operatorname{rad}(T)^2. $$这一步可以像 Theorem 7.3.1 的证明那样直接验证。把它与 (8.52) 比较,并使用 $a+b\le\sqrt{2(a^2+b^2)}$,可得
$$ \|X_{uv}-X_{wz}\|_{\psi_2} \lesssim \|Y_{uv}-Y_{wz}\|_{L^2}. $$于是由 Talagrand 比较不等式(Corollary 8.5.6),证明完成:
$$ \begin{aligned} \mathbb E\sup_{u\in T,v\in S}X_{uv} &\lesssim \mathbb E\sup_{u\in T,v\in S}Y_{uv} \\ &= \mathbb E\sup_{u\in T}\langle g,u\rangle\operatorname{rad}(S) + \mathbb E\sup_{v\in S}\langle h,v\rangle\operatorname{rad}(T) \\ &= w(T)\operatorname{rad}(S)+w(S)\operatorname{rad}(T). \end{aligned} $$在特殊情形 $T=S^{n-1}$ 且 $S=S^{m-1}$ 中,Chevet 不等式给出我们熟悉的算子范数尖锐界:
$$ \mathbb E\|A\|\le CK(\sqrt n+\sqrt m). $$这个结论我们之前在 Section 4.4.2 中用 $\varepsilon$-nets 证明过。但新的方法更加灵活。例如,若把 $T$ 和 $S$ 取成 $\ell^p$ 球,就会得到随机矩阵的 $\|A\|_{p\to q}$ 范数;可以试做 Exercise 8.41。
若 $A$ 是元素独立且服从 $N(0,1)$ 的高斯矩阵,我们甚至可以用尖锐常数 $1$ 证明 Chevet 不等式:
$$ \mathbb E\sup_{x\in T,y\in S}\langle Ax,y\rangle \le w(T)\operatorname{rad}(S)+w(S)\operatorname{rad}(T), $$并且还有一个只差常数因子的反向不等式(Exercise 8.39)。后面在 Section 9.7.1 中,我们会进一步改进高斯 Chevet 不等式。
查看学习笔记:Exercise 8.39 证明作为练习,可以证明 Chevet 不等式的高概率版本(Exercise 8.40)。
8.7 Notes
链式方法的思想已经出现在 Kolmogorov 关于布朗运动连续性定理的证明中;可参见例如 [253, Chapter 1]。Dudley 积分不等式(Theorem 8.1.3)可追溯到 R. Dudley 的工作。我们在 Section 8.1 中的叙述主要遵循 [210, Chapter 11]、[315, Section 1.2] 和 [330, Section 5.3]。
Monte Carlo 方法(第 8.2.1 节)在科学计算中极其流行,尤其是在它与 Markov 链的力量结合起来时;参见例如 [63]。经验过程的丰富理论(见 Definition 8.2.5)在统计学和机器学习中有很多应用;参见 [329, 328, 279, 232]。Lipschitz 大数定律(Theorem 8.2.3)也称为 Wasserstein 大数定律,它大致基于 [330, Example 5.15];另见 [115, Chapter 11]。关于测度输运(第 8.2.3 节)的入门,可参见 [342, 85]。
Section 8.3 中引入的 VC 维数概念源于 V. Vapnik 和 A. Chervonenkis 的基础性工作 [335];现代处理可参见例如 [329, Section 2.6.1]、[210, Section 14.3]、[330, Section 7.2]、[225, Sections 10.2-10.3]、[232, Section 2.2]、[329, Section 2.6]。Pajor Lemma 8.3.7 最早归功于 A. Pajor [268];另见 [128]、[210, Proposition]、[330, Theorem 7.19]、[329, Lemma 2.6.2]。
我们现在称为 Sauer-Shelah 引理(Lemma 8.3.9)的结果,是由 V. Vapnik 与 A. Chervonenkis [335]、N. Sauer [294]、以及 M. Perles 与 S. Shelah [300] 独立证明的。文献中可以找到 Sauer-Shelah 引理的多种证明,例如 [46, Chapter 17]、[225, Sections 10.2-10.3]、[210, Section 14.3]。Sauer-Shelah 引理也有许多变体,参见例如 [158, 312, 313, 15, 337]。增长函数(第 8.3.4 节)最早由 Vapnik-Chervonenkis 引入,他们研究了它的多种性质,包括指数-多项式二分现象(Exercise 8.20)。
Theorem 8.3.13 归功于 R. Dudley [114];参见 [210, Section 14.3]、[329, Theorem 2.6.4]。降维 Lemma 8.3.14 隐含在 Dudley 的证明中;它在 [237] 中被明确表述,并在 [330, Lemma 7.17] 中复现。关于把 VC 理论从 $\{0,1\}$ 推广到一般实值函数类,可参见 [237, 288]、[330, Sections 7.3-7.4]。
自从 V. Vapnik 与 A. Chervonenkis 的基础性工作 [335] 以来,像 Theorem 8.3.15 这样通过 VC 维数控制经验过程的界,一直是统计学习理论的核心主题;参见例如 [232, 32, 329, 288]、[330, Chapter 7]。我们对 Theorem 8.3.15 的展示基于 [330, Corollary 7.18]。虽然这个结果的显式陈述在早期文献中不容易找到,但它可以从 [32, Theorem 6]、[56, Section 5] 推出。
Glivenko-Cantelli 定理(Theorem 8.3.17)是 1933 年的结果 [138, 75],它早于并且部分推动了后来的 VC 理论发展;关于 Glivenko-Cantelli 定理和概率论中的其他一致性结果,可参见 [210, Section 14.2]、[329, 115]。Example 8.3.18 讨论了差异理论中的一个基本问题;关于差异理论的系统处理,见 [224]。
在第 8.4 节中,我们只是接触了统计学习理论的表面。它是概率论、统计学和理论计算机科学交叉处的一个庞大领域。若要更深入地入门,可参见教程 [52, 232] 和书籍 [171, 156, 196]。
第 8.5 节中介绍的泛型链式方法,是 M. Talagrand 自 1985 年以来(在 X. Fernique [124] 更早工作的基础上)提出的用于得到高斯过程尖锐界的方法。我们的叙述基于 Talagrand 的书 [315, 316, 317],这些书非常详细地讨论了泛型链式方法的影响、应用和历史。次高斯过程的上界(Theorem 8.5.2)可见 [315, Theorem 2.2.22];下界,即主控测度定理 8.5.5,可见 [315, Theorem 2.4.1]。Talagrand 比较不等式(Corollary 8.5.6)取自 [315, Theorem 2.4.12]。泛型链式方法的另一种展示可见 [330, Chapter 6]。R. van Handel 最近在 [332, 333] 中给出了主控测度定理的另一种证明。泛型链式方法界的高概率版本(Remark 8.5.4)来自 [316, Theorem 2.2.27] 和 [317, Theorem 2.7.13];S. Dirksen [103] 也用另一种方法证明过它。
第 8.6 节讨论次高斯过程的 Chevet 不等式。对于高斯过程,这一结果可追溯到 S. Chevet [84];常数随后由 Y. Gordon [140] 改进,从而得到我们在 Exercise 8.39 中陈述的结果。关于这个结果的一种讲解可见 [21, Section 9.4]。关于 Chevet 不等式的变体和应用,可见 [321, 7]。
Sudakov 和 Dudley 不等式之间的对数间隙(Exercise 8.5)是最优的;cross-polytope $T=B_1^n$ 可作为例子(见 [317, Exercise 2.5.11])。
次指数 Dudley 不等式(Exercise 8.7)可以很容易地推广到几乎任何类型的尾衰减 [210, Section 11.1]。
局部 Dudley 不等式(Exercise 8.8)是一个在统计学习理论中以及证明随机过程连续性时有用的工具;参见 [330, Section 5.4]。
Exercise 8.10 把 Lipschitz 大数定律推广到更高维。关于它的历史、扩展、改进以及近似最优性,可参见 [85, Chapter 2]。
Exercise 8.21 中得到的界 $O(dk\log k)$ 是最优的 [119]。
增长函数的多项式-指数二分现象(Exercise 8.20),以及把大数定律中的一致收敛(Exercise 8.28)和可学习性(Exercise 8.32)刻画为有限 VC 维数的结论,本质上都可以追溯到 Vapnik 和 Chervonenkis 的原始工作 [335]。
一比特量化(Exercise 8.26)在信号处理和机器学习中已被广泛研究;参见例如 [89, 53, 355, 198, 199, 275, 274, 276, 13, 341, 92, 265, 189, 278, 168, 28, 351, 105, 129, 106, 104, 226, 81, 107]。Exercise 8.26 中的结果来自 [265, Theorem 2.2]。
Exercise 8.27 给出了小球方法的简化版本。小球方法是在 [191] 中提出的,用来放宽随机矩阵理论和机器学习中对分布的假设 [234, 235, 208, 207]。
随机矩阵的 $p\to q$ 范数(Exercise 8.41)最早在 [37] 中针对 $p=2$ 且 $q\in[2,\infty)$ 的有界元素情形被研究;更一般的近期结果可见 [9, 100, 200, 202, 201]。
Exercises
修改 Section 8.1 的链式方法论证,推出 Remark 8.1.6 中的高概率 Dudley 不等式。为此,把 (8.11) 升级为如下高概率版本:
$$ \sup_{t\in T} \bigl( X_{\pi_k(t)}-X_{\pi_{k-1}(t)} \bigr) \le C\varepsilon_{k-1} \left[ \sqrt{\log|T_k|}+z_k \right] \tag{8.53} $$该事件的概率至少为 $1-2e^{-z_k^2}$。选择合适的 $z_k$,使得可以对 (8.12) 中所有项做并集界。
查看学习笔记:Exercise 8.1 证明把 Dudley 不等式与高斯集中(Theorem 7.1.11)结合起来,为高斯过程给出 Remark 8.1.6 的另一种证明。
查看学习笔记:Exercise 8.2 证明在 Theorem 8.1.3 的证明中,我们用 Dudley 积分控制 Dudley 求和。证明二者实际上等价到常数因子:
$$ \sum_{k\in\mathbb Z} 2^{-k} \sqrt{\log\mathcal N(T,d,2^{-k})} \asymp \int_0^\infty \sqrt{\log\mathcal N(T,d,\varepsilon)} \,d\varepsilon. $$ 查看学习笔记:Exercise 8.3 证明令 $e_k$ 为 $\mathbb R^n$ 中的标准基向量,并令
$$ T:= \left\{ \frac{e_k}{\sqrt{1+\log k}}: k=1,\ldots,n \right\}. $$(a) 证明 $w(T)\le C$。
(b) 证明
$$ \int_0^\infty \sqrt{\log\mathcal N(T,d,\varepsilon)} \,d\varepsilon \to\infty \quad\text{as }n\to\infty. $$ 查看学习笔记:Exercise 8.4 证明Sudakov 不等式(Corollary 7.4.2)和 Dudley 不等式(Theorem 8.1.8)对任意有界 $T\subset\mathbb R^n$ 给出
$$ s(T)\lesssim w(T)\lesssim d(T), $$其中
$$ s(T)=\sup_{\varepsilon\gt 0} \varepsilon\sqrt{\log\mathcal N(T,\varepsilon)}, \qquad d(T)=\int_0^\infty \sqrt{\log\mathcal N(T,\varepsilon)} \,d\varepsilon. $$证明这两个界至多差一个对数因子:
$$ s(T)\le d(T)\le C\log(n)\,s(T). $$因此 Sudakov 不等式与 Dudley 不等式都在对数因子意义下尖锐。
查看学习笔记:Exercise 8.5 证明Remark 8.1.7 说明 Dudley 积分的上限可取为 $\operatorname{diam}(T)$。现在改进下限。证明任意有界 $T\subset\mathbb R^n$ 满足
$$ w(T) \le C\int_a^b \sqrt{\log\mathcal N(T,\varepsilon)} \,d\varepsilon, \qquad a=\frac{cw(T)}{\sqrt n}, \quad b=\operatorname{diam}(T). $$ 查看学习笔记:Exercise 8.6 证明设均值为零随机过程 $(X_t)_{t\in T}$ 的增量不是次高斯,而是次指数:
$$ \|X_t-X_s\|_{\psi_1}\le Kd(t,s) \quad\text{for all }t,s\in T. $$证明 Dudley 不等式的如下版本:
$$ \mathbb E\sup_{t\in T}X_t \le CK\int_0^\infty \log\mathcal N(T,d,\varepsilon) \,d\varepsilon. $$ 查看学习笔记:Exercise 8.7 证明令 $(X_t)_{t\in T}$ 是度量空间 $(T,d)$ 上具有次高斯增量的随机过程。给定 $\delta\gt 0$,证明
$$ \mathbb E \sup_{\substack{s,t\in T\\ d(s,t)\le\delta}} |X_s-X_t| \le CK\int_0^\delta \sqrt{\log\mathcal N(T,d,\varepsilon)} \,d\varepsilon. $$ 查看学习笔记:Exercise 8.8 证明对函数类
$$ \mathcal F= \{f:[0,1]\to[0,1]:\|f\|_{\mathrm{Lip}}\le1\}, $$证明对任意 $\varepsilon\in(0,1)$,
$$ \mathcal N(\mathcal F,\|\cdot\|_{L^\infty},\varepsilon) \le e^{C/\varepsilon}. $$ 查看学习笔记:Exercise 8.9 证明把 Theorem 8.2.3 推广到高维。考虑单位立方体 $[0,1]^d$ 上关于 $\|\cdot\|_\infty$ 度量的 $L$-Lipschitz 函数类
$$ \mathcal F= \{f:[0,1]^d\to\mathbb R:\|f\|_{\mathrm{Lip}}\le L\}. $$令 $X,X_1,X_2,\ldots$ 是取值于 $[0,1]^d$ 的 i.i.d. 随机变量。证明
$$ \mathbb E\sup_{f\in\mathcal F} \left| \frac1n\sum_{i=1}^n f(X_i)-\mathbb Ef(X) \right| \le CL \begin{cases} \log(n)n^{-1/2},& d=2,\\ n^{-1/d},& d\gt 2. \end{cases} $$证明时可先假设 $L=1$,并按两步做:(a) 推广覆盖数界为 $\mathcal N(\mathcal F,\|\cdot\|_\infty,\varepsilon)\le e^{C/\varepsilon^d}$;(b) 把 Dudley 界截断为
$$ \mathbb E\sup_{f\in\mathcal F}|X_f| \lesssim \delta+\frac1{\sqrt n} \int_\delta^1 \sqrt{\log\mathcal N(\mathcal F,\|\cdot\|_\infty,\varepsilon)} \,d\varepsilon. \tag{8.54} $$ 查看学习笔记:Exercise 8.10 证明修改 Symmetrization Lemma 6.3.2 的证明,得到经验过程的对称化。令 $\mathcal F$ 是定义域 $\Omega$ 上的布尔函数类,令 $X,X_1,\ldots,X_n$ 是同分布独立样本。证明
$$ \mathbb E\sup_{f\in\mathcal F} \left| \frac1n\sum_{i=1}^n f(X_i)-\mathbb Ef(X) \right| \le 2\mathbb E\sup_{f\in\mathcal F} \left| \frac1n\sum_{i=1}^n \varepsilon_i f(X_i) \right|, $$其中 $\varepsilon_i$ 是独立 Rademacher 随机变量,并且与 $X_i$ 独立。
查看学习笔记:Exercise 8.11 证明令 $\mathcal F$ 是实线上所有形如 $[a,b]\cup[c,d]$ 的集合的指示函数类。证明
$$ \operatorname{vc}(\mathcal F)=4. $$ 查看学习笔记:Exercise 8.12 证明令 $\mathcal F$ 是 $\mathbb R^2$ 中所有圆的指示函数类。证明
$$ \operatorname{vc}(\mathcal F)=3. $$ 查看学习笔记:Exercise 8.13 证明令 $\mathcal F$ 是 $\mathbb R^2$ 中所有闭的坐标轴平行矩形,也就是 $[a,b]\times[c,d]$,的指示函数类。证明
$$ \operatorname{vc}(\mathcal F)=4. $$ 查看学习笔记:Exercise 8.14 证明令 $\mathcal F$ 是 $\mathbb R^2$ 中所有闭的坐标轴平行正方形,也就是 $[a,a+d]\times[b,b+d]$,的指示函数类。证明
$$ \operatorname{vc}(\mathcal F)=3. $$ 查看学习笔记:Exercise 8.15 证明令 $\mathcal F$ 是 $\mathbb R^2$ 中所有凸多边形的指示函数类,不限制顶点个数。证明
$$ \operatorname{vc}(\mathcal F)=\infty. $$ 查看学习笔记:Exercise 8.16 证明证明 Example 8.3.5 的结论:若 $\mathcal F$ 是 $\mathbb R^n$ 中所有半空间的指示函数类,则
$$ \operatorname{vc}(\mathcal F)=n+1. $$ 查看学习笔记:Exercise 8.17 证明证明对任意有限维布尔函数类 $\mathcal F$,
$$ \operatorname{vc}(\mathcal F)\le\dim(\mathcal F), $$其中 $\dim(\mathcal F)$ 表示线性代数意义下的维数,也就是 $\mathcal F$ 中最多能取多少个线性无关函数。
查看学习笔记:Exercise 8.18 证明考虑长度为 $n$ 且至多有 $d$ 个 $1$ 的二进制串集合 $\mathcal F$。证明 Pajor Lemma 8.3.7 与 Sauer-Shelah Lemma 8.3.9 对这个例子都是尖锐的。这个集合也称为 Hamming 球。
查看学习笔记:Exercise 8.19 证明令 $\mathcal F$ 是某个定义域上的布尔函数类。证明增长函数 $\Pi_{\mathcal F}(n)$ 的增长只能是多项式或指数,不存在中间情形。更精确地,证明要么 $\Pi_{\mathcal F}(n)=2^n$ 对所有 $n$ 成立,要么存在多项式 $p(n)$ 使得 $\Pi_{\mathcal F}(n)\le p(n)$ 对所有 $n$ 成立。
查看学习笔记:Exercise 8.20 证明推广 Proposition 8.3.11。令 $\mathcal F_1,\ldots,\mathcal F_k$ 是同一定义域上的布尔函数类。固定任意布尔公式 $\phi:\{0,1\}^k\to\{0,1\}$,并考虑
$$ \mathcal F= \{ \phi(f_1,\ldots,f_k): f_i\in\mathcal F_i \}. $$证明如果每个 $\operatorname{vc}(\mathcal F_i)\le d$,则
$$ \operatorname{vc}(\mathcal F)\le Cdk\log k. $$ 查看学习笔记:Exercise 8.21 证明(a) 令 $\mathcal F$ 与 $\mathcal G$ 是同一定义域上的布尔函数类。证明
$$ \operatorname{vc}(\mathcal F\cup\mathcal G) \le \operatorname{vc}(\mathcal F)+\operatorname{vc}(\mathcal G)+1. $$(b) 给出一个例子,使得 $\mathcal F$ 和 $\mathcal G$ 的 VC dimensions 都为正,并且上面的不等式取等号。
查看学习笔记:Exercise 8.22 证明Theorem 8.3.13 只对 $\varepsilon\in(0,1)$ 表述。问:当 $\varepsilon$ 更大时,覆盖数应满足什么界?
查看学习笔记:Exercise 8.23 证明不用 Pajor Lemma,而是直接使用 Sauer-Shelah Lemma,证明 Theorem 8.3.15 的一个较弱版本:右侧为
$$ C\sqrt{\frac dn}\log\frac{en}{d}, $$其中 $d=\operatorname{vc}(\mathcal F)$。
查看学习笔记:Exercise 8.24 证明学习高维分布很难,但可以用 $O(n)$ 个样本学习所有一维边缘分布。令 $X,X_1,\ldots,X_m$ 是 $\mathbb R^n$ 中的独立同分布随机向量。证明所有一维边缘分布 $\langle X,u\rangle$ 的 CDF 可由经验 CDF 一致估计:
$$ \mathbb E\sup_{u\in S^{n-1},\,t\in\mathbb R} \left| \frac1m\sum_{i=1}^m \mathbf 1_{\{\langle X_i,u\rangle\le t\}} - \mathbb P\{\langle X,u\rangle\le t\} \right| \lesssim \sqrt{\frac nm}. $$ 查看学习笔记:Exercise 8.25 证明量化单位向量 $u\in\mathbb R^n$ 的最简单方法是取每个坐标的 sign,但这可能很不准确。现在考虑更聪明的 random hyperplane quantization。
(a) 令 $A$ 是 $m\times n$ 标准高斯随机矩阵,定义 $\Phi:S^{n-1}\to\{-1,1\}^m$ 为
$$ \Phi(u)=\operatorname{sign}(Au). $$验证对所有 $u,v\in S^{n-1}$,
$$ \mathbb E\frac1m d(\Phi(u),\Phi(v)) = \frac1\pi\rho(u,v), $$其中 $d$ 是 Hamming distance,$\rho$ 是 sphere 上的 geodesic distance。
(b) 使用 VC 大数定律证明下面事件以至少 $0.99$ 的概率发生:
$$ \left| \frac1m d(\Phi(u),\Phi(v)) - \frac1\pi\rho(u,v) \right| \le C\sqrt{\frac nm} \quad\text{for all }u,v\in S^{n-1}. $$ 查看学习笔记:Exercise 8.26 证明大多数随机矩阵结果会假设次高斯尾部或至少有限二阶矩。现在用 VC 理论,特别是 Theorem 8.3.15,证明一个没有矩假设的高瘦随机矩阵可逆性结果。
(a) 令 $\varepsilon,\delta\gt 0$,令 $X$ 是 $\mathbb R^n$ 中的随机向量,满足反集中假设
$$ \mathbb P\{|\langle X,u\rangle|\ge\varepsilon\} \ge \delta \quad \forall u\in S^{n-1}. \tag{8.55} $$令 $A$ 是 $m\times n$ 随机矩阵,其行是 $X$ 的 i.i.d. 副本。证明若 $m\ge C\delta^{-2}n$,则以至少 $0.99$ 的概率,
$$ s_n(A)\ge0.99\,\varepsilon\sqrt{\delta m}. $$(b) 说明 (a) 中的界是最优的。对任意正整数 $m\ge n$ 和正数 $\varepsilon$、$\delta\lt 1/3$,构造一个满足 (8.55) 的随机向量 $X$,使得以它的 i.i.d. 副本为行的随机矩阵 $A$ 满足
$$ \mathbb E s_n(A)\le1.01\,\varepsilon\sqrt{\delta m}. $$ 查看学习笔记:Exercise 8.27 证明Remark 8.3.19 中提到:有限 VC 维数的布尔类是一致 Glivenko-Cantelli 类。证明反向:如果布尔类是一致 Glivenko-Cantelli 类,那么它的 VC 维数必须有限。
查看学习笔记:Exercise 8.28 证明令 $\mathcal F$ 是某集合 $\Omega$ 上的布尔函数类,令 $h$ 是 $\Omega$ 上的一个布尔函数。证明类
$$ \{(f-h)^2:f\in\mathcal F\} $$与 $\mathcal F$ 有相同的 VC 维数。
查看学习笔记:Exercise 8.29 证明第 8.4 节中我们假设标签 $T(X)$ 完全由 $X$ 决定;现实中标签常常只是与 $X$ 相关的随机变量。把学习理论扩展到训练数据
$$ (X_i,Y_i),\qquad i=1,\ldots,n, $$其中 $(X_i,Y_i)$ 是随机对 $(X,Y)$ 的独立副本,$X\in\Omega$ 且 $Y\in\mathbb R$。请把理论推进到 Theorem 8.4.5 对应的层次。
查看学习笔记:Exercise 8.30 证明我们只对布尔类证明了泛化界,但很多数据有实值标签。假设要从随机样本 $X_1,\ldots,X_n\sim\operatorname{Unif}[0,1]$ 上的取值学习一个 Lipschitz 函数 $T:[0,1]\to[0,1]$。考虑假设类
$$ \mathcal F= \{f:[0,1]\to[0,1]:\|f\|_{\mathrm{Lip}}\le1\}, $$并假设 $T\in\mathcal F$。证明第 8.4.2 节的 ERM 算法满足类似 Theorem 8.4.5 的界:
$$ \mathbb E R(f_n^*) \le R(f^*)+\frac C{\sqrt n}. $$ 查看学习笔记:Exercise 8.31 证明Theorem 8.4.5 暗示,要学得好,训练点数至少要和假设类的 VC 维数同阶。证明如果
$$ n\lt \frac12\operatorname{vc}(\mathcal F), $$则存在 $\Omega$ 上的一个分布,使得没有任何算法能够可靠地学习 $\mathcal F$ 中的目标函数。请把“可靠地学习”精确化。
查看学习笔记:Exercise 8.32 证明证明 $\gamma_2$ 泛函被 Dudley 积分控制:
$$ \gamma_2(T,d) \le C\int_0^\infty \sqrt{\log\mathcal N(T,d,\varepsilon)} \,d\varepsilon. $$ 查看学习笔记:Exercise 8.33 证明有时 $\gamma_2$ 泛函 (8.45) 可以显著小于 Dudley 求和 (8.44)。重新考虑 Exercise 8.4 中的例子:
$$ T= \{0\}\cup \left\{ \frac{e_k}{\sqrt{1+\log k}}: k=1,\ldots,n \right\}. $$(a) 构造可容许序列 $(T_k)$,证明 $\gamma_2(T,d)\le C$。
(b) 证明 Dudley 求和无界:
$$ \inf_{(T_k)} \sum_{k=0}^\infty 2^{k/2} \sup_{t\in T}d(t,T_k) \to\infty \quad\text{as }n\to\infty. $$ 查看学习笔记:Exercise 8.34 证明用带有第一步大跳跃的链式方法证明 Remark 8.5.4:
$$ t_0\to\pi_\kappa(t)\to\pi_{\kappa+1}(t)\to\pi_{\kappa+2}(t)\to\cdots\to t, $$其中 $\kappa$ 的选择满足 $u\asymp2^{\kappa/2}$。
查看学习笔记:Exercise 8.35 证明现在处理一般经验过程。令 $\mathcal F$ 是定义域 $\Omega$ 上的实值函数类,令 $X,X_1,\ldots,X_n$ 是独立同分布随机点。设 $d$ 是 $\mathcal F$ 上的度量,并满足
$$ \|f(X)-g(X)\|_{\psi_2} \le d(f,g) \quad\text{for all }f,g\in\mathcal F. $$例如 $d(f,g)=C\|f-g\|_{L^\infty}$ 可用。证明
$$ \mathbb E\sup_{f\in\mathcal F} \left| \frac1n\sum_{i=1}^n f(X_i)-\mathbb Ef(X) \right| \le \frac{C\gamma_2(\mathcal F,d)}{\sqrt n}. $$ 查看学习笔记:Exercise 8.36 证明探索 Corollary 8.5.8 的几个有用变体。令 $(X_x)_{x\in T}$ 是 $T\subset\mathbb R^n$ 上的随机过程,不要求均值为零,并令 $X_0=0$。假设
$$ \|X_x-X_y\|_{\psi_2} \le K\|x-y\|_2 \quad \text{for all }x,y\in T\cup\{0\}. $$(a) 证明期望界:
$$ \mathbb E\sup_{x\in T}|X_x| \le CK\gamma(T). $$(b) 证明高概率界:对任意 $u\ge0$,
$$ \sup_{x\in T}|X_x| \le CK\bigl(w(T)+u\operatorname{rad}(T)\bigr) $$以概率至少 $1-2e^{-u^2}$ 成立。
(c) 推出矩界:
$$ \left( \mathbb E\sup_{x\in T}|X_x|^p \right)^{1/p} \le C\sqrt p\,K\gamma(T). $$ 查看学习笔记:Exercise 8.37 证明Exercise 3.5 在独立次高斯坐标情形给出了 $\ell^p$ 范数期望的最优界。现在去掉独立性。对 $\mathbb R^N$ 中的次高斯随机向量 $X$ 和 $p\in[1,\infty]$,证明
$$ \mathbb E\|X\|_p \le \begin{cases} CK\sqrt p\,N^{1/p},&1\le p\le\log N,\\ CK\sqrt{\log N},&\log N\le p\le\infty, \end{cases} $$其中 $K=\|X\|_{\psi_2}$。
查看学习笔记:Exercise 8.38 证明令 $A$ 是 $m\times n$ 随机矩阵,其元素独立且服从 $N(0,1)$。令 $T\subset\mathbb R^n$、$S\subset\mathbb R^m$ 为任意有界集。
(a) 证明 Theorem 8.6.1 在高斯情形下可取尖锐常数 $1$:
$$ \mathbb E\sup_{x\in T,y\in S}\langle Ax,y\rangle \le w(T)\operatorname{rad}(S)+w(S)\operatorname{rad}(T). $$(b) 证明反向不等式:
$$ \mathbb E\sup_{x\in T,y\in S}\langle Ax,y\rangle \ge c\left[ w(T)\operatorname{rad}(S)+w(S)\operatorname{rad}(T) \right]. $$ 查看学习笔记:Exercise 8.39 证明在 Theorem 8.6.1 的假设下,为
$$ \sup_{x\in T,y\in S}\langle Ax,y\rangle $$证明一个尾界。
查看学习笔记:Exercise 8.40 证明现在可以处理所有 $p,q$,而不只是少数特殊情形。令 $A$ 是 $m\times n$ 随机矩阵,其行 $A_i$ 独立、均值为零且次高斯。取任意 $p,q\in[1,\infty]$。证明
$$ \mathbb E\|A\|_{p\to q} \le CK\left[ r(n,p)w(m,q)+r(m,q')w(n,p') \right], $$其中 $K=\max_i\|A_i\|_{\psi_2}$,$p'$ 与 $q'$ 是 Hölder conjugates,并且
$$ r(n,p)= \begin{cases} 1,&1\le p\le2,\\ n^{1/2-1/p},&2\le p\le\infty, \end{cases} \qquad w(n,p)= \begin{cases} \sqrt p\,n^{1/p},&1\le p\le\log n,\\ \sqrt{\log n},&\log n\le p\le\infty. \end{cases} $$证明这个界对高斯矩阵是紧的:若 $A$ 的元素独立且服从 $N(0,1)$,则
$$ \mathbb E\|A\|_{p\to q} \ge c\left[ r(n,p)w(m,q)+r(m,q')w(n,p') \right]. $$ 查看学习笔记:Exercise 8.41 证明