Abstract:Riemannian Gaussian distributions were initially introduced as basic building blocks for learning models which aim to capture the intrinsic structure of statistical populations of positive-definite matrices (here called covariance matrices). While the potential applications of such models have attracted significant attention, a major obstacle still stands in the way of these applications: there seems to exist no practical method of computing the normalising factors associated with Riemannian Gaussian distributions on spaces of high-dimensional covariance matrices. The present paper shows that this missing method comes from an unexpected new connection with random matrix theory. Its main contribution is to prove that Riemannian Gaussian distributions of real, complex, or quaternion covariance matrices are equivalent to orthogonal, unitary, or symplectic log-normal matrix ensembles. This equivalence yields a highly efficient approximation of the normalising factors, in terms of a rather simple analytic expression. The error due to this approximation decreases like the inverse square of dimension. Numerical experiments are conducted which demonstrate how this new approximation can unlock the difficulties which have impeded applications to real-world datasets of high-dimensional covariance matrices. The paper then turns to Riemannian Gaussian distributions of block-Toeplitz covariance matrices. These are equivalent to yet another kind of random matrix ensembles, here called "acosh-normal" ensembles. Orthogonal and unitary "acosh-normal" ensembles correspond to the cases of block-Toeplitz with Toeplitz blocks, and block-Toeplitz (with general blocks) covariance matrices, respectively.

Lectures on random matrix models. The Riemann-Hilbert approach

Asymptotic analysis of random partitions ∗

Precise Asymptotics for Random Matrices and Random Growth Models

Non-asymptotic theory of random matrices: extreme singular values

A Riemann--Hilbert approach to the perturbation theory for orthogonal polynomials: Applications to numerical linear algebra and random matrix theory

Large n limit of Gaussian random matrices with external source, part I

Universality of a double scaling limit near singular edge points in random matrix models

A path integral approach to sparse non-Hermitian random matrices

Cleaning large correlation matrices: Tools from Random Matrix Theory

Application of Random Matrix Theory in High-Dimensional Statistics

A Riemann-Hilbert problem for biorthogonal polynomials

Multiple orthogonal polynomials in random matrix theory

Lecture Notes on "Random Matrices"

Universality for random matrices with an edge spectrum singularity

A universality theorem for ratios of random characteristic polynomials

Random and Integrable Models in Mathematics and Physics

Random vector and matrix and vector theories: a renormalization group approach

Global Asymptotics of Hermite Polynomials Via Riemann-Hilbert Approach

Precise asymptotics for the spectral radius of a large random matrix

Riemannian statistics meets random matrix theory: towards learning from high-dimensional covariance matrices

A martingale approach to Gaussian fluctuations and laws of iterated logarithm for Ewens-Pitman model