Я пытаюсь провести анализ мощности с помощью пакета python statsmodel.Мне нужно ввести размер эффекта, мощность и альфа, но я не уверен, правильно ли то, что у меня есть для размера эффекта.
Скажем, у меня есть исследование по медицине против плацебо, и в результате продолжительность пребывания в часах.Большинство людей используют стандартизированный размер эффекта, такой как 0,2, 0,5, 0,8, но кажется, что нестандартный размер эффекта больше подходит здесь.Что касается моего исследования, мне интересно только узнать, есть ли у тех, кто принимает лекарства, продолжительность пребывания меньше, чем у плацебо, и для меня - клинически значимый эффект, если он составляет 12 часов или меньше.
SD двух групп вместе составляет 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()