Empirical variogram of the underlying Gaussian fields in the plurigaussian models

Nicolas Desassis (1) , D Renard (1) , Hélène Beucher (1)

1 Équipe Géostatistique

The plurigaussian model is particularly suited to describe categorical regional-ized variables. Starting from a simple principle, the thresholding of one or several Gaussian random fields (GRFs) to obtain categories, the plurigaussian model is well adapted for a wide range of situations. By acting on the form of the thresholding rule and/or the threshold values (which can vary along space) and the variograms of the underlying GRFs, one can generate many spatial configurations for the categorical variables. One difficulty arising with the use of this model is to choose variogram model for the underlying GRFs. Indeed, these latter are hidden by the truncation and we only observe the simple and cross-variograms of the category indicators. In this paper, we propose a method based on the pairwise likelihood to estimate the empirical variogram of the GRFs. It provides an exploratory tool in order to choose a suitable model for each GRF and later to estimate its parameters. We illustrate the efficiency of the method with a Monte-Carlo simulation study. The method presented in this paper is implemented in the R package RGeostats.

Mots clés: Plurigaussian models; empirical variography; pairwise likelihood (PL); underlying Gaussian Random Fields (GRFs)

Domaines: Statistiques [math.ST]

Nicolas Desassis, D Renard, Hélène Beucher. Empirical variogram of the underlying Gaussian fields in the plurigaussian models. 2015. ⟨hal-01213962v1⟩