R 计算Cohen'；s f2效应大小（对于多级模型）

直到现在，我还没有找到计算f²Cohen效应大小的R语法。Selya et al.（2012）对此进行了讨论，并给出了SAS代码。还有Stata代码。如果我敢问，这里有人能把SAS或Stata代码转换成R代码吗？如果我的愿望得以实现，我将非常感谢您的帮助。

请尝试上的

pwr

软件包

您有

cohen.ES（）

函数：

library(pwr)
cohen.ES(test = 'f2', size = 'large')

或者根据比例、卡方检验、t检验、方差分析等使用不同的效应大小计算器

例如：

pwr.anova.test(k = 4, n = 20, sig.level = 0.05, power = 0.8)

返回：

Balanced one-way analysis of variance power calculation 

          k = 4
          n = 20
          f = 0.3788014
  sig.level = 0.05
      power = 0.8

NOTE: n is number in each group

---

我也偶然发现了同样的问题，没有找到一个可用的函数来计算多元回归的f2。但是，我在上找到了公式，并编写了自己的代码来计算：

# If we are evaluating the impact of a set of predictors on an outcome, then the f2 formula is:
R2 = # input your R2
f2 = R2/(1 - R2)

# Else, if we are evaluating the impact of one set of predictors above and beyond a second set of predictors (or covariates), then the f2 formula is:
R2AB = # input your R2AB
R2A = # input your R2A
f2 = (R2AB - R2A)/(1 - R2AB)

---

但是，如果您有多级数据，您可能需要调整此公式。

谢谢您的回答。我想我的问题不够清楚。我引用的两篇文章提到了如何计算多水平回归模型的Cohen's f2。我修改了问题的标题。