Returns the (partial) eta-squared, (partial) omega-squared,
epsilon-squared statistic or Cohen's F for all terms in an anovas.
anova_stats() returns a tidy summary, including all these statistics
and power for each term.
anova_stats(model, digits = 3)A fitted anova-model of class aov or anova. Other
models are coerced to anova.
Amount of digits for returned values.
A data frame with all statistics is returned (excluding confidence intervals).
Levine TR, Hullett CR (2002): Eta Squared, Partial Eta Squared, and Misreporting of Effect Size in Communication Research.
Tippey K, Longnecker MT (2016): An Ad Hoc Method for Computing Pseudo-Effect Size for Mixed Model.
# load sample data
data(efc)
# fit linear model
fit <- aov(
c12hour ~ as.factor(e42dep) + as.factor(c172code) + c160age,
data = efc
)
anova_stats(car::Anova(fit, type = 2))
#> etasq | partial.etasq | omegasq | partial.omegasq | epsilonsq | cohens.f | term | sumsq | df | meansq | statistic | p.value | power
#> ----------------------------------------------------------------------------------------------------------------------------------------------------------
#> 0.212 | 0.224 | 0.209 | 0.221 | 0.209 | 0.537 | as.factor(e42dep) | 4.265e+05 | 3 | 1.422e+05 | 80.299 | < .001 | 1.000
#> 0.004 | 0.005 | 0.002 | 0.003 | 0.002 | 0.071 | as.factor(c172code) | 7352.049 | 2 | 3676.025 | 2.076 | 0.126 | 0.429
#> 0.052 | 0.066 | 0.051 | 0.065 | 0.051 | 0.267 | c160age | 1.052e+05 | 1 | 1.052e+05 | 59.408 | < .001 | 1.000
#> | | | | | | Residuals | 1.476e+06 | 834 | 1770.307 | | |