R语言mutate调用bc.mean.sd逐行计算均值触发告警问题

• dplyr::mutate的默认传参逻辑是将整列向量作为整体传入调用函数，而非逐行传入单个标量值。
• 你之前使用的esc::esc_mean_sd原生支持向量化运算：当传入长度为n的向量参数时，函数会自动逐组匹配入参、逐行计算，返回长度为n的结果列表，此时用$提取的结果长度和原数据行数完全匹配，因此可以正常生成新列。
• estmeansd::bc.mean.sd不支持向量化输入：函数设计上仅接收长度为1的标量作为q1.val、med.val、q3.val、n的入参。当传入长度大于1的整列向量时，函数内部的参数校验、计算逻辑只会提取每个参数向量的第一个元素完成运算，最终仅返回1组计算结果。你看到的警告the condition has length > 1 and only the first element will be used，就是函数内部做参数合法性判断时，传入的判断条件为多元素向量，R自动取第一个元素判定触发的；最终单个计算结果被mutate循环填充到所有行，就出现了所有行返回值和第一行完全一致的问题。

以下两种写法都可以实现逐行计算的需求，可根据自己的编码习惯选择：

方法1：使用rowwise()开启按行运算 #

这是最适配dplyr管道逻辑的写法，计算完成后建议用ungroup()解除行分组，避免后续数据操作受影响。
为避免重复调用函数浪费算力，可以先将单次计算的完整结果存入临时列表列，再一次性提取需要的两个统计量：

library(estmeansd)
library(dplyr)

set.seed(1)
df2 <- df %>%
  rowwise() %>%
  mutate(
    bc_calc = list(bc.mean.sd(q1.val = iqr1, med.val = median, q3.val = iqr2, n = n)),
    est_mean = bc_calc$est.mean,
    est_sd = bc_calc$est.sd
  ) %>%
  select(-bc_calc) %>%
  ungroup()

方法2：使用purrr::pmap做逐行映射 #

通过purrr族的映射函数显式逐行传入参数，逻辑更明确：

library(estmeansd)
library(dplyr)
library(purrr)

set.seed(1)
df2 <- df %>%
  mutate(
    bc_calc = pmap(
      list(q1 = iqr1, med = median, q3 = iqr2, sample_n = n),
      ~bc.mean.sd(q1.val = ..1, med.val = ..2, q3.val = ..3, n = ..4)
    ),
    est_mean = map_dbl(bc_calc, ~.x$est.mean),
    est_sd = map_dbl(bc_calc, ~.x$est.sd)
  ) %>%
  select(-bc_calc)

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