6.3.2. Проверка однородности совокупности дисперсий

 

Для каждого слоя вычисляется несмещенная оценка дисперсии, обозначим эти оценки через μ^*₂(x) , μ^*₂ (y), …, μ^*₂ (w) соответственно. Числа степеней свободы этих оценок

k₁ = п₁–1, k₂ = п₂ – 1, …, k_т = п_w–1.

Гипотеза Н₀ состоит в том, что выборки, по которым определены оценки дисперсии, получены из генеральных совокупностей, обладающих одинаковыми дисперсиями

μ₂(x)= μ₂(y)= … =μ₂(w)= μ₂,

при этом величина дисперсии μ₂ остается неизвестной. Следует выяснить, являются ли величины μ^*₂(x) , μ^*₂ (y), …, μ^*₂ (w) оценками одной и той же _{генеральной дисперсии μ2.}

_{Рассмотрим сначала случай, когда объем выборок по слоям хотя бы частично различается. В такой ситуации применяется критерий однородности Бартлетта. Проверка однородности реализуется в несколько шагов.}

_{Вычисляется усредненная оценка несмещенной дисперсии по всем слоям}

_(6.4)

_{где μ}^*_{2 (i) – несмещенная оценка дисперсии для слоя i.}

_{Рассчитывается значение критерия}

(6.5)

Бартлетт установил, что случайная величина В при условии справедливости нулевой гипотезы распределена приближенно как хи-квадрат с т–1 степенями свободы, если все n_i больше трех. По заданному уровню значимости a, числу степеней свободы т–1 для правосторонней критической области определяется критическое значение с²_кр(m-1;a). Если соблюдается условие

B< с²_кр(m-1;a),

то нет оснований отвергнуть нулевую гипотезу. Если B > с²_кр(m-1;a), то нулевая гипотеза отвергается. Критерий Бартлетта чувствителен к отклонениям распределения от нормального, поэтому к результатам сравнения следует относится осторожно, а при одинаковом объеме всех слоев вместо критерия Бартлетта лучше применять критерий Кочрена (Кохрена).

Итак, если k₁ = k₂ = … = k_т, то применяется критерий Кочрена

(6.6)

где μ_{2, max} – максимальная оценка дисперсии по всем слоям.

Критическая область для критерия Кочрена правосторонняя. Критическую точку   G_кр(k₁,m;) находят по таблице распределения Кочрена, фрагмент которой приведен в табл. Критическая область определяется неравенством G > G_кр(k₁,m;) .

 

6.3.3. Сравнение факторной и остаточной дисперсий

 

_{Пусть все выборки (6.1) характеризуют одну случайную величину Х при различных значениях фактора Φ, т.е. каждый слой соответствует одному количественному или качественному значению фактора. Сравнение дисперсий производится в следующем порядке:}

• _{рассчитывается среднее значение (оценка математического ожидания) по всей совокупности наблюдений}

,

где n=n₁+n₂+…+n_m , а х_ij – j-й элемент i-го слоя;

• вычисляются средние значения для всех слоев (групп) ;

• определяется общая сумма квадратов отклонений наблюдаемых значений от оценки математического ожидания

; (6.7)

• определяется факторная сумма квадратов отклонений средних по слоям от оценки математического ожидания (характеризует рассеяние между слоями)

; (6.8)

• определяется остаточная сумма квадратов отклонений наблюдаемых значений внутри слоя от своей средней

. (6.9)

Величина S_факт характеризует влияние фактора Φ. Это положение можно пояснить следующим образом. Пусть фактор оказывает существенное влияние на величину Х. Тогда результаты наблюдения для одного слоя, вообще говоря, отличаются от результатов, представленных в других слоях. Следовательно, различаются и средние значения по слоям, причем они тем больше отличаются от оценки математического ожидания по всей выборке, чем больше проявляется влияние фактора. Таким образом, сумма квадратов отклонений средних по слоям от общей средней и характеризует влияние фактора (возведение отклонений во вторую степень исключает взаимную компенсацию положительных и отрицательных отклонений).

