在 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 不是事后測試。 這是一個綜合測試。