python - Statsmodels 功效分析 观察次数

Question

我对TTestIndPower类的 statsmodels API 感到困惑。

在plot_power函数中绘制的观察数是样本总数还是每组的样本数？

Answer 1

是每组。

这是我用 statsmodels 生成的 power_plot

from statsmodels.stats.power import TTestIndPower

panalysis = TTestIndPower()

panalysis.plot_power(
    dep_var="nobs",
    nobs=np.arange(5, 101),
    effect_size=np.arange(0.5, 1.5, .2),
    alpha=0.01,
    ax=ax,
)

这是 R 对 0.7 的效应大小 (es) 的确认

> library(pwr)
> pwr.t.test(power = 0.8, sig.level = 0.01, d = 0.7)

 Two-sample t test power calculation 

          n = 49.35336
          d = 0.7
  sig.level = 0.01
      power = 0.8
  alternative = two.sided

 NOTE: n is number in *each* group