É possível aumentar a precisão numerica do R?
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? Se ninguém souber como fazer no R, mas souber de algum pacote que posso utilizar em C++ também seria muito útil. Obrigada, Carolina
Qual a incosistência no resultado? formatC((1.5692524257254742661/(10^41)) + 6.125125476214521152, digits = 50) 2011/4/20 Carolina Pena <carolinapenaa@yahoo.com.br>:
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?
Se ninguém souber como fazer no R, mas souber de algum pacote que posso utilizar em C++ também seria muito útil.
Obrigada,
Carolina _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O
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
participantes (3)
-
Bernardo Rangel Tura -
Carolina Pena -
Henrique Dallazuanna