Skip to content

定义Definition

受限特征值与稀疏可识别性

Restricted eigenvalue condition · Compatibility condition for Lasso · Lasso 兼容条件

在稀疏误差锥上定量比较预测与参数范数,精确定标兼容常数和受限特征值,并处理非零优化容差。

形式陈述 ​

设计矩阵把参数误差 h 送到预测误差 Xh。若不同参数产生几乎相同的预测,单凭拟合就难以辨认系数。受限设计条件只对统计证明实际会遇到的误差方向要求定量可辨,而不要求整个高维参数空间都可逆。

给定实矩阵 X∈Rn×p,其中 n,p≥1,记 Γ=XTX/n。固定 ∅≠S⊆{1,…,p}、s=|S| 和 c≥0,定义

(1)C(S,c)={h∈Rp:‖hSc‖1≤c‖hS‖1}.

这里使用坐标的 ℓ1 与欧氏范数;hS 只保留 S 内坐标。集合是齐次的:把 h 放大不会改变其成员资格。非零成员必有 hS≠0,因此下面分母都非零。

本页把三个常数全部写成平方量,避免把分母或平方的位置混在一起:

(2)κ2(S,c)=infh∈C(S,c)∖{0}hTΓh‖hS‖22,ϕ2(S,c)=infh∈C(S,c)∖{0}shTΓh‖hS‖12,(3)α2(S,c)=infh∈C(S,c)∖{0}hTΓh‖h‖22.

κ 是以支持内欧氏误差为分母的受限特征值常数;ϕ 是兼容常数;α 在本页专指全向量欧氏误差版本。文献可能把其中不同量都简称 RE,使用一个定理时应先核分母,而不是只认名称。记 κ,ϕ,α 为非负平方根。

给定整数 1≤s0≤p,若希望对所有至多 s0 稀疏的真参数统一保证,还须对全部 1≤|S|≤s0 取最小值。本页先给固定 S 的结果;实际不知道真支持,并不会使一个特定 S 的计算自动变成全支持保证。

精确 Lasso 解在噪声事件 ‖XTε/n‖∞≤λ/2 上,由基本不等式落入 C(S,3)。若 ϕ3=ϕ(S,3)>0,则

(4)‖X(β^−β0)‖22n≤9λ2sϕ32,‖β^−β0‖1≤12λsϕ32.

这些结论要求 βSc0=0,目标是半平均平方损失加 λ‖β‖1。若只知道目标差至多 δ>0,不能原样使用式(4);下文给出可直接接纳 δ 的版本。

直觉

无约束的最小特征值遍历所有方向。高维时 p>n,总有某个非零向量落在 ker⁡X,所以全空间最小特征值是零。稀疏分析并不要求每条方向都能被观察到,而是先用惩罚比较证明误差不能主要堆在真支持之外,再在剩余的锥上检查预测能否控制参数。

这套顺序有两项独立责任。基本不等式负责证明“误差属于哪里”,设计条件负责证明“在这片区域内,小预测意味着多小的参数误差”。若把近似解当作精确解,第一项责任已经失效,第二项再漂亮也不能补上。

常数之间能比较什么 ​

对任意 h,Cauchy–Schwarz 不等式给 ‖hS‖1≤s‖hS‖2。同时 ‖hS‖2≤‖h‖2,所以

(5)α2(S,c)≤κ2(S,c)≤ϕ2(S,c).

若全局最小特征值为 λmin(Γ)>0,则 α2≥λmin(Γ)。反过来不成立:受限锥可以完全避开全局零空间。

对固定有限 s,c,锥上还有 ‖h‖22≤(1+c2s)‖hS‖22 和 ‖hS‖12≥‖hS‖22,因此

κ2≤(1+c2s)α2,ϕ2≤sκ2.

所以本页三个固定支持常数是否严格为正是等价的,但其数值以及随 s,n 变化的统一下界并不相同。拿一个 κ 界直接当作同数值的全 ℓ2 界,会丢掉这些因子。文献中更大支持集版本的 RE 还改变了定义域,不能只靠这里的比较替换。

从锥上的不等式走到误差界 ​

令 q=‖Xh‖2/n、a=‖hS‖1、b=‖hSc‖1。精确解满足 q2+λb≤3λa,而兼容条件给 a≤sq/ϕ3。因此

q2≤3λa≤3λsϕ3q.

若 q=0,ϕ3>0 又给 h=0;若 q>0,除以 q 得到式(4)的预测界。再用 a+b≤4a,得到 ‖h‖1≤4sq/ϕ3≤12λs/ϕ32。平方在最后的分母出现两次来源:一次用设计控制 a,一次控制 q。

若改用 κ3>0,同样推理给 ‖hS‖2≤3λs/κ32,其中只控制支持内欧氏误差。若用 α3>0,则直接以 a≤s‖h‖2 得

(6)‖h‖2≤3λsα32.

这也解释了为什么必须先声明使用哪一种常数。

例子与边界

两列相关性可以完整算出 ​

取

Γ=(1rr1),|r|≤1,S={1},c≥1.

锥内非零向量可写成 h=h1(1,t),|t|≤c。支持分母等于 h12,所以

κ2=ϕ2=min|t|≤c(1+2rt+t2)=1−r2.

最小值在 t=−r 取得。全向量分母还要除以 1+t2;最小特征方向 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);实际 q2=1/16、‖h‖1=‖h‖2=1/4。式(4)给 q2≤3/4、‖h‖1≤4,式(6)给 ‖h‖2≤3/2。界很保守,但方向和定标正确。

