Só uma dúvida, antes da eliminação de tais outlies, não seria conveniente verificar se são pontos influentes? Então, como verificaria se são influentes?<br /><em><br />Att.<br />André</em><br /> <br /><br />
<hr style="border-top: 1px solid #ccc;" />
Em 07/08/2012 14:10, <strong>Paulo Justiniano < paulojus@leg.ufpr.br ></strong> escreveu:<br />Os algoritmos (padrao) de quantile() e fivenum() são diferentes paera <br />calculos de quantis.<br />salvo engano<br />A 1a supoe correcoes com interpolacoes lineares entre os valores, em <br />quanto que a 2a trabalha mesmo com os dados apenas<br /><br /><br />On Tue, 7 Aug 2012, Ivan Bezerra Allaman wrote:<br /><br />> Oi Walmes!<br />> <br />> Exatamente! Já haviam me alertando aqui no grupo quanto a isso. O boxplot por default usa 1.5, mais eu sempre altero para 2.5, e<br />> esqueci de alterar no código que enviei. O valor de 2.5 foi um concenso que cheguei junto com o Faria, achamos mais flexível. <br />> <br />> (S,f,P)<br />> Allaman<br />> <br />> <br />> <br />> <br />><br />> \begin{signature}<br />><br />>
<<>>=<br />><br />> Prof. Dr. Ivan Bezerra Allaman<br />><br />> Universidade Estadual de Santa Cruz<br />><br />> Departamento de Ciências Exatas e Tecnológicas<br />><br />> Ilhéus/BA - Brasil<br />><br />> Fone: +55 73 3680-5596<br />><br />> E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com<br />><br />> @<br />><br />> \end{signature}<br />> <br />> <br />> __________________________________________________________________________________________________________
_____________________________<br />> De: Walmes Zeviani <br />> Para: r-br@listas.c3sl.ufpr.br; Ivan Bezerra Allaman <br />> Enviadas: Terça-feira, 7 de Agosto de 2012 11:02<br />> Assunto: Re: [R-br] Eliminar outliers de uma planilha.<br />> <br />> Ivan,<br />> <br />> Ao contrario do que eu também pensava, a boxplot() não usa a quantile() para encontrar os quantis, ela usa a fivenum(). Para grandes<br />> amostras não dá diferença mas para amostras pequenas dá. Você usou 2.5 mas o boxplot() usa 1.5. Eu fiz uma função também, veja<br />> <br />> set.seed(12345)<br />> y <- rbeta(100, 1, 5)<br />> boxplot(y)<br />> abline(h=fivenum(y))<br />> x <- gl(4, 25)<br />> bp <- boxplot(y~x)<br />> bp<br />> bp$out<br />> bp$group<br />> <br />> f <- function(x){<br />> qrtl <- fivenum(x)[c(2,4)]<br />> iqr <- diff(qrtl)<br />> limt <- qrtl+c(-1,1)*1.5*iqr<br />> fi
ndInterval(x, limt)-1<br />> }<br />> <br />> f(y)<br />> <br />> tapply(y, x, f)<br />> <br />> À disposição.<br />> Walmes.<br />> <br />> ==========================================================================<br />> Walmes Marques Zeviani<br />> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)<br />> Departamento de Estatística - Universidade Federal do Paraná<br />> fone: (+55) 41 3361 3573<br />> VoIP: (3361 3600) 1053 1173<br />> e-mail: walmes@ufpr.br<br />> twitter: @walmeszeviani<br />> homepage: <a href="http://www.leg.ufpr.br/~walmes" target="_blank">http://www.leg.ufpr.br/~walmes</a><br />> linux user number: 531218<br />> ==========================================================================<br />> <br />> <br />> <br />><br /><br />