# Bayesian nonstationary and nonparametric covariance estimation for large spatial data

arXiv: https://arxiv.org/abs/2012.05967 · HTML: https://mlsanity.com/paper/2012.05967/

## Key facts

- Cited by 1 papers in the map (data through 14 Oct 2026).
- First posted on arXiv 10 Dec 2020.

## Built on by (most cited first, 1 of 1)

- [Learning non-Gaussian spatial distributions via Bayesian transport maps with parametric shrinkage](https://mlsanity.com/paper/2409.19208/) (2024)

## Abstract (arXiv)

In spatial statistics, it is often assumed that the spatial field of interest is stationary and its covariance has a simple parametric form, but these assumptions are not appropriate in many applications. Given replicate observations of a Gaussian spatial field, we propose nonstationary and nonparametric Bayesian inference on the spatial dependence. Instead of estimating the quadratic (in the number of spatial locations) entries of the covariance matrix, the idea is to infer a near-linear number of nonzero entries in a sparse Cholesky factor of the precision matrix. Our prior assumptions are motivated by recent results on the exponential decay of the entries of this Cholesky factor for Matern-type covariances under a specific ordering scheme. Our methods are highly scalable and parallelizable. We conduct numerical comparisons and apply our methodology to climate-model output, enabling statistical emulation of an expensive physical model.