r - 为什么 exp(log(6)) == exp(log(3)) * 2 在 R 中返回 FALSE？

Question

在 R 中，为什么对于 y > 2 这是错误的？

y <- c(1, 2, 3, 4, 5)
x <- 2*y
exp(log(x)) == exp(log(y)) * 2
[1]  TRUE  TRUE FALSE FALSE FALSE

Answer 1

数值精度。尝试计算差异：

exp(log(x)) - exp(log(y)) * 2

你可以使用类似的东西：

all.equal( exp(log(x)) , exp(log(y)) * 2 )

Answer 2

确切地说，数值精度是原因。试试这个，更简单的计算

1/2+1/3+1/6 #equal to 1    
(1/2+1/3+1/6)-1 # should be 0
[1] -1.110223e-16