在R中从ANOVA获取因变量的估计边际均值

在R中，aov()仅输出方差分析的核心统计量（自由度、平方和、F值等），无法直接给出控制其他变量后的估计边际均值。要实现和SPSS相同的功能，可借助emmeans包完成。

步骤1：安装并加载emmeans包 #

若未安装该包，先执行安装命令：

install.packages("emmeans")

加载包：

library(emmeans)

步骤2：计算自变量ind.v各水平的估计边际均值 #

针对你拟合的my_anova模型，使用emmeans()函数指定目标因子ind.v，函数会自动控制模型中的Control1和Control2变量，输出调整后的均值：

# 计算ind.v各水平的估计边际均值
emm_results <- emmeans(my_anova, ~ ind.v)

# 查看完整结果（含均值、标准误、置信区间）
print(emm_results)

额外功能：两两比较 #

如果需要对ind.v各水平的边际均值进行两两比较，可在结果上调用pairs()函数：

# 执行两两比较并输出显著性结果
pairs(emm_results)

注：你的数据集样本量较小，且模型包含三因素交互项，部分组合可能无观测值，emmeans会基于现有数据进行合理的均值估计。

内容的提问来源于stack exchange，提问作者Naz