Skip to article frontmatterSkip to article content

Quand on dispose de nombreuses observations et que l’on souhaite estimer une moyenne sur un grand volume, le calcul direct de la variance d’estimation peut devenir lourd. Une approche simplifiée consiste à utiliser le principe de combinaison des erreurs élémentaires, qui décompose l’estimation globale en une série d’estimations locales approximativement indépendantes.


i. Grille régulière

Pour une grille régulière, l’estimé global est simplement :

Z=1ni=1nZiZ^* = \frac{1}{n} \sum_{i=1}^n Z_i

Chaque point est supposé représenter un bloc de taille égale. Si la variance d’erreur pour chaque bloc est σe2\sigma_e^2, la variance d’estimation globale est :

σe2(globale)=σe2n\sigma_e^2(\text{globale}) = \frac{\sigma_e^2}{n}

Les erreurs sont approximativement indépendantes, car chaque bloc est estimé avec un seul point.


ii. Échantillonnage aléatoire uniforme (stratifié)

Chaque point représente un bloc positionné aléatoirement dans le domaine. Pour un bloc viv_i, la variance d’estimation est :

σei2=E[(Z(x)Zˉ)2]=1vv(Z(x)Zˉ)2dx=D2(v)\sigma_{e_i}^2 = \mathbb{E}[(Z(x) - \bar{Z})^2] = \frac{1}{v} \int_v (Z(x) - \bar{Z})^2 \, dx = D^2(\cdot \mid v)

Donc la variance d’estimation globale est encore :

σe2=D2(v)n\sigma_e^2 = \frac{D^2(\cdot \mid v)}{n}

iii. Échantillonnage quelconque

Supposons que nous somme en mesure de calculer la variance d’estimation sur un sous domaine viv_i. Ainsi, pour l’ensemble du domaine VV, l’estimé est obtenu en combinant, selon des poids proportionnels aux volumes des sous-domaines, les différents estimés obtenus. On aura donc :

Z=i=1nviVZiZ = \sum_{i=1}^n \frac{v_i}{V} Z_i^*

et la variance d’estimation sera, puisque chaque erreur est considérée indépendante :

σe2=i=1n(viV)2σei2\sigma_e^2 = \sum_{i=1}^n \left( \frac{v_i}{V} \right)^2 \sigma_{e_i}^2

où les variances d’estimation élémentaires sont celles correspondant aux différents sous-domaines que l’on a pu reconnaître.


Remarques


Exemple

Une zone a été estimée directement par krigeage (variance d’estimation 0.36) puis par combinaison de 4 parcelles selon 2 scénarios différents (Fig. 1). Pour l’estimation de chaque parcelle, on n’utilise que les points s’y retrouvant. On calcule les variances d’estimation pour chaque parcelle et on combine le tout suivant le carré des surfaces de chaque parcelle. Les résultats obtenus par subdivision, dans les deux cas, sont quasi-identiques au résultat direct.

Illustration des paramètres du variogramme - effet de pépite, palier et portée.

Figure 1:Illustration des paramètres du variogramme - effet de pépite, palier et portée.