R语言trunc()函数处理多小数位浮点数结果不一致问题咨询

这不是trunc()函数的逻辑问题，本质是R使用的IEEE 754标准64位双精度浮点数存在精度上限——它只能精确表示大约15-17位十进制有效数字。当你输入带有大量小数位的数值字符串时，R会先将其解析为最接近的双精度浮点数，而非完整保留所有小数位，这就导致了看似相似的输入最终存储的实际数值不同，进而让trunc()返回不同结果。

具体验证 #

你可以用print(x, digits=20)查看这些数值实际在R中的存储值：

# 第一个示例
x1 <- 3.99999999999999977799999999999999999999900
print(x1, digits=20)
# [1] 4

x2 <- 3.999999999999999777999999999999999999999000
print(x2, digits=20)
# [1] 3.9999999999999991118215802998747676610925

可以看到，x1被解析后直接等于4.0，所以trunc(x1)返回4；而x2是一个略小于4.0的浮点数，trunc(x2)自然返回3。

补充说明 #

双精度浮点数的最小可区分间隔（机器精度）可以用.Machine$double.eps查看，大约为2.2e-16。当数值和4.0的差距小于这个间隔的一半时，会被舍入为4.0；反之则保留为小于4的数值。你输入的超长小数数值，刚好处于这个舍入阈值的边界附近，所以细微的字符差异就会导致解析后的数值跨过这个阈值，出现trunc()结果不一致的现象。

内容的提问来源于stack exchange，提问作者Bằng Rikimaru