中文

高维噪声矩阵中稀疏子矩阵的检测

统计理论 2013-12-20 v2 统计理论

摘要

我们观测到一个 N×MN\times M 矩阵 Yij=sij+ξijY_{ij}=s_{ij}+\xi_{ij},其中 ξijN(0,1)\xi_{ij}\sim {\mathcal {N}}(0,1)i,ji,j 上独立同分布,且 sijRs_{ij}\in \mathbb {R}。我们检验原假设:对所有 i,ji,jsij=0s_{ij}=0,备择假设为:存在某个大小为 n×mn\times m 的子矩阵,其元素显著,即 sija>0s_{ij}\ge a>0。我们提出了一种检验程序,并计算了渐近检测边界 aa,使得当 MM\to\infty, NN\to\infty, p=n/N0p=n/N\to0, q=m/M0q=m/M\to0 时,最大检验风险趋于 0。我们证明了在某些附加约束下,该边界是渐近锐利极小极大的。我们讨论了与其他检验问题的关系。我们提出了一种能在给定集合内适应未知 (n,m)(n,m) 的检验程序,并计算了自适应锐利速率。我们的检验程序在合成数据上的实现显示出对稀疏、非必为方阵的矩阵具有优良表现。我们将锐利极小极大结果扩展到不同方向:首先,扩展到方差未知的高斯矩阵;其次,扩展到具有指数族分布(非高斯)的随机变量矩阵;最后,扩展到针对高斯元素矩阵的双边备择假设。

关键词

引用

@article{arxiv.1109.0898,
  title  = {Detection of a sparse submatrix of a high-dimensional noisy matrix},
  author = {Cristina Butucea and Yuri I. Ingster},
  journal= {arXiv preprint arXiv:1109.0898},
  year   = {2013}
}

备注

Published in at http://dx.doi.org/10.3150/12-BEJ470 the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm)