Статистический анализ числовых величин (непараметрическая статистика)

Хотя с точки зрения теории математической статистики вероятность совпадения двух элементов выборок равна 0, в реальных выборках экономических данных совпадения встречаются. Так, в рассматриваемых выборках, как видно из табл.1, два раза повторяется величина 2, два раза - величина 7 и три раза - величина 15. В таких случаях говорят о наличии "связанных рангов", а соответствующим совпадающим величинам приписывают среднее арифметическое тех рангов которые они занимают. Так, величины 2 и 2 занимают в объединенной выборке места 3 и 4, поэтому им приписывается ранг . Величины 7 и 7 занимают в объединенной выборке места 8 и 9, поэтому им приписывается ранг . Величины 15, 15 и 15 занимают в объединенной выборке места 13, 14 и 15, поэтому им приписывается ранг .

Следующий шаг - подсчет значения статистики Вилкоксона, т.е. суммы рангов элементов первой выборки

Подсчитаем также сумму рангов элементов второй выборки

Величина может быть использована для контроля вычислений. Дело в том, что суммы рангов элементов первой выборки и второй выборки вместе составляют сумму рангов объединенной выборки, т.е. сумму всех натуральных чисел от 1 до . Следовательно,

В соответствии с ранее проведенными расчетами . Необходимое условие правильности расчетов выполнено. Ясно, что справедливость этого условия не гарантирует правильности расчетов.

Перейдем к расчету статистики Т. Согласно формуле (3)

Следовательно,

Поскольку , то гипотеза однородности принимается на уровне значимости0,05.

Что будет, если поменять выборки местами, вторую назвать первой? Тогда вместо надо рассматривать Имеем

Таким образом, значения статистики критерия отличаются только знаком (можно показать, что это утверждение верно всегда). Поскольку в правиле принятия решения используется только абсолютная величина статистики, то принимаемое решение не зависит от того, какую выборку считаем первой, а какую второй. Для уменьшения объема таблиц принято считать первой выборку меньшего объема.

Продолжим обсуждение критерия Вилкоксона. Правила принятия решений и таблица критических значений для критерия Вилкоксона строятся в предположении справедливости гипотезы полной однородности, описываемой формулой (2). А что будет, если эта гипотеза неверна? Другими словами, какова мощность критерия Вилкоксона?

Пусть объемы выборок достаточно велики, так что можно пользоваться асимптотической нормальностью статистики Вилкоксона. Тогда в соответствии с формулами (1) статистика T будет асимптотически нормальна с параметрами

( 5)

Из формул (5) видно большое значение гипотезы

( 6)

Если эта гипотеза неверна, то, поскольку , справедлива оценка

а потому безгранично растет при росте объемов выборок. В то же время, поскольку

то

( 7)

Следовательно, вероятность отклонения гипотезы H01 , когда она неверна, т.е. мощность критерия Вилкоксона как критерия проверки гипотезы (6), стремится к 1 при возрастании объемов выборок, т.е. критерий Вилкоксона является состоятельным для этой гипотезы при альтернативе

( 8)

Если же гипотеза (6) верна, то статистика T асимптотически нормальна с математическим ожиданием 0 и дисперсией, определяемой формулой

( 9)

Гипотеза (6) является сложной, дисперсия (9), как показывают приводимые ниже примеры, в зависимости от значений и может быть как больше 1, так и меньше 1, но согласно неравенству (7) никогда не превосходит 12.

Приведем пример двух функций распределения и таких, что гипотеза (6) выполнена, а гипотеза (2) - нет. Поскольку

и в случае справедливости гипотезы (2), то для выполнения условия (6) необходимо и достаточно, чтобы

( 10)

а потому естественно в качестве рассмотреть функцию равномерного распределения на интервале (-1 ; 1). Тогда формула (11) переходит в условие

( 11)

Это условие выполняется, если функция является нечетной.

Пример 2. Пусть функции распределения и сосредоточены на интервале , на котором

Тогда

Условие (11) выполнено, поскольку функция является нечетной. Следовательно, Начнем с вычисления

Поскольку

то

С помощью замены переменных получаем, что

В правой части последнего равенства стоят табличные интегралы (см., например, справочник [14, с.71]. Проведя соответствующие вычисления, получаем, что в правой части стоит . Следовательно,

Перейдем к вычислению . Поскольку

то

С помощью замены переменных переходим к табличным интегралам (см., например, справочник [14, с.65]):

Проведя необходимые вычисления, получим, что

Следовательно, для рассматриваемых функций распределения нормированная и центрированная статистика Вилкоксона (см. формулу (4)) асимптотически нормальна с математическим ожиданием 0 и дисперсией (см. формулу (9))

Как легко видеть, дисперсия всегда меньше 1. Это значит, что в рассматриваемом случае гипотеза полной однородности (2) при проверке с помощью критерия Вилкоксона будет приниматься чаще, чем если она на самом деле верна.

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

Обсудим теперь, действительно ли критерий Вилкоксона нацелен на проверку равенства медиан распределений, соответствующих выборкам.

Пример 3. Построим семейство пар функций распределения и таких, что их медианы различны, но для и выполнена гипотеза (6). Пусть распределения сосредоточены на интервале , и на нем а имеет кусочно-линейный график с вершинами в точках . Следовательно,

при ;

на ;

на ;

на ;

при .

Очевидно, что медиана равна , а медиана равна 1/2 .

Согласно соотношению (9) для выполнения гипотезы (6) достаточно определить как функцию из условия

Вычисления дают

Учитывая, что лежит между и 1, не совпадая ни с тем, ни с другим, получаем ограничения на , а именно, Итак, построено искомое семейство пар функций распределения.