如何修改R语言summary函数的舍入规则得到正确的均值舍入结果

问题根源 #

输出偏差由两个内置机制叠加导致：

• R内置的数值舍入采用*银行家舍入（四舍六入五成双）*规则，而非日常使用的“逢五进一”四舍五入：即当舍入位后刚好是5时，会取最近的偶数作为结果，因此round(25.625, 2)默认输出25.62而非25.63
• summary函数的数值展示依赖digits参数控制有效数字长度，默认不会输出多余的小数位，因此没有展示完整的25.625精度，直接按舍入规则显示了两位小数的结果

可实现的调整方案 #

1. 临时提升输出精度查看完整数值 #

直接调用summary时调高digits参数，即可显示完整的均值精确值，避免舍入偏差：

summary(c(1, 12,28,30, 34,25,35, 40), digits = 5)
# 输出结果
#    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
#    1.00   21.75   29.00   25.625   34.25   40.00 

注意：digits参数控制的是有效数字数量，而非固定小数位数，输出会根据数值量级自动调整展示格式

2. 自定义常规四舍五入逻辑 #

如果需要严格实现“逢五进一”的两位小数舍入效果（即把25.625转为25.63），可以通过添加极小浮点数偏移的方式，绕过R默认的银行家舍入规则：

# 自定义四舍五入函数
custom_round <- function(x, digits = 2) {
  scale <- 10^digits
  # 加1e-10的极小偏移，避免浮点数精度误差与银行家舍入的五成双规则
  round(x * scale + 1e-10) / scale
}
# 测试效果
custom_round(25.625, 2)
# 输出 25.63

3. 包装自定义summary函数实现目标输出 #

如果需要每次调用时自动输出调整舍入后的结果，可以基于原生summary写一个轻量包装函数，对输出的所有统计值应用自定义舍入规则：

my_summary <- function(x, digits = 2) {
  # 调用原生summary获取基础统计量
  res <- summary(x)
  # 应用自定义四舍五入逻辑
  scale <- 10^digits
  res <- round(res * scale + 1e-10) / scale
  # 按指定精度格式化输出
  print(res, digits = digits + 1)
}
# 测试效果
my_summary(c(1, 12,28,30, 34,25,35, 40))
# 输出结果
#    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
#    1.00   21.75   29.00   25.63   34.25   40.00 

注意事项 #

• 不建议全局修改R的默认舍入规则，大量内置统计函数依赖原生银行家舍入逻辑，全局调整会引发不可预期的计算偏差
• 自定义舍入时添加的1e-10偏移量是为了适配浮点数存储的精度误差，不会对正常量级的数值计算结果产生干扰

内容的提问来源于stack exchange，提问作者Ruser-lab9