R语言使用tidyverse按分组匹配总计数计算百分比

你之前写的代码存在两处核心逻辑缺失，所以无法得到正确结果：

• 代码中没有定义total变量，直接调用会触发对象不存在的报错
• 没有在Type+Stream_Order的分组维度下，精准提取Group == "all"对应的总计数作为分母，无法匹配到每个分组对应的正确总值

使用dplyr包的分组计算逻辑即可实现，代码可适配任意规模的同结构数据集，只要每个Type+Stream_Order组合下存在唯一一条Group == "all"的总计数记录即可：

library(dplyr)

perc_calc <- perc %>%
  group_by(Type, Stream_Order) %>%
  # 提取当前分组下all对应的总计数，再计算百分比
  mutate(
    total = n[Group == "all"],
    percentage = n / total * 100
    # 如需保留2位小数可替换为下一行代码
    # percentage = round(n / total * 100, 2)
  ) %>%
  ungroup()

针对你给出的测试样例，运行代码后你提到的Type = "Sediment"、Stream_Order = "1"、Group = "stoch"的记录，计算得到的percentage值为22.89，和手动计算结果完全一致。
你可以提前运行以下代码校验数据合规性，避免因缺失all记录导致计算结果出现NA：

# 检查每个维度组合下all记录的数量，返回空表则数据符合计算要求
perc %>%
  group_by(Type, Stream_Order) %>%
  summarise(all_record_num = sum(Group == "all"), .groups = "drop") %>%
  filter(all_record_num != 1)

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