Знаковый тест Уилкоксона

Не путать с критерием суммы рангов Уилкоксона .

— Знаковый ранговый критерий Уилкоксона это непараметрический ранговый критерий для проверки статистических гипотез, используемый либо для проверки местоположения популяции на основе выборки данных, либо для сравнения местоположения двух популяций с использованием двух совпадающих выборок. Одновыборочная версия служит цели, аналогичной цели одновыборочного Стьюдента t -критерия . Для двух совпадающих выборок это тест парных различий -критерию Стьюдента , подобный парному t (также известному как « t -критерий для совпадающих пар» или « t -критерий для зависимых выборок»). Критерий Уилкоксона является хорошей альтернативой t-критерию, когда нельзя предположить нормальное распределение различий между парными особями. Вместо этого он предполагает более слабую гипотезу о том, что распределение этой разницы симметрично относительно центрального значения, и направлен на проверку того, значительно ли это центральное значение отличается от нуля. Критерий Уилкоксона является более мощной альтернативой критерию знаков , поскольку он учитывает величину различий, но требует умеренно сильного предположения о симметрии.

История

Тест назван в честь Фрэнка Уилкоксона (1892–1965), который в одной статье предложил как его, так и критерий суммы рангов для двух независимых выборок. Этот тест был популяризирован Сидни Сигелом (1956) в его влиятельном учебнике по непараметрической статистике. Сигел использовал символ T для статистики теста, и, следовательно, тест иногда называют Уилкоксона T -тестом .

Процедура испытания

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

Для парного выборочного теста данные состоят из выборок. $(X_{1},Y_{1}),\dots ,(X_{n},Y_{n})$. Каждый образец представляет собой пару измерений. В простейшем случае измерения производятся в интервальной шкале . Затем их можно преобразовать в действительные числа , а парный выборочный тест преобразуется в одновыборочный тест путем замены каждой пары чисел. $(X_{i},Y_{i})$ по его разнице $X_{i}-Y_{i}$. В общем, должна быть возможность ранжировать различия между парами. Для этого требуется, чтобы данные были в упорядоченной метрической шкале, типе шкалы, который несет больше информации, чем порядковая шкала, но может иметь меньше информации, чем интервальная шкала.

Данные для одновыборочного теста представляют собой набор образцов действительных чисел. $X_{1},\dots ,X_{n}$. Предположим для простоты, что выборки имеют разные абсолютные значения и ни одна выборка не равна нулю. (Ноли и связи создают некоторые сложности; см. ниже.) Тест выполняется следующим образом:

1. Вычислить $|X_{1}|,\dots ,|X_{n}|$.
2. Сортировать $|X_{1}|,\dots ,|X_{n}|$и используйте этот отсортированный список для присвоения рангов $R_{1},\dots ,R_{n}$: ранг наименьшего наблюдения равен единице, ранг следующего наименьшего наблюдения равен двум и так далее.
3. Позволять $\operatorname {sgn}$ обозначим знаковую функцию : $\operatorname {sgn}(x)=1$ если $x>0$ и $\operatorname {sgn}(x)=-1$ если $x<0$. Статистика теста представляет собой сумму знакового ранга. $T$: $$T=\sum _{i=1}^{N}\operatorname {sgn}(X_{i})R_{i}.$$
4. Создайте $p$-значение путем сравнения $T$ его распределению при нулевой гипотезе.

Ранги определены так, что $R_{i}$ это количество $j$ для чего $|X_{j}|\leq |X_{i}|$. Кроме того, если $\sigma \colon \{1,\dots ,n\}\to \{1,\dots ,n\}$ таков, что $|X_{\sigma (1)}|<\dots <|X_{\sigma (n)}|$, затем $R_{\sigma (i)}=i$ для всех $i$.

