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）将提供替代路线......

在R（WRS2包）中使用Wilcox强大ANOVA进行计划对比

1 个答案: