Метод Холма

Материал из MachineLearning.

(Различия между версиями)
Перейти к: навигация, поиск
(Пример)
Строка 17: Строка 17:
1. <tex>0.005 \cdot (4 - 1 + 1) < 0.05 \qquad \Rightarrow</tex> отклоняем 4-ю нулевую гипотезу.
1. <tex>0.005 \cdot (4 - 1 + 1) < 0.05 \qquad \Rightarrow</tex> отклоняем 4-ю нулевую гипотезу.
2. <tex>0.01 \cdot (4 - 2 + 1) < 0.05 \qquad \Rightarrow</tex> отклоняем 1-ю нулевую гипотезу.
2. <tex>0.01 \cdot (4 - 2 + 1) < 0.05 \qquad \Rightarrow</tex> отклоняем 1-ю нулевую гипотезу.
-
3. <tex>0.03 \cdot (4 - 3 + 1) >= 0.05 \qquad \Rightarrow</tex> принимаем 3-ю и 4-ю нулевую гипотезы.
+
3. <tex>0.03 \cdot (4 - 3 + 1) >= 0.05 \qquad \Rightarrow</tex> принимаем 3-ю и 2-ю нулевую гипотезы.
== Реализации ==
== Реализации ==

Версия 03:52, 28 января 2014

Метод Холма-Бонферрони (также Метод Холма, Поправка Холма-Бонферрони) — один из методов контроля групповой вероятности ошибки (первого рода). Является равномерно более мощным, чем поправка Бонферрони и решает проблему падения мощности при росте числа гипотез.

Содержание

Определение

Пусть p_{(1)}\leq \ldots \leq p_{(m)} — уровни значимости p_i, упорядоченные по неубыванию, H_{(1)}, \ldots, H_{(m)} — соответствующие p_{(i)} гипотезы. Процедура Холма определена следующим образом.

Шаг 1. Если p_{(1)}\geq\frac{\alpha}{m}, принять гипотезы H_{(1)}, \ldots, H_{(m)} и остановиться. Иначе, если p_{(1)}<\frac{\alpha}{m}, отвергнуть гипотезу H_{(1)} и продолжить проверку оставшихся гипотез на уровне значимости \alpha/(m-1).
Шаг 2. Если p_{(2)}\geq\frac{\alpha}{m-1}, принять гипотезы H_{(2)}, \ldots, H_{(m)} и остановиться. Иначе, если p_{(2)}<\frac{\alpha}{m-1}, отвергнуть гипотезу H_{(2)} и продолжить проверку оставшихся гипотез на уровне значимости \alpha/(m-2).
И т.д.

Процедура обеспечивает \operator{FWER}\leq\alpha при любом характере зависимости между p_i.

Альтернативная постановка

При рассмотрении неравенств, деление может быть заменено на умножение, то есть вместо неравенств вида p_{(i)}<\frac{\alpha}{m-i+1} используются неравенства вида p \cdot (m-i+1)<\alpha.

Пример

Рассмотрим проверку 4-х гипотез при \alpha=0.05. Пусть для них получены p-value: 0.01, 0.04, 0.03 and 0.005. Будут проверены следующие неравенства:

 1. 0.005 \cdot (4 - 1 + 1) < 0.05 \qquad \Rightarrow отклоняем 4-ю нулевую гипотезу.
 2. 0.01 \cdot (4 - 2 + 1) < 0.05 \qquad \Rightarrow отклоняем 1-ю нулевую гипотезу.
 3. 0.03 \cdot (4 - 3 + 1) >= 0.05 \qquad \Rightarrow принимаем 3-ю и 2-ю нулевую гипотезы.

Реализации

  • MATLAB: функция multcompare, вычисляющая поправку Бонферрони, не поддерживает, однако, поправку Холма-Бонферрони. Реализация доступна на MATLAB File Exchange
  • R: функция p.adjust (с параметром method="holm") из стандартного пакета stats позволяет получить модифицированные уровни значимости с учетом поправки Холма-Бонферрони.

Ссылки

См. также

Поправка Бонферрони FWER

Личные инструменты