Рассчитать Omega² и частичную Omega² в смешанном Anova

Я анализирую данные в смешанном дизайне, имея 1 между и 1 внутри фактора, используя aov_car() из пакета afex.

Я попытался использовать anova_stats() из пакета sjstats, который отлично работает для внутри- и между проектами, но выдает ошибку:

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

Ошибка на английском языке (переводчик Google):

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
person mailach    schedule 21.05.2019    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Я только что понял, что не так. return в aov_car() должно быть установлено на "aov". Так что это сработало для меня, используя этот код:

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