Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Comme mentionné précédemment, rien ne garantit que le résultat du krigeage d’indicatrices sera une fonction respectant les propriétés d’une fonction de répartititon :

Il faut donc, avant tout calcul, s’assurer que la fonction estimée FKI(x0,z)F_{KI}(x_0, z) respecte ces propriétés, sans quoi des résultats aberrants, comme des probabilités négatives, pourraient apparaître. Pour cela, on effectue des corrections dites “ad hoc”.

La Fig. 1 illustre, de manière exagérée, la violation des propriétés attendues d’une fonction de répartition.

Violation des propriétés d’une fonction de répartition.

Figure 1:Violation des propriétés d’une fonction de répartition.


Correction « ad hoc » de la fonction de répartition estimée

  1. Troncature des bornes :
    Mettre les valeurs négatives de FKI(x,c)F_{KI}(x, c) égales à 0, et celles supérieures à 1 égales à 1.

  2. Correction vers l’avant :
    Soient les seuils ordonnés c1<c2<<cpc_1 < c_2 < \cdots < c_p.
    On définit la fonction corrigée vers l’avant par :

    FKI,avant(x0,c1)=max(0,FKI(x0,c1))F_{KI,\text{avant}}(x_0, c_1) = \max(0, F_{KI}(x_0, c_1))

    Pour i=2,,pi = 2, \ldots, p :

    FKI,avant(x0,ci)=max(FKI,avant(x0,ci1),FKI(x0,ci))F_{KI,\text{avant}}(x_0, c_i) = \max\big(F_{KI,\text{avant}}(x_0, c_{i-1}), F_{KI}(x_0, c_i)\big)
  3. Correction vers l’arrière :
    On définit la fonction corrigée vers l’arrière par :

    FKI,arr(x0,cp)=min(1,FKI(x0,cp))F_{KI,\text{arr}}(x_0, c_p) = \min(1, F_{KI}(x_0, c_p))

    Pour i=p1,,1i = p-1, \ldots, 1 :

    FKI,arr(x0,ci)=min(FKI(x0,ci),FKI,arr(x0,ci+1))F_{KI,\text{arr}}(x_0, c_i) = \min\big(F_{KI}(x_0, c_i), F_{KI,\text{arr}}(x_0, c_{i+1})\big)
  4. Fonction corrigée finale :
    La fonction corrigée est obtenue en faisant la moyenne des deux corrections :

    FKI,corr(x0,ci)=12(FKI,avant(x0,ci)+FKI,arr(x0,ci))F_{KI,\text{corr}}(x_0, c_i) = \frac{1}{2} \big( F_{KI,\text{avant}}(x_0, c_i) + F_{KI,\text{arr}}(x_0, c_i) \big)

Exemple de correction

Le tableau suivant illustre la correction appliquée à une fonction FKIF_{KI} présentant des anomalies :

Seuil ccFKI(x0,c)F_{KI}(x_0,c)FKI,avant(x0,c)F_{KI,\text{avant}}(x_0,c)FKI,arr(x0,c)F_{KI,\text{arr}}(x_0,c)FKI,corr(x0,c)F_{KI,\text{corr}}(x_0,c)
1-0.01 → 0000
20.130.130.130.13
30.240.240.2340.237
40.2380.240.2340.237
50.2340.240.2340.237
60.2370.240.2370.2385
70.530.530.530.53
80.790.790.770.78
90.770.790.770.78
101.02 → 1.01.01.01.0

Exemple visuel

La Fig. 2 illustre l’application des corrections des relations d’ordre présentées dans le tableau précédent. On y observe l’effet de la correction avant et de la correction arrière. Visuellement, la correction préalable consiste à remplacer chaque valeur par le maximum entre celle-ci et la valeur précédente, ce qui garantit une fonction strictement croissante. À l’inverse, la correction arrière part du seuil le plus élevé et remplace chaque valeur par le minimum entre celle-ci et la valeur du seuil situé à un rang supérieur, ce qui produit une fonction strictement décroissante. Les deux corrections définissent une enveloppe dont on prend généralement la moyenne pour obtenir une fonction de répartition corrigée et conforme aux relations d’ordre.

Correction des relations d’ordre de l’exemple du tableau ci-haut.

Figure 2:Correction des relations d’ordre de l’exemple du tableau ci-haut.