R语言中sum计算结果显示为1但sum(vector)==1返回FALSE的问题求助

嘿，我遇到一个特别诡异的问题，想跟大家请教下！

我在R里执行sum(c(.7, .2, .1))的时候，结果显示是1，但当我用sum(c(.7, .2, .1)) == 1做相等判断时，返回的居然是FALSE！类似的情况还有sum(c(.69, .2, .11)) == 1，结果也是FALSE，但sum(c(.65, .2, .15)) == 1却能返回TRUE。

我一开始怀疑是不是我机器的问题？我用的是R 4.3.1版本，系统是macOS Ventura 13.5.2，不管在RStudio控制台还是终端里运行R，都会出现这个情况。

问题原因解析 #

其实这不是你的机器问题，也不是R的bug，这是浮点数精度的固有问题！

计算机是用二进制来存储浮点数的，而很多十进制小数（比如0.7、0.2）转换成二进制后是无限循环的小数，没办法被精确存储。当你把这些数相加时，实际得到的是一个非常接近1，但并不是严格等于1的数值——R在显示的时候做了四舍五入，所以你看到的是1，但底层存储的其实是一个和1有微小差异的数。

你可以用这个命令看看真实的数值：

sprintf("%.20f", sum(c(.7, .2, .1)))

执行后你会发现结果是类似0.99999999999999988898这样的数，自然和1做严格相等判断会返回FALSE。而像0.65、0.2、0.15这些数，它们的二进制表示相加后刚好能精确到1，所以判断会返回TRUE。

解决方案 #

在R里比较浮点数是否相等，不要用==做严格判断，推荐用这两种方式：

1. 使用all.equal()函数，它会自动考虑浮点数的精度误差：

all.equal(sum(c(.7, .2, .1)), 1)

这个命令会返回TRUE，如果差异在可接受的精度范围内。
2. 手动设置一个容忍阈值，判断两个数的差值是否小于这个阈值：

abs(sum(c(.7, .2, .1)) - 1) < 1e-10

这里的1e-10是一个很小的数值，你可以根据实际需求调整。

备注：内容来源于stack exchange，提问作者Nic