English

Outlier-robust additive matrix decomposition

Statistics Theory 2024-01-08 v2 Statistics Theory

Abstract

We study least-squares trace regression when the parameter is the sum of a rr-low-rank matrix and a ss-sparse matrix and a fraction ϵ\epsilon of the labels is corrupted. For subgaussian distributions and feature-dependent noise, we highlight three needed design properties, each one derived from a different process inequality: a "product process inequality", "Chevet's inequality" and a "multiplier process inequality". These properties handle, simultaneously, additive decomposition, label contamination and design-noise interaction. They imply the near-optimality of a tractable estimator with respect to the effective dimensions deff,rd_{eff,r} and deff,sd_{eff,s} of the low-rank and sparse components, ϵ\epsilon and the failure probability δ\delta. The near-optimal rate is r(n,deff,r)+r(n,deff,s)+(1+log(1/δ))/n+ϵlog(1/ϵ)\mathsf{r}(n,d_{eff,r}) + \mathsf{r}(n,d_{eff,s}) + \sqrt{(1+\log(1/\delta))/n} + \epsilon\log(1/\epsilon), where r(n,deff,r)+r(n,deff,s)\mathsf{r}(n,d_{eff,r})+\mathsf{r}(n,d_{eff,s}) is the optimal rate in average with no contamination. Our estimator is adaptive to (s,r,ϵ,δ)(s,r,\epsilon,\delta) and, for fixed absolute constant c>0c>0, it attains the mentioned rate with probability 1δ1-\delta uniformly over all δexp(cn)\delta\ge\exp(-cn). Without matrix decomposition, our analysis also entails optimal bounds for a robust estimator adapted to the noise variance. Our estimators are based on "sorted" versions of Huber's loss. We present simulations matching the theory. In particular, it reveals the superiority of "sorted" Huber's losses over the classical Huber's loss.

Keywords

Cite

@article{arxiv.2310.19136,
  title  = {Outlier-robust additive matrix decomposition},
  author = {Philip Thompson},
  journal= {arXiv preprint arXiv:2310.19136},
  year   = {2024}
}

Comments

This paper studies a broader model but shares content with arXiv:2012.06750 (which will not be further revised). Correction of typos, additional simulations, removal of robust matrix completion. Unlike mentioned in arXiv:2012.06750, (2018) Bellec et all DOES achieve the optimal rate for uncorrupted sparse linear regression (but assuming noise independent of features)

R2 v1 2026-06-28T13:05:16.907Z