如何在R中批量提取各Weed_type的untreated control值并计算防效百分比？

我是R语言新手，现在需要针对每个Weed_type提取Treatment为"untreated control"的predicted.value，再用公式(untreated - predicted.value)/untreated*100计算防效百分比（percentage control）。因为杂草类型很多，现在分步提取再匹配的方法太麻烦，想找更简便的实现方式。

原始数据 #

df <- data.frame(Weed_type = c("weed1", "weed1", "weed1", "weed2", "weed2", "weed2", "weed3", "weed3", "weed3"),
                 Treatment = c("untreated control", "Treatment1", "Treatment2", "untreated control", "Treatment1", "Treatment2", "untreated control", "Treatment1", "Treatment2"),
                 predicted.value = c(23.3, 0.4, 0,  .9, .15, .01, 87, 12,2)
)
df

原繁琐实现步骤 #

1. 提取各杂草类型的对照组值

weed1_c <- df %>% filter(Weed_type == 'weed1' & Treatment == "untreated control")  %>% pull(predicted.value)
weed2_c <- df %>% filter(Weed_type == 'weed2' & Treatment == "untreated control")  %>% pull(predicted.value)
weed3_c <- df %>% filter(Weed_type == 'weed3' & Treatment == "untreated control")  %>% pull(predicted.value)

2. 将对照组值匹配回原数据框

df <- df %>% 
      mutate(untreated = case_when(
        Weed_type=="weed1" ~   weed1_c,
        Weed_type=="weed2" ~   weed2_c,
        Weed_type=="weed3" ~   weed3_c,
      ))

3. 计算防效百分比

df$percentage_control <- (df$untreated  - df$predicted.value) / df$untreated *100
df  

简化实现方案 #

方法一：分组直接计算（推荐） #

利用dplyr的分组功能，在组内直接提取对照组数值并计算，一步完成所有操作：

library(dplyr)

df <- df %>%
  group_by(Weed_type) %>%
  mutate(
    # 提取当前组内对照组的predicted.value
    untreated = predicted.value[Treatment == "untreated control"],
    # 计算防效百分比
    percentage_control = (untreated - predicted.value)/untreated * 100
  ) %>%
  ungroup() # 可选，取消分组恢复普通数据框

df

关键逻辑

• group_by(Weed_type)：按杂草类型分组，后续操作仅在同组内执行
• predicted.value[Treatment == "untreated control"]：在每组中筛选出对照组的数值，自动填充到该组所有行
• 直接在mutate中完成防效计算，无需额外步骤

方法二：提取对照组数据后合并 #

先单独提取对照组数据，再通过left_join合并到原数据框后计算：

library(dplyr)

# 提取对照组数据，重命名列便于合并
control_df <- df %>%
  filter(Treatment == "untreated control") %>%
  select(Weed_type, untreated = predicted.value)

# 合并数据并计算防效百分比
df <- df %>%
  left_join(control_df, by = "Weed_type") %>%
  mutate(percentage_control = (untreated - predicted.value)/untreated * 100)

df

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