[R-br] É possível aumentar a precisão numerica do R?
Bernardo Rangel Tura
tura em centroin.com.br
Quinta Abril 21 06:51:33 BRT 2011
On Wed, 2011-04-20 at 12:35 -0700, Carolina Pena wrote:
> Boa tarde pessoal,
>
> Estou trabalhando com operações que exigem uma precisão numerica muito
> grande.
>
> Por exemplo: (1.5692524257254742661/(10^41)) +
> 6.125125476214521152)
>
> No entanto, o R está gerando resultados inconsistentes. Alguém sabe se
> é possível aumentar a precisão do R?
Carolina,
O que você define como inconsistente?
Observe o código abaixo
rep( (1.5692524257254742661/(10^41) + 6.125125476214521152),32)
[1] 6.125125 6.125125 6.125125 6.125125 6.125125 6.125125 6.125125
6.125125
[9] 6.125125 6.125125 6.125125 6.125125 6.125125 6.125125 6.125125
6.125125
17] 6.125125 6.125125 6.125125 6.125125 6.125125 6.125125 6.125125
6.125125
[25] 6.125125 6.125125 6.125125 6.125125 6.125125 6.125125 6.125125
6.125125
Em todo os casos deu o mesmo valor 6.125125!
Além disso você não deve confundir a representação como armazenamento
observe
> formatC(1.5692524257254742661/(10^41) +
6.125125476214521152,digits=50)
[1] " 6.1251254762145208587753586471080780029296875"
--
[]s
Tura
Mais detalhes sobre a lista de discussão R-br