Scalable Bayesian transport maps for high-dimensional non-Gaussian spatial fields

08/09/2021
by   Matthias Katzfuss, et al.
0

A multivariate distribution can be described by a triangular transport map from the target distribution to a simple reference distribution. We propose Bayesian nonparametric inference on the transport map by modeling its components using Gaussian processes. This enables regularization and accounting for uncertainty in the map estimation, while still resulting in a closed-form and invertible posterior map. We then focus on inferring the distribution of a nonstationary spatial field from a small number of replicates. We develop specific transport-map priors that are highly flexible and are motivated by the behavior of a large class of stochastic processes. Our approach is scalable to high-dimensional fields due to data-dependent sparsity and parallel computations. We also discuss extensions, including Dirichlet process mixtures for marginal non-Gaussianity. We present numerical results to demonstrate the accuracy, scalability, and usefulness of our methods, including statistical emulation of non-Gaussian climate-model output.

READ FULL TEXT

page 2

page 9

page 14

page 18

05/31/2020

Spatial Factor Modeling: A Bayesian Matrix-Normal Approach for Misaligned Data

Multivariate spatially-oriented data sets are prevalent in the environme...
05/31/2019

Greedy inference with layers of lazy maps

We propose a framework for the greedy approximation of high-dimensional ...
07/24/2019

Transport Monte Carlo

In Bayesian inference, transport map is a promising alternative to the c...
11/04/2014

Simple approximate MAP Inference for Dirichlet processes

The Dirichlet process mixture (DPM) is a ubiquitous, flexible Bayesian n...
12/10/2020

Bayesian nonstationary and nonparametric covariance estimation for large spatial data

In spatial statistics, it is often assumed that the spatial field of int...
01/25/2018

A Distributed Framework for the Construction of Transport Maps

The need to reason about uncertainty in large, complex, and multi-modal ...
11/26/2020

A Scalable Partitioned Approach to Model Massive Nonstationary Non-Gaussian Spatial Datasets

Nonstationary non-Gaussian spatial data are common in many disciplines, ...

Code Repositories

BaTraMaSpa

Scalable Bayesian transport maps for high-dimensional non-Gaussian spatial fields


view repo