在 R 中使用 Anova 或 ezANOVA 的 Huynh-Feldt 值有问题吗?

机器算法验证 r 方差分析 重复测量 球形度
2022-04-06 07:27:32

在 R 中对因子之间和因子内(混合设计)运行 ANOVA 时,Huynh-Feldt 校正的输出是否存在问题?我的问题是,当我使用“Anova”或“ezANOVA”函数运行“受试者内(仅)”ANOVA 时,Huynh-Feldt 校正的输出与使用 SPSS 时得到的完全相同。然而,当我在混合设计上运行 ANOVA,并使用“Anova”或“ezANOVA”函数考虑主体效应之间和内部效应时,我的 Greenhouse-Geisser 校正与 SPSS 匹配,但 Huynh-Feldt 校正略有偏差与SPSS的输出相比?这可能是什么原因?我感谢任何人可以提供的任何建议/帮助!

1个回答

如果 Anova 中有错误,ezANOVA 也会出现错误,因为 Mike Lawrence 在 John Fox 的 Anova 包上构建了 ez。而 Anova 包又使用标准的 Ranova命令。R 是正确的。 Peter Dalgaard (pp.3-4) 指出 SAS 和 SPSS 的 HF 计算是错误的,并且已经存在了几十年。这是因为他们使用的是原始公布的公式,而不是 R 使用的更正方程(Lecoutre,1991)。

(我对他们为什么没有修复它的猜测类似于 Excel 仍然存在错误的原因......他们对主要客户有更大的义务与错误兼容而不是正确。修复这个问题也将承认他们错了。虽然,如果它是在 1992 年完成的,那将是合理的。现在它只是开源的一个很好的论据。或者,他们可能不在乎。我想知道一个人是否只是简单地测试它的输出与另一个?)