在混合方差分析中计算 Omega² 和部分 Omega²答案

【问题标题】：Calculate Omega² and partial Omega² in mixed-Anova在混合方差分析中计算 Omega² 和部分 Omega²
【发布时间】：2019-05-21 10:24:03
【问题描述】：

我正在使用 afex 包中的 aov_car() 分析混合设计中的数据，其中 1 介于和 1 之间。

我尝试使用 sjstats 包中的 anova_stats()，它非常适用于设计内部或设计之间，但会引发错误：

Fehler in UseMethod("anova") : 
nicht anwendbare Methode für 'anova' auf Objekt der Klasse "Anova.mlm" angewendet

英文错误（谷歌翻译）：

Error in UseMethod ("anova"):
non-applicable method for 'anova' applied to object of class "Anova.mlm"

表示aov_car()的返回数据类型存在问题

library(sjstats)
library(afex)
data <- data.frame(ID = rep(c(1:24), each = 4),
              list = rep(c(1:4), 24),
              prop = rnorm(96, 0.8, 0.2),
              treatment = rep(c(1:4), each = 24))
data$treatment <- as.factor(data$treatment)
data$list <- as.factor(data$list)

testAov <- aov_car(prop~list*treatment + Error(ID/list), data=data, return = "Anova")
anova_stats(testAov)

如果有人能帮我解决这个问题，我将非常感激！

关于我的系统和 R 版本的信息：

    R version 3.6.0 (2019-04-26)  
    Platform: x86_64-pc-linux-gnu (64-bit)  
    Running under: Arch Linux

【问题讨论】：

标签： r anova mixed-models

【解决方案1】：

我刚刚发现出了什么问题。 aov_car() 中的 return 必须设置为 "aov"。所以它使用这段代码对我有用：

testAov <- aov_car(prop~list*treatment + Error(ID/list), data=data, return = "aov")
anova_stats(testAov)

【讨论】：