Abstract
We investigate a connection between spatial statistics and statistical physics to obtain new covariance functions with direct physical interpretation for spatial random fields. These covariance functions are based on the exponential Boltzmann-Gibbs representation and use an energy functional to represent interactions between the values of the random field at different points in space. This formulation results in closed-form generalized covariance functions, which display infinite variance in Euclidean spaces of dimension larger than one. We propose regularization schemes in real and reciprocal (spectral) space that lead to well-behaved covariance structures. The real-space regularization parameter allows a continuous interpolation between the Boltzmann-Gibbs covariance and the exponential covariance. We also propose discretized approximations on regular grids, and we show that they represent reparametrized versions of the well-known Besag and Leroux lattice models. We then discuss parameter estimation and spatial prediction for the regularized Boltzmann-Gibbs covariance model in two dimensions. We recommend using the pairwise difference likelihood that combines satisfactory estimation performance and good scalability with many observation points. The predictive performance of the regularized covariance function is assessed by means of cross-validation statistics. Irregularly-spaced samples from the Walker Lake dataset are used, and spatial prediction is conducted by means of ordinary kriging. The regularized Boltzmann-Gibbs covariance yields improved predictive performance compared to the exponential covariance model.
Highlights
The main ideas of statistical field theory have been used to develop Boltzmann-Gibbs random fields with a new class of covariance structure [24, 28, 25], including cross-covariance models for multivariate data [29]
We investigate theoretical properties and statistical techniques for a form of Gaussian Boltzmann-Gibbs models whose energy functional H is related to the spatial variation of the realizations of the field x(s) and of its gradient function ∇x(s)
{X(s, ω) : s ∈ Rd, ω ∈ Ω} is a scalar, real-valued spatial random field defined on a probability space (Ω, F, P ), where Ω is the sample space, F is a σ-field of subspaces of Ω, and P is a probability measure [49]
Summary
We consider spatial models defined over a domain D which we choose as a compact subset of Rd, where s ∈ Rd denotes the location vector. The boundary ∂D of D is supposed to be a piecewise smooth (i.e., infinitely differentiable) manifold of dimension d − 1. {X(s, ω) : s ∈ Rd, ω ∈ Ω} is a scalar, real-valued spatial random field defined on a probability space (Ω, F , P ), where Ω is the sample space, F is a σ-field of subspaces of Ω, and P is a probability measure [49]. The realizations of the random field X(s; ω) for ω ∈ Ω are deterministic functions denoted by x(s). We will consider below mx(s) = 0 unless otherwise specified
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.