Коррекция Бонферрони
Реализуем множественное тестирование гипотез с помощью коррекции Бонферрони, которую мы рассмотрели на слайдах. Для этого воспользуйтесь импортированной функцией multipletests().
Используйте уровень значимости для одиночного теста .05 и посмотрите, как коррекция Бонферрони влияет на подготовленный список p-значений.
Это упражнение является частью курса
Практика статистических вопросов для собеседований на Python
Инструкции к упражнению
- Вычислите список скорректированных по Бонферрони p-значений с помощью импортированной функции
multipletests(). - Выведите результаты множественного тестирования гипотез, хранящиеся в индексе 0 переменной
p_adjusted. - Выведите сами p-значения, хранящиеся в индексе 1 переменной
p_adjusted.
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
from statsmodels.sandbox.stats.multicomp import multipletests
pvals = [.01, .05, .10, .50, .99]
# Create a list of the adjusted p-values
p_adjusted = multipletests(____, alpha=____, method='bonferroni')
# Print the resulting conclusions
print(____)
# Print the adjusted p-values themselves
print(____)