在 R 中使用鲁棒方差分析进行事后比较的问题答案

【问题标题】：Problem with post hoc comparison using Robust ANOVA in R在 R 中使用鲁棒方差分析进行事后比较的问题
【发布时间】：2021-03-19 22:43:17
【问题描述】：

我正在尝试使用 R 中的 WRS2 包进行稳健的 ANOVA，包括事后测试。我的数据结构如下： 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

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

我很乐意得到回复！

最好的

【问题讨论】：

【解决方案1】：

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

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

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

【讨论】：