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.

Fast estimation of the look-elsewhere effect using Gaussian random fields

Fast Gaussian particle filtering algorithm

The profile likelihood ratio and the look elsewhere effect in high energy physics

On fast Generation of Cosmological Random Fields

Optimal and fast detection of spatial clusters with scan statistics

Fast sampling of parameterised Gaussian random fields

Memory parameter estimation for long range dependent random fields

Linear approximation to the statistical significance autocovariance matrix in the asymptotic regime

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

Simultaneous Inference for Non-Stationary Random Fields, with Application to Gridded Data Analysis

Non-Parametric Approximations for Anisotropy Estimation in Two-dimensional Differentiable Gaussian Random Fields

Simulation of Gaussian random field in a ball

High-Frequency Tail Index Estimation by Nearly Tight Frames

Adaptive Bayesian estimation using a Gaussian random field with inverse Gamma bandwidth

Fast likelihood evaluation using meshfree approximations for reconstructing compact binary sources

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

Fast marginalization algorithm for optimizing gravitational wave detection, parameter estimation and sky localization

Applying Simulation-Based Inference to Spectral and Spatial Information from the Galactic Center Gamma-Ray Excess

Generalized Steppingstone Sampling: Efficient marginal likelihood estimation in gravitational wave analysis of Pulsar Timing Array data

Sampling Gaussian stationary random fields: A stochastic realization approach

An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach