Skip to Main content Skip to Navigation
Journal articles

Efficiently estimating some common geostatistical models by ‘energy–variance matching’ or its randomized ‘conditional–mean’ versions

Didier A. Girard 1
1 IPS - Inférence Processus Stochastiques
LJK - Laboratoire Jean Kuntzmann
Abstract : We consider the problem of fitting an isotropic zero-mean stationary Gaussian field model to (possibly noisy) observations, when the model belongs to the Matérn family with known regularity index ν > 0 , or to the spherical family. For estimating the correlation range (also called “decorrelation length”) and the variance of the field, two simple estimating functions based on the so-called “conditional Gaussian Gibbs-energy mean” (CGEM) and the empirical variance (EV) were recently introduced. This article presents an extensive Monte Carlo simulation study for problems with around a thousand observations and settings including large, moderate, and even “small”, correlation ranges. The known observation sites are either on a 2D grid (including a case of “very non-uniform” grid spacings) or randomly uniformly distributed on a simple 2D region. Some experiments for a grid with missing values are also analyzed.
This study empirically demonstrates that, for all the (possibly random) uniform designs, the statistical efficiency CGEM−EV compared to exact maximum likelihood (ML) is globally very satisfactory (except a degradation for the very extremal ranges in some contexts) provided the signal-to-noise ratio (SNR) is strong enough and ν too large, this SNR restriction being alleviated as ν decreases. For the “very non-uniform” design, a simple weighting of EV restores this efficiency. In the less favorable cases, the statistical loss remains in fact acceptable: e.g. for the largest considered index (ν = 3/2 ) and a “not strong enough” SNR, it may happen (in fact only for large ranges) that CGEM-EV almost doubles the mean squared error for the range parameter or for the widely used combination of the two parameters known as microergodic-parameter. Furthermore an important conclusion for computational efficiency is that the use of the natural fast randomized-trace version of CGEM-EV does not significantly degrade this statistical efficiency.
Complete list of metadatas

Cited literature [29 references]  Display  Hide  Download
Contributor : Didier A. Girard <>
Submitted on : Friday, December 14, 2018 - 3:09:39 PM
Last modification on : Friday, November 20, 2020 - 10:52:01 AM
Long-term archiving on: : Friday, March 15, 2019 - 4:15:42 PM


Files produced by the author(s)




Didier A. Girard. Efficiently estimating some common geostatistical models by ‘energy–variance matching’ or its randomized ‘conditional–mean’ versions. Spatial Statistics, Elsevier, 2017, 21, Part A, pp.1-26. ⟨10.1016/j.spasta.2017.01.001⟩. ⟨hal-00515832v4⟩



Record views


Files downloads