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

by   Helmut Harbrecht, et al.

Centered Gaussian random fields (GRFs) indexed by compacta such as smooth, bounded Euclidean domains or smooth, compact and orientable manifolds are 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 are known 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 a novel compressive algorithm 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.


Data sparse multilevel covariance estimation in optimal complexity

We consider the ℋ^2-formatted compression and computational estimation o...

On the Numerical Approximation of the Karhunen-Loève Expansion for Random Fields with Random Discrete Data

Many physical and mathematical models involve random fields in their inp...

Scalable Parallel Factorizations of SDD Matrices and Efficient Sampling for Gaussian Graphical Models

Motivated by a sampling problem basic to computational statistical infer...

Equivalence of measures and asymptotically optimal linear prediction for Gaussian random fields with fractional-order covariance operators

We consider Gaussian measures μ, μ̃ on a separable Hilbert space, with f...

Adaptive Regularization for Weight Matrices

Algorithms for learning distributions over weight-vectors, such as AROW ...

Computing eigenvalues of semi-infinite quasi-Toeplitz matrices

A quasi-Toeplitz (QT) matrix is a semi-infinite matrix of the form A=T(a...

Bi-Parametric Operator Preconditioning

We extend the general operator preconditioning framework [R. Hiptmair, C...

Please sign up or login with your details

Forgot password? Click here to reset