现在rmanova改为调用稳健的混合设计 ANOVA 函数tsplit，不幸的是 WRS2 包尚未包含相应的事后测试。（mcp2atm仅用于t2way）。如果您需要进行事后测试，您现在必须使用原始 WRS 包。

这个包包含bwmcp变量A（之间）、B（内部）的多重比较以及这两者之间的交互作用的函数。如果您只对其中一个感兴趣，您可以分别使用bwamcpfor A 和bwbmcpfor Bbwimcp进行交互。这些函数的用法并不那么简单：

# load the latest WRS package
> source("https://dornsife.usc.edu/assets/sites/239/docs/Rallfun-v38.txt")

# convert your data for the function from long format into list format. 

> data <- bw2list(yourdf, 4, c(7,8))

# In this example "4" stands for the column where your independent groups are located.
# c(7,8) indicates the position of the columns with the within groups data

现在我们已经准备好数据，假设我们的主体间变量有 3 个水平（组），而我们的主体内变量（如上所示）有 2 个水平。该函数将如下所示：

> bwmcp(3,2, data, tr = 0.2) # uses by default 20% trimmed means

您可以使用 WRS2 包并运行以下事后测试：sppbb（在效果 p 内）、sppba（在效果之间）和 sppbi（交互效果）。

您可以在此处阅读更多相关信息：https ://link.springer.com/article/10.3758/s13428-019-01246-w