A fast and slightly robust covariance estimator
Abstract
Let Z = \Z1, …, Zn\ i.i.d. P ⊂ Rd from a distribution P with mean zero and covariance . Given a dataset X such that dham(X, Z) ≤ n, we are interested in finding an efficient estimator that achieves err(, ) := \|-12-12 - I\| op ≤ 1/2. We focus on the low contamination regime = o(1/d). In this regime, prior work required either (d3/2) samples or runtime that is exponential in d. We present an algorithm that, for subgaussian data, has near-linear sample complexity n = (d) and runtime O((n+d)ω + 12), where ω is the matrix multiplication exponent. We also show that this algorithm works for heavy-tailed data with near-linear sample complexity, but in a smaller regime of . Concurrent to our work, Diakonikolas et al. [2024] give Sum-of-Squares estimators that achieve similar sample complexity but with large polynomial runtime.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.