在rstatix函数anova_test中，如何获得ges和pes作为效果大小？

Question

它说，这里，其中一个可以得到一个或另一个或两者。即使在设置effect.size = c("ges", "pes")时，我也能够将每一个单独地得到，但这两者都不能结合在一起。相反，我只得到"pes“。当我使用自己的数据和使用来自hangover包的{WRS2}数据集时，我也遇到了同样的问题。对于宿醉数据，我的代码是：

anova_test(data = hangover, dv = symptoms, wid = id, between = group, within = time, effect.size = c("ges", "pes"))

我将非常感谢您的帮助！

Answer 1

这似乎是文档或代码中的一个错误。文档确实表示您可以同时指定两者，但是相关代码使用if / else只返回其中的一个：

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
}

但是，正如这里所讨论的，并被@phen全方位提到，rstatix目前计算的ges是错误的。

Answer 2

您可以使用effect.size =“两者”来获得两者。

但是，您可能希望了解此包的github页面中提出的有关如何计算effect.size的问题：

https://github.com/kassambara/rstatix/issues/132