R中MANOVA分析如何获取每个DV的partial eta squared效应量

默认调用effectsize::eta_squared()或heplots::etasq()仅返回全局MANOVA效应量，是因为函数默认参数适配多元效应量输出场景，调整参数即可直接得到每个因变量（DV）对应的偏η²（np2），无需逐个拟合单因素ANOVA模型，结果也和eta_squared文档中「mlm/maov类模型会单独计算每个响应变量效应量」的描述完全匹配。

具体操作步骤 #

• 第一步：拟合manova模型，示例代码如下：

# 基于iris数据集的示例MANOVA拟合
manova_mod <- manova(cbind(Sepal.Length, Sepal.Width, Petal.Length, Petal.Width) ~ Species, data = iris)

• 第二步：调用effectsize::eta_squared()时显式设置multivariate = FALSE，即可输出各DV对应的偏η²：

# partial = TRUE指定输出偏η²，multivariate = FALSE指定按单个DV拆分计算
effectsize::eta_squared(manova_mod, partial = TRUE, multivariate = FALSE)

• 若使用heplots包的etasq函数，添加anova = TRUE参数即可得到对应结果：

heplots::etasq(manova_mod, anova = TRUE)

上述方法输出的各DV偏η²结果与SPSS中勾选「Estimates of effect size」得到的结果完全一致，直接复用已拟合的MANOVA模型，计算效率更高。

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