Hélène Guis a écrit :Bonjour
Merci pour ces réponses. Si j'ai bien compris:
1) Je peux extraire les résidus de ma régression logistique avec la fonction "residuals". Mais quel type de résidus faut-il choisir? (déviance, pearson, working, response, partial)?
2) Une fois que l'on a extrait les résidus, quel fonction appliquer? Luc Anselin précise dans un de ses docs de travail (Pages 10 et 11 de An intro to spatial regression analysis in R :
www.sal.uiuc.edu/stuff/stuff-sum/pdf/spdepintro.pdf) que ce n'est pas pareil d'appliquer la fonction lm.morantest directement (fonction qui prend comme argument un modèle linaire x.lm) que d'appliquer la fonction moran.test sur des résidus qu'on aurait extrait par ailleurs. Il dit que la 2eme façon de faire (avec moran.test sur des résidus) est incorrecte. Je n'ai pas compris pourquoi cela n'était pas pareil. Pouvez-vous m'éclairer sur ce point?
J'imagine qu'il faut tenir compte du fait que les résidus sont estimés et non observés, ce qui doit réduire le nombre de ddl et diminuer la puissance du test. Juste une hypothèse...
Je ne pense pas qu'il existe de méthode programmée, voire de théorie, pour faire le test de Moran sur les résidus d'un modèle linéaire généralisé. En regardant rapidement sur le net, je suis tombé sur ce
working paper qui semble confirmer ce que je pense et propose une solution pour des modèles log-linéaires (Poisson). Ils utilisent alors soit les résidus de Pearson, soit les résidus de la déviance pour effectuer le test de Moran. Mais c'est un papier de recherche, non publié etc ==> prudence.
Par ailleurs, il y a quelques outils permettant d'évaluer l'autocorrélation ou le variogramme pour des modèles de Poisson (cf dans geoRglm). Tu peux aussi essayer de fouiller dans le package DCluster. Si c'est le cas, tu peux reformuler ta régression logistique en régression de Poisson et utiliser une des méthodes dispo.
Il faudrait en fait poser cette question sur un forum spécialisé en géostats. Voir la "task view" sur les méthodes spatiales avec R à
http://cran.r-project.org/src/contrib/Views/Spatial.html, et la liste R-SIG-Geo
Je ne sais donc plus que faire s'il ne faut pas utiliser moran.test sur des résidus, et si la fonction lm.morantest ne s'applique que pour des modèles linéaires et non logistiques...
Merci encore
Hélène
Désolé de de pas être plus informatif, mais le plus sage est de poser la question à des spécialistes. Je te conseille de passer par la liste R-SIG-Geo.
Merci de nous tenir au courant si tu résouds le problème. L'idéal serait que tu postes alors le code avec un extrait des données dans la section "Fichiers", rubrique "Scripts et fonctions R". Cela pourra alors servir à tout le monde.
Bon courage,
Renaud