Какова гипотеза NULL для взаимодействия в двухстороннем ANOVA?


20

Допустим, у нас есть два фактора (A и B), каждый из которых имеет два уровня (A1, A2 и B1, B2) и переменную отклика (y).

При выполнении двухстороннего ANOVA типа:

y~A+B+A*B

Мы проверяем три нулевые гипотезы:

  1. Нет разницы в средствах фактора А
  2. Там нет разницы в средствах фактора B
  3. Нет взаимодействия между факторами А и В

Когда записано, первые две гипотезы легко сформулировать (для 1 это )H0:μA1=μA2

Но как сформулировать гипотезу 3?

редактировать : и как это будет сформулировано для случая более двух уровней?

Благодарю.


3
У меня нет репутации, чтобы позволить мне редактировать, но я думаю, что вы хотите (или если вы хотите двойной индекс) [упс, у него есть автоматически отправил сообщение, что: или ] μ A 1H0=μA1=μA2μA1H_0 = \mu_{A1}=\mu_{A2}\mu_{A_1}
Бен Болкер

1
Упс, не видел, что вы используете заглавные буквы для обозначения имени фактора и их уровней - исправьте это (следуя нотации @Ben).
ХЛ

Ответы:


18

Я думаю, что важно четко разделить гипотезу и соответствующий ей тест. Для следующего я предполагаю сбалансированную схему CRF- между субъектами (равные размеры ячеек, запись Кирка: полностью рандомизированный факториальный дизайн).pq

i j A k B 1 i n 1 j p 1 k q Y i j k = μ j k + ϵ i ( j k ) ,Yijk - это наблюдение при обработке фактора и при обработке фактора с , и . МодельijAkB1in1jp1kqYijk=μjk+ϵi(jk),ϵi(jk)N(0,σϵ2)

Дизайн:  B1BkBq A1μ11μ1kμ1qμ1.Ajμj1μjkμjqμj.Apμp1μpkμpqμp. μ.1μ.kμ.qμ

j k ϵ i ( j k ) i ( ) j k iμjk - ожидаемое значение в ячейке , - ошибка, связанная с измерением человека в этой ячейке. Запись указывает, что индексы являются фиксированными для любого данного человека потому что этот человек наблюдается только в одном состоянии. Несколько определений для эффектов:jkϵi(jk)i()jki

jAμj.=1qk=1qμjk (среднее ожидаемое значение для обработки фактора )jA

kBμ.k=1pj=1pμjk (среднее ожидаемое значение для обработки фактора )kB

j A p j = 1 α j = 0αj=μj.μ (эффект обработки фактора , )jAj=1pαj=0

k B q k = 1 β k = 0βk=μ.kμ (эффект обработки фактора , )kBk=1qβk=0

j A k B p j = 1 ( α β ) j k = 0(αβ)jk=μjk(μ+αj+βk)=μjkμj.μ.k+μ
(эффект взаимодействия для комбинации обработки фактора с обработкой фактора ,jAkBj=1p(αβ)jk=0k=1q(αβ)jk=0)

j A k B p j = 1 α ( k ) j = 0αj(k)=μjkμ.k
(условный основной эффект для обработки фактора рамках фиксированной обработки фактора ,jAkBj=1pαj(k)=01qk=1qαj(k)=αjj,k)

k B j A q k = 1 β ( j ) k = 0βk(j)=μjkμj.
(условный основной эффект для обработки фактора рамках фиксированной обработки фактора ,kBjAk=1qβk(j)=01pj=1pβk(j)=βkj,k)

С этими определениями модель также может быть записана как: Yijk=μ+αj+βk+(αβ)jk+ϵi(jk)

Это позволяет нам выразить нулевую гипотезу отсутствия взаимодействия несколькими эквивалентными способами:

  1. 0 μ j k = μ + α j + β kH0I:jk(αβ)jk2=0
    (все отдельные члены взаимодействия равны , так что . Это означает, что эффекты лечения обоих факторов - как определено выше - везде аддитивны.)0μjk=μ+αj+βkj,k

  2. H0I:αj(k)αj(k)=0jk,k(kk)
    (all conditional main effects for any treatment j of factor A are the same, and therefore equal αj. This is essentially Dason's answer.)

  3. H0I:βk(j)βk(j)=0j,jk(jj)
    (all conditional main effects for any treatment k of factor B are the same, and therefore equal βk.)

  4. H0I : в диаграмме, которая показывает ожидаемые значения с уровнями фактора на оси и уровнями фактора нарисованными в виде отдельных линий, различных линий параллельны.μjkAxBq


1
Действительно впечатляющий ответ Каракала - спасибо.
Тал Галили

9

Взаимодействие говорит нам, что уровни фактора А оказывают различное влияние в зависимости от того, какой уровень фактора Б вы применяете. Таким образом, мы можем проверить это с помощью линейного контраста. Пусть C = (A1B1 - A1B2) - (A2B1 - A2B2), где A1B1 обозначает среднее значение для группы, получившей A1 и B1, и так далее. Итак, здесь мы смотрим на A1B1 - A1B2, который является эффектом, который оказывает фактор B, когда мы применяем A1. Если нет взаимодействия, это должно быть то же самое, что эффект, который имеет B, когда мы применяем A2: A2B1 - A2B2. Если они одинаковы, то их разница должна быть 0, чтобы мы могли использовать тесты:

H0:C=0vs.HA:C0.


1
Спасибо, Дейсон, это помогло. Кроме того, после прочтения вашего ответа мне вдруг стало ясно, что я не совсем уверен, как это обобщается, если у нас есть больше факторов. Не могли бы вы посоветовать? Еще раз спасибо. Тал
Тал Галили

2
Вы можете проверить несколько контрастов одновременно. Так, например, если у A было три уровня, а у B было 2, мы могли бы использовать два контраста: C1 = (A1B1 - A2B1) - (A2B1 - A2B2) и C2 = (A2B1 - A2B2) - (A3B1 - A3B2) и использовать 2 Тест на степень свободы для одновременной проверки, если C1 = C2 = 0. Также интересно отметить, что C2 мог бы быть одинаково (A1B1 - A1B2) - (A3B1 - A3B2), и мы придумали бы то же самое.
Дейсон

Привет @Dason: у тебя, похоже, несколько аккаунтов. Не могли бы вы заполнить форму по адресу stats.stackexchange.com/contact и попросить объединить их? Это упростит вам использование этого сайта (и даст вам общую чистую репутацию обеих учетных записей).
whuber
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.