到目前为止,我还没有找到用于计算f²Cohen效果大小的R语法。这是来自Selya等人的paper。 (2012)谈论它并提供SAS代码。还有Stata代码here。如果我敢问,这里的任何人都可以将SAS或Stata代码转换为R代码吗?如果我的愿望得到批准,我真的很感激帮助。
答案 0 :(得分:0)
尝试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
答案 1 :(得分:0)
我也偶然发现了一个相同的问题,并且没有找到可用的函数来计算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)
但是,如果您有多级数据,则可能需要调整此公式。