是否有python（scipy）函数来确定获得目标功率所需的参数？

Question

在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包中有类似的功能吗？

Answer 1

马特得到所需n（每组）的答案几乎是正确的，但有一个小错误。

给定d（均值差）、s（标准差）、sig（显著性水平，通常为0.05）和幂（通常为0.80），计算每组观察数的公式为：

n= (2s^2 * ((z_(sig/2) + z_power)^2) / (d^2)

从他的公式中可以看出

n = s * ((zp + z)**2) / (d**2)

“s”部分是错误的。再现r的功能的正确函数是：

def sample_power_difftest(d, s, power=0.8, sig=0.05):
    z = norm.isf([sig/2]) 
    zp = -1 * norm.isf([power])
    n = (2*(s**2)) * ((zp + z)**2) / (d**2)
    return int(round(n[0]))

希望这有帮助。