Skip to content

定理Theorem

Lasso 基本不等式与预测误差

Lasso basic inequality · Lasso prediction error bound

从带优化容差的目标比较推导预测界与松弛锥,分清噪声事件、稀疏结构、设计条件和可计算证书各自的作用。

形式陈述 ​

一个 Lasso 候选可以已经把训练目标解得很准,却仍离真实系数很远。要把优化保证转成统计保证,先把观测拆成信号与噪声,再问惩罚能否压住噪声与每列设计的相关性。本页给出这条比较链的起点。

采用固定设计的线性回归模型

Y=Xβ0+ε,X∈Rn×p,n,p≥1.

不另加截距;若需要未惩罚截距,应先明确处理方式,再对被惩罚的设计使用以下结论。向量的 ℓ1 范数是各坐标绝对值之和,ℓ∞ 范数是最大绝对值。固定 λ>0,定义

Qλ(β)=‖Y−Xβ‖222n+λ‖β‖1.

设算法返回可测候选 β~,并有经过认证的容差 δ≥0,使

(1)Qλ(β~)≤minβQλ(β)+δ.

目标连续且 λ‖β‖1 在无穷远趋于无穷,故最小值存在。Lasso 对偶间隙采用不除以 n 的损失;将这里的目标乘以 n,对应旧页惩罚 nλ。若那里的可行间隙为 G,这里可取 δ=G/n,不能直接把 G 当作 δ。

记

h=β~−β0,q2=‖Xh‖22n,W=XTεn.

q2 是这些设计点上的均值预测误差;它不是系数误差,也还不是任意新输入分布上的风险。由式(1)和真参数作为比较器,得到基本不等式

(2)12q2+λ(‖β~‖1−‖β0‖1)≤WTh+δ.

在噪声事件

(3)Eλ={‖W‖∞≤λ/2}

上,不需要任何满秩或受限特征值条件,就有

(4)q2+λ‖β~‖1≤3λ‖β0‖1+2δ.

若真支持 S={j:βj0≠0},记 hS 为只保留 S 内坐标的向量,Sc 为其补集,则同一事件上还有更精细的结论

(5)q2+λ‖hSc‖1≤3λ‖hS‖1+2δ.

式(5)在 δ=0 时给出锥约束 ‖hSc‖1≤3‖hS‖1;在 δ>0 时只给出带截距的松弛约束。后续使用受限设计条件时,必须保留这一区别。

直觉

展开平方后,噪声的影响只剩 WTh:算法若沿某列改变系数,能利用多少噪声,取决于噪声与该列的相关性。惩罚在每个坐标上都收取 λ 倍绝对值,因此选择 λ 压住全部坐标相关性,就能让结构惩罚承担这项随机扰动。

稀疏性的作用是让惩罚差具有方向。真支持以外原来都是零,新增一个系数一定增加绝对值惩罚;真支持以内,系数向零移动却可能节省惩罚。式(5)比较的正是这两份费用。它并未说算法知道 S,只是分析时用真实支持把误差分成两部分。

优化容差 δ 是算法尚未排除的额外预算。只要容差非零,候选就可能拿其中一部分预算在真支持外放入小系数。因此“目标非常接近最优”和“误差严格处于同一个齐次锥”不能互换。

两次三角不等式完成证明 ​

因为 Y−Xβ~=ε−Xh,目标比较展开为

‖ε−Xh‖22−‖ε‖222n+λ(‖β~‖1−‖β0‖1)≤δ.

平方差等于 q2/2−WTh,这就证明式(2)。再用$\ell_\infty$–$\ell_1$ Hölder 不等式

WTh≤‖W‖∞‖h‖1≤λ2‖h‖1.

为得到不依赖稀疏性的式(4),只需用 ‖h‖1≤‖β~‖1+‖β0‖1,移项并乘以二。为得到式(5),则用

‖β0‖1−‖β0+h‖1≤‖hS‖1−‖hSc‖1,

以及 ‖h‖1=‖hS‖1+‖hSc‖1。支持内的系数由三角不等式控制,支持外则因为 βSc0=0 而得到精确的惩罚增加。证明至此完全是确定性代数;概率只用来说明事件(3)多常发生。

例子与边界

相同的预测,可以对应完全不同的坐标 ​

设两列完全相同且 ‖X1‖22/n=1,取 Y=X1、β0=(1,0)T、无噪声、λ=1/4。目标仅通过 s=β1+β2 使用拟合值。对非负系数,绝对值惩罚也只取决于 s,所以全部

β~=(a,3/4−a)T,0≤a≤3/4

