Kriging Beyond Gaussianity: Fixed-Domain Asymptotic Optimality
Abstract
Under squared-error loss, the optimal predictor of a non-Gaussian random field is its conditional mean, which may be difficult to compute. We consider a broad class of non-Gaussian random fields obtained by transforming finitely many independent Gaussian fields sharing a common correlation function, and study when kriging entails no first-order loss in mean squared prediction error (MSPE) under fixed-domain asymptotics. We compare the observed-data conditional mean, kriging with the observed covariance, kriging with the latent correlation, and a latent-information oracle. These predictors satisfy a finite-sample MSPE ordering. Under conditions linking transformation regularity to the high-frequency decay of the latent spectral density, all four MSPEs are asymptotically equivalent. For radially nonincreasing spectral densities with polynomial tails, the required regularity order is independent of both the spatial dimension and the number of latent fields. Exponential correlation also yields finite-sample bounds. Examples include Weibull, skew-Gaussian, and copula-based fields, while numerical comparisons illustrate finite-sample differences. Counterexamples show why covariance structure and first-order transformation regularity alone are insufficient. Finally, standardized prediction errors converge to a Gaussian scale mixture, which is standard normal exactly when the transformation's gradient norm is almost surely constant.
Related
- On the Information Bottleneck for Ergodic Finite-State Markov Chains: Stationary Approximation and Asymptotic Optimality
- Translate, then decide: asymptotic optimality of latent-space control in bistable systems
- Exact Bayes Regret and Asymptotic Optimality in High-Dimensional Gaussian Bandits
- Finite-sample nonparametric mean tests: Leave-one-out duality and asymptotic optimality
- Pollak's Minimax Quickest Change Detection: Non-Asymptotic Optimality
BibTeX
@misc{bevilacqua2026kriging,
title = {{Kriging Beyond Gaussianity: Fixed-Domain Asymptotic Optimality}},
author = {Moreno Bevilacqua and Hector Araya},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.24115},
eprint = {2610.24115},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.24115},
}
RIS
TY - UNPB
TI - Kriging Beyond Gaussianity: Fixed-Domain Asymptotic Optimality
AU - Moreno Bevilacqua
AU - Hector Araya
PY - 2026
DA - 2026/10/16
PB - arXiv
DO - 10.48550/arxiv.2610.24115
UR - https://arxiv.org/abs/2610.24115
AB - Under squared-error loss, the optimal predictor of a non-Gaussian random field is its conditional mean, which may be difficult to compute. We consider a broad class of non-Gaussian random fields obtained by transforming finitely many independent Gaussian fields sharing a common correlation function, and study when kriging entails no first-order loss in mean squared prediction error (MSPE) under fixed-domain asymptotics. We compare the observed-data conditional mean, kriging with the observed covariance, kriging with the latent correlation, and a latent-information oracle. These predictors satisfy a finite-sample MSPE ordering. Under conditions linking transformation regularity to the high-frequency decay of the latent spectral density, all four MSPEs are asymptotically equivalent. For radially nonincreasing spectral densities with polynomial tails, the required regularity order is independent of both the spatial dimension and the number of latent fields. Exponential correlation also yields finite-sample bounds. Examples include Weibull, skew-Gaussian, and copula-based fields, while numerical comparisons illustrate finite-sample differences. Counterexamples show why covariance structure and first-order transformation regularity alone are insufficient. Finally, standardized prediction errors converge to a Gaussian scale mixture, which is standard normal exactly when the transformation's gradient norm is almost surely constant.
ER -