Implémentation en R d'une alternative ("CGEM-EV") au maximum de vraisemblance pour des champs de type Matérn

Après une introduction à l'alternative CGEM-EV qui permet d'estimer les paramètres de champs aléatoires Gaussiens de type Matérn, nous présentons son implémentation en R et commentons l'exécution du script ScriptForDemo.R qui exploite 2 fichiers source R annexés. Cette implémentation résout les systèmes linéaires invoqués, par une classique méthode itérative du gradient conjugué associée à un produit de matrice-vecteur par FFT. Elle permet donc de refaire des expérimentations comme celles de hal.archives-ouvertes.fr/hal-00515832/en/, mais pour des données de plus grande taille.