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