为何R语言signif函数对1.275取3位有效数字未向上取整？

核心原因是二进制浮点数的精度限制：十进制的1.275无法被二进制浮点数精确存储，R中实际存储的是一个略小于1.275的近似值（执行sprintf("%.20f", 1.275)可以看到真实值是1.27499999999999982236）。signif函数基于这个近似值做舍入时，会判定它更接近1.27，因此返回1.27而非1.28。

R内置的round函数默认采用银行家舍入（四舍六入五成双）规则，但它是按指定小数位数舍入的。如果需要针对有效数字实现该规则，可以通过以下方式：

方法1：自定义函数（基础版） #

先计算对应有效数字所需的小数位数，再用round处理，最后用signif规范格式：

bankers_signif <- function(x, digits = 3) {
  if (x == 0) return(0)
  # 计算需要保留的小数位数
  decimals <- digits - 1 - floor(log10(abs(x)))
  # 用round执行银行家舍入
  rounded_val <- round(x, digits = decimals)
  # 确保有效数字位数，去除末尾冗余零
  signif(rounded_val, digits = digits)
}

测试示例：

• bankers_signif(1.265, 3) → 返回1.26（符合你的要求）
• bankers_signif(1.275, 3) → 返回1.28

方法2：精确十进制处理（避免浮点数误差） #

如果要彻底规避浮点数近似带来的影响，可以使用decimal包处理精确十进制数值：

# 先安装包：install.packages("decimal")
library(decimal)

bankers_signif_exact <- function(x, digits = 3) {
  x_dec <- decimal(x)
  exp <- floor(log10(abs(x_dec)))
  decimals <- digits - 1 - exp
  # 精确执行银行家舍入
  rounded_dec <- round(x_dec, digits = decimals)
  # 转换回数值并规范有效数字
  as.numeric(signif(rounded_dec, digits = digits))
}

测试示例：

• bankers_signif_exact(1.275, 3) → 返回1.28
• bankers_signif_exact(1.265, 3) → 返回1.26

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