Gridding and Parameter Expansion for Scalable Latent Gaussian Models of Spatial Multivariate Data

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Peruzzi, Michele, Banerjee, Sudipto, Dunson, David B., Finley, Andrew O.
Natura: Preprint
Pubblicazione: 2021
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866915216737959936
author Peruzzi, Michele
Banerjee, Sudipto
Dunson, David B.
Finley, Andrew O.
author_facet Peruzzi, Michele
Banerjee, Sudipto
Dunson, David B.
Finley, Andrew O.
contents Scalable spatial GPs for massive datasets can be built via sparse Directed Acyclic Graphs (DAGs) where a small number of directed edges is sufficient to flexibly characterize spatial dependence. The DAG can be used to devise fast algorithms for posterior sampling of the latent process, but these may exhibit pathological behavior in estimating covariance parameters. In this article, we introduce gridding and parameter expansion methods to improve the practical performance of MCMC algorithms in terms of effective sample size per unit time (ESS/s). Gridding is a model-based strategy that reduces the number of expensive operations necessary during MCMC on irregularly spaced data. Parameter expansion reduces dependence in posterior samples in spatial regression for high resolution data. These two strategies lead to computational gains in the big data settings on which we focus. We consider popular constructions of univariate spatial processes based on Matérn covariance functions and multivariate coregionalization models for Gaussian outcomes in extensive analyses of synthetic datasets comparing with alternative methods. We demonstrate effectiveness of our proposed methods in a forestry application using remotely sensed data from NASA's Goddard LiDAR, Hyper-Spectral, and Thermal imager (G-LiHT).
format Preprint
id arxiv_https___arxiv_org_abs_2101_03579
institution arXiv
publishDate 2021
record_format arxiv
spellingShingle Gridding and Parameter Expansion for Scalable Latent Gaussian Models of Spatial Multivariate Data
Peruzzi, Michele
Banerjee, Sudipto
Dunson, David B.
Finley, Andrew O.
Methodology
Applications
Computation
Scalable spatial GPs for massive datasets can be built via sparse Directed Acyclic Graphs (DAGs) where a small number of directed edges is sufficient to flexibly characterize spatial dependence. The DAG can be used to devise fast algorithms for posterior sampling of the latent process, but these may exhibit pathological behavior in estimating covariance parameters. In this article, we introduce gridding and parameter expansion methods to improve the practical performance of MCMC algorithms in terms of effective sample size per unit time (ESS/s). Gridding is a model-based strategy that reduces the number of expensive operations necessary during MCMC on irregularly spaced data. Parameter expansion reduces dependence in posterior samples in spatial regression for high resolution data. These two strategies lead to computational gains in the big data settings on which we focus. We consider popular constructions of univariate spatial processes based on Matérn covariance functions and multivariate coregionalization models for Gaussian outcomes in extensive analyses of synthetic datasets comparing with alternative methods. We demonstrate effectiveness of our proposed methods in a forestry application using remotely sensed data from NASA's Goddard LiDAR, Hyper-Spectral, and Thermal imager (G-LiHT).
title Gridding and Parameter Expansion for Scalable Latent Gaussian Models of Spatial Multivariate Data
topic Methodology
Applications
Computation
url https://arxiv.org/abs/2101.03579