用于计算Cohen的f2效果大小的R代码（适用于多级模型）

Question

到目前为止，我还没有找到用于计算f²Cohen效果大小的R语法。这是来自Selya等人的paper。 （2012）谈论它并提供SAS代码。还有Stata代码here。如果我敢问，这里的任何人都可以将SAS或Stata代码转换为R代码吗？如果我的愿望得到批准，我真的很感激帮助。

Answer 1

尝试CRAN上的pwr包：

https://cran.r-project.org/web/packages/pwr/index.html

您拥有cohen.ES()功能：

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

或者根据您是否使用比例，卡方检验，t检验，ANOVA等进行不同的效果大小计算器

举个例子：

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

Answer 2

我也偶然发现了一个相同的问题，并且没有找到可用的函数来计算f2的多元回归。但是，我在Quick-R Power Analysis page上找到了公式，并编写了自己的代码来对其进行计算：

# 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)

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