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