使用R Tidyverse计算分组变量最优组均值与其他组重复的差值

原代码问题说明 #

你之前的代码错误出在分组逻辑：按Diversity和Replicate分组后，每个分组仅包含1条optimal处理的记录，此时计算的mean(Soil_N[Soilwater=="optimal"]) 本质还是单个重复的optimal值，没有用到整个optimal组的均值。

修改后可运行代码 #

library(tidyverse)

# 原始数据集
df <- data.frame(Soilwater = c("optimal", "optimal", "optimal", "optimal", "optimal", 
                               "40", "40", "40", "40", "40", 
                               "30","30","30","30","30", 
                               "20", "20","20","20","20",
                               "10","10","10","10","10", 
                               "optimal", "optimal", "optimal", "optimal", "optimal", 
                               "40", "40", "40", "40", "40", 
                               "30","30","30","30","30", 
                               "20", "20","20","20","20",
                               "10","10","10","10","10"), 
                 Diversity = c("High","High","High","High","High","High","High","High","High","High",   
                               "High","High","High","High","High","High","High","High","High","High",
                               "High","High","High","High","High", 
                               "Low", "Low", "Low","Low","Low","Low","Low","Low","Low","Low",
                               "Low","Low","Low","Low","Low","Low","Low","Low","Low","Low",
                               "Low","Low","Low","Low","Low"),
                 Soil_N = c(50,45, 49, 48, 49, 69, 68, 69, 70, 67, 79, 78, 79, 78, 77, 89, 89, 87, 88, 89, 99, 98, 97, 98, 98, 120,    
                            121,    121,    120,    122,    134,    131,    132,    134,    131,    145,    148,    149,    147,    
                            148,    159,    159,    157,    156,    157,    169,    167,    167,    168,    164))

# 核心计算代码
df_result <- df %>%
  # 按多样性分组，分别计算高、低多样性组下optimal处理的平均土壤氮含量
  group_by(Diversity) %>%
  mutate(opt_mean_N = mean(Soil_N[Soilwater == "optimal"])) %>%
  # 可选：生成重复编号，方便核对单个重复的计算结果，不需要可以删除该步骤
  group_by(Soilwater, Diversity) %>%
  mutate(Replicate = row_number()) %>%
  # 计算非标准化差值和标准化差值
  mutate(
    unstandard_diff = opt_mean_N - Soil_N,
    standard_diff = (opt_mean_N - Soil_N)/opt_mean_N
  ) %>%
  ungroup()

# 可选：按土壤含水量和多样性分组，汇总差值的平均结果
df_summary <- df_result %>%
  group_by(Soilwater, Diversity) %>%
  summarise(
    mean_unstd_diff = mean(unstandard_diff),
    mean_std_diff = mean(standard_diff),
    .groups = "drop"
  )

代码说明 #

• 分组逻辑符合试验设计：高、低多样性组的optimal均值分开计算，同组内所有行共享对应optimal均值，避免了重复级别的计算偏差
• 可直接输出单重复级别的差值结果，也可通过附加的汇总步骤得到处理水平的统计结果
• 全程使用tidyverse语法，无循环语句，符合你的需求

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