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

by   Matthias Katzfuss, et al.

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.


page 2

page 9

page 14

page 18


Transport map unadjusted Langevin algorithms

Langevin dynamics are widely used in sampling high-dimensional, non-Gaus...

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

Multivariate spatially-oriented data sets are prevalent in the environme...

Limit theorems for semidiscrete optimal transport maps

We study statistical inference for the optimal transport (OT) map (also ...

Simple approximate MAP Inference for Dirichlet processes

The Dirichlet process mixture (DPM) is a ubiquitous, flexible Bayesian n...

Greedy inference with layers of lazy maps

We propose a framework for the greedy approximation of high-dimensional ...

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...

Code Repositories


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

view repo

Please sign up or login with your details

Forgot password? Click here to reset