R语言处理无限分数时的数值精度问题：均值偏差求和异常

首先你遇到的不是R的计算异常，而是双精度浮点数的精度限制导致的舍入误差：

• 理论上一组数据的偏差之和确实为0，但R默认使用64位双精度浮点数存储数值，这类格式只能精确表示有限的十进制小数（比如整数、分母是2的幂的分数）。像20/7这种无限循环小数，无法被双精度浮点数精确存储，在计算x - mean(x)时会产生微小的精度损失，最终求和得到的是一个极其接近0的极小值（-4.440892e-16就是科学计数法表示的接近0的数）。

关于你问的核心问题：

• R默认不支持无限长数字的精确计算，但它有专门的扩展包可以实现高精度计算：
  • gmp包：专注于高精度整数运算
  • Rmpfr包：支持任意精度的浮点数运算

比如用Rmpfr重新计算你的例子，就能得到精确的0：

library(Rmpfr)
# 用128位精度定义数值向量
x <- mpfr(c(1,2,2,3,3,4,5), precBits = 128)
sum(x - mean(x))

日常分析中，这种极小的舍入误差完全可以忽略，你可以用all.equal(sum(x - mean(x)), 0)来判断结果是否近似为0，而不用纠结严格等于0。

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