形式陈述
设计矩阵把参数误差 h 送到预测误差 X h 。若不同参数产生几乎相同的预测,单凭拟合就难以辨认系数。受限设计条件只对统计证明实际会遇到的误差方向要求定量可辨,而不要求整个高维参数空间都可逆。
给定实矩阵 理路 矩阵 Matrix 以有限行列集合为索引、取值于半环,并以中间指标求和定义乘法的函数。 X ∈ R n × p ,其中 n , p ≥ 1 ,记 Γ = X T X / n 。固定 ∅ ≠ S ⊆ { 1 , … , p } 、s = | S | 和 c ≥ 0 ,定义
(1) C ( S , c ) = { h ∈ R p : ‖ h S c ‖ 1 ≤ c ‖ h S ‖ 1 } . 这里使用坐标的 ℓ 1 与欧氏范数 理路 赋范向量空间 Normed vector space 带满足正定、齐次与三角不等式范数的向量空间。 ;h S 只保留 S 内坐标。集合是齐次的:把 h 放大不会改变其成员资格。非零成员必有 h S ≠ 0 ,因此下面分母都非零。
本页把三个常数全部写成平方量 ,避免把分母或平方的位置混在一起:
(2) κ 2 ( S , c ) = inf h ∈ C ( S , c ) ∖ { 0 } h T Γ h ‖ h S ‖ 2 2 , ϕ 2 ( S , c ) = inf h ∈ C ( S , c ) ∖ { 0 } s h T Γ h ‖ h S ‖ 1 2 , (3) α 2 ( S , c ) = inf h ∈ C ( S , c ) ∖ { 0 } h T Γ h ‖ h ‖ 2 2 . κ 是以支持内欧氏误差为分母的受限特征值常数;ϕ 是兼容常数;α 在本页专指全向量欧氏误差版本。文献可能把其中不同量都简称 RE,使用一个定理时应先核分母,而不是只认名称。记 κ , ϕ , α 为非负平方根。
给定整数 1 ≤ s 0 ≤ p ,若希望对所有至多 s 0 稀疏的真参数统一保证,还须对全部 1 ≤ | S | ≤ s 0 取最小值。本页先给固定 S 的结果;实际不知道真支持,并不会使一个特定 S 的计算自动变成全支持保证。
精确 Lasso 解在噪声事件 ‖ X T ε / n ‖ ∞ ≤ λ / 2 上,由基本不等式 理路 Lasso 基本不等式与预测误差 Lasso basic inequality · Lasso prediction error bound 从带优化容差的目标比较推导预测界与松弛锥,分清噪声事件、稀疏结构、设计条件和可计算证书各自的作用。 落入 C ( S , 3 ) 。若 ϕ 3 = ϕ ( S , 3 ) > 0 ,则
(4) ‖ X ( β ^ − β 0 ) ‖ 2 2 n ≤ 9 λ 2 s ϕ 3 2 , ‖ β ^ − β 0 ‖ 1 ≤ 12 λ s ϕ 3 2 . 这些结论要求 β S c 0 = 0 ,目标是半平均平方损失加 λ ‖ β ‖ 1 。若只知道目标差至多 δ > 0 ,不能原样使用式(4);下文给出可直接接纳 δ 的版本。
直觉
无约束的最小特征值遍历所有方向。高维时 p > n ,总有某个非零向量落在 ker X ,所以全空间最小特征值是零。稀疏分析并不要求每条方向都能被观察到,而是先用惩罚比较证明误差不能主要堆在真支持之外,再在剩余的锥上检查预测能否控制参数。
这套顺序有两项独立责任。基本不等式负责证明“误差属于哪里”,设计条件负责证明“在这片区域内,小预测意味着多小的参数误差”。若把近似解当作精确解,第一项责任已经失效,第二项再漂亮也不能补上。
常数之间能比较什么
对任意 h ,Cauchy–Schwarz 不等式 理路 Cauchy–Schwarz 不等式 Cauchy–Schwarz inequality · 柯西–施瓦茨不等式 内积的绝对值不超过两向量范数之积,且等号精确刻画线性相关。 给
‖ h S ‖ 1 ≤ s ‖ h S ‖ 2 。同时 ‖ h S ‖ 2 ≤ ‖ h ‖ 2 ,所以
(5) α 2 ( S , c ) ≤ κ 2 ( S , c ) ≤ ϕ 2 ( S , c ) . 若全局最小特征值为 λ min ( Γ ) > 0 ,则
α 2 ≥ λ min ( Γ ) 。反过来不成立:受限锥可以完全避开全局零空间。
对固定有限 s , c ,锥上还有
‖ h ‖ 2 2 ≤ ( 1 + c 2 s ) ‖ h S ‖ 2 2 和
‖ h S ‖ 1 2 ≥ ‖ h S ‖ 2 2 ,因此
κ 2 ≤ ( 1 + c 2 s ) α 2 , ϕ 2 ≤ s κ 2 . 所以本页三个固定支持常数是否严格为正是等价的,但其数值以及随 s , n 变化的统一下界并不相同。拿一个 κ 界直接当作同数值的全 ℓ 2 界,会丢掉这些因子。文献中更大支持集版本的 RE 还改变了定义域,不能只靠这里的比较替换。
从锥上的不等式走到误差界
令 q = ‖ X h ‖ 2 / n 、a = ‖ h S ‖ 1 、b = ‖ h S c ‖ 1 。精确解满足
q 2 + λ b ≤ 3 λ a ,而兼容条件给
a ≤ s q / ϕ 3 。因此
q 2 ≤ 3 λ a ≤ 3 λ s ϕ 3 q . 若 q = 0 ,ϕ 3 > 0 又给 h = 0 ;若 q > 0 ,除以 q 得到式(4)的预测界。再用 a + b ≤ 4 a ,得到
‖ h ‖ 1 ≤ 4 s q / ϕ 3 ≤ 12 λ s / ϕ 3 2 。平方在最后的分母出现两次来源:一次用设计控制 a ,一次控制 q 。
若改用 κ 3 > 0 ,同样推理给
‖ h S ‖ 2 ≤ 3 λ s / κ 3 2 ,其中只控制支持内欧氏误差。若用 α 3 > 0 ,则直接以
a ≤ s ‖ h ‖ 2 得
(6) ‖ h ‖ 2 ≤ 3 λ s α 3 2 . 这也解释了为什么必须先声明使用哪一种常数。
例子与边界
两列相关性可以完整算出
取
Γ = ( 1 r r 1 ) , | r | ≤ 1 , S = { 1 } , c ≥ 1. 锥内非零向量可写成 h = h 1 ( 1 , t ) ,| t | ≤ c 。支持分母等于 h 1 2 ,所以
κ 2 = ϕ 2 = min | t | ≤ c ( 1 + 2 r t + t 2 ) = 1 − r 2 . 最小值在 t = − r 取得。全向量分母还要除以 1 + t 2 ;最小特征方向 t = − sign ( r ) 属于锥,因此
α 2 = 1 − | r | . 当 r = 1 / 2 时,κ 2 = ϕ 2 = 3 / 4 ,α 2 = 1 / 2 。取无噪声 β 0 = ( 1 , 0 ) 、λ = 1 / 4 ,唯一 Lasso 解为 ( 3 / 4 , 0 ) ;实际 q 2 = 1 / 16 、‖ h ‖ 1 = ‖ h ‖ 2 = 1 / 4 。式(4)给 q 2 ≤ 3 / 4 、‖ h ‖ 1 ≤ 4 ,式(6)给 ‖ h ‖ 2 ≤ 3 / 2 。界很保守,但方向和定标正确。
当 r = 1 ,h = ( 1 , − 1 ) 既在锥内又满足 X h = 0 ,三个常数全为零。此时两种一稀疏真参数 ( 1 , 0 ) 和 ( 0 , 1 ) 产生相同响应分布。不是把优化再做精确一些就能恢复哪一列真正非零。
全局奇异不一定破坏指定支持
考虑三个单位尺度列,其 Gram 矩阵为
Γ = ( 1 0 0 0 1 1 0 1 1 ) , S = { 1 } . 第二、第三列重复,故 λ min ( Γ ) = 0 。但任意 h 都满足
h T Γ h = h 1 2 + ( h 2 + h 3 ) 2 ≥ h 1 2 ,取 h = ( 1 , 0 , 0 ) 又达到等号,所以对任意有限 c 都有 κ 2 = ϕ 2 = 1 。第一列的稀疏信号仍可得到式(4)的控制。
这并非“该矩阵对所有一稀疏信号都好”。若改成 S = { 2 } 、c ≥ 1 ,向量 ( 0 , 1 , − 1 ) 就在锥中,受限常数变成零。固定支持保证与统一稀疏保证的量词在这里有实际区别。
可识别性、等距性与支持选择的边界
若对每个至多 s 0 的支持都有 κ ( S , c ) > 0 ,且 c ≥ 1 ,就不可能有非零的至多 2 s 0 稀疏向量 v 满足 X v = 0 。证明是把其支持分成大小至多 s 0 的两部分,选 ℓ 1 质量较大的一半作为 S ,于是 v ∈ C ( S , 1 ) ⊆ C ( S , c ) ,与正的 κ 矛盾。因此两个至多 s 0 稀疏参数不能产生相同均值。
受限等距性质则要求对每个至多 k 稀疏的 v 有
( 1 − d k ) ‖ v ‖ 2 2 ≤ ‖ X v ‖ 2 2 / n ≤ ( 1 + d k ) ‖ v ‖ 2 2 。它同时要求上下界,而且检验的是严格稀疏向量;本页的锥允许许多小的非零坐标。适当阶数和常数的等距条件可以推出 Lasso 所需的受限界,但不能仅凭两个名称都含“受限”就视为同一条件。
即便 Γ 全局正定,也未必让原始 Lasso 选中真支持。支持恢复 理路 Lasso 支持恢复与不可表示条件 Lasso support recovery · Lasso sign consistency · Lasso irrepresentable condition 用活动集上的显式候选和非活动集的严格对偶余量证明符号恢复,区分可识别性、参数小误差、变量筛选和坐标推断。 给出一个最小特征值为 1 − 3 2 / 5 > 0 的三列设计,Lasso 却持续选择额外变量。受限界说明误差可以变小,支持选择还要防止零坐标上残留任何非零值。
推论与应用
近似解:把松弛量与锥宽一起处理
假定真支持非空,噪声事件与基本不等式成立,候选的归一化目标容差为 δ ≥ 0 。现在有
(7) q 2 + λ b ≤ 3 λ a + 2 δ . 令 ϕ 4 = ϕ ( S , 4 ) > 0 。分两种情况,而不是直接宣称 h ∈ C ( S , 3 ) :
若 a ≥ 2 δ / λ ,则 b ≤ 3 a + 2 δ / λ ≤ 4 a ,所以 h ∈ C ( S , 4 ) 。又 q 2 ≤ 4 λ a ,兼容条件给 q ≤ 4 λ s / ϕ 4 ,继而 a + b ≤ 5 a ≤ 20 λ s / ϕ 4 2
若 a < 2 δ / λ ,式(7)直接给 q 2 < 8 δ 、b < 8 δ / λ 、a + b < 10 δ / λ ,无需再调用受限条件
两种情况合并得到
(8) q 2 ≤ max { 16 λ 2 s ϕ 4 2 , 8 δ } , ‖ h ‖ 1 ≤ max { 20 λ s ϕ 4 2 , 10 δ λ } . 第一种情况下如果 q = 0 ,正的兼容常数先给 h = 0 ,无须除以零。若 S = ∅ ,不用定义式(2)的常数;基本不等式直接给
q 2 + λ ‖ h ‖ 1 ≤ 2 δ 。
例如正交设计 Γ = I 2 、S = { 1 } 、λ = 0.1 、δ = 0.001 ,有 ϕ 4 2 = 1 。式(8)给 q 2 ≤ 0.16 、‖ h ‖ 1 ≤ 2 ;较小的优化项分别为 0.008 和 0.1 。这不是预测误差必然为 0.16 ,而是此组条件可认证的上界。随着 δ 降低,统计项不会跟着消失。
怎样接到坐标推断
去偏 Lasso 理路 去偏 Lasso 与坐标置信区间 Debiased Lasso · de-biased Lasso · de-sparsified Lasso · desparsified Lasso · 去稀疏 Lasso 用近似逆矩阵修正 Lasso 的残差得分,将坐标误差拆成高斯主项与可控制的乘积余项,并据此构造置信区间。 需要初始总误差 L = ‖ β ~ − β 0 ‖ 1 的控制。本页式(8)提供一种有条件的上界;把它乘上该页可计算的逆矩阵行缺陷 δ j ,再与 σ v j / n 比较,就能判断余项是否足够小。这里优化容差 δ 与逆行缺陷 δ j 是不同量。
若列尺度有界、次高斯噪声使 λ 取 log p / n 量级,并且相关兼容常数有统一正下界,精确解的预测界是 s log p / n 量级,ℓ 1 界是 s log p / n 量级。式(8)说明近似解还须把 δ 控制在相容尺度,不能只报告迭代次数。验证统一设计条件本身可能很困难;这些公式并不宣称已提供高维矩阵上廉价的全支持认证算法。
无噪声欠定解码需要再指定恢复规则
本页的稀疏可识别性排除两个短支持参数产生同一观测,却尚未说明一个可计算的规则会返回哪一个可行点。基追踪的稀疏恢复证书 理路 基追踪的稀疏恢复证书 Basis pursuit recovery certificate · RIP certificate for sparse recovery · 基追踪与受限等距恢复 从无噪声欠定观测建立 ℓ1 解码器,用完整分块证明和可精算的单纯形矩阵分开认证优化最优性与统一稀疏恢复。 给出 min { ‖ z ‖ 1 : A z = y } 、常数一误差锥和受限奇异值的完整分块证明;其七乘八单纯形矩阵可以对所有一稀疏输入统一认证。另一份二乘三矩阵的所有二列都独立,甚至 δ 2 < 1 ,基追踪却偏爱目标 2 / 3 的错误向量。这里的锥 RE、严格稀疏下界和指定解码器的成功因此仍须分别核对;有噪声 Lasso 的常数三锥及本页全部误差界不因此改变。
参考资料