当 r=1,h=(1,−1) 既在锥内又满足 Xh=0,三个常数全为零。此时两种一稀疏真参数 (1,0) 和 (0,1) 产生相同响应分布。不是把优化再做精确一些就能恢复哪一列真正非零。

全局奇异不一定破坏指定支持 ​

考虑三个单位尺度列,其 Gram 矩阵为

Γ=(100011011),S={1}.

第二、第三列重复,故 λmin(Γ)=0。但任意 h 都满足 hTΓh=h12+(h2+h3)2≥h12,取 h=(1,0,0) 又达到等号,所以对任意有限 c 都有 κ2=ϕ2=1。第一列的稀疏信号仍可得到式(4)的控制。

这并非“该矩阵对所有一稀疏信号都好”。若改成 S={2}、c≥1,向量 (0,1,−1) 就在锥中,受限常数变成零。固定支持保证与统一稀疏保证的量词在这里有实际区别。

可识别性、等距性与支持选择的边界 ​

若对每个至多 s0 的支持都有 κ(S,c)>0,且 c≥1,就不可能有非零的至多 2s0 稀疏向量 v 满足 Xv=0。证明是把其支持分成大小至多 s0 的两部分,选 ℓ1 质量较大的一半作为 S,于是 v∈C(S,1)⊆C(S,c),与正的 κ 矛盾。因此两个至多 s0 稀疏参数不能产生相同均值。

受限等距性质则要求对每个至多 k 稀疏的 v 有 (1−dk)‖v‖22≤‖Xv‖22/n≤(1+dk)‖v‖22。它同时要求上下界,而且检验的是严格稀疏向量;本页的锥允许许多小的非零坐标。适当阶数和常数的等距条件可以推出 Lasso 所需的受限界,但不能仅凭两个名称都含“受限”就视为同一条件。

即便 Γ 全局正定,也未必让原始 Lasso 选中真支持。支持恢复给出一个最小特征值为 1−32/5>0 的三列设计,Lasso 却持续选择额外变量。受限界说明误差可以变小,支持选择还要防止零坐标上残留任何非零值。

推论与应用

近似解:把松弛量与锥宽一起处理 ​

假定真支持非空,噪声事件与基本不等式成立,候选的归一化目标容差为 δ≥0。现在有

(7)q2+λb≤3λa+2δ.

令 ϕ4=ϕ(S,4)>0。分两种情况,而不是直接宣称 h∈C(S,3):

  • 若 a≥2δ/λ,则 b≤3a+2δ/λ≤4a,所以 h∈C(S,4)。又 q2≤4λa,兼容条件给 q≤4λs/ϕ4,继而 a+b≤5a≤20λs/ϕ42
  • 若 a<2δ/λ,式(7)直接给 q2<8δ、b<8δ/λ、a+b<10δ/λ,无需再调用受限条件

两种情况合并得到

(8)q2≤max{16λ2sϕ42,8δ},‖h‖1≤max{20λsϕ42,10δλ}.

第一种情况下如果 q=0,正的兼容常数先给 h=0,无须除以零。若 S=∅,不用定义式(2)的常数;基本不等式直接给 q2+λ‖h‖1≤2δ。

例如正交设计 Γ=I2、S={1}、λ=0.1、δ=0.001,有 ϕ42=1。式(8)给 q2≤0.16、‖h‖1≤2;较小的优化项分别为 0.008 和 0.1。这不是预测误差必然为 0.16,而是此组条件可认证的上界。随着 δ 降低,统计项不会跟着消失。

怎样接到坐标推断 ​

去偏 Lasso需要初始总误差 L=‖β~−β0‖1 的控制。本页式(8)提供一种有条件的上界;把它乘上该页可计算的逆矩阵行缺陷 δj,再与 σvj/n 比较,就能判断余项是否足够小。这里优化容差 δ 与逆行缺陷 δj 是不同量。

若列尺度有界、次高斯噪声使 λ 取 log⁡p/n 量级,并且相关兼容常数有统一正下界,精确解的预测界是 slog⁡p/n 量级,ℓ1 界是 slog⁡p/n 量级。式(8)说明近似解还须把 δ 控制在相容尺度,不能只报告迭代次数。验证统一设计条件本身可能很困难;这些公式并不宣称已提供高维矩阵上廉价的全支持认证算法。

无噪声欠定解码需要再指定恢复规则 ​

本页的稀疏可识别性排除两个短支持参数产生同一观测,却尚未说明一个可计算的规则会返回哪一个可行点。基追踪的稀疏恢复证书给出 min{‖z‖1:Az=y}、常数一误差锥和受限奇异值的完整分块证明;其七乘八单纯形矩阵可以对所有一稀疏输入统一认证。另一份二乘三矩阵的所有二列都独立,甚至 δ2<1,基追踪却偏爱目标 2/3 的错误向量。这里的锥 RE、严格稀疏下界和指定解码器的成功因此仍须分别核对;有噪声 Lasso 的常数三锥及本页全部误差界不因此改变。

参考资料
  • Peter J. Bickel, Ya'acov Ritov and Alexandre B. Tsybakov, Simultaneous Analysis of Lasso and Dantzig Selector, Annals of Statistics 37(4), 2009,§3 的 RE 定义及 §7 的误差结论。原文主要常数以支持内欧氏范数为分母;本文分别标明三种 convention,并重新推导所用常数。
  • 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 的常数定义、§2.1 的误差推导、§2.2 的条件关系以及 §11 的噪声情形。式(8)是本文从带容差基本不等式作二分得到的版本,不是把原文精确解定理直接换一个符号。
关系图谱16 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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