使用statsmodel进行功率分析

时间:2019-04-22 13:57:40

标签: python statistics statsmodels

我正在尝试使用python statsmodel软件包进行功率分析。我需要输入效果大小,功效和Alpha,但是我不确定效果大小是否正确。

说我对药物与安慰剂进行了一项研究,结果是停留时间(以小时为单位)。大多数人使用标准化效果大小,例如0.2、0.5、0.8,但似乎非标准化效果大小在这里更合适。对于我的研究,我只想看看接受药物治疗的患者的住院时间是否比安慰剂短,并且对我而言,如果住院时间小于或等于12小时,则对临床意义重大。

两组的标准差为34.9小时。

那么效果大小将等于12 / 34.9吗?

from matplotlib import pyplot
from statsmodels.stats.power import TTestIndPower
# parameters for power analysis
my_effect_size = 12/34.9
effect_sizes = array([my_effect_size, 0.2, 0.5, 0.8])
sample_sizes = array(range(5, 150))
# calculate power curves from multiple power analyses
analysis = TTestIndPower()
analysis.plot_power(dep_var='nobs', nobs=sample_sizes, effect_size=effect_sizes)
pyplot.show()

0 个答案:

没有答案
相关问题