[R-br] trabalhando com números pequenos

salah salah3.1416 em gmail.com
Segunda Dezembro 12 13:59:49 BRST 2016


Caro  Luiz


O Cesar Rabak já lhe deu uma excelente dica de leitura
Compreenda qual o limite de representação numérica da sua máquina
Experimente os limites da sua máquina
Estude as bibliotecas gmp e Rmpfr

https://en.wikipedia.org/wiki/Floating_point

options(digits = 22)

.Machine

##
as.numeric("14472334024676221")
# [1] 14472334024676220

## se sua máquina representasse esse número, não seria necessário usar a 
biblioteca Rmpfr
as.numeric("-0.999999999999999999999999998975999999999999914663505297831951000646311908486498")
# [1] -1

## epsilon da maquina
eps = 1
while(eps + 1 > 1){eps = eps/2}
(eps = eps*2)

##
1 - 0.00000000000000001
1 - 0.0000000000000001

p <- mpfr(c(500,600,700)^(-10), precBits=255) - mpfr(1,precBits=255)
class(p)
str(p)

bons estudos!


Mais detalhes sobre a lista de discussão R-br