形式陈述
一个 Lasso 候选可以已经把训练目标解得很准,却仍离真实系数很远。要把优化保证转成统计保证,先把观测拆成信号与噪声,再问惩罚能否压住噪声与每列设计的相关性。本页给出这条比较链的起点。
采用固定设计的线性回归模型 理路 线性回归统计模型 Linear regression model · Linear model 响应的条件均值由设计变量对未知系数线性表示,并显式规定误差结构的统计模型。
Y = X β 0 + ε , X ∈ R n × p , n , p ≥ 1. 不另加截距;若需要未惩罚截距,应先明确处理方式,再对被惩罚的设计使用以下结论。向量的 ℓ 1 范数是各坐标绝对值之和,ℓ ∞ 范数是最大绝对值。固定 λ > 0 ,定义
Q λ ( β ) = ‖ Y − X β ‖ 2 2 2 n + λ ‖ β ‖ 1 . 设算法返回可测候选 β ~ ,并有经过认证的容差 δ ≥ 0 ,使
(1) Q λ ( β ~ ) ≤ min β Q λ ( β ) + δ . 目标连续且 λ ‖ β ‖ 1 在无穷远趋于无穷,故最小值存在。Lasso 对偶间隙 理路 Lasso 的最优性与对偶间隙 Lasso optimality conditions · Lasso duality gap · Lasso primal-dual certificate 从残差相关性核验 Lasso 的零与非零坐标,并用可行对偶值认证剩余优化误差。 采用不除以 n 的损失;将这里的目标乘以 n ,对应旧页惩罚 n λ 。若那里的可行间隙为 G ,这里可取 δ = G / n ,不能直接把 G 当作 δ 。
记
h = β ~ − β 0 , q 2 = ‖ X h ‖ 2 2 n , W = X T ε n . q 2 是这些设计点上的均值预测误差;它不是系数误差,也还不是任意新输入分布上的风险。由式(1)和真参数作为比较器,得到基本不等式
(2) 1 2 q 2 + λ ( ‖ β ~ ‖ 1 − ‖ β 0 ‖ 1 ) ≤ W T h + δ . 在噪声事件
(3) E λ = { ‖ W ‖ ∞ ≤ λ / 2 } 上,不需要任何满秩或受限特征值条件,就有
(4) q 2 + λ ‖ β ~ ‖ 1 ≤ 3 λ ‖ β 0 ‖ 1 + 2 δ . 若真支持 S = { j : β j 0 ≠ 0 } ,记 h S 为只保留 S 内坐标的向量,S c 为其补集,则同一事件上还有更精细的结论
(5) q 2 + λ ‖ h S c ‖ 1 ≤ 3 λ ‖ h S ‖ 1 + 2 δ . 式(5)在 δ = 0 时给出锥约束 ‖ h S c ‖ 1 ≤ 3 ‖ h S ‖ 1 ;在 δ > 0 时只给出带截距的松弛约束。后续使用受限设计条件 理路 受限特征值与稀疏可识别性 Restricted eigenvalue condition · Compatibility condition for Lasso · Lasso 兼容条件 在稀疏误差锥上定量比较预测与参数范数,精确定标兼容常数和受限特征值,并处理非零优化容差。 时,必须保留这一区别。
直觉
展开平方后,噪声的影响只剩 W T h :算法若沿某列改变系数,能利用多少噪声,取决于噪声与该列的相关性。惩罚在每个坐标上都收取 λ 倍绝对值,因此选择 λ 压住全部坐标相关性,就能让结构惩罚承担这项随机扰动。
稀疏性的作用是让惩罚差具有方向。真支持以外原来都是零,新增一个系数一定增加绝对值惩罚;真支持以内,系数向零移动却可能节省惩罚。式(5)比较的正是这两份费用。它并未说算法知道 S ,只是分析时用真实支持把误差分成两部分。
优化容差 δ 是算法尚未排除的额外预算。只要容差非零,候选就可能拿其中一部分预算在真支持外放入小系数。因此“目标非常接近最优”和“误差严格处于同一个齐次锥”不能互换。
两次三角不等式完成证明
因为 Y − X β ~ = ε − X h ,目标比较展开为
‖ ε − X h ‖ 2 2 − ‖ ε ‖ 2 2 2 n + λ ( ‖ β ~ ‖ 1 − ‖ β 0 ‖ 1 ) ≤ δ . 平方差等于 q 2 / 2 − W T h ,这就证明式(2)。再用$\ell_\infty$–$\ell_1$ Hölder 不等式 理路 Hölder 不等式 Hölder's inequality 共轭指数下函数乘积的 L¹ 范数由各自的 Lᵖ 范数乘积控制。
W T h ≤ ‖ W ‖ ∞ ‖ h ‖ 1 ≤ λ 2 ‖ h ‖ 1 . 为得到不依赖稀疏性的式(4),只需用
‖ h ‖ 1 ≤ ‖ β ~ ‖ 1 + ‖ β 0 ‖ 1 ,移项并乘以二。为得到式(5),则用
‖ β 0 ‖ 1 − ‖ β 0 + h ‖ 1 ≤ ‖ h S ‖ 1 − ‖ h S c ‖ 1 , 以及 ‖ h ‖ 1 = ‖ h S ‖ 1 + ‖ h S c ‖ 1 。支持内的系数由三角不等式控制,支持外则因为 β S c 0 = 0 而得到精确的惩罚增加。证明至此完全是确定性代数;概率只用来说明事件(3)多常发生。
例子与边界
相同的预测,可以对应完全不同的坐标
设两列完全相同且 ‖ X 1 ‖ 2 2 / n = 1 ,取 Y = X 1 、β 0 = ( 1 , 0 ) T 、无噪声、λ = 1 / 4 。目标仅通过 s = β 1 + β 2 使用拟合值。对非负系数,绝对值惩罚也只取决于 s ,所以全部
β ~ = ( a , 3 / 4 − a ) T , 0 ≤ a ≤ 3 / 4 都是精确最优解。它们的预测误差全部是 q 2 = ( 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 不能省掉松弛量
取 Γ = X T X / n = I 2 、Y = X ( 1 , 0 ) T 、λ = 1 / 10 。精确解是 β ^ = ( 9 / 10 , 0 ) T 。候选
β ~ = ( 9 / 10 , 2 / 5 ) T 的目标差为
δ = 1 2 ( 2 / 5 ) 2 + 1 10 2 5 = 3 25 . 此时 h = ( − 1 / 10 , 2 / 5 ) T ,所以 ‖ h S c ‖ 1 = 2 / 5 > 3 / 10 = 3 ‖ h S ‖ 1 ,不在精确锥内。令 0 < λ < 1 趋于零,并取候选 ( 1 − λ , 4 λ ) ,目标差变为 12 λ 2 → 0 ,同样的锥违背仍存在。绝对 gap 很小不会自动恢复精确解的几何条件。
训练设计上的预测,何时能迁移
若未来输入 x ∗ 与拟合资料独立,且二阶矩矩阵为 Σ ∗ = E [ x ∗ x ∗ T ] ,给定候选后的未来均值预测误差是 h T Σ ∗ h 。这里重新记 Γ = X T X / n 。若存在有限的 C ≥ 0 满足矩阵序关系 Σ ∗ ⪯ C Γ ,便得到
h T Σ ∗ h ≤ C h T Γ h = C q 2 . 没有这样的覆盖条件,式(4)只回答原设计上的问题。上面的重复列训练资料只看到系数之和,新输入 ( 1 , 0 ) 却单独读取第一坐标;在这个新输入上,训练不可识别的方向会重新出现。新增响应噪声还须另计,不能把均值误差当成完整预测损失。
推论与应用
用尾部条件选择惩罚
假设给定 X 后,噪声 ε i 相互独立、均值为零,且是共同尺度上界为 σ > 0 的次高斯随机变量 理路 次高斯随机变量 Sub-Gaussian random variable · 次高斯尺度 · Subgaussian variance proxy 用全实数上的中心化指数矩上界定义次高斯尺度,计算独立加权平均的尾界,并区分方差、尺度代理与条件次高斯假设。 :
E [ e t ε i ∣ X ] ≤ e σ 2 t 2 / 2 。再假定每列 ‖ X j ‖ 2 2 / n ≤ 1 。线性组合的指数矩相乘,对 u > 0 给
Pr ( | W j | > u ∣ X ) ≤ 2 exp ( − n u 2 2 σ 2 ) . 对 p 列作并集界 理路 并集界 Union bound · Boole 不等式 多个坏事件中至少一个发生的概率,不超过各事件概率之和。 ,不要求不同 W j 独立。给定 0 < a < 1 ,选择
(6) λ = 2 σ 2 log ( 2 p / a ) n 便以至少 1 − a 的条件概率保证事件(3),从而同时保证式(4)、(5)。若列范数上界为 n L ,式(6)再乘 L ;未经归一化便沿用原数值 λ 会改变结论。只有有限方差时,不能直接使用这条指数尾公式。
例如 n = 100 , p = 10 , a = 0.05 , σ = 1 ,式(6)给 λ ≈ 0.692327 。若已知 ‖ β 0 ‖ 1 ≤ 2 且 δ ≤ 0.01 ,式(4)给 q 2 ≤ 6 λ + 0.02 ≈ 4.173964 。这份界可能保守,但每个条件和尺度都能核对;“看起来噪声不大”不能替代式(3)的概率保证。
近似稀疏与计算预算
对任意事先指定的坐标集 T ,同样分解得到
(7) q 2 + λ ‖ h T c ‖ 1 ≤ 3 λ ‖ h T ‖ 1 + 4 λ ‖ β T c 0 ‖ 1 + 2 δ . 这里 ‖ β T c 0 ‖ 1 是被忽略的真实尾部,来自
‖ β 0 ‖ 1 − ‖ β 0 + h ‖ 1 ≤ ‖ h T ‖ 1 − ‖ h T c ‖ 1 + 2 ‖ β T c 0 ‖ 1 。它与优化误差是两项不同的费用;继续迭代只能减少后者。
在真稀疏情形,受限特征值与兼容常数 理路 受限特征值与稀疏可识别性 Restricted eigenvalue condition · Compatibility condition for Lasso · Lasso 兼容条件 在稀疏误差锥上定量比较预测与参数范数,精确定标兼容常数和受限特征值,并处理非零优化容差。 把式(5)进一步变为 λ 2 | S | 量级的预测界和 λ | S | 量级的参数界,并给出非零 δ 时如何放大锥。支持恢复 理路 Lasso 支持恢复与不可表示条件 Lasso support recovery · Lasso sign consistency · Lasso irrepresentable condition 用活动集上的显式候选和非活动集的严格对偶余量证明符号恢复,区分可识别性、参数小误差、变量筛选和坐标推断。 还需检查最小信号和非活动列的得分余量。反方向看,小对偶间隙能验证式(1),却不能计算含未知真参数的 q 、S 或噪声事件;它是统计证明的一个输入,不是其余条件的替代品。
完整的同设计、重复列与有限迭代比较见正则化学习保证终点练习 。
参考资料