使用R中的Wilcox robust ANOVA进行计划对比（WRS2包）
问题描述：
我使用WRS2对我的数据集运行双向ANOVA和posthoc测试。然而在对比之下，我没有找到我需要的比较。是否有可能计划自己的对比？使用R中的Wilcox robust ANOVA进行计划对比（WRS2包）
答
不幸的是，它看起来像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）会提供替代路线...