在R（WRS2包）中使用Wilcox強大ANOVA進行計划對比

Question

我使用WRS2在我的數據集上運行雙向ANOVA和posthoc測試。 然而在對比中，我找不到我需要的比較。 是否可以規划自己的對比？

Answer 1

不幸的是，看起來對比在mcp2atm函數中是硬編碼的。 一些來自該函數內容的摘錄（您可以通過在命令提示符下鍵入mcp2atm來查看：您還可以將函數dump() mcp2atm到外部文件）。

...
temp <- con2way(J, K)
conA <- temp$conA
conB <- temp$conB
conAB <- temp$conAB
if (!op) {  ## note: op is hardcoded to FALSE earlier in the function
    Factor.A <- lincon1(x, con = conA, tr = tr, alpha = alpha)
    Factor.B <- lincon1(x, con = conB, tr = tr, alpha = alpha)
    Factor.AB <- lincon1(x, con = conAB, tr = tr, alpha = alpha)
}
## ... stuff for generating contrast names ...
contrasts <- as.data.frame(cbind(conA, conB, conAB))
colnames(contrasts) <- c(dnamesA, dnamesB, dnamesAB)
rownames(contrasts) <- colnames(dataWide)
...

con2way和lincon1是內部WRS2函數（你可以使用WRS2::con2way和WRS2:::lincon1看到它們）。 （如果你真的想深入研究源代碼，下載源代碼包並解壓縮它可能是有意義的。）

要進行自定義對比，您可能需要與WRS2 （ maintainer("WRS2") ，或者可能嘗試使用R- maintainer("WRS2")的問題跟蹤器 ...）我不知道是否還有其他任何強大的統計機器（參見強大的統計任務視圖 ）將提供替代路線......