使用R中的Wilcox robust ANOVA進行計劃對比（WRS2包）

Question

我使用WRS2對我的數據集運行雙向ANOVA和posthoc測試。然而在對比之下，我沒有找到我需要的比較。是否有可能計劃自己的對比？

Answer 1

不幸的是，它看起來像mcp2atm函數中的對比是硬編碼的。一些摘自該功能的內容（你可以通過在命令提示符處輸入mcp2atm來看到：你也可以通過dump()函數訪問外部文件）。

... 
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")，或者可能試試issues tracker on R-forge ......）我不知道是否有其他強大的統計機器（見robust stats task view）會提供替代路線...