R语言中为何两仅首数不同的数值表达式判断为相等？

这不是溢出错误，而是双精度浮点数的精度丢失问题，核心原因如下：

• R中的numeric类型本质是IEEE 754标准的双精度浮点数，这类浮点数仅用53位二进制位存储有效数字，换算为十进制约为15-17位有效数字。
• 你用到的大数76561197960265728本身已经是17位十进制数，接近双精度浮点数的有效位数上限。当它加上13532728或13532740时，两个加数的差异（12）相对于大数的量级来说，已经小到超出双精度浮点数的分辨能力。
• 双精度浮点数在这个量级上的**最小可分辨单位（ULP）**远大于12。也就是说，这两个相加结果会被舍入到同一个可精确表示的双精度数值，因此等式判断结果为TRUE。

可以用R代码验证这个精度限制：

# 查看该大数的下一个可表示的双精度数
next_val <- nextafter(76561197960265728, Inf)
# 计算最小可分辨单位（ULP）
ulp <- next_val - 76561197960265728
print(ulp) # 结果远大于12，说明无法区分相差12的两个数

简单总结：双精度浮点数能表示的数值范围极大（对应.Machine$double.xmax），但有效数字的位数是有限的，当数值大到一定程度，小的数值差异会被精度限制“抹平”。

内容的提问来源于stack exchange，提问作者Gustavo Gomes