Допустим, у нас есть два фактора (A и B), каждый из которых имеет два уровня (A1, A2 и B1, B2) и переменную отклика (y).
При выполнении двухстороннего ANOVA типа:
y~A+B+A*B
Мы проверяем три нулевые гипотезы:
- Нет разницы в средствах фактора А
- Там нет разницы в средствах фактора B
- Нет взаимодействия между факторами А и В
Когда записано, первые две гипотезы легко сформулировать (для 1 это )
Но как сформулировать гипотезу 3?
редактировать : и как это будет сформулировано для случая более двух уровней?
Благодарю.
hypothesis-testing
anova
Таль Галили
источник
источник
H_0 = \mu_{A1}=\mu_{A2}
\mu_{A_1}
Ответы:
Я думаю, что важно четко разделить гипотезу и соответствующий ей тест. Для следующего я предполагаю сбалансированную схему 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 - это наблюдение при обработке фактора и при обработке фактора с , и . Модельi j A k B 1≤i≤n 1≤j≤p 1≤k≤q Yя J K= μJ K+ ϵя ( j k ),εя ( j k )∼ N( 0 , σ2ε)
Дизайн: А 1...A j...А п Б 1μ11...μJ 1...μр 1μ0,1.....................B kμ1 к...μJ K...μр кμ, К.....................B qμ1 кв...μJ Q...μр дμ, Q μ1....μДж ....μстр .μ
j k ϵ i ( j k ) i ( ) j k iμJ K - ожидаемое значение в ячейке , - ошибка, связанная с измерением человека в этой ячейке. Запись указывает, что индексы являются фиксированными для любого данного человека потому что этот человек наблюдается только в одном состоянии. Несколько определений для эффектов:J K ϵi(jk) i () jk i
jAμj.=1q∑qk=1μjk (среднее ожидаемое значение для обработки фактора )j A
kBμ.k=1p∑pj=1μjk (среднее ожидаемое значение для обработки фактора )k B
j A ∑ p j = 1 α j = 0αJ= μДж .- μ (эффект обработки фактора , )J A ΣпJ = 1αJ= 0
k B ∑ q k = 1 β k = 0βК= μ, К- μ (эффект обработки фактора , )К В ΣQк = 1βК= 0
j A k B ∑ p j = 1 ( α β ) j k = 0( α β)J K= μJ K- ( μ + αJ+ βК) = μJ K- μДж .- μ, К+ μ J A К В ΣпJ = 1( α β)J K= 0∧ΣQк = 1( α β)J K= 0 )
(эффект взаимодействия для комбинации обработки фактора с обработкой фактора ,
j A k B ∑ p j = 1 α ( k ) j = 0α( к )J= μJ K- μ, К J A К В ΣпJ = 1α( к )J= 0∧1QΣQк = 1α( к )J= αJ∀J , K )
(условный основной эффект для обработки фактора рамках фиксированной обработки фактора ,
k B j A ∑ q k = 1 β ( j ) k = 0β( J )К= μJ K- μДж . К В J A ΣQк = 1β( J )К= 0∧1пΣпJ = 1β( J )К= βК∀J , K )
(условный основной эффект для обработки фактора рамках фиксированной обработки фактора ,
С этими определениями модель также может быть записана как:Yя J K= μ + αJ+ βК+ ( α β)J K+ ϵя ( j k )
Это позволяет нам выразить нулевую гипотезу отсутствия взаимодействия несколькими эквивалентными способами:
0 μ j k = μ + α j + β kЧАС0я: ∑JΣК( α β)2J K= 0 0 μJ K= μ + αJ+ βК∀ j , k
(все отдельные члены взаимодействия равны , так что . Это означает, что эффекты лечения обоих факторов - как определено выше - везде аддитивны.)
(all conditional main effects for any treatment
(all conditional main effects for any treatment
источник
Взаимодействие говорит нам, что уровни фактора А оказывают различное влияние в зависимости от того, какой уровень фактора Б вы применяете. Таким образом, мы можем проверить это с помощью линейного контраста. Пусть C = (A1B1 - A1B2) - (A2B1 - A2B2), где A1B1 обозначает среднее значение для группы, получившей A1 и B1, и так далее. Итак, здесь мы смотрим на A1B1 - A1B2, который является эффектом, который оказывает фактор B, когда мы применяем A1. Если нет взаимодействия, это должно быть то же самое, что эффект, который имеет B, когда мы применяем A2: A2B1 - A2B2. Если они одинаковы, то их разница должна быть 0, чтобы мы могли использовать тесты:
источник