INLA-MRA: A Bayesian method for large spatiotemporal datasets

04/21/2020
by   Luc Villandré, et al.
0

Large spatiotemporal datasets are a challenge for conventional Bayesian models, because of the cubic computational complexity of the algorithms for obtaining the Cholesky decomposition of the covariance matrix in the multivariate normal density. Moreover, standard numerical algorithms for posterior estimation, such as Markov Chain Monte Carlo (MCMC), are intractable in this context, as they require thousands, if not millions, of costly likelihood evaluations. To overcome those limitations, we propose INLA-MRA, a method that mixes an estimation algorithm inspired by INLA (Integrated Nested Laplace Approximation), and a model taking advantage of the sparse covariance structure produced by the Multi-Resolution Approximation (MRA) approach. INLA-MRA extends MRA to spatiotemporal data, while also facilitating the approximation of the hyperparameter marginal posterior distributions. We apply INLA-MRA to large MODIS Level 3 Land Surface Temperature (LST) datasets, sampled between May 18 and May 31, 2012 in the western part of the state of Maharashtra, India. We find that INLA-MRA can produce realistic prediction surfaces over regions where concentrated missingness, caused by sizable cloud cover, is observed. Through a validation analysis, we also find that predictions tend to be very accurate.

READ FULL TEXT

page 4

page 18

research
11/02/2019

Bayesian model averaging with the integrated nested Laplace approximation

The integrated nested Laplace approximation (INLA) for Bayesian inferenc...
research
04/07/2021

Laplace-aided variational inference for differential equation models

Ordinary differential equation (ODE) model whose regression curves are a...
research
03/04/2021

On MCMC for variationally sparse Gaussian processes: A pseudo-marginal approach

Gaussian processes (GPs) are frequently used in machine learning and sta...
research
09/16/2021

Statistical Inference for Bayesian Risk Minimization via Exponentially Tilted Empirical Likelihood

The celebrated Bernstein von-Mises theorem ensures that credible regions...
research
06/30/2020

Spatiotemporal Multi-Resolution Approximations for Analyzing Global Environmental Data

Technological developments and open data policies have made large, globa...
research
04/21/2020

Bayesian Optimization of Hyperparameters when the Marginal Likelihood is Estimated by MCMC

Bayesian models often involve a small set of hyperparameters determined ...
research
03/04/2020

Bayesian System ID: Optimal management of parameter, model, and measurement uncertainty

We evaluate the robustness of a probabilistic formulation of system iden...

Please sign up or login with your details

Forgot password? Click here to reset