Average-Case Integrality Gap for Non-Negative Principal Component Analysis
Average-Case Integrality Gap for Non-Negative Principal Component Analysis
复制标题
非负主成分分析的平均情况完整性差距
DOI:
--
复制
发表时间:
2020
期刊:
影响因子:
--
通讯作者:
Alexander S. Wein
中科院分区:
文献类型:
--
作者:
A. Bandeira;Dmitriy Kunisky;Alexander S. Wein
Montanari and Richard (2015) asked whether a natural semidefinite programming (SDP) relaxation can effectively optimize $\mathbf{x}^{\top}\mathbf{W} \mathbf{x}$ over $\|\mathbf{x}\| = 1$ with $x_i \geq 0$ for all coordinates $i$, where $\mathbf{W} \in \mathbb{R}^{n \times n}$ is drawn from the Gaussian orthogonal ensemble (GOE) or a spiked matrix model. In small numerical experiments, this SDP appears to be tight for the GOE, producing a rank-one optimal matrix solution aligned with the optimal vector $\mathbf{x}$. We prove, however, that as $n \to \infty$ the SDP is not tight, and certifies an upper bound asymptotically no better than the simple spectral bound $\lambda_{\max}(\mathbf{W})$ on this objective function. We also provide evidence, using tools from recent literature on hypothesis testing with low-degree polynomials, that no subexponential-time certification algorithm can improve on this behavior. Finally, we present further numerical experiments estimating how large $n$ would need to be before this limiting behavior becomes evident, providing a cautionary example against extrapolating asymptotics of SDPs in high dimension from their efficacy in small "laptop scale" computations.
DOI:
10.1214/18-aos1697
发表时间:
2015-09
期刊:
The Annals of Statistics
影响因子:
--
作者:
I. Johnstone;A. Onatski
通讯作者:
I. Johnstone;A. Onatski