R语言中根据标准差小数位数对均值进行取整

需求：计算数据框指定列的标准差与均值，先将标准差取整至1位有效数字并去除末尾零，再将均值取整至与该标准差相同的小数位数。

现有代码：

library("tidyverse")

# 创建示例数据
data <- data.frame(name = c("A", "B", "C"), 
                   value = c(0.1, 0.11, 0.111,
                             5.0, 0.0003, 0.00002,
                             0.5, 0.13, 0.113))

# 计算标准差（保留1位有效数字）、均值，去重
data <- data %>%
  mutate(sd = signif(sd(value), 1), mean = mean(value), .by = name) %>%
  select(- value) %>%
  distinct()

# 去除标准差末尾零
data$sd <- data$sd %>% 
  str_remove("0+$") 

解决方案 #

核心是先确定每个标准差对应的小数位数，再用该位数对均值取整。以下提供两种可靠实现方式：

方法1：基于数值计算确定小数位数（推荐） #

通过对数运算直接计算有效数字的位置，避免字符串处理的格式问题：

library("tidyverse")

data <- data.frame(name = c("A", "B", "C"), 
                   value = c(0.1, 0.11, 0.111,
                             5.0, 0.0003, 0.00002,
                             0.5, 0.13, 0.113))

data_processed <- data %>%
  mutate(
    # 计算保留1位有效数字的标准差
    sd_raw = signif(sd(value), 1),
    # 计算原始均值
    mean = mean(value),
    # 推导标准差需要的小数位数：大于等于1的数小数位为0，小于1的数取-log10后的整数部分
    decimal_places = ifelse(sd_raw >= 1, 0, -floor(log10(sd_raw)))
  ) %>%
  # 按name去重，保留唯一分组结果
  distinct(name, .keep_all = TRUE) %>%
  # 对均值取整，同时格式化标准差去除末尾零和多余小数点
  mutate(
    mean_rounded = round(mean, decimal_places),
    sd = format(sd_raw, nsmall = decimal_places) %>% 
      str_remove("0+$") %>% 
      str_remove("\\.$") # 处理类似"3."的情况
  ) %>%
  select(name, sd, mean_rounded)

print(data_processed)

方法2：基于字符串提取小数位数 #

如果偏好字符串处理，可直接从格式化后的标准差中提取小数部分长度：

library("tidyverse")

data <- data.frame(name = c("A", "B", "C"), 
                   value = c(0.1, 0.11, 0.111,
                             5.0, 0.0003, 0.00002,
                             0.5, 0.13, 0.113))

data_processed <- data %>%
  mutate(sd = signif(sd(value), 1), mean = mean(value), .by = name) %>%
  select(-value) %>%
  distinct() %>%
  mutate(
    # 将标准差转为字符串，拆分整数与小数部分
    sd_str = as.character(sd),
    # 提取小数位数：有小数点则取小数部分长度，否则为0
    decimal_places = case_when(
      str_detect(sd_str, "\\.") ~ nchar(str_split(sd_str, "\\.")[[1]][2]),
      TRUE ~ 0
    ),
    # 按小数位数取整均值
    mean_rounded = round(mean, decimal_places),
    # 清理标准差格式：去除末尾零和多余小数点
    sd = sd_str %>% str_remove("0+$") %>% str_remove("\\.$")
  ) %>%
  select(name, sd, mean_rounded)

print(data_processed)

输出示例 #

两种方法都会得到类似如下的结果：

name    sd mean_rounded
1    A   0.1          0.1
2    B     3            2
3    C   0.2          0.2

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