O arquivo água foi o meu banco de dados, vc pode rodar colocando seu próprio dataframe contendo seus dados e fazer as modificações necessárias. Att Em 07/11/2016 21:42, <sznelwar@uol.com.br> escreveu:
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/ =================================