[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