WRS2 包中使用 bwtrim 的简单方差分析会引发错误

Question

出于教育目的，我尝试使用 R 中的 WRS2 包进行简单稳健的方差分析。我想在两种不同条件（A 和 B）下比较变量 var1：

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)
bwtrim(formula = var1 ~ condition, id = id, data = df)

不幸的是，代码导致了一个我无法解释的错误：

Error in '[.data.frame'(mf, , 3) : undefined columns selected

。我的代码或对 bwtrim 函数的理解有什么问题？

Answer 1

根据文档，

bwtrim

是

专为一个受试者间变量和一个受试者内变量设计

您的公式右侧需要两个分类变量。

如果您正在使用 WRS2 软件包寻找简单的单向重复测量方差分析来进行修剪均值，请尝试

pairdepb

pairdepb(df$var1, df$condition, df$id)
# Call:
# pairdepb(y = df$var1, groups = df$condition, blocks = df$id)
# 
#         psihat  ci.lower ci.upper    test     crit   sig
# A vs. B   1.93 -23.19144 27.05144 1.60483 20.88889 FALSE