Calibrated Bayesian inference for random fields on large irregular domains using the debiased spatial Whittle likelihood

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Goodwin, Thomas, Guillaumin, Arthur, Quiroz, Matias, Villani, Mattias, Kohn, Robert
Natura: Preprint
Pubblicazione: 2025
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866916766291066880
author Goodwin, Thomas
Guillaumin, Arthur
Quiroz, Matias
Villani, Mattias
Kohn, Robert
author_facet Goodwin, Thomas
Guillaumin, Arthur
Quiroz, Matias
Villani, Mattias
Kohn, Robert
contents Bayesian inference for stationary random fields is computationally demanding. Whittle-type likelihoods in the frequency domain based on the fast Fourier Transform (FFT) have several appealing features: i) low computational complexity of only $\mathcal{O}(n \log n)$, where $n$ is the number of spatial locations, ii) robustness to assumptions of the data-generating process, iii) ability to handle missing data and irregularly spaced domains, and iv) flexibility in modelling the covariance function via the spectral density directly in the spectral domain. It is well known, however, that the Whittle likelihood suffers from bias and low efficiency for spatial data. The debiased Whittle likelihood is a recently proposed alternative with better frequentist properties. We propose a methodology for Bayesian inference for stationary random fields using the debiased spatial Whittle likelihood, with an adjustment from the composite likelihood literature. The adjustment is shown to give a well-calibrated Bayesian posterior as measured by coverage properties of credible sets, without sacrificing the quasi-linear computation time. We apply the method to simulated data and two real datasets.
format Preprint
id arxiv_https___arxiv_org_abs_2505_23330
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Calibrated Bayesian inference for random fields on large irregular domains using the debiased spatial Whittle likelihood
Goodwin, Thomas
Guillaumin, Arthur
Quiroz, Matias
Villani, Mattias
Kohn, Robert
Methodology
Bayesian inference for stationary random fields is computationally demanding. Whittle-type likelihoods in the frequency domain based on the fast Fourier Transform (FFT) have several appealing features: i) low computational complexity of only $\mathcal{O}(n \log n)$, where $n$ is the number of spatial locations, ii) robustness to assumptions of the data-generating process, iii) ability to handle missing data and irregularly spaced domains, and iv) flexibility in modelling the covariance function via the spectral density directly in the spectral domain. It is well known, however, that the Whittle likelihood suffers from bias and low efficiency for spatial data. The debiased Whittle likelihood is a recently proposed alternative with better frequentist properties. We propose a methodology for Bayesian inference for stationary random fields using the debiased spatial Whittle likelihood, with an adjustment from the composite likelihood literature. The adjustment is shown to give a well-calibrated Bayesian posterior as measured by coverage properties of credible sets, without sacrificing the quasi-linear computation time. We apply the method to simulated data and two real datasets.
title Calibrated Bayesian inference for random fields on large irregular domains using the debiased spatial Whittle likelihood
topic Methodology
url https://arxiv.org/abs/2505.23330