Знаковая сумма $T$ тесно связан с двумя другими статистическими данными испытаний. Сумма положительного ранга $T^{+}$ и сумма отрицательного ранга $T^{-}$ определяются $${\begin{aligned}T^{+}&=\sum _{1\leq i\leq n,\ X_{i}>0}R_{i},\\T^{-}&=\sum _{1\leq i\leq n,\ X_{i}<0}R_{i}.\end{aligned}}$$Потому что $T^{+}+T^{-}$ равно сумме всех рангов, что $1+2+\dots +n=n(n+1)/2$, эти три статистики связаны соотношением: $${\begin{aligned}T^{+}&={\frac {n(n+1)}{2}}-T^{-}={\frac {n(n+1)}{4}}+{\frac {T}{2}},\\T^{-}&={\frac {n(n+1)}{2}}-T^{+}={\frac {n(n+1)}{4}}-{\frac {T}{2}},\\T&=T^{+}-T^{-}=2T^{+}-{\frac {n(n+1)}{2}}={\frac {n(n+1)}{2}}-2T^{-}.\end{aligned}}$$Потому что $T$, $T^{+}$, и $T^{-}$ несут одну и ту же информацию, любой из них может быть использован в качестве тестовой статистики.

Сумма положительного ранга и сумма отрицательного ранга имеют альтернативные интерпретации, которые полезны для теории, лежащей в основе теста. Определите среднее значение Уолша $W_{ij}$ быть ${\tfrac {1}{2}}(X_{i}+X_{j})$. Затем: $${\begin{aligned}T^{+}=\#\{W_{ij}>0\colon 1\leq i\leq j\leq n\},\\T^{-}=\#\{W_{ij}<0\colon 1\leq i\leq j\leq n\}.\end{aligned}}$$

Нулевые и альтернативные гипотезы

Одновыборочный тест

Одновыборочный знаково-ранговый критерий Уилкоксона можно использовать для проверки того, поступают ли данные из симметричной совокупности с указанным центром (который соответствует медиане , среднему и псевдомедиане ). Если известен населенный пункт, его можно использовать для проверки симметричности данных относительно его центра.

Чтобы формально объяснить нулевую и альтернативную гипотезы, предположим, что данные состоят из независимых и одинаково распределенных выборок из распределения $F$. Если $F$ можно считать симметричным, то нулевая и альтернативная гипотезы таковы:

Нулевая гипотеза H ₀

$F$ симметричен относительно $\mu =0$.

Односторонняя альтернативная гипотеза H ₁

$F$ симметричен относительно $\mu <0$.

Односторонняя альтернативная гипотеза H ₂

$F$ симметричен относительно $\mu >0$.

Двусторонняя альтернативная гипотеза H ₃

$F$ симметричен относительно $\mu \neq 0$.

Если вдобавок $\Pr(X=\mu )=0$, затем $\mu$ является медианой $F$. Если эта медиана уникальна, то критерий суммы знаковых рангов Уилкоксона становится проверкой местоположения медианы. Когда среднее значение $F$ определено, то среднее значение равно $\mu$, и тест также является проверкой местоположения среднего значения.

Ограничение на симметричность альтернативного распределения является весьма строгим, но для односторонних тестов его можно ослабить. Скажи это $F$ стохастически меньше , чем распределение, симметричное относительно нуля , если $F$-распределенная случайная величина $X$ удовлетворяет $\Pr(X<-x)\geq \Pr(X>x)$ для всех $x\geq 0$. Сходным образом, $F$ стохастически больше , чем распределение, симметричное относительно нуля , если $\Pr(X<-x)\leq \Pr(X>x)$ для всех $x\geq 0$. Тогда критерий суммы знаковых рангов Уилкоксона также можно использовать для следующих нулевых и альтернативных гипотез:

Нулевая гипотеза H ₀

$F$ симметричен относительно $\mu =0$.

Односторонняя альтернативная гипотеза H ₁

$F$ стохастически меньше, чем распределение, симметричное относительно нуля.

Односторонняя альтернативная гипотеза H ₂

$F$ стохастически больше, чем распределение, симметричное относительно нуля.

Гипотезу о том, что данные являются IID, можно ослабить. Каждая точка данных может быть взята из другого распределения, при условии, что все распределения считаются непрерывными и симметричными относительно общей точки. $\mu _{0}$. Точки данных не обязаны быть независимыми, пока условное распределение каждого наблюдения с учетом остальных симметрично относительно $\mu _{0}$.

Тест парных данных

Поскольку тест парных данных возникает в результате анализа парных разностей, его нулевые и альтернативные гипотезы могут быть получены из гипотез теста одной выборки. В каждом случае они становятся утверждениями о поведении различий. $X_{i}-Y_{i}$.

Позволять $F(x,y)$ быть совместным кумулятивным распределением пар $(X_{i},Y_{i})$. В этом случае нулевая и альтернативная гипотезы таковы:

Нулевая гипотеза H ₀

Наблюдения $X_{i}-Y_{i}$ симметричны относительно $\mu =0$.

Односторонняя альтернативная гипотеза H ₁

Наблюдения $X_{i}-Y_{i}$ симметричны относительно $\mu <0$.

Односторонняя альтернативная гипотеза H ₂

Наблюдения $X_{i}-Y_{i}$ симметричны относительно $\mu >0$.

Двусторонняя альтернативная гипотеза H ₃

Наблюдения $X_{i}-Y_{i}$ симметричны относительно $\mu \neq 0$.

Их также можно выразить более непосредственно через исходные пары:

Нулевая гипотеза H ₀

Наблюдения $(X_{i},Y_{i})$ взаимозаменяемы что , а это означает, $(X_{i},Y_{i})$ и $(Y_{i},X_{i})$ имеют одинаковое распределение. Эквивалентно, $F(x,y)=F(y,x)$.

Односторонняя альтернативная гипотеза H ₁

Для некоторых $\mu <0$, пары $(X_{i},Y_{i})$ и $(Y_{i}+\mu ,X_{i}-\mu )$ имеют одинаковое распределение.

Односторонняя альтернативная гипотеза H ₂

Для некоторых $\mu >0$, пары $(X_{i},Y_{i})$ и $(Y_{i}+\mu ,X_{i}-\mu )$ имеют одинаковое распределение.

Двусторонняя альтернативная гипотеза H ₃

Для некоторых $\mu \neq 0$, пары $(X_{i},Y_{i})$ и $(Y_{i}+\mu ,X_{i}-\mu )$ имеют одинаковое распределение.

Нулевая гипотеза об обменности может возникнуть в результате парного эксперимента с экспериментальной и контрольной группой. Рандомизация лечения и контроля внутри каждой пары делает наблюдения взаимозаменяемыми. Для сменного дистрибутива $X_{i}-Y_{i}$ имеет то же распределение, что и $Y_{i}-X_{i}$, и, следовательно, при нулевой гипотезе распределение симметрично относительно нуля.

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

Нулевая гипотеза H ₀

Наблюдения $(X_{i},Y_{i})$ являются взаимозаменяемыми.

Односторонняя альтернативная гипотеза H ₁

Различия $X_{i}-Y_{i}$ стохастически меньше, чем распределение, симметричное относительно нуля, то есть для каждого $x\geq 0$, $Pr(X_{i}<Y_{i}-x)\geq \Pr(X_{i}>Y_{i}+x)$.

Односторонняя альтернативная гипотеза H ₂

Различия $X_{i}-Y_{i}$ стохастически больше, чем распределение, симметричное относительно нуля, то есть для каждого $x\geq 0$, $Pr(X_{i}<Y_{i}-x)\leq \Pr(X_{i}>Y_{i}+x)$.

Нули и связи

В реальных данных иногда бывает, что есть выборка $X_{i}$ что равно нулю или паре $(X_{i},Y_{i})$ с $X_{i}=Y_{i}$. Также может случиться, что имеются связанные образцы. Это означает, что для некоторых $i\neq j$, у нас есть $X_{i}=X_{j}$ (в случае одной выборки) или $X_{i}-Y_{i}=X_{j}-Y_{j}$ (в случае парной выборки). Это особенно характерно для дискретных данных. В этом случае процедура тестирования, определенная выше, обычно не определена, поскольку не существует способа однозначно ранжировать данные. (Единственное исключение – если имеется один образец $X_{i}$ который равен нулю и не имеет других нулей или связей.) По этой причине необходимо изменить статистику теста.

Нули

В оригинальной статье Уилкоксона не затрагивался вопрос наблюдений (или, в случае парной выборки, различий), равных нулю. Однако в более поздних исследованиях он рекомендовал убрать нули из выборки. Затем к полученным данным можно было применить стандартный критерий знакового ранга, если не было никаких связей. Теперь это называется процедурой сокращенной выборки.

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

0, 2, 3, 4, 6, 7, 8, 9, 11, 14, 15, 17, −18.

Процедура сокращенной выборки удаляет ноль. Оставшимся данным присваиваются знаковые ранги:

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, −12.

Это имеет одностороннее p -значение $55/2^{12}$, и, следовательно, выборка не является достоверно положительной ни на каком уровне значимости $\alpha <55/2^{12}\approx 0.0134$. Пратт утверждает, что можно было бы ожидать, что уменьшение количества наблюдений определенно не приведет к тому, что данные будут выглядеть более позитивными. Однако если нулевое наблюдение уменьшается на величину меньше 2 или если все наблюдения уменьшаются на величину меньше 1, то ранги со знаком становятся:

−1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, −13.

Это имеет одностороннее p -значение $109/2^{13}$. Следовательно, выборка будет считаться значимо положительной на любом уровне значимости. $\alpha >109/2^{13}\approx 0.0133$. Парадокс заключается в том, что если $\alpha$ находится между $109/2^{13}$ и $55/2^{12}$, то уменьшение незначительной выборки приводит к тому, что она становится значительно положительной .

Поэтому Пратт предложил процедуру нулевого знакового ранга. Эта процедура включает нули при ранжировании выборок. Однако он исключает их из тестовой статистики или, что то же самое, определяет $\operatorname {sgn}(0)=0$. Пратт доказал, что процедура нулевого знакового ранга имеет несколько желательных свойств, которых нет в процедуре сокращенной выборки:

1. Увеличение наблюдаемых значений не делает существенно положительную выборку незначительной, а незначительная выборка не делает существенно отрицательной.
2. Если распределение наблюдений симметрично, то значения $\mu$ которые тест не отклоняет из интервала.
3. Выборка является значимо положительной, незначимой или существенно отрицательной тогда и только тогда, когда это так, когда нулям присвоены произвольные ненулевые знаки, тогда и только тогда, когда это так, когда нули заменяются ненулевыми значениями, которые меньше по абсолютной величине, чем любое ненулевое наблюдение.
4. Для фиксированного порога значимости $\alpha$, а для теста, который рандомизирован для получения точного уровня $\alpha$, вероятность назвать набор наблюдений значимо положительным (соответственно существенно отрицательным) является неубывающей (соответственно невозрастающей) функцией наблюдений.

Пратт отмечает, что когда процедура нулевого знакового ранга объединяется с процедурой среднего ранга для разрешения связей, полученный тест является последовательным тестом против альтернативной гипотезы, согласно которой для всех $i\neq j$, $\Pr(X_{i}+X_{j}>0)$ и $\Pr(X_{i}+X_{j}<0)$ отличаются по крайней мере на фиксированную константу, независящую от $i$ и $j$.

Процедура нулевого знакового ранга имеет тот недостаток, что при появлении нулей нулевое распределение тестовой статистики меняется, поэтому таблицы p -значений больше нельзя использовать.

Когда данные представлены по шкале Лайкерта с равноотстоящими друг от друга категориями, процедура нулевого знакового ранга с большей вероятностью сохранит частоту ошибок типа I, чем процедура сокращенной выборки.

С точки зрения статистической эффективности не существует идеального правила обработки нулей. Коновер нашел примеры нулевых и альтернативных гипотез, которые показывают, что ни один из методов Уилкоксона и Пратта не является однозначно лучшим, чем другой. При сравнении дискретного равномерного распределения с распределением, в котором вероятности линейно возрастают слева направо, метод Пратта превосходит метод Уилкоксона. При тестировании биномиального распределения с центром в нуле, чтобы увидеть, является ли параметр каждого испытания Бернулли ${\tfrac {1}{2}}$, метод Уилкоксона превосходит метод Пратта.

Галстуки

Когда данные не имеют связей, ранги $R_{i}$ используются для расчета статистики теста. При наличии связей звания не определяются. Есть два основных подхода к решению этой проблемы.

Наиболее распространенная процедура обработки ничьих, первоначально рекомендованная Уилкоксоном, называется процедурой среднего ранга или процедурой среднего ранга. номера от 1 до n Эта процедура присваивает наблюдениям , причем два наблюдения получают одинаковый номер тогда и только тогда, когда они имеют одинаковое абсолютное значение. Эти числа условно называются рангами, хотя множество этих чисел не равно $\{1,\dots ,n\}$ (кроме случаев отсутствия связей). Ранг, присвоенный наблюдению, представляет собой среднее из возможных рангов, которые оно могло бы иметь, если бы связи были нарушены всеми возможными способами. После присвоения рангов статистика теста рассчитывается так же, как обычно.

Например, предположим, что наблюдения удовлетворяют $$|X_{3}|<|X_{2}|=|X_{5}|<|X_{6}|<|X_{1}|=|X_{4}|=|X_{7}|.$$В этом случае, $X_{3}$ присвоено 1 ранг, $X_{2}$ и $X_{5}$ присваиваются звания $(2+3)/2=2.5$, $X_{6}$ присвоен ранг 4, и $X_{1}$, $X_{4}$, и $X_{7}$ присваиваются звания $(5+6+7)/3=6$. Формально предположим, что существует набор наблюдений, имеющих одинаковое абсолютное значение. $v$, что $k-1$ наблюдения имеют абсолютное значение меньше, чем $v$, и это $\ell$ наблюдения имеют абсолютное значение меньше или равное $v$. Если связи между наблюдениями с абсолютным значением $v$ были нарушены, то эти наблюдения заняли бы ряды $k$ через $\ell$. Таким образом, процедура среднего ранга присваивает им ранг $(k+\ell )/2$.

При использовании процедуры среднего ранга нулевое распределение отличается при наличии связей. Процедура среднего ранга также имеет некоторые недостатки, аналогичные недостаткам процедуры сокращенной выборки для нулей. Вполне возможно, что выборка может быть оценена как значительно положительная с помощью процедуры среднего ранга; но увеличение некоторых значений с целью разрыва связей или разрыв связей каким-либо образом приводит к получению выборки, которую тест считает незначительной. Однако увеличение всех наблюдаемых значений на одинаковую величину не может превратить существенно положительный результат в незначимый, а незначительный - в существенно отрицательный. При этом, если наблюдения распределены симметрично, то значения $\mu$ которые тест не отклоняет из интервала.

Другим распространенным вариантом разрешения ничьих является процедура разрешения ничьих. В процедуре разрешения конфликтов наблюдениям присваиваются отдельные ранги в наборе $\{1,\dots ,n\}$. Ранг, присвоенный наблюдению, зависит от его абсолютного значения и правила разрешения ничьей. Наблюдениям с меньшими абсолютными значениями всегда присваиваются меньшие ранги, как и в стандартном тесте суммы рангов. Правило разрешения конфликтов используется для присвоения рангов наблюдениям с одинаковым абсолютным значением. Одним из преимуществ правил разрешения конфликтов является то, что они позволяют использовать стандартные таблицы для вычисления p -значений.

Случайный тай-брейк разрывает ничью случайным образом. При случайном выявлении связей нулевое распределение такое же, как и при отсутствии связей, но результат теста зависит не только от данных, но и от дополнительных случайных выборов. Усреднение рангов по возможным случайным выборам приводит к процедуре среднего ранга. Можно также сообщить о вероятности отклонения всех случайных выборов. Преимущество случайного определения тай-брейка состоит в том, что вероятность того, что выборка будет признана значимо положительной, не уменьшается при увеличении числа наблюдений. Консервативный подход разрывает связи в пользу нулевой гипотезы. При выполнении одностороннего теста, при котором отрицательные значения $T$ имеют тенденцию быть более значимыми, связи разрываются путем присвоения более низких рангов отрицательным наблюдениям и более высоких - положительным. Когда тест дает положительные значения $T$ существенны, связи разрываются в другую сторону, и когда большие абсолютные значения $T$ значимы, связи разрываются, чтобы сделать $|T|$ как можно меньше. Пратт отмечает, что, когда ничья вероятны, консервативная процедура разрешения ничьей «предположительно имеет низкую эффективность, поскольку она равносильна разрыву всех связей в пользу нулевой гипотезы».

Процедура среднего ранга может не совпадать с процедурой определения ничьей. Пратт приводит следующий пример. Предположим, что наблюдения:

1, 1, 1, 1, 2, 3, −4.

Процедура среднего ранга присваивает им знаковые ранги.

2.5, 2.5, 2.5, 2.5, 5, 6, −7.

Этот образец значительно положителен на одностороннем уровне. $\alpha =14/2^{7}$. С другой стороны, любое правило тай-брейка будет присваивать ранги

1, 2, 3, 4, 5, 6, −7.

На том же одностороннем уровне $\alpha =14/2^{7}$, это не существенно.

Два других варианта решения ничьей основаны на усреднении результатов тай-брейка. В методе средней статистики тестовая статистика $T$ рассчитывается для всех возможных способов разрыва ничьих, а окончательная статистика представляет собой среднее значение статистики разрыва ничьих. В средней вероятности методе значение p вычисляется для каждого возможного способа разрыва связей, а окончательное значение p представляет собой среднее значение p -значений, связанных с ничьей.

Вычисление нулевого распределения

Вычисление значений p требует знания распределения $T$ при нулевой гипотезе. Для этого распределения не существует закрытой формулы. Однако для небольших значений $n$, распределение может быть вычислено точно. При нулевой гипотезе, что данные симметричны относительно нуля, каждый $X_{i}$ с такой же вероятностью будет положительным, как и отрицательным. Поэтому вероятность того, что $T=t$ при нулевой гипотезе равно числу комбинаций знаков, которые дают $T=t$ разделить на количество возможных комбинаций знаков $2^{n}$. Это можно использовать для расчета точного распределения $T$ при нулевой гипотезе.

Вычисление распределения $T$ рассмотрение всех возможностей требует вычислений $2^{n}$ суммы, что невыполнимо для всех, кроме самых маленьких $n$. Однако существует эффективная рекурсия для распределения $T^{+}$. Определять $u_{n}(t^{+})$ — количество комбинаций знаков, для которых $T^{+}=t^{+}$. Это равно числу подмножеств $\{1,\dots ,n\}$ какая сумма $t^{+}$. Базовые случаи рекурсии: $u_{0}(0)=1$, $u_{0}(t^{+})=0$ для всех $t^{+}\neq 0$, и $u_{n}(t^{+})=0$ для всех $t<0$ или $t>n(n+1)/2$. Рекурсивная формула $$u_{n}(t^{+})=u_{n-1}(t^{+})+u_{n-1}(t^{+}-n).$$Формула верна, поскольку каждое подмножество $\{1,\dots ,n\}$ что в сумме равно $t^{+}$ либо не содержит $n$, и в этом случае это также подмножество $\{1,\dots ,n-1\}$, или оно содержит $n$, в этом случае удаление $n$ из подмножества производит подмножество $\{1,\dots ,n-1\}$ что в сумме равно $t^{+}-n$. Согласно нулевой гипотезе, функция массы вероятности $T^{+}$ удовлетворяет $\Pr(T^{+}=t^{+})=u_{n}(t^{+})/2^{n}$. Функция $u_{n}$ тесно связана с целочисленной функцией распределения .

Если $p_{n}(t^{+})$ это вероятность того, что $T^{+}=t^{+}$ при нулевой гипотезе, когда существуют $n$ образцы, затем $p_{n}(t^{+})$ удовлетворяет аналогичной рекурсии: $$2p_{n}(t^{+})=p_{n-1}(t^{+})+p_{n-1}(t^{+}-n)$$с аналогичными граничными условиями. Существует также рекурсивная формула для кумулятивной функции распределения. $\Pr(T^{+}\leq t^{+})$.

Для очень больших $n$, даже приведенная выше рекурсия слишком медленная. В этом случае нулевое распределение может быть аппроксимировано. Нулевые распределения $T$, $T^{+}$, и $T^{-}$ асимптотически нормальны со средними и дисперсиями: $${\begin{aligned}\mathbf {E} [T^{+}]&=\mathbf {E} [T^{-}]={\frac {n(n+1)}{4}},\\\mathbf {E} [T]&=0,\\\operatorname {Var} (T^{+})&=\operatorname {Var} (T^{-})={\frac {n(n+1)(2n+1)}{24}},\\\operatorname {Var} (T)&={\frac {n(n+1)(2n+1)}{6}}.\end{aligned}}$$

Лучшие приближения можно получить, используя расширения Эджворта. Использование разложения Эджворта четвертого порядка показывает, что: $$\Pr(T^{+}\leq k)\approx \Phi (t)+\phi (t){\Big (}{\frac {3n^{2}+3n-1}{10n(n+1)(2n+1)}}{\Big )}(t^{3}-3t),$$где $$t={\frac {k+{\tfrac {1}{2}}-{\frac {n(n+1)}{4}}}{\sqrt {\frac {n(n+1)(2n+1)}{24}}}}.$$Техническая основа этих разложений весьма сложна, поскольку обычные разложения Эджворта применимы к суммам непрерывных случайных величин IID, в то время как $T^{+}$ представляет собой сумму неодинаково распределенных дискретных случайных величин. Однако окончательный результат заключается в том, что приведенное выше расширение имеет ошибку $O(n^{-3/2})$, как обычное расширение Эджворта четвертого порядка.

Производящая функция момента $T$ имеет точную формулу: $$M(t)={\frac {1}{2^{n}}}\prod _{j=1}^{n}(1+e^{jt}).$$

Когда присутствуют нули и используется процедура обнуления знакового ранга, или когда присутствуют связи и используется процедура среднего ранга, нулевое распределение $T$ изменения. Кюртон получил нормальное приближение для этой ситуации. Предположим, что исходное количество наблюдений было $n$ и количество нулей было $z$. Коррекция связи $$c=\sum t^{3}-t,$$где сумма по всем размерам $t$ каждой группы связанных наблюдений. Ожидание $T$ по-прежнему равно нулю, в то время как ожидание $T^{+}$ является $$\mathbf {E} [T^{+}]={\frac {n(n+1)}{4}}-{\frac {z(z+1)}{4}}.$$Если $$\sigma ^{2}={\frac {n(n+1)(2n+1)-z(z+1)(2z+1)-c/2}{6}},$$затем $${\begin{aligned}\operatorname {Var} (T)&=\sigma ^{2},\\\operatorname {Var} (T^{+})&=\sigma ^{2}/4.\end{aligned}}$$

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

Уилкоксон первоначально определил статистику суммы рангов Уилкоксона как $\min(T^{+},T^{-})$. Ранние авторы, такие как Сигел последовал за Уилкоксоном. Это подходит для двусторонних проверок гипотез, но не может использоваться для односторонних проверок.

Вместо присвоения рангов от 1 до n также можно присваивать ранги от 0 до $n-1$. Это так называемые модифицированные ранги . Модифицированная сумма знакового ранга $T_{0}$, модифицированная сумма положительного ранга $T_{0}^{+}$и модифицированная сумма отрицательного ранга $T_{0}^{-}$ определяются аналогично $T$, $T^{+}$, и $T^{-}$ но с измененными званиями вместо обычных званий. Вероятность того, что сумма двух независимых $F$-распределенная случайная величина положительна, может быть оценена как $2T_{0}^{+}/(n(n-1))$. Когда рассмотрение ограничивается непрерывными распределениями, это несмещенная оценка минимальной дисперсии $p_{2}$.

Пример

$i$ $x_{2,i}$ $x_{1,i}$ $x_{2,i}-x_{1,i}$ $\operatorname {sgn}$ ${\text{abs}}$ 1 125 110 1 15 2 115 122  –1 7 3 130 125 1 5 4 140 120 1 20 5 140 140   0 6 115 124  –1 9 7 140 123 1 17 8 125 137  –1 12 9 140 135 1 5 10 135 145  –1 10

упорядочить по абсолютной разнице

$i$ $x_{2,i}$ $x_{1,i}$ $x_{2,i}-x_{1,i}$ $\operatorname {sgn}$ ${\text{abs}}$ $R_{i}$ $\operatorname {sgn} \cdot R_{i}$ 5 140 140   0     3 130 125 1 5 1.5 1.5 9 140 135 1 5 1.5 1.5 2 115 122  –1 7 3  –3 6 115 124  –1 9 4  –4 10 135 145  –1 10 5  –5 8 125 137  –1 12 6  –6 1 125 110 1 15 7 7 7 140 123 1 17 8 8 4 140 120 1 20 9 9

$\operatorname {sgn}$ это знаковая функция , ${\text{abs}}$ является абсолютной величиной , и $R_{i}$ это звание . Обратите внимание, что пары 3 и 9 связаны по абсолютному значению. Им будут присвоены ранги 1 и 2, поэтому каждый получает среднее значение этих рангов — 1,5.

$W=1.5+1.5-3-4-5-6+7+8+9=9$

$|W|<W_{\operatorname {crit} (\alpha =0.05,\ 9{\text{, two-sided}})}=15$

$\therefore {\text{failed to reject }}H_{0}$ что медиана парных разностей отлична от нуля.

The $p$-значение для этого результата $0.6113$

Размер эффекта

Основная статья: Манн-Уитни_U_test § Ранг-бисериальная_корреляция

Чтобы вычислить величину эффекта для теста знакового ранга, можно использовать ранг-бисериальную корреляцию .

Если сообщается тестовая статистика T , ранговая корреляция r равна тестовой статистике T, разделенной на общую сумму рангов S , или r = T / S . В приведенном выше примере статистика теста равна T = 9. Размер выборки, равный 9, имеет общую сумму рангов S = (1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 9) = 45. Следовательно, ранговая корреляция равна 9/45, поэтому r = 0,20.

тестовая статистика T , эквивалентным способом вычисления ранговой корреляции является использование разницы в пропорциях между двумя суммами рангов, что представляет собой формулу простой разности Керби (2014). Если сообщается Продолжая текущий пример, размер выборки равен 9, поэтому общая сумма рангов равна 45. T — меньшая из двух сумм рангов, поэтому T равна 3 + 4 + 5 + 6 = 18. Только из этой информации можно сделать вывод, что Оставшуюся сумму рангов можно вычислить, поскольку это общая сумма S минус T или, в данном случае, 45 - 18 = 27. Далее, две пропорции суммы рангов: 27/45 = 60% и 18/45 = 40%. Наконец, ранговая корреляция представляет собой разницу между двумя пропорциями (0,60 минус 0,40), следовательно, r = 0,20.

Реализации программного обеспечения

• R включает реализацию теста как wilcox.test(x,y, paired=TRUE), где x и y — векторы одинаковой длины.
• ALGLIB включает реализацию знакового рангового теста Уилкоксона на C++, C#, Delphi, Visual Basic и т. д.
• GNU Octave реализует различные односторонние и двусторонние версии теста в wilcoxon_test функция.
• SciPy включает в себя реализацию знакового рангового теста Уилкоксона на Python.
• Accord.NET включает реализацию знакового рангового теста Уилкоксона на C# для приложений .NET.
• MATLAB реализует этот тест, используя «тест суммы рангов Уилкоксона», как [p,h] = signrank(x,y) также возвращает логическое значение, указывающее решение теста. Результат h = 1 указывает на отклонение нулевой гипотезы, а h = 0 указывает на неспособность отвергнуть нулевую гипотезу на уровне значимости 5%.
• Пакет Julia HypothesisTests включает в себя тест знакового ранга Уилкоксона, как value(SignedRankTest(x, y)).
• SAS PROC UNIVARIATE включает ранговый тест со знаком Уилкоксона в заголовках рамок «Тесты на определение местоположения» как «Знаковый ранг». Несмотря на то, что эта процедура рассчитывает S-статистику, а не W-статистику, полученное значение p все равно можно использовать для этого теста.

См. также

• Критерий Манна – Уитни – Вилкоксона
• Знаковый тест

Внешние ссылки

• Знако-ранговый критерий Уилкоксона в R
• Пример использования знаково-рангового критерия Уилкоксона
• Онлайн-версия теста
• Таблица критических значений знаково-рангового критерия Уилкоксона
• Краткое руководство психолога-экспериментатора Карла Л. Вюнша - Непараметрические оценки величины эффекта (авторские права Карла Л. Вюнша, 2015 г.)
• Керби, Д.С. (2014). Простая формула разницы: подход к обучению непараметрической корреляции. Комплексная психология , том 3, статья 1. doi:10.2466/11.IT.3.1. ссылка на статью