MLLGSTSep 19, 2023

Posterior Contraction Rates for Matérn Gaussian Processes on Riemannian Manifolds

arXiv:2309.10918v315 citationsh-index: 32
Originality Incremental advance
AI Analysis

This work addresses the problem of uncertainty quantification in geometric machine learning for researchers and practitioners, providing foundational theoretical insights but is incremental as it builds on prior computational tools and analyses.

The paper tackles the theoretical performance comparison between intrinsic Matérn Gaussian processes on Riemannian manifolds and extrinsic Euclidean processes, proving optimal contraction rates for both and showing they coincide when smoothness parameters are matched, with empirical illustrations indicating intrinsic processes can achieve better practical performance.

Gaussian processes are used in many machine learning applications that rely on uncertainty quantification. Recently, computational tools for working with these models in geometric settings, such as when inputs lie on a Riemannian manifold, have been developed. This raises the question: can these intrinsic models be shown theoretically to lead to better performance, compared to simply embedding all relevant quantities into $\mathbb{R}^d$ and using the restriction of an ordinary Euclidean Gaussian process? To study this, we prove optimal contraction rates for intrinsic Matérn Gaussian processes defined on compact Riemannian manifolds. We also prove analogous rates for extrinsic processes using trace and extension theorems between manifold and ambient Sobolev spaces: somewhat surprisingly, the rates obtained turn out to coincide with those of the intrinsic processes, provided that their smoothness parameters are matched appropriately. We illustrate these rates empirically on a number of examples, which, mirroring prior work, show that intrinsic processes can achieve better performance in practice. Therefore, our work shows that finer-grained analyses are needed to distinguish between different levels of data-efficiency of geometric Gaussian processes, particularly in settings which involve small data set sizes and non-asymptotic behavior.

Code Implementations1 repo
Foundations

The foundational work for this paper's niche, ranked by how specifically the neighbourhood builds on it — not by global fame.

Your Notes