Наблюдения внутри одного слоя различаются из-за воздействия случайных причин. Именно сумма квадратов отклонений наблюдаемых значений в каждом слое от среднего значения в слое и характеризует воздействие этих причин, т.е. величина S_ост отражает суммарное влияние случайных причин на значение величины Х.

Величина S_общ, как сумма квадратов отклонений конкретных значений от среднего значения, характеризует суммарное влияние фактора и случайных причин. Можно показать, что

S_общ = S_ост + S_факт,

_{тогда для вычисления остаточной суммы квадратов можно воспользоваться более простым соотношением}

_{Sост = Sобщ – Sфакт.}

_{Разделив суммы квадратов отклонений на соответствующее число степеней свободы, получим оценки общей, факторной и остаточной дисперсий:}

(6.10)

_{Если средние значения случайной величины, вычисленные по отдельным выборкам одинаковы, то оценки факторной и остаточной дисперсий являются несмещенными оценками генеральной дисперсии и различаются несущественно. Тогда сопоставление оценок этих дисперсий по критерию Р. Фишера}

F = μ_2факт /μ_2ост

должно показать, что нулевую гипотезу о равенстве факторной и остаточной дисперсий отвергнуть нет оснований. Если μ_2факт < μ_2ост,  то нет необходимости прибегать к вычислению критерия Р. Фишера – из неравенства сразу следует вывод о выполнении нулевой гипотезы. Итак, из справедливости гипотезы о равенстве средних величин по группам следует соблюдение гипотезы о равенстве факторной и остаточной дисперсий.

Если нулевая гипотеза о равенстве средних величин по слоям является ложной, то с увеличением расхождения между слоями возрастает оценка факторной дисперсии, а вместе с ней и величина критерия F = μ_2факт /μ_2ост. В результате значение F превысит критическое значение, и гипотеза о равенстве дисперсий будет отвергнута.

Рассуждая от противного, можно доказать справедливость утверждений: из справедливости (ложности) гипотезы о дисперсиях следует истинность (ложность) гипотезы о математических ожиданиях. Таким образом, вместо проверки нулевой гипотезы Н₀ о равенстве средних значений для совокупности выборок следует проверить гипотезу о равенстве факторной и остаточной дисперсий.

Пример 6.2. Методом дисперсионного анализа при уровне значимости 0,05 проверить нулевую гипотезу о равенстве средних значений по слоям, применительно к результатам наблюдений, табл. 6.1. Предполагается, что выборки принадлежат нормальному распределению, а каждый слой соответствует некоторому значению фактора Φ.

Решение. Необходимо проверить однородность дисперсий, а затем непосредственно провести дисперсионный анализ. Проверим гипотезу об однородности дисперсий. Для этого вычислим:

• оценки математического ожидания по слоям (групповые средние)

m_{гр 1} = 263,93; m_{гр 2} = 262,95; m_{гр 3} = 265,32;

• несмещенные оценки дисперсии по слоям

m^*₂ (1) = 29,79; m^*₂ (2) = 54,20; m^*₂ (3) = 34,61;

• усредненную оценку несмещенной дисперсии по всем слоям

m^*₂ =(29,79×5+54,20×6+34,61×6)/17 = 40,11.

• значение критерия Бартлетта

В=a/c=0,56/1,08=0,52,

где  а=2,303(17lg 40,11 – (5 lg 29,79+6 lg 54,20+6 lg 34,61)) = 0,56;

 с = 1+(1/5+1/6+1/6 –1/17)/[3(3–1)]=1,08.

Критическое значение хи-квадрат для правосторонней области

c²_кр(2; 0,05)=6,0.

Поскольку величина В меньше c²_кр(2; 0,05), отвергнуть нулевую гипотезу об однородности дисперсий нет оснований.

_{Дисперсионный анализ предусматривает вычисление:}

• _{суммы квадратов}

S_общ = 701,65; S_факт = 19,81; S_ост = 681,84;

• оценок дисперсий

μ^*_{2 общ}=701,65/19=36,93; μ^*_{2 факт}=19,81/2=9,91; μ^*_{2 ост}=681,84/17=40,10.

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