<div><span style="font-size: medium; color: #0000ff;">Eu tentei rodar, e diz que não existe o arquivo água!</span></div>
<div> </div>
<div>Eu trabalho assim</div>
<div>#---------------------Analise de Cook's distance -------------------------------------------</div>
<div>install.packages(sfsmisc); library(sfsmisc)<br />
<div>analise<-lm(CONSUMO~factor(GEST)*factor(MANEJO),data=agua)<br />
<div>n<-length(agua$CONSUMO)  # número de observações</div>
<div>n.plot(cooks.distance(analise),seq(1:n),cex=.5,nam=agua$ANIMAL)</div>
<div>criterio<-4/analise$df.residual</div>
<div>abline(v=criterio)</div>
<div>#--------------------------------Eliminacao de outiliers--------------------------------------</div>
<div>ifelse(cooks.distance(analise)>criterio,1,0)->agua$cook</div>
<div>analisecook<-lm(CONSUMO~factor(MANEJO)*factor(GEST),data=agua,subset=(cook==0 & GEST!=0 )) #eliminando outliers.<br />==================================<br />Fernando Souza<br />Zootecnista, DSc. Produção e Alimentação Animal<br />celular: (31)99796-8781 (Vivo) / (31)97358-4685 (Tim)<br />e-mail:nandodesouza@gmail.com<br />Lattes: <a href="http://lattes.cnpq.br/6519538815038307" target="_blank">http://lattes.cnpq.br/6519538815038307</a><br />blog: <a href="https://producaoanimalcomr.wordpress.com/" target="_blank">https://producaoanimalcomr.wordpress.com/</a><br />================================= <br /><br /></div>
</div>
</div>