感谢 @hcp4715 和 @李晓煦 两位老师之前的精彩回答!
这个问题是我几个月前提的,当时还很少用R来做传统意义上的方差分析,所以比较想知道“如何使用R做方差分析、简单效应检验、事后多重比较”。
时隔多月,在已有回答和很多其他参考资料的基础上,我在自己的R包(bruceR)里面新增了两个函数:MANOVA和 EMMEANS —— 有了它们,大家就可以用“全宇宙最简单的代码”进行几乎任何类型的方差分析(ANOVA)、简单效应检验(simple-effect analyses)、事后多重比较(post-hoc multiple comparisons)!
具体内容我已在一篇知乎专栏文章中做了介绍:包寒吴霜:bruceR包重磅更新:方差分析、简单效应检验、事后多重比较zhuanlan.zhihu.com
在本回答中,我将展示 MANOVA和 EMMEANS这两个R函数的强大和简约之处。方差分析总流程
首先,这两个函数完整覆盖了方差分析的几乎所有分析流程。大家在上图中已经看到:方差分析涉及主效应和交互作用,如果交互作用显著,我们需要继续做简单效应检验(注:简单效应检验 ≠事后检验或事后比较,两者是独立的两件事!)。而对于因素的水平数≥3的情况,在检验完主效应(或者简单主效应)之后,还需要做事后比较(可以是两两配对比较、顺序依次比较、正交多项式比较等等),这其中可能会涉及到多重比较校正的问题(如Bonferroni校正、Tukey's HSD、FDR校正)。
例如,我们有一个“2×3两因素被试间设计”的数据