Minimum mean square distance estimation of a subspace

01/18/2011
by   Olivier Besson, et al.
0

We consider the problem of subspace estimation in a Bayesian setting. Since we are operating in the Grassmann manifold, the usual approach which consists of minimizing the mean square error (MSE) between the true subspace U and its estimate Û may not be adequate as the MSE is not the natural metric in the Grassmann manifold. As an alternative, we propose to carry out subspace estimation by minimizing the mean square distance (MSD) between U and its estimate, where the considered distance is a natural metric in the Grassmann manifold, viz. the distance between the projection matrices. We show that the resulting estimator is no longer the posterior mean of U but entails computing the principal eigenvectors of the posterior mean of U U^T. Derivation of the MMSD estimator is carried out in a few illustrative examples including a linear Gaussian model for the data and a Bingham or von Mises Fisher prior distribution for U. In all scenarios, posterior distributions are derived and the MMSD estimator is obtained either analytically or implemented via a Markov chain Monte Carlo simulation method. The method is shown to provide accurate estimates even when the number of samples is lower than the dimension of U. An application to hyperspectral imagery is finally investigated.

READ FULL TEXT
research
02/07/2020

Explicit Mean-Square Error Bounds for Monte-Carlo and Linear Stochastic Approximation

This paper concerns error bounds for recursive equations subject to Mark...
research
06/01/2016

Hyperspectral Subspace Identification Using SURE

Remote sensing hyperspectral sensors collect large volumes of high dimen...
research
10/05/2013

Moments and Root-Mean-Square Error of the Bayesian MMSE Estimator of Classification Error in the Gaussian Model

The most important aspect of any classifier is its error rate, because t...
research
10/01/2013

Joint Bayesian estimation of close subspaces from noisy measurements

In this letter, we consider two sets of observations defined as subspace...
research
09/28/2020

Replica Analysis of the Linear Model with Markov or Hidden Markov Signal Priors

This paper estimates free energy, average mutual information, and minimu...
research
05/06/2021

Random density matrices: analytical results for mean root fidelity and mean square Bures distance

Bures distance holds a special place among various distance measures due...

Please sign up or login with your details

Forgot password? Click here to reset