It says here that one can get one or the other or both. I have been able to get each one separately but not both together even when I set effect.size = c("ges", "pes")
. Instead, I only get "pes". I have the same problem when I use my own data and when I use the hangover
dataset from the {WRS2}
package. For the hangover data, my code is:
anova_test(data = hangover, dv = symptoms, wid = id, between = group, within = time, effect.size = c("ges", "pes"))
I would be very grateful for your help!
This seems to be an error in either the docs or the code. The docs do say you can specify both, but the relevant code uses if / else to return just one or the other:
add_anova_effect_size <- function(res.anova.summary, effect.size = "ges", observed = NULL){
ss.exists <- "SSn" %in% colnames(res.anova.summary$ANOVA)
if(!ss.exists){
return(res.anova.summary)
}
if("pes" %in% effect.size){
res.anova.summary <- res.anova.summary %>%
add_partial_eta_squared()
}
else {
res.anova.summary <- res.anova.summary %>%
add_generalized_eta_squared(observed)
}
res.anova.summary
}
But also, as discussed here and mentioned by @Phenomniverse, rstatix currently calculates ges incorrectly.