<p dir="ltr">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.<br>
Att</p>
<div class="gmail_extra"><br><div class="gmail_quote">Em 07/11/2016 21:42,  <<a href="mailto:sznelwar@uol.com.br">sznelwar@uol.com.br</a>> escreveu:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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 ------------------------------<wbr>-------------</div>
<div>install.packages(sfsmisc); library(sfsmisc)<br>
<div>analise<-lm(CONSUMO~factor(<wbr>GEST)*factor(MANEJO),data=<wbr>agua)<br>
<div>n<-length(agua$CONSUMO)  # número de observações</div>
<div>n.plot(cooks.distance(analise)<wbr>,seq(1:n),cex=.5,nam=agua$<wbr>ANIMAL)</div>
<div>criterio<-4/analise$df.<wbr>residual</div>
<div>abline(v=criterio)</div>
<div>#-----------------------------<wbr>---Eliminacao de outiliers---------------------<wbr>-----------------</div>
<div>ifelse(cooks.distance(analise)<wbr>>criterio,1,0)->agua$cook</div>
<div>analisecook<-lm(CONSUMO~<wbr>factor(MANEJO)*factor(GEST),<wbr>data=agua,subset=(cook==0 & GEST!=0 )) #eliminando outliers.<br>==============================<wbr>====<br>Fernando Souza<br>Zootecnista, DSc. Produção e <wbr>Alimentação Animal<br>celular: (31)99796-8781 (Vivo)<wbr> / (31)97358-4685 (Tim)<br><a href="mailto:e-mail%3Anandodesouza@gmail.com" target="_blank">e-mail:nandodesouza@gmail.com</a><br>Lattes: <a href="http://lattes.cnpq.br/6519538815038307" target="_blank">http://lattes.cnpq.br/<wbr>6519538815038307</a><br>blog: <a href="https://producaoanimalcomr.wordpress.com/" target="_blank">https://<wbr>producaoanimalcomr.wordpress.<wbr>com/</a><br>==============================<wbr>=== <br><br></div>
</div>
</div></blockquote></div></div>