都是精确最优解。它们的预测误差全部是 q2=(1−3/4)2=1/16;取 a=0 时,参数误差却是 h=(−1,3/4)T,其 ℓ1 范数为 7/4,支持也选到了另一列。

式(5)没有失效:左侧是 1/16+(1/4)(3/4)=1/4,右侧是 3(1/4)(1)=3/4。它控制了一项本来就很小的预测误差,却没有把不可区分的两列分开。这正是进一步引入设计条件的原因。

非零 gap 不能省掉松弛量 ​

取 Γ=XTX/n=I2、Y=X(1,0)T、λ=1/10。精确解是 β^=(9/10,0)T。候选 β~=(9/10,2/5)T 的目标差为

δ=12(2/5)2+11025=325.

此时 h=(−1/10,2/5)T,所以 ‖hSc‖1=2/5>3/10=3‖hS‖1,不在精确锥内。令 0<λ<1 趋于零,并取候选 (1−λ,4λ),目标差变为 12λ2→0,同样的锥违背仍存在。绝对 gap 很小不会自动恢复精确解的几何条件。

训练设计上的预测,何时能迁移 ​

若未来输入 x∗ 与拟合资料独立,且二阶矩矩阵为 Σ∗=E[x∗x∗T],给定候选后的未来均值预测误差是 hTΣ∗h。这里重新记 Γ=XTX/n。若存在有限的 C≥0 满足矩阵序关系 Σ∗⪯CΓ,便得到

hTΣ∗h≤ChTΓh=Cq2.

没有这样的覆盖条件,式(4)只回答原设计上的问题。上面的重复列训练资料只看到系数之和,新输入 (1,0) 却单独读取第一坐标;在这个新输入上,训练不可识别的方向会重新出现。新增响应噪声还须另计,不能把均值误差当成完整预测损失。

推论与应用

用尾部条件选择惩罚 ​

假设给定 X 后,噪声 εi 相互独立、均值为零,且是共同尺度上界为 σ>0 的次高斯随机变量: E[etεi∣X]≤eσ2t2/2。再假定每列 ‖Xj‖22/n≤1。线性组合的指数矩相乘,对 u>0 给

Pr(|Wj|>u∣X)≤2exp(−nu22σ2).

对 p 列作并集界,不要求不同 Wj 独立。给定 0<a<1,选择

(6)λ=2σ2log⁡(2p/a)n

便以至少 1−a 的条件概率保证事件(3),从而同时保证式(4)、(5)。若列范数上界为 nL,式(6)再乘 L;未经归一化便沿用原数值 λ 会改变结论。只有有限方差时,不能直接使用这条指数尾公式。

例如 n=100,p=10,a=0.05,σ=1,式(6)给 λ≈0.692327。若已知 ‖β0‖1≤2 且 δ≤0.01,式(4)给 q2≤6λ+0.02≈4.173964。这份界可能保守,但每个条件和尺度都能核对;“看起来噪声不大”不能替代式(3)的概率保证。

近似稀疏与计算预算 ​

对任意事先指定的坐标集 T,同样分解得到

(7)q2+λ‖hTc‖1≤3λ‖hT‖1+4λ‖βTc0‖1+2δ.

这里 ‖βTc0‖1 是被忽略的真实尾部,来自 ‖β0‖1−‖β0+h‖1≤‖hT‖1−‖hTc‖1+2‖βTc0‖1。它与优化误差是两项不同的费用;继续迭代只能减少后者。

在真稀疏情形,受限特征值与兼容常数把式(5)进一步变为 λ2|S| 量级的预测界和 λ|S| 量级的参数界,并给出非零 δ 时如何放大锥。支持恢复还需检查最小信号和非活动列的得分余量。反方向看,小对偶间隙能验证式(1),却不能计算含未知真参数的 q、S 或噪声事件;它是统计证明的一个输入,不是其余条件的替代品。

完整的同设计、重复列与有限迭代比较见正则化学习保证终点练习。

参考资料
  • Peter J. Bickel, Ya'acov Ritov and Alexandre B. Tsybakov, Simultaneous Analysis of Lasso and Dantzig Selector, Annals of Statistics 37(4), 2009,§3 与 Appendix B:噪声相关性、基本比较和受限设计。该文采用加权列范数与另一目标定标;本文统一为半平均平方损失,式(1)中的加性容差和式(7)均在正文直接推导。
  • Sara van de Geer and Peter Bühlmann, On the Conditions Used to Prove Oracle Results for the Lasso, Electronic Journal of Statistics 3, 2009,§§2、11:从基本不等式到兼容条件和有噪声情形。
关系图谱21 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系