Eu tentei rodar, e diz que não existe o arquivo água!

 

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/
=================================