在R中基于多分组变量计算指定行百分比差异报错求助

我有一个包含多分组变量的数据集，想要计算以form为"t25"的行作为对照组的百分比差异，同时保留所有分组变量，但运行代码时出现错误：

percent_diff <- df %>%
  group_by(days,density,batch,form) %>%
  mutate(pct_diff = 100 * (value - value[form=="t25"]/value[form=="t25"]))

报错信息：

"Caused by error: ! pct_diff must be size 1, not 0."

期望form为"t25"的行的pct_diff值为0，另外两种form的行显示对应的差异值。示例数据如下：

value<-runif(n=36, min=50, max=100) 
days<-c("d1","d1","d1","d1","d1","d1","d1","d1","d1","d1","d1","d1","d2","d2","d2","d2","d2","d2","d2","d2","d2","d2","d2","d2","d3","d3","d3","d3","d3","d3","d3","d3","d3","d3","d3","d3") 

density<-c("1K","1K","1K","2K","2K","2K","1K","1K","1K","2K","2K","2K","1K","1K","1K","2K","2K","2K","1K","1K","1K","2K","2K","2K","1K","1K","1K","2K","2K","2K","1K","1K","1K","2K","2K","2K") 

form<-c("t25","6w","24w","t25","6w","24w","t25","6w","24w","t25","6w","24w","t25","6w","24w","t25","6w","24w","t25","6w","24w","t25","6w","24w","t25","6w","24w","t25","6w","24w","t25","6w","24w","t25","6w","24w") 

batch<-c("r1","r1","r1","r1","r1","r1","r2","r2","r2","r2","r2","r2","r1","r1","r1","r1","r1","r1","r2","r2","r2","r2","r2","r2","r1","r1","r1","r1","r1","r1","r2","r2","r2","r2","r2","r2")

df<-data.frame(days,density,form,batch,value)

错误原因 #

1. 分组逻辑错误：你把form加入了group_by参数，导致每个分组内只有单一form类型。处理非"t25"的分组时，value[form=="t25"]会返回空向量，触发长度不匹配的错误。
2. 公式计算错误：原公式的括号位置错误，正确的百分比差异计算应该是(当前值 - 对照值)/对照值 * 100，而非原写法里的当前值 - (对照值/对照值)。

修正后的代码 #

library(dplyr)

percent_diff <- df %>%
  # 仅按days、density、batch分组，确保每组包含所有form类型
  group_by(days, density, batch) %>%
  # 以同组内form=="t25"的value为对照计算百分比差异
  mutate(pct_diff = 100 * (value - value[form == "t25"]) / value[form == "t25"]) %>%
  # 重新按所有原始变量分组，满足保留所有分组变量的需求
  group_by(days, density, batch, form)

逻辑说明 #

• 分组时排除form，让每个分组包含同一days、density、batch下的所有form数据，才能正确提取"t25"的对照值。
• 公式(value - value[form == "t25"]) / value[form == "t25"] * 100：
  • 对于form=="t25"的行，value与对照值相等，计算结果为0，符合预期。
  • 对于其他form的行，自动计算与同组对照值的百分比差异。
• 最后重新按所有变量分组，确保保留全部原始分组维度。

验证结果 #

运行修正后的代码后，form=="t25"的行中pct_diff列全部为0；6w和24w的行则显示对应的百分比差异值，所有原始分组变量均被保留。

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