Mathematics > Statistics Theory
[Submitted on 8 Apr 2013 (this version), latest version 4 Nov 2014 (v4)]
Title:Minimax Risk of Matrix Denoising by Singular Value Thresholding
View PDFAbstract:An unknown $m$ by $n$ matrix $X_0$ is to be estimated from noisy measurements $Y = X_0 + Z$, where the noise matrix $Z$ has i.i.d Gaussian entries. A popular matrix denoising scheme solves the nuclear norm penalization problem $\min_X || Y - X ||_F^2/2 + \lambda ||X||_* $, where $ ||X||_*$ denotes the nuclear norm (sum of singular values). This is the analog, for matrices, of $\ell_1$ penalization in the vector case. It has been empirically observed that, if $X_0$ has low rank, it may be recovered quite accurately from the noisy measurement $Y$. In a proportional growth framework where the rank $r_n$, number of rows $m_n$ and number of columns $n$ all tend to $\infty$ proportionally to each other ($ r_n/m_n -> \rho$, $m_n/n ->\beta$), we evaluate the asymptotic minimax MSE $M(\rho, \beta) = \lim_{m_n,n \goto \infty} \inf_\lambda \sup_{rank(X) \leq r_n} MSE(X,\hat{X}_\lambda)$ Our formulas involve incomplete moments of the quarter- and semi-circle laws ($\beta = 1$, square case) and the Marčenko-Pastur law ($\beta < 1$, non square case). We also show that any least-favorable matrix $X_0$ has norm "at infinity". The nuclear norm penalization problem is solved by applying soft thresholding to the singular values of $Y$. We also derive the minimax threshold, namely the value $\lambda^*(\rho)$ which is the optimal place to threshold the singular values. All these results are obtained for general (non square, non symmetric) real matrices. Comparable results are obtained for square symmetric nonnegative- definite matrices.
Submission history
From: Matan Gavish [view email][v1] Mon, 8 Apr 2013 01:39:01 UTC (140 KB)
[v2] Thu, 11 Apr 2013 04:55:32 UTC (141 KB)
[v3] Mon, 15 Apr 2013 02:13:10 UTC (142 KB)
[v4] Tue, 4 Nov 2014 12:45:39 UTC (224 KB)
Current browse context:
math.ST
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
Connected Papers (What is Connected Papers?)
CORE Recommender (What is CORE?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.