Acho que não fui claro o suficiente: os desvios deveriam, de fato, ser iguais (os números estão à mesma distância de suas respectivas médias), porém o R não retrata assim (independentemente da quantidade de dígitos) e eu fiquei surpreso.<br>
<br>Na minha época, sei que havia coisinhas que evitávamos, como elevar números muito grandes ao quadrado para somá-los e depois dividi-los, quando era possível contornar o problema aplicando a divisão previamente, mas nas máquinas de hoje eu achei que não era necessário preocupar-se com tais acontecimentos. <br>
<br>Alguma outra pista?<br><br>C.<br><br><div class="gmail_quote">Em 14 de abril de 2011 17:13, Gustavo Henrique de Carvalho <span dir="ltr"><<a href="mailto:gustavo.bio@gmail.com">gustavo.bio@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><a href="http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f" target="_blank">http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f</a><br>

<br>
> print(y$dp, digits = 13)<br>
[1] 0.005773502691896 0.005773502691897<br>
<br>
2011/4/14 Benilton Carvalho <<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>>:<br>
<div><div></div><div class="h5">> Cezar,<br>
><br>
> o que vc esta' procurando e' o seguinte:<br>
><br>
> all.equal(y$dp - refer, rep(0, 2))<br>
><br>
> b<br>
><br>
> 2011/4/14 Cézar Freitas <<a href="mailto:cezarcamelo@gmail.com">cezarcamelo@gmail.com</a>>:<br>
>> Olá.<br>
>><br>
>> Bem, escrevo apenas para ilustrar algo que me surpreendeu. Estava plotando<br>
>> gráficos e os colorindo segundo seu desvio padrão estava ou não acima de um<br>
>> ponto de corte. Reproduzo abaixo a parte do código que interessa:<br>
>><br>
>> y=as.data.frame(matrix(c(4.96,5.02,4.97,5.03,4.96,5.03), ncol=3))<br>
>> y$dp=apply(y,1,sd)<br>
>><br>
>>> y<br>
>>     V1   V2   V3          dp<br>
>> 1 4.96 4.97 4.96 0.005773503<br>
>> 2 5.02 5.03 5.03 0.005773503<br>
>><br>
>> refer=y$dp[1]<br>
>><br>
>> y$dp<=refer<br>
>><br>
>> [1]  TRUE FALSE<br>
>><br>
>> Sei que o R é extremamente preciso e até podemos alterar sua mantissa nos<br>
>> cálculos, mas isso não é rotineiro para cálculos tão simples. O pessoal do<br>
>> cálculo numérico saberia dizer o que evitar para não cair nessas armadilhas<br>
>> no futuro?<br>
>><br>
>> Abraços,<br>
>> Cézar Freitas<br>
>><br>
>> _______________________________________________<br>
>> R-br mailing list<br>
>> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
>> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
>><br>
>><br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
><br>
</div></div></blockquote></div><br>