r - 在 R WSR2 中使用鲁棒方差分析进行事后比较

Question

我正在尝试使用 R 中的 WRS2 包进行稳健的 ANOVA，包括事后测试。我的真实数据结构是这样的：1 个因变量，1 个因变量（具有 2 个不同条件）和 1 个因子内（具有两个条件）。我有 110 个科目。

直到 bwtrim 它工作正常但是当我尝试使用这行代码检查事后：

WRS2::sppbb（公式 = 值 ~ 系统 * 机动，id = ID，数据 = df）

我收到以下错误：

combn(levels(mf[, fixvar]), 2) 中的错误：n < m

我不知道该怎么做，我尝试调试，但似乎找不到错误。我检查了我的数据，一切都完成了——所以没有丢失数据。此外，主题等的索引是正确的。在这个问题上，我真的很感激帮助。

我很乐意得到回应！

最好的

一个可重现的数据集是这样的：

df <- read.table(header=TRUE, text="
ID System Maneuver value  
1 A E 3
1 A R 2  
2 A E 3
2 A R 2   
3 A E 2
3 A R 6  
4 B E 5  
4 B R 3 
5 B E 6 
5 B R 6 
6 B E 4
6 B R 5
7 A E 3
7 A R 2  
8 A E 3
8 A R 2   
9 A E 2
9 A R 6  
10 B E 5  
10 B R 3 
11 B E 6 
11 B R 6 
12 B E 4
12 B R 5
13 A E 3
13 A R 2  
14 A E 3
14 A R 2   
15 A E 2
15 A R 6  
16 B E 5  
16 B R 3 
17 B E 6 
17 B R 6 
18 B E 4
18 B R 5
")
WRS2::bwtrim(formula = value ~ System * Maneuver, id = ID, data = df)

WRS2::sppbb(formula = value ~ System * Maneuver, id = ID, data = df)

标签： r