Abstract:The distribution of centered Gaussian random fields (GRFs) indexed by compacta such as smooth, bounded Euclidean domains or smooth, compact and orientable manifolds is determined by their covariance operators. We consider centered GRFs given as variational solutions to coloring operator equations driven by spatial white noise, with an elliptic self-adjoint pseudodifferential coloring operator from the Hörmander class. This includes the Matérn class of GRFs as a special case. Using biorthogonal multiresolution analyses on the manifold, we prove that the precision and covariance operators, respectively, may be identified with bi-infinite matrices and finite sections may be diagonally preconditioned rendering the condition number independent of the dimension p of this section. We prove that a tapering strategy by thresholding applied on finite sections of the bi-infinite precision and covariance matrices results in optimally numerically sparse approximations. That is, asymptotically only linearly many nonzero matrix entries are sufficient to approximate the original section of the bi-infinite covariance or precision matrix using this tapering strategy to arbitrary precision. The locations of these nonzero matrix entries can be determined a priori. The tapered covariance or precision matrices may also be optimally diagonally preconditioned. Analysis of the relative size of the entries of the tapered covariance matrices motivates novel, multilevel Monte Carlo (MLMC) oracles for covariance estimation, in sample complexity that scales log-linearly with respect to the number p of parameters. In addition, we propose and analyze novel compressive algorithms for simulating and kriging of GRFs. The complexity (work and memory vs. accuracy) of these three algorithms scales near-optimally in terms of the number of parameters p of the sample-wise approximation of the GRF in Sobolev scales.

Random action of compact Lie groups and minimax estimation of a mean pattern

Precise Error Analysis of Regularized M-estimators in High-dimensions

Learning extremal graphical structures in high dimensions

Wasserstein-based Minimax Estimation of Dependence in Multivariate Regularly Varying Extremes

Geometric planted matchings beyond the Gaussian model

Noisy recovery from random linear observations: Sharp minimax rates under elliptical constraints

Kernel regression, minimax rates and effective dimensionality: beyond the regular case

Mostly contracting random maps

Multilevel approximation of Gaussian random fields: Covariance compression, estimation, and spatial prediction

Random fields of multivariate test statistics, with applications to shape analysis

On the Concentration of the Minimizers of Empirical Risks

Minimax rates for latent position estimation in the generalized random dot product graph

Adaptive estimation of a quadratic functional by model selection

Diffeomorphisms Groups and Pattern Matching in Image Analysis

Annealed quantitative estimates for the quadratic 2D-discrete random matching problem

New bounds for normal approximation on product spaces with applications to monochromatic edges, random sums and an infinite de Jong CLT

Statistical analysis of random objects via metric measure Laplacians

Spurious minimizers in non uniform Fourier sampling optimization

Smoothness Estimation for Whittle-Matérn Processes on Closed Riemannian Manifolds

Adaptive minimax estimation for discretely observed Lévy processes

Likelihood landscape and maximum likelihood estimation for the discrete orbit recovery model