[Résolu]RMSE et IDW (Pondération Inverse à la Distance)

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Tibo Neudin
Messages : 11
Enregistré le : 01 Avr 2016, 15:25

[Résolu]RMSE et IDW (Pondération Inverse à la Distance)

Messagepar Tibo Neudin » 09 Juin 2018, 12:48

Bonjour à tous,

J'ai une question autour du choix de paramètre de puissance sur une IDW. Précisons que j'ai données très simple, un comptage et des données GPS (qui je transforme en Lambert 93).

En effet c'est normalement assez subjectif comme choix. Je ne sais donc que rarement quoi prendre comme paramètre.

J'ai donc voulu faire une IDW avec un paramètre de puissance à 2 (le standard si je ne dis pas de bêtise dans le package gstat) qui me donne un résultat acceptable (comme map j'entennds).

Ensuite je fais cherche le RMSE:

Code : Tout sélectionner

IDW.out <- vector(length = length(Carto_L93))
for (i in 1:length(Carto_L93)) {
    IDW.out[i] <- idw(Count ~ 1, Carto_L93[-i,], Carto_L93[i,], idp=2.0)$var1.pred
}
###Plot
OP <- par(pty="s", mar=c(4,3,0,0))
plot(IDW.out ~ Carto_L93$Count, asp=1, xlab="Observed", ylab="Count", pch=16,
     col=rgb(0,0,0,0.5))
abline(lm(IDW.out ~ Carto_L93$Count), col="red", lw=2,lty=2)
abline(0,1)
###RMSE: 318.8981
sqrt( sum((IDW.out - Carto_L93$Count)^2) / length(Carto_L93))
par(OP)


Ce qui me donne donc un RMSE de 318.8981...
Lorsque je refais une formule avec ce nouveau paramètre de puissance. Je me retrouve avec une carte dénuée d'intérêts.
Voila ma formule:

Code : Tout sélectionner

lzn.IDW <- idw(log(Count)~1,Carto_L93, grille_L93, idp=318.8981)

A votre avis. Est-ce que j'ai fait une erreur?
Que choisir comme paramètre de puissance?
Est-ce que mes données ne sont en faite pas super adaptées?

Par avance merci.
Bonne weekend à tous.
Tibo

Tibo Neudin
Messages : 11
Enregistré le : 01 Avr 2016, 15:25

Re: RMSE et IDW (Pondération Inverse à la Distance)

Messagepar Tibo Neudin » 13 Juin 2018, 17:17

Hello,

Je vais me répondre tout seul.
Au cas où cette info peut aider une autre personne.

Après avoir discuté avec un professionnel, il semble que la plupart du temps, on ne prenne pas en considération se paramètre.
La plupart du temps, on prend 2 (la valeur par défaut).

A titre d'info, la méthode est bonne, et reste applicable. Mais si le RMSE sort quelque chose d'étrange, on se satisfera du paramètre par défaut.

J'espère avoir aider quelqu'un :)
Cordialement,
Tibo


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité