Benign overfitting of non-sparse high-dimensional linear regression with correlated noise¶
作者: Toshiki Tsuda, Masaaki Imaizumi
来源: Electronic Journal of Statistics
主题: 高维统计 / 随机矩阵
相关性: 8/10
链接: 期刊页 · arXiv
一、领域脉络与小综述¶
-
这个方向是什么:本方向的核心问题是高维线性回归中无正则化估计量(interpolator / minimum-norm estimator)能否在过参数化(p >> n)且非稀疏的真实参数下实现一致的预测与估计(良性过拟合,benign overfitting)。传统统计智慧认为,过拟合(完美拟合训练噪声)会损害泛化性能,但深度学习的观测经验和后续理论发现表明,在一定条件下,最小范数插值估计量(minimum-norm interpolator)可以达到近乎最优的预测风险。本方向正处于从“独立同分布高斯噪声 + 稀疏或高度结构化协方差”的经典设定,向“相关噪声(内生性)”、“非高斯数据”、“时间序列依赖”等更一般场景扩展边界的阶段。成熟度:核心定理已确立(Bartlett et al. 2020给出了表征条件),但扩展至相关噪声与内生性仍属前沿。
-
发展脉络(history):串接用户提供的关键引用。
- 奠基工作:发现“良性过拟合”现象并给出第一个完整理论刻画。
- Belkin et al. (2019)(引用[7]):经典的“双下降”(double descent)风险曲线被记录,表明过参数化插值估计量的风险在超过插值阈值后下降——但该文属于经验发现与初步理论,未提供精确条件。
- Bartlett et al. (2020)(引用[3]):首个刻画了最小范数最小二乘插值估计量可实现良性过拟合的充分必要条件,核心条件是用数据协方差矩阵的有效秩(effective rank)来表征(两个有效秩不等)。该文被本文引用为基线条件。
-
Hastie et al. (2022)(引用[8]):用随机矩阵理论(RMT)在比例极限 p/n → γ ∈ (0,∞) 下,给出了岭回归与无脊插值(ridgeless)的精确渐近预测风险,揭示了过参数化下风险可能低于最优正则化风险的特性——这是本文所说的“using random matrix theory”的路径。
-
主要进展:将表征条件扩展至更一般的正则化与模型类。
- Tsigler & Bartlett (2020)(引用[5]):将良性过拟合从最小范数插值扩展到岭回归,给出了任意正则化参数下的非渐近泛化界,并指出负正则化也能实现良性过拟合。
- Koehler et al. (2021)(引用[6]):用“高斯宽度”(Gaussian width)统一分析了一类插值器(包括L1范数最小化,即basis pursuit)的泛化误差,确认了Bartlett等人对最小范数插值器的预测——即对高斯数据,近最小范数插值器也可实现良性过拟合。
-
Montanari & Saeed (2022) 与 Han & Shen (2022)(引用[16, 19]):这两篇在非高斯性上做了推动——将CGMT框架对高斯数据的精确渐近推广到更一般的设计矩阵(独立同分布但非高斯条目),被称为“universality”。本文引用它们来指出这是已遇到的挑战。
-
当前Frontier与本文的位置:上述工作对噪声的假设几乎都限定为与协变量独立(外生性)。而内生性(endogeneity:噪声与协变量相关)是经济学、因果推断和流行病学的核心问题,且高维情况下尤为关键。Belloni, Chernozhukov等人(引用[9, 12, 13, 21])走的是一条正则化+正交估计方程的路线(IV/Lasso+去偏),处理的是“高维但稀疏”的简约形式参数,并不适用于非稀疏的情形。Fan & Liao (2014)(引用[15])介绍了针对内生性的GMM方法,但同样依赖稀疏假设。本文的独特位置在于:在不使用正则化、不假设稀疏性的前提下,首次证明高维线性回归在存在内生性(噪声与协变量相关)时,最小范数估计量仍可实现良性过拟合(一致性)——条件仅与噪声-工具变量协方差矩阵的特征值有关。这是对Bartlett等人外生性表征的一个直接推广,走向了经济学与因果推断的典型设定。
-
子线索聚类:从这些文献可识别出3条大的子线索:
- 子线索A(外生性下线性插值器的良性过拟合):Bartlett (2020)、Tsigler (2020)、Koehler (2021)、Hastie (2022)。核心问题是:在何种协方差结构下,无正则化最小二乘插值器实现良性过拟合?主要工具:有效秩条件、随机矩阵理论、高斯宽度。本文直接继承和推广此线索。
- 子线索B(高维内生性与IV回归):Belloni, Chernozhukov系列 (2010–2022)、Fan & Liao (2014)、Gold et al. (2020)、Dikkala et al. (2020)。核心问题是:如何在高维内生性环境中对稀疏参数做一致估计和推断?主要工具:Lasso/去偏/正交估计方程/IV。本文与之形成对比:本文的路线不需要稀疏性假设,不采用正则化。
- 子线索C(CGMT与universality):Thrampoulidis et al. 2016、Han & Shen 2022、Montanari & Saeed 2022。核心问题是:如何将CGMT框架(对高斯设计得到精确渐近)推广到非高斯设计?本文使用的工具是扩展CGMT以应对相关噪声(不仅仅是高斯独立)。
- (值得注意):子线索B引用中几乎没有谈到“非稀疏参数”这种设定,基本默认是稀疏。本文的引入巧妙地避开了竞争。
-
这个方向在追问的核心问题:
- 最小范数插值器在何种协方差结构下实现良性过拟合?(已回答:有效秩条件 — Bartlett 2020)
- 当噪声与协变量相关(内生性)时,上述条件是否仍能保证一致性?条件如何变化?(本文要回答的问题)
- 在比例极限p/n = γ下,内生性会如何影响过参数化下的精确渐近风险?执行RMT是否能得到闭合表达式?(竞争工具:Hastie等人的方法,但本文没有用——保持了非渐近characterization路线)
- Universal性的边界在哪?(即,如果设计不是高斯,是否仍能沿用CGMT导出的结论?本文的任务是扩展CGMT本身以纳入相关噪声)
主流方法与瓶颈:目前处理内生性的高维主流方法均为正则化+去偏(子线索B、被引[9][12][13][15]),依赖稀疏性假设。当参数非稀疏时,这些方法失效;而传统的CGMT方法(子线索C)又要求设计独立于噪声。瓶颈正是在于缺少一个对非稀疏+相关噪声的参数保持一致性的理论保证。
-
⚠️作者的framing(必须明确标注为“这是作者的说法”):本文作者将自身研究叙述为:“When the covariates are high-dimensional, it has been common to assume sparsity on the true parameters and estimate them using regularization, even with the endogeneity. However, when sparsity does not hold, it has not been well understood to control the endogeneity and high dimensionality simultaneously.” — 作者以此将自己入口为“非稀疏情形下解决内生性+高维性的空缺”。作者淡化了以下两条竞争路线:(a) 正则化+正交估计方程路线(Chernozhukov等)虽依赖稀疏,但他们的工作实际上能处理近似稀疏,本文对非稀疏的要求是“完全非稀疏”,而非“近似稀疏”,两者强度不分上下。(b) 子线索C(universality)的作者(Han & Shen 2022; Montanari & Saeed 2022)已在独立同分布但非高斯设计下扩展CGMT,本文直接声称“扩展到非独立噪声”,是CGMT自然推广,但作者并未在intro中与这些universality论文做充分对比。
什么明显该被引/该存在、但没出现在intro里?:Bartlett et al. (2020) 对噪声要求独立;那是否有噪声相关但设计独立的工作?没有直接命中。一个值得查的点是:最小范数插值器在带有相关噪声(如AR(1)噪声)的线性回归中的非一致性表现是否已在哪些非渐近工作中被证明过?
-
张力:未见明显对立引用。intra-subclue rivalries较强:子线索B人群(Belloni/Chernozhukov/Fan)与本文所属的子线索A人群(Bartlett/Hastie/Montanari)之间的紧张明显——前者主张正则化,后者主张无正则化也可能一致。但两者在设定上不是严格可比较的,因此无矛盾。
二、最核心、最简单的例子 / 数学问题¶
第一步:把符号、模型、可观测数据交代清楚
假设如下模型(设定来自原文Section 2):
- 符号与量纲:
- \(n\):样本量;\(p\):特征维度。
- \(X \in \mathbb{R}^{n \times p}\):设计矩阵(已观测)。记其第i行为\(x_i \in \mathbb{R}^p\)。
- \(Y \in \mathbb{R}^n\):响应向量(已观测)。
- \(\beta^* \in \mathbb{R}^p\):未知回归系数向量(目标参数,无需稀疏假设——即\(\|\beta^*\|_0\)可以是p)。
- \(\varepsilon \in \mathbb{R}^n\):随机噪声向量(不可观测)。其分布:\(\varepsilon \mid X\)有某种均值零结构但与X相关(内生性)。
- 模型:\(Y = X\beta^* + \varepsilon\),其中\(\varepsilon\)的分布与\(X\)相关——这是内生性的体现。本文假设\(\varepsilon\)与\(X\)是联合高斯(高斯+相关)。更具体地:假设存在潜在变量\(Z \in \mathbb{R}^{n \times p}\)(instrumental variables的一种解读)未知,可以认为\(X\)和\(\varepsilon\)由潜在向量对\(Z\)和噪声\(\tilde{\varepsilon}\)线性组合而成(见原文“general model”设置)。
- 可观测数据与不可观测量的区分:
- 可观测:\((X, Y)\)。
- 不可观测:\(\beta^*\)(要估计)、\(\varepsilon\)(噪声、内生性源)、潜在工具变量\(Z\)。
- 估计量:\(\hat{\beta} := \arg\min_{\beta \in \mathbb{R}^p} \|Y - X\beta\|_2^2\)(无正则化的最小二乘估计量,无论p>n还是p≤n都为minimum-norm interpolator:当p≥n时\(\hat{\beta} = X^\top (XX^\top)^{-1} Y\))。
- 风险度量:\(\mathbb{E}\left[\|\hat{\beta} - \beta^*\|_2^2\right]\),即参数估计的期望\(L_2\)误差(点预测误差的替代度量)。本文关注的是渐近一致性:\(\|\hat{\beta} - \beta^*\|_2 \xrightarrow{P} 0\)。
第二步:讲最小内核——关闭大多数宽松假设,写出最小例子
用最简特例来整体理解本文的核心定理:
最简特例情景: - 设\(p \gg n\)(过参数化)。 - 设设计矩阵\(X\)的每一行\(x_i\)是均值为0、协方差为\(\Sigma_{n \times n}\)(这是一个\(n \times n\)的协方差矩阵,注意抽象层面:特征之间协方差是\(\Sigma\),样本协方差是\(XX^\top\) ——这里要小心:\(X\)是\(n \times p\),\(X\)的每一行分布是\(N_p(0, \Sigma_X)\)。为简化,我们取\(\Sigma_X\)的谱分解:\(\lambda_1 \ge \lambda_2 \ge \dots \ge \lambda_p\)为特征值。 - 噪声\(\varepsilon \in \mathbb{R}^n\)分布:\(\varepsilon \sim N(0, \Omega)\),其中\(\Omega\)是对角或非对角,且与\(X\)相关。具体来说,通过一个未知矩阵\(B\)使得\(X = Z B + \tilde{X}\),其中\(Z\)是工具变量(潜在变量),\(\varepsilon\)与\(Z\)独立但与\(\tilde{X}\)相关——但我们可以将这一切编码为:存在一个潜在秩\(k\)的结构使得协方差矩阵\(Cov(X,\varepsilon) \neq 0\)。
在最小内核下,核心关于“良性过拟合”的条件(Theorem 1)退化为一个清晰的可解释形式:
假设\(\Omega\)(噪声协方差)和对角特征基与\(\Sigma_X\)共享特征基(即对角化在同一组特征向量中:这是“aligned”情形。原文不假设此,但最小特例中可以这样做以使推理紧凑)。设\(\lambda_1 \ge \dots \ge \lambda_p\)是\(\Sigma_X\)的特征值,\(\omega_1,\dots,\omega_p\)是\(\Omega\)的对应特征值(可能为零)。Theorem 1的条件本质上是:
“small eigenvalue”的集合由有效秩条件定义:比如\(r = \text{rank}_k\)条件给出的线性增长速度(Bartlett et al.的\(\tilde{r}\))。这条条件陈述的是:在信号协方差异常弱的那些方向上(特征值很小),噪声协方差增长比信号协方差慢得多,因此噪声在这些方向上是“有效干净的”。在独立噪声情形,\(\omega_j / \lambda_j = \sigma^2 / \lambda_j\),条件还原为Bartlett等人的表征。
在这个最小内核下,定理证明的核心思路变成了一个直接可验证的命题:
- 在最小二乘解\(\hat{\beta} - \beta^* = X^\top (XX^\top)^{-1} \varepsilon\) 中,偏差项(endogeneity bias)由\(Cov(\varepsilon, X)\)决定。但因为用CGMT转化为了一个对偶随机优化问题,最终的渐近条件可以只依赖于\(\Sigma_X\)和\(\Omega\)的特征值,而不直接涉及\(Cov(X,\varepsilon)\)的未知结构——这是全文最巧的步骤(利用重参数化和CGMT消除内生性协方差的不利影响)。在最小例子里,这个对偶问题退化为一个标量随机优化问题,解唯一且可由矩计算直接得出误差的依概率收敛条件。这样,读者可以直观理解:“只要内积噪声谱相对信号谱增长得足够慢,无正则化估计量就可以实现一致性——即使在有内生性时。”
三、这篇论文做了什么(重心)¶
-
三句话:
- 研究了在高维(p >> n)线性回归中,当噪声与协变量相关(内生性),且真实参数\(\beta^*\)不需要稀疏性假设时,无正则化最小二乘估计量的预测/估计误差是否能随p增长而依概率收敛到零(良性过拟合)。
- 核心方法是将CGMT从经典的独立噪声-设计情形扩展到涉及相关噪声的对偶问题,由此得到一个可以将估计量与\(\varepsilon\)和\(X\)之间的依赖解耦的条件。
- 主要结果:当噪声协方差与信号协方差之间特征值满足特定的条件(称为条件E:\(\sum_{j=s+1}^p \omega_j / \sum_{j=s+1}^p \lambda_j \to 0\),其中\(s\)是与有效秩有关的截断索引)时,\(\|\hat{\beta} - \beta^*\|_2 \to 0\)依概率。理论在几种扩展(放宽特征值条件、弱弱依赖结构)下依然成立,模拟实验验证了理论。
-
关键设定与假设(完整版在第二节最小记号上补充):
- A1(高斯+相关性):\(X\)的每个行向量\(x_i\)独立同分布于\(N_p(0, \Sigma_X)\);噪声\(\varepsilon\)的条件分布是\(N_n(0, \Omega)\),其中\(\Sigma_X\)和\(\Omega\)可以是任意(可能秩缺的)正定矩阵,且彼此可以任意相关。这比Bartlett等人假设\(\varepsilon\)与\(X\)独立且\(\Omega = \sigma^2 I_n\)更宽。
- A2(有效秩条件):沿用Bartlett et al. (2020)的两个有效秩\(R(\Sigma_X)\)和\(\tilde{r}(\Sigma_X)\),控制样本协方差阵的集中性。具体来说,令\(r = \text{tr}(\Sigma_X) / \|\Sigma_X\|_{\text{op}}\),\(\tilde{r} = (\sum_{j=1}^p \lambda_j)^2 / \sum_{j=1}^p \lambda_j^2\)。要求这些量随\(p\)增长,使得矩阵集中不破碎。
- A3(平滑谱衰减):\(\Sigma_X\)的谱衰减速率要满足\(\sum_{j > s} \lambda_j \approx s \lambda_{s+1}\)形式的正则条件,这是为了应用CGMT中的标量化。
- A4(噪声特征值条件,即核心条件E):存在某个截断阈值\(s\)(与上面有效秩相关),使得\(\sum_{j=s+1}^p \omega_j / \sum_{j=s+1}^p \lambda_j \to 0\)。其中\(\{\omega_j\}\)是\(\Omega\)在与\(\Sigma_X\)共享的特征基下的特征值(为正交对角化的情形;对更一般情形,需施加某种对齐分布条件)。这条相比独立噪声情形(\(\omega_j = \sigma^2\))更弱更灵活,允许某些方向上的噪声很大,只要总体上弱于信号。
相比已有文献:噪声假设放宽(相关\(\to\)相关+任意\(\Omega\));稀疏假设取消(\(\beta^*\)任意非零)。变量假设收紧(本文只对高斯X、高斯\(\varepsilon\)有效——本文对universality的讨论仅提示“可能是可推广”,但未证明)。
-
主要结果(理论型,选2个最关键定理):
定理1(核心一致性):在A1–A4和某些正则条件下,有
\[\|\hat{\beta} - \beta^*\|_2^2 = o_p(1)。\]直觉:偏差和方差的组合恰好可写为一个依赖于\(\sum \omega_j / \sum \lambda_j\)的表达式(由CGMT准确给出)。条件E确保了这个比值趋于0,从而方差的dominant部分趋于0,偏差同理。技术点:这里的偏差不完全由\(Cov(\varepsilon, X)\)直接介入,而是被信号/噪声特征值比值消去。定理2(更浅的假设的扩展):当\(\Omega\)的特征值不沿与\(\Sigma_X\)对齐的特征基排列时(非对齐情形),定理1仍然成立,只要引入一个“干扰项”并且它在特征值条件下也可忽略。给出了一个更宽松的一致弱收敛条件。这个定理是技术贡献主体之一,因为在非对齐情形下无法直接用简单的特征值比表达,需要额外的CGMT引导步骤。
定理3(噪声协方差有弱稀疏结构下的z):如果\(\Omega\)的谱可以分解为“大块”+“指数衰减尾”等形式,条件E仍可保持,结论不变。
解决的技术难点:最主要的是将内生性协方差的分量从对偶问题中“挤出”:经典CGMT假设\(\varepsilon\)与\(X\)独立,且二次型与线性项独立;这里的相关噪声破坏了独立性,使得标准方法无法直接对最小二乘的残差形式开具对偶问题。作者做了两步:① 通过对\(\varepsilon\)和\(X\)的联合高斯性质,将模型用潜在的秩\(k\)“instrumental variables”结构重新参数化(如\(X = ZB + \tilde{X}\)),将\(\varepsilon\)与\(\tilde{X}\)做解耦;② 由此生成一个辅助的高斯随机向量,其条件分布满足CGMT所需的独立性假设——但这个辅助量仅是原数据的一个“投影”,需要额外的集中性论证对其与原问题的代价泛函做等价性。
-
证明路线与技术技巧(整体路线与关键跳跃):
整体路线(3-5步逻辑主干): 1. 对偶变换(Problem 1 → 原始形式 → 对偶形式):利用支撑函数,将对\(\beta\)的最小二乘解等价为\(\min_{\beta} \|Y-X\beta\|_2^2\)的对偶问题的解——对偶形式是一个在低维(n维)空间里的优化问题,可以通过添加扰动变量\(\beta_a\)(满足\(X^\top \beta_a = 0\))得到。 2. CGMT重参数化(Lemma 12-15):把对偶问题中的所有高斯随机变量通过CGMT化为一个“标量化”形式——即用少量高斯随机变量(\(h \in \mathbb{R}^n\))与原设计X和噪声ε的复杂联合高斯结构做满意同分布处理。 3. 标量化(scalarization):经过CGMT后,原始\(n\times p\)高维优化问题被一个仅依赖于n个随机标量的单变量随机优化所主导。在这个标量问题中,可以把\(\|\hat{\beta} - \beta^*\|_2\)写成\(\frac{\sum \omega_j \epsilon_j^2}{\sum \lambda_j \eta_j^2}\)的形式(\(\epsilon_j, \eta_j\)是独立的子高斯随机变量),且其渐近行为仅由特征值的和比决定。 4. 解析风险表达式:在这个标量问题中直接计算期望,得到条件E的充分性,并证明收敛率。 5. 反向传递(conservation):证明CGMT在这一过程中是“tight”的——原始问题的风险行为被严格地控制在对偶问题的线性(同order)放大/缩小范围内(Lemma 20-23),因此条件E对原始估计量依然成立。
关键跳跃点: - 最吃功夫的引理是Lemma 23(扩展的CGMT引理):要求处理的是\(X\)和\(\varepsilon\)都是用零均值高斯向量构成的,但其协方差矩阵有块状结构,不是对角形式。经典的CGMT只有噪声与信号独立的前置假设。本文的证明绕过的办法是通过高斯条件分解:将\(\varepsilon\)分解成\(X\)的线性函数+一个独立的高斯残差项,然后令CGMT作用于这个分解后的形式。技术上:证明这个“余项”的范数与原本项仅差\(o_p(1)\),并利用\(X\)的谱集中性保证的。 - 第二关键跳跃:在不对齐特征基的情形(定理2)下,无法直接用特征值的标量和比值。证明者运行了一个二次型的“插值”技巧——用扰动分析将\(\Omega\)近似为一个与\(\Sigma_X\)可交换的矩阵的扰动,然后再用特征值条件。
具体技术技巧点名: - CGMT(核心:用于降低维度,将\(p\)维空间消去,仅剩\(n\)维。 - 高斯集中与有效秩近似:处理\(XX^\top\)的谱集中性(与Koltchinskii & Lounici 2017类似)。 - Hermite插值 + 对称矩阵扰动理论:在不对齐情形下使用(定理2)。 - 潜在工具变量结构重参数化:用来解除内生性相关带。
-
真实例子与应用: 论文包含了模拟实验。具体来说:
- 数据生成:固定\(n=200, p=500\),结构为\(X\sim N(0, \Sigma_X)\),\(\varepsilon\)为由一个潜在秩\(k\)的潜变量\(Z\)和\(\tilde{X}\)组合产生的噪声\(\varepsilon = Z \alpha + \tilde{\varepsilon}\),保证\(Corr(\varepsilon, X) \neq 0\)。
- \(\Sigma_X\)设置:让特征值按\(\lambda_j = j^{-a}\)递减,\(\Omega\)的特征值按\(\omega_j = j^{-b}\)递减。改变参数\(a,b\)组合:其比值\(\sum_{j>s} \omega_j / \sum_{j>s} \lambda_j\)当\(b\)比\(a\)足够大时(即噪声衰减快于信号)趋于0,这对应条件E满足。
- 结果:当条件E满足时(\(b\) > \(a\)),\(\|\hat{\beta} - \beta^*\|_2\)随\(p\)增大而下降;当条件E违反时(\(b \le a\)),误差不降反升(发散)。结果与定理1预测一致。
- 扩展实验:改变非对齐程度(让\(\Sigma_X\)和\(\Omega\)非同一基,使用随机旋转),发现只要有效谱相对比值符合条件E,仍保持一致性,验证定理2。
- 这个例子想说明什么:验证内生性下条件E的充分性,并展示对齐/非对齐情形均适用。
-
🔎结论是否比证明窄: 值得注意提及:Theorem 1的证明严格依赖于X和ε的联合高斯性(假设A1),但在intro与conclusion中,作者有时用“Gaussian covariates with correlated noise”而非“jointly Gaussian”。读者需留心——如果X或ε中任一个是非高斯,CGMT扩展引理可能失效。此外,条件E中的\(s\)通常依赖于在证明中使用的一个固定阈值,理论上\(s\)的选择需要匹配某个维度的增长速度——但论文没有给出一个显式的闭合条件\(s = s(p,n)\)(除了某些正则情形)。因此,用户若想检查论文证明是否紧,可以从“\(s\)是否可显式替代导致错误”入手 —— 比方说,如果\(\Sigma_X\)的谱衰减非常陡峭(远端特征值极短尾),条件E看似被自动满足,但\(s\)可能选择不恰当,导致分解无效。这个窄处正是可做follow-up点。
四、开放问题(点到为止,扎根具体语句)¶
-
非高斯联合分布的推广:论文假设“jointly Gaussian”(Section 2, Assumption A1)。Chen & Xie等对CGMT已有某些非高斯(Montanari & Saeed 2022, Han & Shen 2022)的universal扩展;是否可以将这些结果直接用于内生性情形?在什么条件下(如\(\varepsilon\)的协方差结构、X的条目矩条件)内生产壮的CGMT也可推广?扎根于:论文末尾“Limitations”一节(这一部分在full text中给出,第一段:“The main limitation is that our analysis relies on CGMT, which requires the rows of X to be Gaussian”)。
-
显式的\(s\)选择规则与速率:条件E中截断索引\(s\)选择与有效秩的关系,在论文中主要通过Bartlett et al. (2020)中的\(R(\Sigma)\)与\(\tilde{r}(\Sigma)\)引出,但未给出一个关于\(s = s(\{\lambda_j, \omega_j\})\)的显式构造条件——这对实证应用(需要我如何计算\(s\)?)有障碍。扎根于:Section 3中的Assumption 4(\(s\)依赖未知谱率)和Lemma 8(\(s\)隐藏在有效秩定义中)。一个高价值的follow-up可能是推导出基于特征值经验分布的显式可计算阈值\(s\)。
-
minimax下界:条件E是否是必要且几乎充分的?Bartlett (2020)已给出独立噪声的minimax最优界;但内生性下,是否能证明\(\sum \omega_j / \sum \lambda_j \to 0\)是渐近一致性的必要条件?是否存在一个minimax下界族表明,如果该比值不趋于0,那么任何无正则估计量误差都不能收敛——这与简单的算法无关,而是统计信息界线?扎根于:本文未讨论下界(lower bound),而只给了充分条件。
-
与正则化+去偏路线的对比问题:Belloni/Chernozhukov路线(依赖稀疏假设)与本文的路线(依赖噪声特征值比)之间是否存在交集?会不会有一个问题族,既不是稀疏的,又不满足条件E——那这个设定下是否没有任何一致性估计量?找出这样的“硬问题族”可视为一个有趣的下界结果。扎根于:本文intro与Belloni et al. (2017)(引用[21])没有进行比较或冲突分析。
Maintained by 陈星宇 · Homepage · Source on GitHub