Просмотр

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

У вас нет прав на выполнение этой операции по следующим причинам:

• Запрошенное действие могут выполнять только участники из групп emailconfirmed, Администраторы.
• Вы должны подтвердить ваш адрес электронной почты перед правкой страниц. Пожалуйста, введите и подтвердите ваш адрес эл. почты на странице настроек.

Вы можете просмотреть и скопировать исходный текст этой страницы:

{{TOCright}} '''Теорема Неймана — Пирсона''' (также '''лемма Неймана — Пирсона''') — фундаментальный результат [[Проверка статистических гипотез|теории проверки статистических гипотез]], определяющий наиболее мощный критерий заданного уровня значимости для различения двух простых статистических гипотез. Критерий сравнивает [[Отношение правдоподобия|отношение правдоподобия]] альтернативной и нулевой моделей с порогом: нулевая гипотеза отвергается для наблюдений, которые относительно более вероятны при альтернативе.<ref name="lehmann">{{книга |автор=Lehmann E. L., Romano J. P. |заглавие=Testing Statistical Hypotheses |место=New York |издательство=Springer |год=2005 |издание=3rd ed. |isbn=978-0-387-98864-1}}</ref> Теорема формализует компромисс между ошибками первого и второго рода. Если допустимая вероятность ложного отклонения нулевой гипотезы заранее ограничена, то отношение правдоподобия задаёт область отклонения, для которой вероятность обнаружить данную альтернативу максимальна. Результат был получен Ежи Нейманом и Эгоном Пирсоном и опубликован в 1933 году.<ref name="neyman_pearson">{{статья |автор=Neyman J., Pearson E. S. |заглавие=IX. On the Problem of the Most Efficient Tests of Statistical Hypotheses |издание=Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character |год=1933 |том=231 |страницы=289–337 |doi=10.1098/rsta.1933.0009}}</ref> == Постановка задачи == Пусть наблюдение <tex>X</tex> принимает значения в измеримом пространстве. Рассматриваются две простые гипотезы :<tex>H_0:\;X\sim P_0,</tex> :<tex>H_1:\;X\sim P_1.</tex> Гипотеза называется простой, если она полностью задаёт распределение наблюдения. При каждой из гипотез не остаётся неизвестных параметров, влияющих на распределение. Составная гипотеза, напротив, соответствует семейству распределений. Статистический критерий можно задать измеримой функцией :<tex>\varphi(x)\in[0,1].</tex> Значение <tex>\varphi(x)</tex> интерпретируется как условная вероятность отклонить <tex>H_0</tex> после получения результата <tex>x</tex>. Обычный, нерандомизированный критерий принимает только значения <tex>0</tex> и <tex>1</tex>. Промежуточные значения описывают [[Рандомизированный статистический критерий|рандомизированный критерий]]: окончательное решение дополнительно использует независимый случайный механизм. Вероятность ошибки первого рода и мощность критерия при заданной альтернативе равны соответственно :<tex>\alpha(\varphi)={\rm E}_0\varphi(X),</tex> :<tex>\pi_1(\varphi)={\rm E}_1\varphi(X).</tex> Здесь <tex>{\rm E}_0</tex> и <tex>{\rm E}_1</tex> обозначают математические ожидания при распределениях <tex>P_0</tex> и <tex>P_1</tex>. Ошибка первого рода состоит в отклонении истинной нулевой гипотезы. Ошибка второго рода имеет вероятность <tex>1-\pi_1(\varphi)</tex> и состоит в сохранении <tex>H_0</tex>, когда верна <tex>H_1</tex>. Критерий имеет уровень значимости не выше <tex>\alpha</tex>, если :<tex>{\rm E}_0\varphi(X)\leq\alpha.</tex> Требуется найти среди всех таких критериев критерий с наибольшей мощностью. Иными словами, вероятность ошибки первого рода служит ограниченным ресурсом, а вероятность обнаружения конкретной альтернативы — оптимизируемой величиной. == Формулировка теоремы == Предположим, что распределения <tex>P_0</tex> и <tex>P_1</tex> имеют плотности <tex>p_0</tex> и <tex>p_1</tex> относительно общей меры <tex>\mu</tex>. Плотностью здесь может быть как плотность относительно меры Лебега, так и функция вероятностей относительно считающей меры. Для пары вероятностных мер общую доминирующую меру всегда можно выбрать как <tex>P_0+P_1</tex>. Для любого уровня <tex>0<\alpha<1</tex> существуют число <tex>k\geq0</tex> и вероятность рандомизации <tex>0\leq\gamma\leq1</tex>, при которых критерий <tex>\varphi^*</tex> имеет следующий вид:<ref name="lehmann" /> На множестве <tex>p_1(x)>k p_0(x)</tex> полагается :<tex>\varphi^*(x)=1.</tex> На множестве <tex>p_1(x)<k p_0(x)</tex> полагается :<tex>\varphi^*(x)=0.</tex> На границе <tex>p_1(x)=k p_0(x)</tex> полагается :<tex>\varphi^*(x)=\gamma,</tex> причём <tex>k</tex> и <tex>\gamma</tex> можно выбрать так, чтобы :<tex>{\rm E}_0\varphi^*(X)=\alpha.</tex> Тогда для любого другого критерия <tex>\varphi</tex> уровня не выше <tex>\alpha</tex> выполняется :<tex>{\rm E}_1\varphi(X)\leq{\rm E}_1\varphi^*(X).</tex> Следовательно, <tex>\varphi^*</tex> является ''наиболее мощным критерием уровня <tex>\alpha</tex>'' для проверки <tex>H_0</tex> против <tex>H_1</tex>. При естественных условиях верно и обратное утверждение: всякий наиболее мощный критерий должен иметь указанную структуру, за исключением множеств нулевой вероятности и возможного выбора решения на границе.<ref name="neyman_pearson" /> Если <tex>p_0(x)>0</tex>, правило удобно записывать через отношение правдоподобия :<tex>L(x)={p_1(x)\over p_0(x)}.</tex> Большие значения <tex>L(x)</tex> свидетельствуют в пользу <tex>H_1</tex>, поэтому нулевая гипотеза отвергается при <tex>L(x)>k</tex>. На множестве, где <tex>p_0(x)=0</tex>, но <tex>p_1(x)>0</tex>, наблюдение невозможно при нулевой гипотезе и критерий всегда отвергает её. Точки, в которых обе плотности равны нулю, на вероятности ошибок не влияют. В литературе отношением правдоподобия иногда называют обратную величину <tex>p_0(x)/p_1(x)</tex>. При таком соглашении нулевая гипотеза отвергается при малых, а не при больших значениях статистики. Различие относится только к направлению записи правила. == Обоснование оптимальности == Доказательство опирается на то, что критерий <tex>\varphi^*</tex> принимает решение в соответствии со знаком разности <tex>p_1-kp_0</tex>. Для любого другого критерия <tex>\varphi</tex> поточечно выполняется :<tex>\bigl(\varphi^*(x)-\varphi(x)\bigr)\bigl(p_1(x)-k p_0(x)\bigr)\geq0.</tex> Там, где первая плотность в сравнении больше, критерий <tex>\varphi^*</tex> принимает максимально возможное значение; там, где она меньше, — минимально возможное. На границе второй множитель равен нулю. Интегрирование по мере <tex>\mu</tex> даёт :<tex>{\rm E}_1\bigl(\varphi^*(X)-\varphi(X)\bigr)-k{\rm E}_0\bigl(\varphi^*(X)-\varphi(X)\bigr)\geq0.</tex> Если уровень <tex>\varphi</tex> не превосходит <tex>\alpha</tex>, а размер <tex>\varphi^*</tex> равен <tex>\alpha</tex>, то :<tex>{\rm E}_0\bigl(\varphi^*(X)-\varphi(X)\bigr)\geq0.</tex> Поскольку <tex>k\geq0</tex>, отсюда следует :<tex>{\rm E}_1\varphi^*(X)-{\rm E}_1\varphi(X)\geq0.</tex> Таким образом, никакое иное перераспределение области отклонения при том же ограничении на ошибку первого рода не может увеличить мощность против заданного распределения <tex>P_1</tex>. Доказательство одновременно объясняет вид оптимального правила: в область отклонения в первую очередь включаются наблюдения с наибольшей относительной правдоподобностью альтернативы. Если множество <tex>p_1=kp_0</tex> имеет нулевую вероятность и нет иных участков, на которых перестановка решений не меняет обе вероятности, наиболее мощный критерий единственен почти наверное. При положительной вероятности границы возможны эквивалентные способы рандомизации. == Выбор порога и рандомизация == Порог определяется распределением отношения правдоподобия при нулевой гипотезе. Для критерия точного уровня должно выполняться :<tex>P_0\{L(X)>k\}+\gamma P_0\{L(X)=k\}=\alpha.</tex> Сначала выбирается порог <tex>k</tex>, при котором вероятность строгого превышения не больше <tex>\alpha</tex>, а вероятность нестрогого превышения не меньше <tex>\alpha</tex>. Затем значение <tex>\gamma</tex> заполняет возможный разрыв. Построение критерия эквивалентно упорядочению исходов по величине <tex>L</tex> и включению исходов с наибольшими значениями до исчерпания допустимой вероятности ошибки первого рода. Для непрерывного распределения <tex>L(X)</tex> при <tex>H_0</tex> вероятность равенства порогу обычно равна нулю, и критерий можно сделать нерандомизированным. В дискретных моделях распределение статистики имеет скачки. Тогда критерий без рандомизации нередко достигает лишь уровня, строго меньшего <tex>\alpha</tex>; рандомизация позволяет получить точный уровень и максимальную мощность в полном классе критериев. В приложениях случайное решение может быть нежелательно, поэтому часто выбирают консервативный нерандомизированный критерий и указывают его фактический размер. === Дискретный пример === Пусть имеется одно наблюдение <tex>X</tex> из [[Распределение Бернулли|распределения Бернулли]], а гипотезы имеют вид :<tex>H_0:\;P_0\{X=1\}={1\over5},\qquad H_1:\;P_1\{X=1\}={3\over5}.</tex> Отношение правдоподобия равно <tex>3</tex> при <tex>X=1</tex> и <tex>1/2</tex> при <tex>X=0</tex>. Поэтому исход <tex>X=1</tex> должен входить в область отклонения первым. При требуемом уровне <tex>\alpha=1/10</tex> безусловное отклонение при <tex>X=1</tex> дало бы размер <tex>1/5</tex>, превышающий ограничение. Оптимальный критерий отвергает <tex>H_0</tex> при этом исходе с вероятностью <tex>1/2</tex> и никогда не отвергает её при <tex>X=0</tex>. Его размер равен <tex>1/10</tex>, а мощность — <tex>3/10</tex>. == Пример: проверка среднего нормального распределения == Пусть <tex>X_1,\ldots,X_n</tex> — независимые наблюдения из [[Нормальное распределение|нормального распределения]] со средним <tex>\mu</tex> и известной дисперсией <tex>\sigma^2</tex>. Рассматриваются простые гипотезы :<tex>H_0:\;\mu=\mu_0,\qquad H_1:\;\mu=\mu_1,</tex> где <tex>\mu_1>\mu_0</tex>. Логарифм отношения правдоподобия равен :<tex>\log L(X_1,\ldots,X_n)={\mu_1-\mu_0\over\sigma^2}\sum_{i=1}^n X_i-{n(\mu_1^2-\mu_0^2)\over2\sigma^2}.</tex> При фиксированных параметрах это выражение строго возрастает с выборочным средним <tex>\overline X</tex>. Поэтому критерий Неймана — Пирсона отвергает нулевую гипотезу при больших значениях выборочного среднего. Порог уровня <tex>\alpha</tex> имеет вид :<tex>\overline X>\mu_0+{\sigma\over\sqrt n}z_{1-\alpha},</tex> где <tex>z_{1-\alpha}</tex> — квантиль стандартного нормального распределения, определяемый равенством :<tex>\Phi(z_{1-\alpha})=1-\alpha.</tex> Мощность этого критерия при указанной альтернативе составляет :<tex>\pi(\mu_1)=1-\Phi\left(z_{1-\alpha}-{\sqrt n(\mu_1-\mu_0)\over\sigma}\right).</tex> Формула показывает, что мощность возрастает с объёмом выборки и стандартизованным расстоянием между средними. Критическая граница не зависит от конкретного значения <tex>\mu_1</tex>, если <tex>\mu_1>\mu_0</tex>. Поэтому тот же критерий является равномерно наиболее мощным для односторонней составной альтернативы <tex>\mu>\mu_0</tex>.<ref name="casella_berger">{{книга |автор=Casella G., Berger R. L. |заглавие=Statistical Inference |место=Pacific Grove |издательство=Duxbury Press |год=2002 |издание=2nd ed. |isbn=978-0-534-24312-8}}</ref> Это свойство связано с монотонностью отношения правдоподобия и не является общим для составных гипотез. == Составные гипотезы и связанные критерии == === Равномерно наиболее мощные критерии === Теорема непосредственно сравнивает только два полностью заданных распределения. Если альтернатива содержит несколько распределений, наиболее мощные критерии для разных её точек могут иметь разные области отклонения. Критерий называется [[Равномерно наиболее мощный критерий|равномерно наиболее мощным]], если при соблюдении ограничения на размер он обладает не меньшей мощностью, чем любой конкурирующий критерий, одновременно для каждого распределения из альтернативы. Равномерно наиболее мощный критерий существует не всегда. В однопараметрических семействах с монотонным отношением правдоподобия критерии Неймана — Пирсона для односторонних альтернатив часто имеют общую критическую область. Это приводит к теореме Карлина — Рубина и объясняет оптимальность многих односторонних критериев для экспоненциальных семейств.<ref name="lehmann" /> Для двусторонних альтернатив единая область, оптимальная во всех направлениях, обычно отсутствует; тогда вводятся дополнительные требования, например несмещённость или инвариантность критерия. === Обобщённый критерий отношения правдоподобия === Для параметрических составных гипотез широко применяется [[Критерий отношения правдоподобия|обобщённый критерий отношения правдоподобия]]. Если нулевая гипотеза задаётся множеством параметров <tex>\Theta_0</tex>, а полная модель — множеством <tex>\Theta</tex>, его статистика обычно записывается как :<tex>\lambda(x)={\sup_{\theta\in\Theta_0}p_\theta(x)\over\sup_{\theta\in\Theta}p_\theta(x)}.</tex> Нулевая гипотеза отвергается при малых значениях <tex>\lambda</tex>. Если обе гипотезы простые, это правило сводится к пороговому сравнению обычного отношения правдоподобия. Однако для составных гипотез оптимальность обобщённого критерия не следует из теоремы Неймана — Пирсона. Его обоснование может опираться на специальные свойства модели, инвариантность или асимптотическое распределение статистики.<ref name="casella_berger" /> При наличии мешающих параметров ограничение на ошибку первого рода должно выполняться для всех допустимых значений этих параметров либо в явно заданном усреднённом смысле. Простая подстановка оценок вместо неизвестных величин не сохраняет точный уровень автоматически. == Статистическая интерпретация == === Уровень, мощность и наблюдаемое решение === Уровень значимости — характеристика процедуры при многократном повторении эксперимента под нулевой гипотезой, а не вероятность истинности <tex>H_0</tex> после получения данных. Теорема не приписывает гипотезам апостериорных вероятностей. Она отвечает на условный вопрос: какая процедура лучше всего обнаруживает конкретную альтернативу при заранее ограниченной частоте ложных срабатываний. [[P-значение|Достигаемый уровень значимости]] позволяет представить семейство вложенных критериев, но не заменяет альтернативную модель в утверждении теоремы. Два критерия одинакового уровня могут иметь существенно различную мощность; для её сравнения необходимо указать распределение при альтернативе. === Байесовское правило решений === В [[Байесовская статистика|байесовской постановке]] оптимальное решение для двух простых моделей также часто имеет вид порогового сравнения отношения правдоподобия. Причина порога иная: он определяется априорными вероятностями гипотез и потерями от различных решений. В постановке Неймана — Пирсона априорные вероятности не требуются, а порог выбирается из ограничения на ошибку первого рода. Совпадение математической формы не означает тождества критериев оптимальности. == Связь с классификацией и машинным обучением == Двоичную [[Классификация|классификацию]] можно представить как различение двух распределений признаков. Пусть <tex>p_0(x)</tex> и <tex>p_1(x)</tex> — условные плотности признакового описания для отрицательного и положительного классов. Для решающего правила <tex>\delta</tex> доля ложноположительных решений и чувствительность в генеральной совокупности равны :<tex>{\rm FPR}(\delta)=P_0\{\delta(X)=1\},</tex> :<tex>{\rm TPR}(\delta)=P_1\{\delta(X)=1\}.</tex> Теорема Неймана — Пирсона утверждает, что пороговое правило по <tex>p_1(x)/p_0(x)</tex> максимизирует <tex>{\rm TPR}</tex> среди всех правил с заданным верхним ограничением на <tex>{\rm FPR}</tex>. При изменении порога такие правила, с возможной рандомизацией на границе, образуют верхнюю границу достижимой [[ROC-кривая|ROC-кривой]] для фиксированной пары распределений.<ref name="neyman_pearson" /> Эта интерпретация применяется в задачах обнаружения сигналов, медицинской диагностики, фильтрации событий и выявления аномалий, где стоимость ложного срабатывания задаёт явное ограничение. В практическом [[Машинное обучение|машинном обучении]] плотности классов обычно неизвестны. Отношение правдоподобия заменяют оценкой по обучающей выборке либо используют классификатор, оценка которого монотонно связана с этим отношением. Порог затем калибруют по отдельным данным. Теорема описывает оптимальное правило для истинных распределений, но не гарантирует точного контроля частоты ошибок после оценивания модели на конечной выборке. Для статистически обоснованного контроля необходимы независимая проверочная выборка, учёт неопределённости оценивания или специальные методы конечновыборочной калибровки. Результат также разделяет качество ранжирования и выбор рабочего порога. Обученная оценочная функция может хорошо упорядочивать объекты, однако требуемый уровень ложноположительных решений достигается только после калибровки порога относительно распределения отрицательного класса. == Ограничения и типичные ошибки применения == * '''Простота гипотез.''' Точная оптимальность относится к одной фиксированной нулевой и одной фиксированной альтернативной модели. Максимизация правдоподобия по множествам параметров не превращает составную задачу в прямое следствие теоремы. * '''Зависимость от модели.''' Неверно заданные плотности могут дать область отклонения, не обладающую ни заявленным уровнем, ни наибольшей мощностью для фактических данных. * '''Смешение размера и мощности.''' Ограничение вероятности ошибки первого рода не означает одновременной минимизации обеих ошибок. Оптимальность формулируется при фиксированном ограничении и для конкретной альтернативы. * '''Обратное отношение.''' При статистике <tex>p_1/p_0</tex> отклонению соответствуют большие значения, а при <tex>p_0/p_1</tex> — малые. Использование неверного направления меняет критическую область. * '''Игнорирование дискретности.''' В дискретной модели точный заданный уровень может требовать рандомизации. Нерандомизированный критерий следует характеризовать его фактическим размером. * '''Настройка по тем же данным.''' Выбор модели, статистики или порога после просмотра проверочных данных способен увеличить действительную вероятность ложного отклонения. Уровень сохраняется только для процедуры, случайность которой полностью учтена при калибровке. * '''Множественные проверки.''' Теорема рассматривает одно решение между двумя гипотезами и сама по себе не контролирует семейную вероятность ошибки или долю ложных открытий при большом числе проверок. Теорема Неймана — Пирсона служит эталоном оптимальности статистического решения: при полностью заданных моделях она даёт исчерпывающий ответ, но при переходе к оценённым, составным или адаптивно выбранным моделям требуются дополнительные аргументы. == См. также == * [[Проверка статистических гипотез]] * [[Отношение правдоподобия]] * [[Критерий отношения правдоподобия]] * [[Мощность статистического критерия]] * [[Равномерно наиболее мощный критерий]] * [[Уровень значимости]] * [[Ошибка первого рода]] * [[ROC-кривая]] == Примечания == {{примечания}} == Литература == * {{статья |автор=Neyman J., Pearson E. S. |заглавие=IX. On the Problem of the Most Efficient Tests of Statistical Hypotheses |издание=Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character |год=1933 |том=231 |страницы=289–337 |doi=10.1098/rsta.1933.0009}} * {{книга |автор=Lehmann E. L., Romano J. P. |заглавие=Testing Statistical Hypotheses |место=New York |издательство=Springer |год=2005 |издание=3rd ed. |isbn=978-0-387-98864-1}} * {{книга |автор=Casella G., Berger R. L. |заглавие=Statistical Inference |место=Pacific Grove |издательство=Duxbury Press |год=2002 |издание=2nd ed. |isbn=978-0-534-24312-8}}

• Шаблон:TOCright
• Шаблон:Книга
• Шаблон:Переводы
• Шаблон:Примечания
• Шаблон:Статья

Возврат к странице Теорема Неймана-Пирсона.