19,205 research outputs found
OptShrink: An algorithm for improved low-rank signal matrix denoising by optimal, data-driven singular value shrinkage
The truncated singular value decomposition (SVD) of the measurement matrix is
the optimal solution to the_representation_ problem of how to best approximate
a noisy measurement matrix using a low-rank matrix. Here, we consider the
(unobservable)_denoising_ problem of how to best approximate a low-rank signal
matrix buried in noise by optimal (re)weighting of the singular vectors of the
measurement matrix. We exploit recent results from random matrix theory to
exactly characterize the large matrix limit of the optimal weighting
coefficients and show that they can be computed directly from data for a large
class of noise models that includes the i.i.d. Gaussian noise case.
Our analysis brings into sharp focus the shrinkage-and-thresholding form of
the optimal weights, the non-convex nature of the associated shrinkage function
(on the singular values) and explains why matrix regularization via singular
value thresholding with convex penalty functions (such as the nuclear norm)
will always be suboptimal. We validate our theoretical predictions with
numerical simulations, develop an implementable algorithm (OptShrink) that
realizes the predicted performance gains and show how our methods can be used
to improve estimation in the setting where the measured matrix has missing
entries.Comment: Published version. The algorithm can be downloaded from
http://www.eecs.umich.edu/~rajnrao/optshrin
An early warning indicator for atmospheric blocking events using transfer operators
The existence of persistent midlatitude atmospheric flow regimes with
time-scales larger than 5-10 days and indications of preferred transitions
between them motivates to develop early warning indicators for such regime
transitions. In this paper, we use a hemispheric barotropic model together with
estimates of transfer operators on a reduced phase space to develop an early
warning indicator of the zonal to blocked flow transition in this model. It is
shown that, the spectrum of the transfer operators can be used to study the
slow dynamics of the flow as well as the non-Markovian character of the
reduction. The slowest motions are thereby found to have time scales of three
to six weeks and to be associated with meta-stable regimes (and their
transitions) which can be detected as almost-invariant sets of the transfer
operator. From the energy budget of the model, we are able to explain the
meta-stability of the regimes and the existence of preferred transition paths.
Even though the model is highly simplified, the skill of the early warning
indicator is promising, suggesting that the transfer operator approach can be
used in parallel to an operational deterministic model for stochastic
prediction or to assess forecast uncertainty
Optimal Estimation and Rank Detection for Sparse Spiked Covariance Matrices
This paper considers sparse spiked covariance matrix models in the
high-dimensional setting and studies the minimax estimation of the covariance
matrix and the principal subspace as well as the minimax rank detection. The
optimal rate of convergence for estimating the spiked covariance matrix under
the spectral norm is established, which requires significantly different
techniques from those for estimating other structured covariance matrices such
as bandable or sparse covariance matrices. We also establish the minimax rate
under the spectral norm for estimating the principal subspace, the primary
object of interest in principal component analysis. In addition, the optimal
rate for the rank detection boundary is obtained. This result also resolves the
gap in a recent paper by Berthet and Rigollet [1] where the special case of
rank one is considered
- …