Множественная проверка гипотез

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

Перейти к: навигация, поиск
Зависимость вероятности наличия ложных отклонений от мощности семейства гипотез.
Зависимость вероятности наличия ложных отклонений от мощности семейства гипотез.

Множественная проверка гипотез — класс задач прикладной статистики, возникающий при необходимости одновременного построения семейства статистических выводов. В каждое такое утверждение, представляющее собой вывод об отвержении гипотезы при проверке статистических гипотез, изначально закладывается возможность ошибки (ложного отклонения гипотезы, ошибки первого рода). Вероятность такой ошибки ограничивается неким достаточно малым значением — уровнем значимости \alpha (например, часто берут \alpha=0.05), выражающим степень неуверенности в справедливости заключения. При одновременном построении m заключений это приводит к тому, что вероятность совершения хотя бы одной такой ошибки оказывается ограничена сверху величиной вида 1-(1-\alpha)^m, которая является недопустимо большой уже для достаточно малых m, то есть, с достаточно большой вероятностью рассматриваемое множество выводов содержит ошибочные. Для устранения этого эффекта, называемого эффектом множественных сравнений, разработан ряд статистических техник. Процедура принятия или отклонения гипотез перестраивается таким образом, чтобы некая величина, являющаяся многомерным обобщением ошибки первого рода, была ограничена сверху наперёд заданным числом.

Содержание

Примеры

  • В 1950 году Джозеф Райн проводил эксперименты по выявлению людей с экстрасенсорными способностями. Каждый испытуемый должен был угадать цвет каждой карты в последовательности из 10 карт. В таком опыте вероятность случайно угадать 9 или 10 карт из 10 составляет 11*2-10 ≈ 0.01 — достаточно редкое событие, чтобы признать его маловероятным, поэтому Райн считал его свидетельством того, что испытуемый — экстрасенс. Однако, поскольку в опыте участвовало 1000 человек, вероятность того, что хотя бы один из них случайно угадает 9 или 10 карт составляла 1-(1-11*2-10)1000 ≈ 0.9998. Действительно, из 1000 испытуемых 12 человек угадали 9 из 10 карт, двое — все 10 карт. В последующих экспериментах ни один из них не подтвердил свои способности.
  • Пусть инвестор выбирает, в какие из множества инвестиционных фондов вложить свои деньги. Для оценки инвестиционной привлекательности каждого фонда он вынужден использовать статистический критерий, чтобы отделить фонды, получившие в прошлом прибыль за счёт правильной стратегии управления, от тех, которые получили прибыль только в силу влияния случайных факторов. Чем больше множество фондов, из которых инвестор выбирает, тем больше вероятность того, что какие-то из них покажутся ему потенциально прибыльными, даже если их успех имел не систематический характер.
  • Пусть сравниваются способности нового и уже существующего лекарств влиять на ряд симптомов некоторого заболевания, причём, поскольку интерес представляют выводы отдельно по каждому симптому, проверяются покомпонентные гипотезы превосходства нового лекарства над старым. Чем больше список рассматриваемых симптомов, тем более вероятно, что новый препарат окажется лучше старого хотя бы в отношении одного из симптомов.
  • Пусть проверяется безопасность нового лекарства путём проверки наличия большого количества побочных эффектов. Чем больше побочных эффектов рассматривается, тем более вероятно, что хотя бы один из них будет признан связанным с употреблением лекарства.

Обозначения

Пусть H = \{H_1,\ldots,H_m\} — множество нулевых гипотез, проверяемых против альтернатив общего вида H_{Ai} = \bar{H}_i, \: i=1,\ldots,m. Если нулевая гипотеза верна, будем писать H_i=0, и H_i=1 в противном случае.

За p =\{p_1,\ldots,p_m\} обозначим множество пи-величин, полученных при проверке соответствующих гипотез H_i подходящими статистическими критериями.

Пусть M_0=\{i:\:H_i=0\} и M_1=\{i:\:H_i=1\} — неизвестные множества индексов верных и неверных нулевых гипотез, m_0=\left|M_0\right|, m_1=\left|M_1\right|, \left|M_0\cup M_1\right|=m. Количество отклонённых нулевых гипотез R и количество принятых W = m-R — наблюдаемые случайные величины, в то время как величины S,T,U и V из приведённой ниже таблицы являются ненаблюдаемыми.

Число принятых
гипотез
Число отвергнутых
гипотез
Всего
Число верных
гипотез
U V m_0
Число неверных
гипотез
T S m_1
Всего W R m

Задача состоит в том, чтобы выбрать метод, допускающий минимальное число ложных отклонений гипотез V и ложных принятий T. Как и в одномерном случае, подход заключается в том, чтобы, ограничив меру ошибки первого рода малой величиной \alpha, выбрать процедуру, минимизирующую меру ошибки второго рода, т.е., максимизирующую мощность.

Типы семейств гипотез

Иерархические, неиерархические

Меры ошибки первого рода при множественной проверке

Существует большое количество величин, обобщающих определение ошибки первого рода на задачу множественной проверки гипотез. Наиболее популярны следующие:

  • FWER — familywise error rate, групповая вероятность ошибки (первого рода). \operator{FWER}=\operator{P}\left(V\geq 1\right)
  • k-FWER, gFWER — generalized familywise error rate, обобщённая групповая вероятность ошибки (первого рода). k-\operator{FWER}=\operator{P}\left(V\geq k\right)
  • FDR — false discovery rate, средняя доля ложных отклонений гипотез (среди всех отклонений). \operator{FDR}= \mathrm{E}\!\left [\frac{V}{R} \cdot I \left(R>0\right)\right ]
  • FDP — false discovery proportion, доля ложных отклонений гипотез (среди всех отклонений). \operator{FDP}=\frac{V}{R}\cdot I \left(R>0\right)

Связь между мерами ошибки первого рода

  • \frac{\mathrm{E\left[V\right]}}{m}\leq\operator{FDR}\leq\operator{FWER}\leq\mathrm{E\left[V\right]
  • При m=m_0 (все рассматриваемые гипотезы верны) \operator{FWER}=\operator{FDR}
  • \operator{FDR}= \mathrm{E}\!\left[ \operator{FDP}\right ]

Литература

  1. E. L. Lehmann, J. P. Romano: Testing statistical hypotheses. Springer, New York, 2005 (3rd edition), 786 pages. Chapter 9: Multiple testing and simultaneous inference.

Ссылки

Примечания