我需要找到进行卡方检验所需的样本量。我的信息是alpha(0.05)和power(0.8)。我没有效果大小。如何找到所需的样本量。
答案 0 :(得分:1)
首先,加载pwr库。如果你想找到N,你必须估计效果大小w,并从你的列联表中知道df(row-1)(col-1)。因此,对于第一组中0.3的比例和第二组中0.5的比例(ES函数将从这些中给出效果大小),2x2列联表(df =(2-1)(2- 1)),所需显着性水平为0.95,所需功效为0.8:
> library(pwr)
> pwr.chisq.test(w=ES.w1(0.3,0.5),df=1,sig.level=0.05,power=0.80)
Chi squared power calculation
w = 0.3651484
N = 58.86645
df = 1
sig.level = 0.05
power = 0.8
NOTE: N is the number of observations
如果您输入数字并收到如下错误:
Error in uniroot(function(N) eval(p.body) - power, c(1 + 0.0000000001, :
f() values at end points not of opposite sign
......这里的其他地方有几个答案可以帮助你完成这个任务;见https://stackoverflow.com/search?q=Error+in+uniroot