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

research
02/14/2023

Transport map unadjusted Langevin algorithms

Langevin dynamics are widely used in sampling high-dimensional, non-Gaus...
research
05/31/2020

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

Multivariate spatially-oriented data sets are prevalent in the environme...
research
03/17/2023

Limit theorems for semidiscrete optimal transport maps

We study statistical inference for the optimal transport (OT) map (also ...
research
11/04/2014

Simple approximate MAP Inference for Dirichlet processes

The Dirichlet process mixture (DPM) is a ubiquitous, flexible Bayesian n...
research
05/31/2019

Greedy inference with layers of lazy maps

We propose a framework for the greedy approximation of high-dimensional ...
research
06/16/2022

Bayesian conjugacy in probit, tobit, multinomial probit and extensions: A review and new results

A broad class of models that routinely appear in several fields can be e...
research
07/30/2020

A Vecchia Approximation for High-Dimensional Gaussian Cumulative Distribution Functions Arising from Spatial Data

We introduce an approach to quickly and accurately approximate the cumul...

Please sign up or login with your details

Forgot password? Click here to reset