python-3.x - 计算一个样本单侧测试的功效

Question

我有下面的代码，我用它来尝试计算在给定样本量下检测与下面一样大的效应量的统计能力。我正在计算一个样本测试的功效，备择假设是样本均值更大。无论我将 alpha 参数更改为什么，我似乎都得到了 1.0 的幂，这看起来很大且不正确。我在下面做错了吗？如果是这样，有人可以指出我需要改变什么吗？如果是正确的，有人可以解释为什么吗？

价值观：

mean of population:
abar=35.95678241834582

mean of sample:
bbar=42.434961073666734

standard deviation of population:
astd=36.668227035470316

observations in sample:
nb=13933

代码：

from statsmodels.stats.power import TTestIndPower

power_analysis = TTestIndPower()

effct_size = (bbar-abar)/astd

power = power_analysis.solve_power(effect_size = effct_size,
                                   nobs1=nb,
                                         power=None,
                                         alpha = 0.05,
                                  alternative='larger')

power

输出：

1.0

标签： python-3.xstatsmodelshypothesis-test