Graphe d'influence des points aberrantes sur les résidus.

Bonjour ,

Je voudrais obtenir le graphe des résidus suivant (fichier ci joint) j'essayé avec cet code mais ça marche pas et j'arrive pas a comprendre l'erreur malgré que R n'affiche aucun message d'erreur !  si quelqu'un parmi vous a une idée merci de me signaler .

Code:

X11()
op <- par(fg='white', bg='black',col='white', col.axis='white', col.lab='white',col.main='white',col.sub='white')
cd <- cooks.distance(Modfinj1)
m <- max(cd)
plot(GLM.1$res ~ GLM.1$fitted.values,cex=1+5*cd/m,col=heat.colors(100)[ceiling(100*cd/m)],pch=16)
abline(h=0,lty=3)
abline(h=c(-2,0,2),lty=c(2,1,2),col="green")
par(op)


Merci beaucoup d'avance.

Bonjour,

Peut être qu'une idée de tes données "Modfinj1" et de ce que tu obtiens serait un plus pour t'aider.

Cordialement,

Khalid

Voilà le Modfinj1 juste un nom pour mon modèle :

Code:

Modfinj1<-lm(DV_max~X1+X2+X3+X4+X5,data=MergeTestP)

.

Ok ... donc je voulais dire peux tu faire un head() ou quelque chose comme ça pour qu'on puisse tester et voir ce que ton code fait?

Idem as tu un png de ce que tu obtiens?

Désolé j'ai pa compris l'utilité de head() , et voilà ce
que j'obtient ci joint.

Le problème est au niveau de cette option " cex=1+5*cd/m,col=heat.colors(100)[ceiling(100*cd/m)]" car quand je le supprime ça marche mais c'est pas ce que cherche comme j'ai déjà montrer au début (l'influence de point):

Code:

X11()
op <- par(fg='white', bg='black',col='white', col.axis='white', col.lab='white',col.main='white',col.sub='white')
cd <- cooks.distance(Modfinj1)
m <- max(cd)
plot(Modfinj1$res ~ Modfinj1$fitted.values,pch=16,ylim=c(-2.5,2.5))
abline(h=0,lty=3)
abline(h=c(-2,0,2),lty=c(2,1,2),col="green")
par(op)


!

Ce code marche très bien chez moi :

Code:

mydata <- data.frame(x=rchisq(100,5),y=rnorm(100))
glm1 <- glm(y~x,data=mydata)
par(fg='white', bg='black',col='white', col.axis='white', col.lab='white',col.main='white',col.sub='white')
cd <- cooks.distance(glm1)
m <- max(cd)
plot(glm1$res ~ glm1$fitted.values,cex=1+5*cd/m,col=heat.colors(100)[ceiling(100*cd/m)],pch=16,ylim=c(-2.5,2.5))
abline(h=0,lty=3)
abline(h=c(-2,0,2),lty=c(2,1,2),col="green")

Essaye d'enlever le ylim ou d'étendre les limites.

Nik

Bon merci beaucoup , enfin ça marche mais j'ai changé le logiciel avec Rstudio !

problème de package peut être .