Skip to content

定理Theorem

Lasso 支持恢复与不可表示条件

Lasso support recovery · Lasso sign consistency · Lasso irrepresentable condition

用活动集上的显式候选和非活动集的严格对偶余量证明符号恢复,区分可识别性、参数小误差、变量筛选和坐标推断。

形式陈述 ​

支持恢复问的是:哪些真实系数恰好非零。它比均值预测准确更苛刻,也不由参数误差趋零自动推出;一个本来为零的系数即使估成 10−8,在未阈值化的支持里仍算被选中。本页给出一条可以逐项核查的有限样本充分条件,并完整说明其失效方式。

在固定设计的线性模型 Y=Xβ0+ε 中,记

Γ=XTX/n,W=XTε/n,S={j:βj0≠0},z=sign(βS0).

这里 X∈Rn×p、n,p≥1,假定 1≤|S|<p,ΓSS 正定,λ>0,并采用半平均平方损失加 λ‖β‖1。下标 SS 表示取对应行列的子矩阵。空支持可用零解 KKT 单独检查;全支持没有下面的非活动集条件。

先构造一个只在真支持上非零的候选

(1)bS=βS0+ΓSS−1(WS−λz),bSc=0.

这不是可实施的未知支持搜索算法,而是证明候选真解的工具。若同时满足

(2)sign(bS)=z,(3)‖WSc−ΓScSΓSS−1WS+λΓScSΓSS−1z‖∞<λ,

则 b 是唯一 Lasso 最优解,且其支持和符号与 β0 完全相同。式(3)取严格不等式,是本页保证所有解都排除非活动列并获得唯一性的充分余量;把它换成非严格式仍能配合式(2)验证某个最优解,却不能原样引用下面的唯一性证明。

一种更容易解释的充分条件把式(3)拆开。若存在 0<η≤1 使

(4)‖ΓScSΓSS−1z‖∞≤1−η,(5)‖WSc−ΓScSΓSS−1WS‖∞<ηλ,

并且最小非零信号满足

(6)βmin:=minj∈S|βj0|>‖ΓSS−1(WS−λz)‖∞,

那么式(2)、(3)成立。式(4)称为针对真实符号的不可表示条件,式(5)控制活动列无法解释的噪声相关性,式(6)防止信号被收缩或噪声翻过零点。这三项控制不同的失败来源。

直觉

先假定支持与符号都猜对了。活动坐标上的 ℓ1 惩罚便有固定斜率 λz,所以可以解一个普通线性系统得到式(1)。但它只让活动列达到平衡;还必须检查剩余列有没有超过阈值、要求进入模型的相关性。这正是式(3)的任务。

不可表示条件检查一种“借道进入”的机制。活动系数为了支付惩罚而收缩,留下沿活动列的残差。如果某个非活动列能与这些残差过强地对齐,它就会被 Lasso 选进来,即使没有任何观测噪声。条件依赖真实符号 z,而不只是两两相关性的绝对值。

最小信号条件的作用不同。即便所有列彼此正交,一个幅度小于 λ 的真实系数也会被软阈值压成零。良好设计不能替代足够的信号强度,足够强的信号也不能替代非活动列的对偶余量。

证明:从候选到唯一真解 ​

Lasso 的 KKT 条件在本页定标下是 XT(Y−Xb)/n=λs,其中非零坐标 sj=sign(bj),零坐标允许 |sj|≤1。活动部分由式(1)直接得到

WS+ΓSS(βS0−bS)=λz.

非活动部分则恰好是式(3)左端向量。因此式(2)、(3)让候选满足全部 KKT,得到全局最优性。

为了证明唯一性,先说明所有最优解都有相同的拟合值。若两个最优解的 Xb 不同,平方损失在这两个拟合值的中点严格小于端点平均,而 ℓ1 惩罚在中点至多等于端点平均,于是中点目标严格更小,矛盾。故所有最优解的残差与残差得分都相同。

式(3)使每个非活动得分严格小于 λ。任何最优解若在这个坐标非零,KKT 就要求得分绝对值恰等于 λ,矛盾。因此所有最优解都只用 S。最后 ΓSS 正定意味着 XS 满列秩;相同拟合值只能来自相同的活动系数,唯一性得证。

例子与边界

一个相关设计,三项条件都能核算 ​

取

Γ=(11/21/21),β0=(1,0)T,λ=1/4,W=(1/20,1/10)T.

这里 S={1}、z=1。活动候选是 b1=1+1/20−1/4=4/5,仍为正。不可表示量为 1/2,可取 η=1/2;投影后噪声是

W2−12W1=110−140=340<18=ηλ.

非活动得分等于 3/40+1/8=1/5<1/4。因此唯一解为 (4/5,0),支持恢复正确。它仍有系数收缩误差 (−1/5,0),说明支持正确不等于参数已经等于真值。该 W 可由满列秩设计上的某个噪声向量实现;这项确定性检查本身未声称该噪声事件具有任何特定概率。

