在R中有一个非常有用的函数,它有助于确定双边t检验的参数,以获得目标统计功率。
函数名为power.prop.test
。
http://stat.ethz.ch/R-manual/R-patched/library/stats/html/power.prop.test.html
您可以使用以下方法调用它:
power.prop.test(p1 = .50, p2 = .75, power = .90)
它会告诉你获得这种能量所需的样本量。这在确定测试样本量时非常有用。
scipy包中有类似的功能吗?
Tags:
一些基本的功率计算现在可以在statsmodels中获得
http://statsmodels.sourceforge.net/devel/stats.html#power-and-sample-size-calculationshttp://jpktd.blogspot.ca/2013/03/statistical-power-in-statsmodels.html
博客文章还没有考虑到statsmodels代码的最新更改。另外,我还没有决定要提供多少包装器函数,因为许多功率计算只是减少到基本分布。
在R统计中
使用R的
pwr
包我已经使用下面的n公式和scipy.stats中的逆生存函数
norm.isf
成功地复制了这个函数马特得到所需n(每组)的答案几乎是正确的,但有一个小错误。
给定d(均值差)、s(标准差)、sig(显著性水平,通常为0.05)和幂(通常为0.80),计算每组观察数的公式为:
从他的公式中可以看出
“s”部分是错误的。再现r的功能的正确函数是:
希望这有帮助。
相关问题 更多 >
编程相关推荐