A pairwise likelihood approach for the empirical estimation of the underlying variograms in the plurigaussian models

Nicolas Desassis (1) , Didier Renard (1) , Hélène Beucher (1) , Sylvain Petiteau (1) , Xavier Freulon (1)

1 Équipe Géostatistique

The plurigaussian model is particularly suited to describe categorical regionalized variables. Starting from a simple principle, the thresh- olding 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 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 semiparametric 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 param- eters. 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: underlying Gaussian Random Fields (GRFs); Plurigaussian models; empirical variography; pairwise likelihood (PL)

Domaines: Statistiques [math.ST]

Nicolas Desassis, Didier Renard, Hélène Beucher, Sylvain Petiteau, Xavier Freulon. A pairwise likelihood approach for the empirical estimation of the underlying variograms in the plurigaussian models. 2015. ⟨hal-01213962v2⟩