在R中使用Wilcox robust ANOVA进行计划对比(WRS2软件包)

Question

我使用WRS2对我的数据集进行了双向方差分析和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)会提供另一种方法……