我尝试在 R 中执行稳健的混合方差分析,因为违反了同方差性、正态分布和协方差相等的假设。我在主题因素(群体)之间有一个,在主题因素(情况)内有一个。我总共有两个问题。
问题1:如何获取对应的效果大小?
mixedanova <- WRS2::bwtrim(formula = value ~ group*situation, id = id, data = data.long)
输出:
价值 | df1 | df2 | p 值 | |
---|---|---|---|---|
组 | 162.9427 | 2 | 33.9744 | 0 |
情况 | 136.4997 | 1 | 38.5618 | 0 |
组:情况 | 61.4333 | 2 | 2040.28 | 0 |
问题 2:如何解释事后测试?(为什么值大于 1?)(为什么当我更改 nboot 时输出没有改变?)
sppba(formula = value ~ group*situation, id = id, data = data.long, est = "mom", nboot = 500)
Test statistics:
Test whether the corrresponding population parameters are the same:
p-value: 0
组 | 估计 |
---|---|
1-2 | 0.8178253 |
1-3 | 1.4685165 |
2-3 | 0.6506912 |
sppbb(formula = value ~ group * situation, id = id, data = data.long, est = "mom", nboot = 500)
Test statistics:
Test whether the corrresponding population parameters are the same:
p-value: 0
估计 |
|
---|---|
情况1-情况2 | 0.2575347 |
sppbi(formula = value ~ group * situation, id = id, data = data.long, est = "mom", nboot = 500)
Test statistics:
Test whether the corrresponding population parameters are the same:
p-value: 0
标题 1 | 标题 2 |
---|---|
情况1-情况2 1-2 | -0.58068783 |
情况1-情况2 1-3 | 0.07747543 |
情况1-情况2 2-3 | 0.65816327 |
提前感谢您的帮助!!
我已通读 Mair & Wilcox (2020) 的论文“使用 WRS2 的稳健统计方法”,并尝试改变函数 bwtrim、sppba、sppbb 和 sppbi,但不幸的是没有成功。
你有答案了吗?我有完全相同的问题。将不胜感激您的帮助。
许多问候
丽莎