如何在dplyr的reframe()函数中添加行标签列

需求背景 #

使用dplyr::reframe()对分组数据计算分位数时，每组会生成对应分位数的多行结果，但缺少明确的行标签（如lowCI、median、hiCI）。希望无需提前知道结果行数，直接在流程内完成标签添加。

示例数据集：

d <- data.frame(group = rep(letters[1:3],each=30),
                condition = rep(1:3,times=30),
                score = sample(1:10,size=90,replace=T))

解决方案1：在reframe()内直接生成标签列 #

这是最简洁的方案，利用reframe()在分组内会自动绑定同长度向量的特性，直接在reframe()中同时定义标签列和分位数列：

d %>%
  group_by(group, condition) %>%
  reframe(
    quantiles = c("lowCI", "median", "hiCI"),
    sumScore = quantile(score, probs = c(0.025, 0.5, 0.975))
  )

效果：每组会生成3行数据，quantiles列与sumScore列一一对应，无需后续额外处理。

解决方案2：利用enframe转换带命名的分位数结果 #

如果需要保留分位数的概率命名（如2.5%）再转换为自定义标签，可以结合tibble::enframe()：

library(tibble)
library(dplyr)

d %>%
  group_by(group, condition) %>%
  reframe(
    quantile(score, probs = c(0.025, 0.5, 0.975)) %>%
      enframe(name = "quantiles", value = "sumScore") %>%
      mutate(quantiles = case_when(
        quantiles == "2.5%" ~ "lowCI",
        quantiles == "50%" ~ "median",
        quantiles == "97.5%" ~ "hiCI"
      ))
  )

效果：先将quantile()返回的带命名向量转换为两列，再通过mutate()替换为自定义标签，适合需要处理分位数命名的场景。

解决方案3：用rep_len自动适配结果行数 #

如果偏好后续添加列的方式，可以使用rep_len()替代rep()，它会自动循环向量直到匹配目标长度，无需提前知道结果行数：

d %>%
  group_by(group, condition) %>%
  reframe(sumScore = quantile(score, probs = c(0.025, 0.5, 0.975))) %>%
  add_column(quantiles = rep_len(c("lowCI", "median", "hiCI"), length.out = nrow(.)))

效果：rep_len()会根据当前数据框的行数（nrow(.)）自动重复标签向量，避免了手动指定重复次数的麻烦。

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