Table of Contents
Fetching ...

Near-Optimal Tensor PCA via Normalized Stochastic Gradient Ascent with Overparameterization

Shihong Ding, Yihong Gu, Yuanshi Liu, Cong Fang

TL;DR

This work tackles tensor PCA in the spiked tensor model by proposing a normalized stochastic gradient ascent method with overparameterization (NSGA). It demonstrates that, without global or spectral initialization, the planted vector $v_*$ can be recovered with constant probability whenever $N\lambda^2\ge\widetilde{Ω}(d^{\lceil k/2\rceil})$ for the odd/even cases, achieving near-optimal sample complexity akin to $d^{k/2}$ thresholds known from SoS relaxations. The analysis reveals two-phase learning dynamics: an initial alignment phase that pushes toward the rank-1 signal and a secondary estimation phase where the dynamics approximate a one-dimensional SGD with controllable noise, benefiting from overparameterization to avoid trapping and improve generalization. This provides the first solid theoretical evidence that overparameterization can yield statistical advantages beyond exact parameterization in nonconvex tensor problems and points to potential extensions to other homogeneous models and online learning settings.

Abstract

We study the Order-$k$ ($k \geq 4$) spiked tensor model for the tensor principal component analysis (PCA) problem: given $N$ i.i.d. observations of a $k$-th order tensor generated from the model $\mathbf{T} = λ\cdot v_*^{\otimes k} + \mathbf{E}$, where $λ> 0$ is the signal-to-noise ratio (SNR), $v_*$ is a unit vector, and $\mathbf{E}$ is a random noise tensor, the goal is to recover the planted vector $v_*$. We propose a normalized stochastic gradient ascent (NSGA) method with overparameterization for solving the tensor PCA problem. Without any global (or spectral) initialization step, the proposed algorithm successfully recovers the signal $v_*$ when $Nλ^2 \geq \widetildeΩ(d^{\lceil k/2 \rceil})$, thereby breaking the previous conjecture that (stochastic) gradient methods require at least $Ω(d^{k-1})$ samples for recovery. For even $k$, the $\widetildeΩ(d^{k/2})$ threshold coincides with the optimal threshold under computational constraints, attained by sum-of-squares relaxations and related algorithms. Theoretical analysis demonstrates that the overparameterized stochastic gradient method not only establishes a significant initial optimization advantage during the early learning phase but also achieves strong generalization guarantees. This work provides the first evidence that overparameterization improves statistical performance relative to exact parameterization that is solved via standard continuous optimization.

Near-Optimal Tensor PCA via Normalized Stochastic Gradient Ascent with Overparameterization

TL;DR

This work tackles tensor PCA in the spiked tensor model by proposing a normalized stochastic gradient ascent method with overparameterization (NSGA). It demonstrates that, without global or spectral initialization, the planted vector can be recovered with constant probability whenever for the odd/even cases, achieving near-optimal sample complexity akin to thresholds known from SoS relaxations. The analysis reveals two-phase learning dynamics: an initial alignment phase that pushes toward the rank-1 signal and a secondary estimation phase where the dynamics approximate a one-dimensional SGD with controllable noise, benefiting from overparameterization to avoid trapping and improve generalization. This provides the first solid theoretical evidence that overparameterization can yield statistical advantages beyond exact parameterization in nonconvex tensor problems and points to potential extensions to other homogeneous models and online learning settings.

Abstract

We study the Order- () spiked tensor model for the tensor principal component analysis (PCA) problem: given i.i.d. observations of a -th order tensor generated from the model , where is the signal-to-noise ratio (SNR), is a unit vector, and is a random noise tensor, the goal is to recover the planted vector . We propose a normalized stochastic gradient ascent (NSGA) method with overparameterization for solving the tensor PCA problem. Without any global (or spectral) initialization step, the proposed algorithm successfully recovers the signal when , thereby breaking the previous conjecture that (stochastic) gradient methods require at least samples for recovery. For even , the threshold coincides with the optimal threshold under computational constraints, attained by sum-of-squares relaxations and related algorithms. Theoretical analysis demonstrates that the overparameterized stochastic gradient method not only establishes a significant initial optimization advantage during the early learning phase but also achieves strong generalization guarantees. This work provides the first evidence that overparameterization improves statistical performance relative to exact parameterization that is solved via standard continuous optimization.
Paper Structure (16 sections, 30 theorems, 179 equations, 2 figures, 2 algorithms)

This paper contains 16 sections, 30 theorems, 179 equations, 2 figures, 2 algorithms.

Key Result

Corollary 3.1

Suppose Assumption ass-2 holds. For any unit vector $u\in\mathbb{R}^{d^k}$, the following second-order moment condition holds for $\left\langle u,\mathop{\mathrm{flat}}\limits(\mathbf{E})\right\rangle$

Figures (2)

  • Figure 1: The performance of the algorithms for tensor PCA ($k$ even). Combining overparameterization, our algorithm elevates the recovery threshold of gradient-based methods from the light blue region to the green striped region.
  • Figure 2: The figure shows the convergence behavior among several gradient-based methods, where the horizontal axis represents the number of iterations (or sample size) $N$, and the vertical axis represents the squared estimation error $\left\|\hat{v}^{(N)}-v_*\right\|^2$. The signal vector $v_*$ is sampled uniformly from the unit sphere. In the comparison, "NSGA (Matrix)" denotes our proposed method (corresponding to Algorithm \ref{['SGA']} when the order $k$ is even, and Algorithm \ref{['k-odd']} when $k$ is odd). "SGA (Vector)" represents the standard stochastic gradient ascent under exact parameterization, while "SGA_PS (Vector)" represents the exact parameterization of SGA with projection onto the sphere arous2020online. Finally, "SGA_Acc (Vector)" corresponds to the accelerated variant of SGA lan2020first.

Theorems & Definitions (59)

  • Corollary 3.1
  • Theorem 4.1
  • Corollary 4.1
  • Theorem 4.2
  • proof
  • Corollary 4.2
  • Corollary 4.3
  • Corollary 4.4
  • Remark 4.1
  • Lemma A.1
  • ...and 49 more