Eu trabalho assim

#---------------------Analise de Cook's distance -------------------------------------------

install.packages(sfsmisc); library(sfsmisc)

analise<-lm(CONSUMO~factor(GEST)*factor(MANEJO),data=agua)

n<-length(agua$CONSUMO)  # número de observações

n.plot(cooks.distance(analise),seq(1:n),cex=.5,nam=agua$ANIMAL)

criterio<-4/analise$df.residual

abline(v=criterio)

#--------------------------------Eliminacao de outiliers--------------------------------------

ifelse(cooks.distance(analise)>criterio,1,0)->agua$cook

analisecook<-lm(CONSUMO~factor(MANEJO)*factor(GEST),data=agua,subset=(cook==0 & GEST!=0 )) #eliminando outliers.
==================================
Fernando Souza
Zootecnista, DSc. Produção e Alimentação Animal
celular: (31)99796-8781 (Vivo) / (31)97358-4685 (Tim)
e-mail:nandodesouza@gmail.com
Lattes: http://lattes.cnpq.br/6519538815038307
blog: https://producaoanimalcomr.wordpress.com/
=================================

Em Sex, Out 14, 2016 em 3:52 , Walmes Zeviani via R-br <r-br@listas.c3sl.ufpr.br> escreveu:

Além da distância de Cook, você tem mais opções de medidas de influência com a inflence.measures(). Dê uma olhada aqui para ver exemplos http://leg.ufpr.br/~walmes/cursoR/mgest/1medidas-influen.html. Eu gosto de usar o DFits como medida.

À disposição.

Walmes.

​

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.