Efficient methods for Gaussian Markov random fields under sparse linear constraints

06/03/2021
by   David Bolin, et al.
0

Methods for inference and simulation of linearly constrained Gaussian Markov Random Fields (GMRF) are computationally prohibitive when the number of constraints is large. In some cases, such as for intrinsic GMRFs, they may even be unfeasible. We propose a new class of methods to overcome these challenges in the common case of sparse constraints, where one has a large number of constraints and each only involves a few elements. Our methods rely on a basis transformation into blocks of constrained versus non-constrained subspaces, and we show that the methods greatly outperform existing alternatives in terms of computational cost. By combining the proposed methods with the stochastic partial differential equation approach for Gaussian random fields, we also show how to formulate Gaussian process regression with linear constraints in a GMRF setting to reduce computational cost. This is illustrated in two applications with simulated data.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
01/30/2022

The Circular Matern Covariance Function and its Link to Markov Random Fields on the Circle

The link between Gaussian random fields and Markov random fields is well...
research
04/24/2023

Faster estimation of the Knorr-Held Type IV space-time model

In this paper we study the type IV Knorr Held space time models. Such mo...
research
08/06/2013

Invariances of random fields paths, with applications in Gaussian Process Regression

We study pathwise invariances of centred random fields that can be contr...
research
02/18/2018

Spatial modelling with R-INLA: A review

Coming up with Bayesian models for spatial data is easy, but performing ...
research
06/21/2022

Efficient Inference of Spatially-varying Gaussian Markov Random Fields with Applications in Gene Regulatory Networks

In this paper, we study the problem of inferring spatially-varying Gauss...
research
04/17/2018

Bayesian model-data synthesis with an application to global Glacio-Isostatic Adjustment

We introduce a framework for updating large scale geospatial processes u...
research
05/17/2018

Efficient simulation of Gaussian Markov random fields by Chebyshev polynomial approximation

This paper presents an algorithm to simulate Gaussian random vectors who...

Please sign up or login with your details

Forgot password? Click here to reset