Covariance-based rational approximations of fractional SPDEs for computationally efficient Bayesian inference

09/10/2022
by   Zhen Xiong, et al.
0

The stochastic partial differential equation (SPDE) approach is widely used for modeling large spatial datasets. It is based on representing a Gaussian random field u on ℝ^d as the solution of an elliptic SPDE L^β u = 𝒲 where L is a second-order differential operator, 2β (belongs to natural number starting from 1) is a positive parameter that controls the smoothness of u and 𝒲 is Gaussian white noise. A few approaches have been suggested in the literature to extend the approach to allow for any smoothness parameter satisfying β>d/4. Even though those approaches work well for simulating SPDEs with general smoothness, they are less suitable for Bayesian inference since they do not provide approximations which are Gaussian Markov random fields (GMRFs) as in the original SPDE approach. We address this issue by proposing a new method based on approximating the covariance operator L^-2β of the Gaussian field u by a finite element method combined with a rational approximation of the fractional power. This results in a numerically stable GMRF approximation which can be combined with the integrated nested Laplace approximation (INLA) method for fast Bayesian inference. A rigorous convergence analysis of the method is performed and the accuracy of the method is investigated with simulated data. Finally, we illustrate the approach and corresponding implementation in the R package rSPDE via an application to precipitation data which is analyzed by combining the rSPDE package with the R-INLA software for full Bayesian inference.

READ FULL TEXT

page 13

page 15

research
11/12/2017

The SPDE approach for Gaussian random fields with general smoothness

A popular approach for modeling and inference in spatial statistics is t...
research
11/24/2022

Numerical Approximation of Gaussian random fields on Closed Surfaces

We consider the numerical approximation of Gaussian random fields on clo...
research
04/13/2019

Regularity and convergence analysis in Sobolev and Hölder spaces for generalized Whittle-Matérn fields

We analyze several Galerkin approximations of a Gaussian random field ZD...
research
03/27/2023

Integrated Nested Laplace Approximations for Large-Scale Spatial-Temporal Bayesian Modeling

Bayesian inference tasks continue to pose a computational challenge. Thi...
research
04/10/2022

Parallelized integrated nested Laplace approximations for fast Bayesian inference

There is a growing demand for performing larger-scale Bayesian inference...
research
11/26/2019

Accounting for spatial varying sampling effort due to accessibility in Citizen Science data: A case study of moose in Norway

Citizen Scientists together with an increasing access to technology prov...
research
02/08/2023

Regularity and numerical approximation of fractional elliptic differential equations on compact metric graphs

The fractional differential equation L^β u = f posed on a compact metric...

Please sign up or login with your details

Forgot password? Click here to reset