La méthode consiste à combiner Régression Linéaire Multiple (RLM) et krigeage. Les RLM, indépendemment pour chacune des 21 années et chacun des 4 prédictands T7, Pan, MJJ et EQ, sont calculées en fonction de trois prédicteurs: altitude, latitude et longitude. À titre d’exemple, une illustration est fournie pour EQ 1991 (fig. a). L’apprentissage se fait sur k-1 stations météorologiques. Une validation croisée est effectuée en cachant de manière itérative chacune des k stations. Les Root-Mean-Square Errors (RMSE), ou racines carrées des erreurs moyennes, vont de 0,4 à 0,8°C pour les températures de juillet (T7), de 93 à 157 mm pour les précipitations annuelles (Pan), de 0,25 à 0,53 mm pour MJJ, et enfin de 2,3 à 4,6°C·mm-1 pour EQ. Les coefficients de détermination (R2) vont de 0,3 à 0,7 selon les années et les indices. Les résidus des RLM sont ensuite krigés (fig. b). Les valeurs positives (négatives) indiquent que la régression a sous(sur)-estimé la valeur climatique. Enfin, la somme des estimations issues de la RLM et des résidus krigés est calculée (fig. c).