[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