两个取值都在 的变量,可能一个几乎总为零,另一个经常在两端跳动。只看取值范围,Hoeffding 给二者同样的界;如果已知前者方差很小,Bernstein 还能利用这份信息。
本页是 矩阵 Bernstein理路矩阵 Bernstein 不等式Matrix Bernstein inequality · 矩阵伯恩斯坦不等式用逐项谱范数与矩阵方差控制独立中心化 Hermitian 矩阵和,借助 Lieb 凹性完成非交换指数矩证明,并计算随机图邻接矩阵的统一方向误差。 的 情形,公式并不新。独立入口的作用是让只学过概率的读者完成证明、反解样本量,而不必先经过矩阵谱理论与 Lieb 凹性。把实数 看作 Hermitian 矩阵时,谱范数就是 ,矩阵方差就是 。若用较大的 ,右端随方差参数增加而增大,因此仍合法。矩阵页中同时控制全部方向的责任仍留在原页。
形式陈述
两份预算必须同时成立
设整数 ,实变量 相互 独立理路独立性Statistical independence从概率表理解独立性,区分两两、相互和条件独立,并用可计算反例澄清零协方差与条件均值的限度。,。假设存在确定的 ,使 几乎处处,并有确定的 方差理路方差Variance随机变量相对其均值的平方偏差期望,也是最佳常数平方预测的剩余误差。预算
当 时,对任意 ,
与 可以是保守上界,但必须由模型或已证事实支持。若 ,各 几乎处处,任何正阈值的尾概率都为零;无需把零方差代进后面的优化参数。若原本 ,也直接使用这个退化结论。
控制单次可能的冲击, 控制累计平方波动。两者有不同单位: 与 相同, 带平方单位,分母中 才能相加。
对原始观测 ,中心化后 可以统一取 。通常不能自动取 :例如 Bernoulli 的一次取值 离均值为 ,超过半个区间宽度。
直觉
证明:为什么幅度与方差一起进入指数矩
先看一个中心化变量 ,。对整数 ,
后一个不等式在 取等,其后每增加一项,阶乘至少再乘 。利用指数级数和 ,当 时有
有界,指数级数可以用确定的可积常数 控制后逐项取期望。再用 ,得到
把每项的 矩母函数理路矩母函数Moment-generating function · MGF在存在邻域内以 E[e^{tX}] 编码随机变量各阶矩的函数。相乘,独立性将对数指数矩相加。因此对 和 ,
Chernoff 方法理路Chernoff 方法与 Chernoff 界Chernoff method · Chernoff bounds从指数矩与 Markov 不等式推导尾界,给出独立 Bernoulli 和的乘法形式、KL 形式及适用条件。便给出
取合法参数 ;因 ,它严格小于 。此时 ,代回指数得到
这只是一个方便而有效的参数选择,不需要声称它精确最小化前面整个有理函数。对 重复同一证明,再用 并集界理路并集界Union bound · Boole 不等式多个坏事件中至少一个发生的概率,不超过各事件概率之和。合并两侧,即得式 (1)。
从尾界算出可实施的保证
令 ,。反解 ,正根为
更好记但稍松的充分阈值为
式 (4) 不小于式 (3),因为 。这里的线性系数为 ;不要与另一些直接由更细 MGF 优化得到的界的 版本混用。
对于 IID 观测,共同均值 ,,已知 ,取 得
因此一个充分样本量是
若确定 ,直接用一个样本的退化结论,比继续使用式 (5) 更合适。所有这些结论针对预先固定的样本数;反复检查区间并在喜欢的时刻停止,需要额外的全时间校准。
例子与边界
罕见故障:同一个任务节省多少样本
设每次检测独立,故障指标 。未知的 满足事先给定的模型限制 。目标是把故障率的绝对误差控制在 内,失败概率不超过 。
范围是 ,所以中心化幅度统一取 。由于 在 上递增,可以取 。依次计算
所以取 已满足式 (5)。只使用范围的 Hoeffding 界则要求
差异来自“故障率至多百分之一”所提供的有效方差预算。这个限制必须独立于此次观测成立;看了少量零故障后自行宣称 ,不能获得同样保证。没有这份先验限制时,可用 经验 Bernstein理路经验 Bernstein 均值界Empirical Bernstein mean bound · 经验方差均值界对已知有界范围的 IID 观测,以经过校准的样本方差构造固定样本均值区间,解释标准差估计误差、零样本方差与重复查看的边界。 支付方差估计本身的不确定性。
这一算例的 仍然未知,所以估计任务没有消失。上面比较的是充分样本量,不是精确最小样本量,也不是把置信区间长度当成真实误差。
方差区与单次冲击区
当 时,分母 ,尾界不超过 。这时偏差主要按平方尺度受控。
当 时,分母不超过 ,得到 。此时幅度预算限制了单次冲击,总体上界呈线性指数尺度。两段是在解释同一个公式的主导项,并不是声称真实分布的尾在 处突然改变形状。
ψ₁ 变量理路次指数尾随机变量(ψ₁)Sub-exponential random variable (Orlicz) · Psi-1 random variable · ψ₁ 次指数尾用绝对值指数矩定义 ψ₁ 尺度,推导局部 MGF 与两尺度和界,并以指数分布和 Pareto 尾辨析同名的重尾卷积类。也有类似的两尺度和界,但它未必有有限的逐项幅度 。在那里必须用局部指数矩参数替代这里的硬边界,不能把从样本中看见的最大值当成总体幅度上界。
缺少哪一项,会坏在哪里
只有有限方差而没有幅度或合适指数矩,不能得到式 (1)。一个方差固定、偶尔极大的变量,仍可令样本平均以远大于所套指数界的概率失准;MoM理路分组均值中位数估计Median-of-means estimator · MoM 均值估计将独立样本分成不相交的组,先平均再取中位数,以有限方差取得依赖置信度的均值误差保证,并明确组数、余数和尺度条件。通过改变估计量处理这种情况。
即使每项都满足同样的幅度与方差限制,依赖也会破坏和的保证。若所有项都是同一个公平符号 ,则 ,。错误地把它们当独立会设置 ,在 时给出 。实际和的方差为 ,独立乘法步骤已经失败。
另一方面,Bernstein 不总比 Hoeffding 更小。若只知道 ,使用最坏方差 ,式 (5) 多出的线性项会使它比 Hoeffding 的充分样本量更保守。拥有更多形式不意味着每次都应该选同一个公式;应比较在当前假设下合法的界。
推论与应用
复算与迁移
- 对 个独立中心化项,已知各项方差至多 、幅度至多 。取 ,分别用式 (3)、(4) 计算总和阈值,为什么较松式仍然有效?
- 改用两条生产线的检测数据:A线有 个 Bernoulli 观测,已知 ;B线有 个 Bernoulli 观测,已知 。两组内及组间的全部观测相互独立,但不要求 。把这 个指标直接平均,估计的目标是什么?写出中心化总和、方差预算,并在 下给平均误差的半径。能否沿用上一例 的共同方差预算?
核对第一题:,精确反解阈值约为 ,方便式约为 ;后者对平方根作了合法放松。
第二题的平均值估计
不是自动估计两条生产线故障率的等权平均 。写 为两组指标,则中心化总和为
每项中心化幅度可取 。两个参数上界都小于 ,方差 在这些区间内递增;独立性使方差相加,得到
现在直接用独立总和的式 (3),不调用 IID 均值简式。取 ,
所以平均误差半径 ,满足 ;用方便式 (4) 则给稍大的半径 。保证无需两条生产线有相同故障率。
不能沿用旧的共同预算 :例如允许的 已使总方差等于 。同时,也不能只把估计目标的名字改成共同的 ;分组采样比例决定了这里的目标权重。若想估计两条线的等权平均,可以改用 ,再按不同权重重新计算幅度与方差预算。
放进自然滤过,得到可预测方差的特例
本页独立中心化变量令 ,并采用自然滤过 ,则 ,且条件方差为确定的 。因而标量 Freedman 不等式理路标量 Freedman 不等式Scalar Freedman inequality · 可预测方差鞅尾界用有界鞅差及可预测条件方差控制预算内的全时间越界,证明指数过程并给出适应门控与分层方差预算的实际计算。在固定时刻、确定方差预算下包含本页尾界;原有 还允许对 与 分别使用 Freedman,再取并集界,恢复本页双侧式前的因子2。原有标量证明仍提供不必先学鞅的直接入口。若滤过预先揭示其他有关未来尺度的信息,独立增量之间的关系本身并不足以保证条件方差还是确定数,需要重新检查对整个过去的信息独立性。
参考资料