I've conducted a Welch Anova from summary statistics with the following valid results:
# Welch's ANOVA (One-way analysis of means)
data: group and scores
F value df1 df2 p-value
5.35 2 16.83 0.015912 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
summary statistics
group size mean sd
<chr> <int> <dbl> <dbl>
A 10 77.3 11.324016
B 10 89.3 4.808557
C 10 84.7 5.292552
underlying data
df <- data.frame(group = rep(c('A', 'B', 'C'), each = 10),
score = c(64, 66, 68, 75, 78, 94, 98, 79, 71, 80,
91, 92, 93, 85, 87, 84, 82, 88, 95, 96,
79, 78, 88, 94, 92, 85, 83, 85, 82, 81))
#summary
group size mean sd
<chr> <int> <dbl> <dbl>
A 10 77.3 11.324016
B 10 89.3 4.808557
C 10 84.7 5.292552
I found some references on how to calculate omega squared: https://peterstatistics.com/Packages/python-docs/effect_sizes/eff_size_omega_sq.html
Since i dont have calculated any sum of squares, which formula is appropriate to retrieve omega squared? The reference states, all formulas give the same result. But somehow the results using the formulas from Kirk 1996 and Caroll and Nordholm 1975 differ from the formula provided by Hays 1973 (and Albers and Lakens 2018). Any advice?