使用R语言WRS2包的bwtrim执行稳健ANOVA时出现错误求助

问题场景 #

使用WRS2包的bwtrim函数对配对设计数据（每个被试分A、B两个条件）做稳健ANOVA，运行代码后触发错误：Error in '[.data.frame'(mf, , 3) : undefined columns selected

错误原因 #

bwtrim的id参数仅用于重复测量次数≥3次的重复测量设计，而你的数据是配对设计（每个被试仅2次观测），此时不需要指定id参数。指定该参数后，函数会错误地尝试解析数据结构，导致找不到对应列的报错。

修正代码 #

df <- data.frame(
  id = c(1, 1, 2, 2, 3, 3, 5, 5),
  condition = c("A", "B", "A", "B", "A", "B", "A", "B"),
  var1 = c(1.70, 0.80, 0.80, 1.78, 5.61, 1.53, 7.64, 3.92),
  var2 = c(0.85, 2.41, 1.50, 1.3, 0.45, 1.53, 0.00, 1.34)
)
df$condition <- as.factor(df$condition)
library(WRS2)
# 配对设计无需指定id参数
bwtrim(formula = var1 ~ condition, data = df)

额外说明 #

若你的数据是重复测量设计（每个被试有3次及以上重复观测），才需要添加id = id参数来指定被试标识。

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