在 R 中使用稳健方差分析进行事后比较的问题

Question

我正在尝试在 R 中使用 WRS2 package 进行稳健的方差分析，包括事后测试。 我的数据结构是这样的：1个因变量（res）2个自变量：（组）之间-（bla）1个ID变量引用主题（id）

我的代码如下所示：

fileHB=read.csv2("trial_II.csv");
fileHB$Hb=as.numeric(fileHB$Hb);

m1 <- WRS2::bwtrim(res ~ Group*Bla,id=id, data = fileHB)

到目前为止，它工作正常，但是当我尝试使用这行代码检查事后：

sppba(res ~ Group*Bla,id,data = fileHB)

我收到以下错误：

error in combn(levels(mf[, fixvar]), 2) : n < m

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

我很乐意得到回应！

最好的

Answer 1

前段时间我遇到了同样的错误，尝试将组转换为因素。 对我来说，tjos 解决了我的问题

fileHB$Group <- as.factor(fileHB$Group)

如果不是，另一种可能性是您缺少一些数据，因此您的组的大小不相等。 然后对我来说，事后没有工作，我不得不进行估算，但这取决于您的缺失案例如何，哪种策略更适合您。

Answer 2

嗨，“sppba” function 不是事后测试。 这是一个综合测试。