形式陈述
支持恢复问的是:哪些真实系数恰好非零。它比均值预测准确更苛刻,也不由参数误差趋零自动推出;一个本来为零的系数即使估成 10 − 8 ,在未阈值化的支持里仍算被选中。本页给出一条可以逐项核查的有限样本充分条件,并完整说明其失效方式。
在固定设计的线性模型 理路 线性回归统计模型 Linear regression model · Linear model 响应的条件均值由设计变量对未知系数线性表示,并显式规定误差结构的统计模型。 Y = X β 0 + ε 中,记
Γ = X T X / n , W = X T ε / n , S = { j : β j 0 ≠ 0 } , z = sign ( β S 0 ) . 这里 X ∈ R n × p 、n , p ≥ 1 ,假定 1 ≤ | S | < p ,Γ S S 正定,λ > 0 ,并采用半平均平方损失加 λ ‖ β ‖ 1 。下标 S S 表示取对应行列的子矩阵。空支持可用零解 KKT 单独检查;全支持没有下面的非活动集条件。
先构造一个只在真支持上非零的候选
(1) b S = β S 0 + Γ S S − 1 ( W S − λ z ) , b S c = 0. 这不是可实施的未知支持搜索算法,而是证明候选真解的工具。若同时满足
(2) sign ( b S ) = z , (3) ‖ W S c − Γ S c S Γ S S − 1 W S + λ Γ S c S Γ S S − 1 z ‖ ∞ < λ , 则 b 是唯一 Lasso 最优解,且其支持和符号与 β 0 完全相同。式(3)取严格不等式,是本页保证所有解都排除非活动列并获得唯一性的充分余量;把它换成非严格式仍能配合式(2)验证某个最优解,却不能原样引用下面的唯一性证明。
一种更容易解释的充分条件把式(3)拆开。若存在 0 < η ≤ 1 使
(4) ‖ Γ S c S Γ S S − 1 z ‖ ∞ ≤ 1 − η , (5) ‖ W S c − Γ S c S Γ S S − 1 W S ‖ ∞ < η λ , 并且最小非零信号满足
(6) β min := min j ∈ S | β j 0 | > ‖ Γ S S − 1 ( W S − λ z ) ‖ ∞ , 那么式(2)、(3)成立。式(4)称为针对真实符号的不可表示条件 ,式(5)控制活动列无法解释的噪声相关性,式(6)防止信号被收缩或噪声翻过零点。这三项控制不同的失败来源。
直觉
先假定支持与符号都猜对了。活动坐标上的 ℓ 1 惩罚便有固定斜率 λ z ,所以可以解一个普通线性系统得到式(1)。但它只让活动列达到平衡;还必须检查剩余列有没有超过阈值、要求进入模型的相关性。这正是式(3)的任务。
不可表示条件检查一种“借道进入”的机制。活动系数为了支付惩罚而收缩,留下沿活动列的残差。如果某个非活动列能与这些残差过强地对齐,它就会被 Lasso 选进来,即使没有任何观测噪声。条件依赖真实符号 z ,而不只是两两相关性的绝对值。
最小信号条件的作用不同。即便所有列彼此正交,一个幅度小于 λ 的真实系数也会被软阈值压成零。良好设计不能替代足够的信号强度,足够强的信号也不能替代非活动列的对偶余量。
证明:从候选到唯一真解
Lasso 的 KKT 条件 理路 Lasso 的最优性与对偶间隙 Lasso optimality conditions · Lasso duality gap · Lasso primal-dual certificate 从残差相关性核验 Lasso 的零与非零坐标,并用可行对偶值认证剩余优化误差。 在本页定标下是
X T ( Y − X b ) / n = λ s ,其中非零坐标 s j = sign ( b j ) ,零坐标允许 | s j | ≤ 1 。活动部分由式(1)直接得到
W S + Γ S S ( β S 0 − b S ) = λ z . 非活动部分则恰好是式(3)左端向量。因此式(2)、(3)让候选满足全部 KKT,得到全局最优性。
为了证明唯一性,先说明所有最优解都有相同的拟合值。若两个最优解的 X b 不同,平方损失在这两个拟合值的中点严格小于端点平均,而 ℓ 1 惩罚在中点至多等于端点平均,于是中点目标严格更小,矛盾。故所有最优解的残差与残差得分都相同。
式(3)使每个非活动得分严格小于 λ 。任何最优解若在这个坐标非零,KKT 就要求得分绝对值恰等于 λ ,矛盾。因此所有最优解都只用 S 。最后 Γ S S 正定意味着 X S 满列秩;相同拟合值只能来自相同的活动系数,唯一性得证。
例子与边界
一个相关设计,三项条件都能核算
取
Γ = ( 1 1 / 2 1 / 2 1 ) , β 0 = ( 1 , 0 ) T , λ = 1 / 4 , W = ( 1 / 20 , 1 / 10 ) T . 这里 S = { 1 } 、z = 1 。活动候选是
b 1 = 1 + 1 / 20 − 1 / 4 = 4 / 5 ,仍为正。不可表示量为 1 / 2 ,可取 η = 1 / 2 ;投影后噪声是
W 2 − 1 2 W 1 = 1 10 − 1 40 = 3 40 < 1 8 = η λ . 非活动得分等于 3 / 40 + 1 / 8 = 1 / 5 < 1 / 4 。因此唯一解为 ( 4 / 5 , 0 ) ,支持恢复正确。它仍有系数收缩误差 ( − 1 / 5 , 0 ) ,说明支持正确不等于参数已经等于真值。该 W 可由满列秩设计上的某个噪声向量实现;这项确定性检查本身未声称该噪声事件具有任何特定概率。
若同一设计无噪声,则唯一解变为 ( 3 / 4 , 0 ) 。这与受限设计页的预测、参数计算使用同一组定标;是否恢复支持是额外核查出来的事实。
全局正定、预测趋准,支持仍可一直错误
取三个列的 Gram 矩阵
Γ = ( 1 0 3 / 5 0 1 3 / 5 3 / 5 3 / 5 1 ) , β 0 = ( 1 , 1 , 0 ) T , W = 0. 其特征值为 1 , 1 + 3 2 / 5 , 1 − 3 2 / 5 ,全部为正,所以全空间以及受限锥 理路 受限特征值与稀疏可识别性 Restricted eigenvalue condition · Compatibility condition for Lasso · Lasso 兼容条件 在稀疏误差锥上定量比较预测与参数范数,精确定标兼容常数和受限特征值,并处理非零优化容差。 上的参数方向都可辨认。可是对真支持 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。实际唯一解为
β ^ = ( 6 7 , 6 7 , 1 14 ) T . 验证很短:有
Γ − 1 1 = ( 10 / 7 , 10 / 7 , − 5 / 7 ) T ,
所以 β ^ = β 0 − λ Γ − 1 1 的三个坐标全为正,且残差得分恰为 λ 1 。正定性给出唯一性。其均值预测误差为
q 2 = ( β ^ − β 0 ) T Γ ( β ^ − β 0 ) = λ 2 1 T Γ − 1 1 = 3 140 . 更一般地,对每个 0 < λ < 7 / 10 ,同一形式都让第三坐标等于 5 λ / 7 > 0 。当 λ ↓ 0 ,预测和参数误差都趋零,但原始 Lasso 支持始终是三个变量。这给出了“小误差不推出恰好为零”的具体序列。
正交设计也会漏掉弱信号
若 Γ = I p 且无噪声,Lasso 是逐坐标软阈值:
β ^ j = sign ( β j 0 ) ( | β j 0 | − λ ) + 。
此时不可表示量为零,非活动得分也为零;但若某个非零系数满足 | β j 0 | ≤ λ ,它仍被删掉。边界等号也返回零,所以保证保留真支持需要严格的信号余量。
重复列则有更根本的障碍。若 X 1 = X 2 ,真参数 ( 1 , 0 ) 与 ( 0 , 1 ) 的全部观测分布相同,任何算法都不能在两个模型上同时以趋一概率恢复各自不同的支持。这个障碍来自实验设计中的不可识别性,不是 Lasso 独有的缺点。
推论与应用
从确定性余量到概率保证
式(1)–(3)是固定噪声向量下的充分条件。要写“以至少 1 − a 概率恢复”,还须在已声明的噪声模型下证明式(5)、(6)共同成立的概率;不能只看到 ‖ W ‖ ∞ ≤ λ / 2 ,就认为投影后噪声也满足同一个界。矩阵 Γ S S − 1 和 Γ S c S Γ S S − 1 可能放大原始噪声。
例如用矩阵行绝对值和范数 ‖ B ‖ ∞ → ∞ = max i ∑ j | B i j | ,可得
‖ Γ S S − 1 ( W S − λ z ) ‖ ∞ ≤ ‖ Γ S S − 1 ‖ ∞ → ∞ ‖ W S ‖ ∞ + λ ‖ Γ S S − 1 z ‖ ∞ . 把一个已证明的噪声上界代入右端,才得到可使用的 beta-min 要求。该上界可能保守,但暴露了活动设计病态会把噪声和收缩同时放大的原因。
优化近似与阈值化支持
小优化 gap 不保证候选的零坐标恰好为零。即便精确解唯一且非活动 KKT 余量严格,给某个非活动坐标加入任意趋零的小数,也得到目标差趋零但支持错误的候选。因此本页精确 KKT 定理不能把 β ^ 不加说明地换成数值迭代中任何近似向量。
有误差界时,可以明确改变报告规则。若某估计量满足
‖ β ~ − β 0 ‖ ∞ ≤ r ,且 β min > 2 r ,定义
S ^ r = { j : | β ~ j | > r } . 真零坐标的估计绝对值至多 r ,不会被保留;真非零坐标的估计绝对值至少 β min − r > r ,所以 S ^ r = S 。一个 ℓ 1 上界也能充当保守的 r ,因为 ‖ h ‖ ∞ ≤ ‖ h ‖ 1 。这里保证的是明确阈值化后的支持,不能把它写成原始 Lasso 自动支持恢复。
支持恢复与坐标推断是两个终点
去偏 Lasso 理路 去偏 Lasso 与坐标置信区间 Debiased Lasso · de-biased Lasso · de-sparsified Lasso · desparsified Lasso · 去稀疏 Lasso 用近似逆矩阵修正 Lasso 的残差得分,将坐标误差拆成高斯主项与可控制的乘积余项,并据此构造置信区间。 给预先指定坐标构造置信区间,其修正后系数通常不再稀疏;它的余项条件可以在未恢复完整支持时成立。相反,即使这次选对了支持,也没有自动得到无偏系数或校准的标准误。若用同一资料先筛选再报告区间,还要考虑选择后推断 理路 选择后推断与截断正态校准 Post-selection inference · Selective inference · 选择性推断 在同一观测先被筛选再用于推断时,按选择事件重算零分布,并由截断正态枢轴得到条件有效的检验和区间。 的目标与条件分布。
实际报告可以分别回答三件事:原设计的均值预测误差多大,系数估计误差有什么条件界,非零名单是否有经过验证的恢复保证。每件事的证据应对应其自己的定理。
参考资料