Geometrical and statistical properties of M-estimates of scatter on Grassmann manifolds
We consider data from the Grassmann manifold G(m,r) of all vector subspaces of dimension r of Rm, and focus on the Grassmannian statistical model which is of common use in signal processing and statistics. Canonical Grassmannian distributions GΣ on G(m,r) are indexed by parameters Σ from the manifold M=Possym1(m) of positive definite symmetric matrices of determinant 1. Robust M-estimates of scatter (GE) for general probability measures P on G(m,r) are studied. Such estimators are defined to be the maximizers of the Grassmannian log-likelihood −ℓP(Σ) as function of Σ. One of the novel features of this work is a strong use of the fact that M is a CAT(0) space with known visual boundary at infinity ∂M. We also recall that the sample space G(m,r) is a part of ∂M, show the distributions GΣ are SL(m,R)--quasi-invariant, and that ℓP(Σ) is a weighted Busemann function. Let Pn=(δU1+⋯+δUn)/n be the empirical probability measure for n-samples of random i.i.d. subspaces Ui∈G(m,r) of common distribution P, whose support spans Rm. For Σn and ΣP the GEs of Pn and P, we show the almost sure convergence of Σn towards Σ as n→∞ using methods from geometry, and provide a central limit theorem for the rescaled process Cn=tr(ΣP−1Σn)mg−1Σng−1, where Σ=gg with g∈SL(m,R) the unique symmetric positive-definite square root of Σ.