Scalable Spatio-Temporal Smoothing via Hierarchical Sparse Cholesky Decomposition

07/19/2022
by   Marcin Jurek, et al.
0

We propose an approximation to the forward-filter-backward-sampler (FFBS) algorithm for large-scale spatio-temporal smoothing. FFBS is commonly used in Bayesian statistics when working with linear Gaussian state-space models, but it requires inverting covariance matrices which have the size of the latent state vector. The computational burden associated with this operation effectively prohibits its applications in high-dimensional settings. We propose a scalable spatio-temporal FFBS approach based on the hierarchical Vecchia approximation of Gaussian processes, which has been previously successfully used in spatial statistics. On simulated and real data, our approach outperformed a low-rank FFBS approximation.

READ FULL TEXT

page 7

page 14

page 18

research
06/30/2020

Hierarchical sparse Cholesky decomposition with applications to high-dimensional spatio-temporal filtering

Spatial statistics often involves Cholesky decomposition of covariance m...
research
03/16/2022

Scalable marginalization of latent variables for correlated data

Marginalization of latent variables or nuisance parameters is a fundamen...
research
12/29/2021

Correlation-based sparse inverse Cholesky factorization for fast Gaussian-process inference

Gaussian processes are widely used as priors for unknown functions in st...
research
03/05/2018

Banded Spatio-Temporal Autoregressions

We propose a new class of spatio-temporal models with unknown and banded...
research
04/07/2022

Neural Implicit Flow: a mesh-agnostic dimensionality reduction paradigm of spatio-temporal data

High-dimensional spatio-temporal dynamics can often be encoded in a low-...
research
08/19/2022

Exploring seismic hazard in the Groningen gas field using adaptive kernel smoothing and inhomogeneous summary statistics

The discovery of gas in Groningen in 1959 has been a massive boon to the...

Please sign up or login with your details

Forgot password? Click here to reset