r - 在 R(WRS2 包)中使用 Wilcox 稳健方差分析进行计划对比

Question

我使用 WRS2 对我的数据集运行双向方差分析和事后测试。然而在对比中我找不到我需要的比较。是否可以计划我自己的对比？

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 )会提供替代路线......