Я работаю с некоторым исследовательским пространственным анализом в R с использованием пакета spdep.
Я наткнулся на вариант настройки p- значений локальных показателей пространственной ассоциации (LISA), рассчитанных с помощью localmoran
функции. Согласно документам он нацелен на:
... корректировка значения вероятности для нескольких тестов.
Далее в документации p.adjustSP
я читал, что доступны следующие варианты:
Методы корректировки включают поправку Бонферрони («bonferroni»), в которой значения p умножаются на количество сравнений. Четыре менее консервативных исправления также включены Holm (1979) («holm»), Hochberg (1988) («hochberg»), Hommel (1988) («hommel») и Benjamini & Hochberg (1995). ('"fdr"') соответственно. Опция сквозного доступа («none») также включена.
Первые четыре метода разработаны, чтобы дать сильный контроль над частотой ошибок по семейным обстоятельствам. Кажется, нет причин использовать немодифицированную поправку Бонферрони, потому что в ней преобладает метод Холма, который также действителен при произвольных допущениях.
Методы Хохберга и Хоммеля действительны, когда тесты гипотез независимы или когда они неотрицательно связаны (Sarkar, 1998; Sarkar and Chang, 1997). Метод Хоммеля является более мощным, чем метод Хохберга, но разница обычно невелика, и р-значения Хохберга вычисляются быстрее.
Методы «BH» (он же «fdr») и «BY» Бенджамини, Хохберга и Екутиэли контролируют уровень ложных открытий, ожидаемую долю ложных открытий среди отвергнутых гипотез. Частота ложных обнаружений является менее строгим условием, чем частота семейных ошибок, поэтому эти методы более эффективны, чем другие.
Пара вопросов, которые появились:
- Проще говоря - какова цель этой корректировки?
- Нужно ли использовать такие исправления?
- Если да - как выбрать из доступных вариантов?