Skip to Content

§1. Уточнение модели: от пар к приращениям

Методы этой главы предназначены для выявления неоднородности реализаций выборок X1,…,XnX_1, \dots, X_n и Y1,…,YnY_1, \dots, Y_n одинакового размера, которые нельзя считать независимыми между собой (§2 главы 14).

Приращения

Прежде всего уточним статистическую модель из §1 главы 14 применительно к этой ситуации. Вычислим приращения Zi=Yi−XiZ_i = Y_i - X_i, i=1,…,ni = 1, \dots, n, и разложим каждое на две части:

Zi=θ+εi,Z_i = \theta + \varepsilon_i,

где θ\theta — интересующий нас эффект воздействия, систематический сдвиг (будем считать его положительным), а εi\varepsilon_i — случайная ошибка, в которую входит влияние неучтённых факторов на ZiZ_i.

По кирпичикам: YiY_i и XiX_i одного человека тесно связаны — у гипертоника давление высокое и до, и после. Но в разности Yi−XiY_i - X_i «личные особенности» человека в основном сокращаются, и остаётся то, что сделало лекарство, плюс шум.

Смысл по-пацански: не сравнивай всех «до» со всеми «после» — сравни каждого с самим собой. Тогда разница между людьми перестаёт мешать.

Допущение Д6

В дополнение к допущениям Д1–Д3 из §1 главы 14 предположим, что выполнено условие

  • Д6. Случайные величины ε1,…,εn\varepsilon_1, \dots, \varepsilon_n независимы и имеют непрерывные (вообще говоря, разные) распределения такие, что
P(εi≤0)=P(εi≥0)=1/2,i=1,…,n.\mathbf P(\varepsilon_i \le 0) = \mathbf P(\varepsilon_i \ge 0) = 1/2, \qquad i = 1, \dots, n.

Это означает, что медианы распределений ошибок εi\varepsilon_i равны нулю (§2 главы 7). Разные распределения — это нормально: у одного человека шум больше, у другого меньше. Важно лишь, что шум в среднем по вероятности не сдвигает ни вверх, ни вниз.

Почему допущений главы 14 мало

Замечание 1. Предположения Д1–Д3 из §1 главы 14 не обеспечивают одинаковой распределённости εi\varepsilon_i и даже их независимости. Пусть X1X_1 и X2X_2 независимы и распределены по закону N(0,1)\mathcal N(0, 1). Положим Y1=X1+X2Y_1 = X_1 + X_2 и Y2=X1−X2Y_2 = X_1 - X_2. Нетрудно проверить, что Y1Y_1 и Y2Y_2 распределены по закону N(0,2)\mathcal N(0, 2) и независимы, так как cov(Y1,Y2)=0\mathrm{cov}(Y_1, Y_2) = 0 (П9). Но

Z1=Y1−X1=X2∼N(0,1),Z2=Y2−X2=X1−2X2∼N(0,5).Z_1 = Y_1 - X_1 = X_2 \sim \mathcal N(0, 1), \qquad Z_2 = Y_2 - X_2 = X_1 - 2X_2 \sim \mathcal N(0, 5).

Отсюда ε1=Z1−θ∼N(−θ,1)\varepsilon_1 = Z_1 - \theta \sim \mathcal N(-\theta, 1), а ε2=Z2−θ∼N(−θ,5)\varepsilon_2 = Z_2 - \theta \sim \mathcal N(-\theta, 5). Мало того, что у ε1\varepsilon_1 и ε2\varepsilon_2 разные распределения, они ещё и зависимы, то есть нарушается условие Д6:

cov(ε1,ε2)=cov(X2,X1)−2 cov(X2,X2)=−2≠0.\mathrm{cov}(\varepsilon_1, \varepsilon_2) = \mathrm{cov}(X_2, X_1) - 2\,\mathrm{cov}(X_2, X_2) = -2 \ne 0.

Смысл по-пацански: то, что каждая кучка сама по себе «правильная», ещё ничего не говорит о разностях. Условия надо ставить прямо на приращения.

Гипотеза и альтернатива

Итак, пусть выполнено Д6. Рассмотрим задачу проверки гипотезы H0′ ⁣:θ=0H_0'\colon \theta = 0 против альтернативы H3′ ⁣:θ>0H_3'\colon \theta > 0. Штрих указывает, что гипотеза и сдвиговая альтернатива задаются не для пары законов (F,G)(F, G), как в §1 главы 14, а для распределений приращений ZiZ_i. Для решения используем критерии знаков (§2) и знаковых рангов Уилкоксона (§3). Их обобщения для многомерных данных — в §2 главы 23.

Проверь себя

Тест

У 20 водителей измерили время торможения выспавшимися и после бессонной ночи. Что нужно проверять?
Обновлено