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

Question

在R中，为什么y>2为false？

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