若同一设计无噪声,则唯一解变为 (3/4,0)。这与受限设计页的预测、参数计算使用同一组定标;是否恢复支持是额外核查出来的事实。

全局正定、预测趋准,支持仍可一直错误 ​

取三个列的 Gram 矩阵

Γ=(103/5013/53/53/51),β0=(1,1,0)T,W=0.

其特征值为 1,1+32/5,1−32/5,全部为正,所以全空间以及受限锥上的参数方向都可辨认。可是对真支持 S={1,2},不可表示量是 (3/5,3/5)(1,1)T=6/5>1。

在 λ=1/10 时,只用真支持的候选为 (9/10,9/10,0),第三列的残差得分是 3/25=0.12>λ,违反 KKT。实际唯一解为

β^=(67,67,114)T.

验证很短:有 Γ−11=(10/7,10/7,−5/7)T, 所以 β^=β0−λΓ−11 的三个坐标全为正,且残差得分恰为 λ1。正定性给出唯一性。其均值预测误差为

q2=(β^−β0)TΓ(β^−β0)=λ21TΓ−11=3140.

更一般地,对每个 0<λ<7/10,同一形式都让第三坐标等于 5λ/7>0。当 λ↓0,预测和参数误差都趋零,但原始 Lasso 支持始终是三个变量。这给出了“小误差不推出恰好为零”的具体序列。

正交设计也会漏掉弱信号 ​

若 Γ=Ip 且无噪声,Lasso 是逐坐标软阈值: β^j=sign(βj0)(|βj0|−λ)+。 此时不可表示量为零,非活动得分也为零;但若某个非零系数满足 |βj0|≤λ,它仍被删掉。边界等号也返回零,所以保证保留真支持需要严格的信号余量。

重复列则有更根本的障碍。若 X1=X2,真参数 (1,0) 与 (0,1) 的全部观测分布相同,任何算法都不能在两个模型上同时以趋一概率恢复各自不同的支持。这个障碍来自实验设计中的不可识别性,不是 Lasso 独有的缺点。

推论与应用

从确定性余量到概率保证 ​

式(1)–(3)是固定噪声向量下的充分条件。要写“以至少 1−a 概率恢复”,还须在已声明的噪声模型下证明式(5)、(6)共同成立的概率;不能只看到 ‖W‖∞≤λ/2,就认为投影后噪声也满足同一个界。矩阵 ΓSS−1 和 ΓScSΓSS−1 可能放大原始噪声。

例如用矩阵行绝对值和范数 ‖B‖∞→∞=maxi∑j|Bij|,可得

‖ΓSS−1(WS−λz)‖∞≤‖ΓSS−1‖∞→∞‖WS‖∞+λ‖ΓSS−1z‖∞.

把一个已证明的噪声上界代入右端,才得到可使用的 beta-min 要求。该上界可能保守,但暴露了活动设计病态会把噪声和收缩同时放大的原因。

优化近似与阈值化支持 ​

小优化 gap 不保证候选的零坐标恰好为零。即便精确解唯一且非活动 KKT 余量严格,给某个非活动坐标加入任意趋零的小数,也得到目标差趋零但支持错误的候选。因此本页精确 KKT 定理不能把 β^ 不加说明地换成数值迭代中任何近似向量。

有误差界时,可以明确改变报告规则。若某估计量满足 ‖β~−β0‖∞≤r,且 βmin>2r,定义

S^r={j:|β~j|>r}.

真零坐标的估计绝对值至多 r,不会被保留;真非零坐标的估计绝对值至少 βmin−r>r,所以 S^r=S。一个 ℓ1 上界也能充当保守的 r,因为 ‖h‖∞≤‖h‖1。这里保证的是明确阈值化后的支持,不能把它写成原始 Lasso 自动支持恢复。

支持恢复与坐标推断是两个终点 ​

去偏 Lasso给预先指定坐标构造置信区间,其修正后系数通常不再稀疏;它的余项条件可以在未恢复完整支持时成立。相反,即使这次选对了支持,也没有自动得到无偏系数或校准的标准误。若用同一资料先筛选再报告区间,还要考虑选择后推断的目标与条件分布。

实际报告可以分别回答三件事:原设计的均值预测误差多大,系数估计误差有什么条件界,非零名单是否有经过验证的恢复保证。每件事的证据应对应其自己的定理。

参考资料
  • Peng Zhao and Bin Yu, On Model Selection Consistency of Lasso, JMLR 7, 2006, pp. 2541–2563,§2、Proposition 1 与 Appendix A 的 KKT 证明。原文目标为未归一化平方和加 λold‖β‖1,对应本文 λold=2nλ;本文将其机制写成固定噪声下的严格余量证书,并单独证明唯一性。
  • 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.1、6–8:预测条件与变量选择条件的区别。本文三列反例的解与误差均在正文直接核算。
关系图谱18 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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