Quelle est l'hypothèse NULL pour l'interaction dans une ANOVA bidirectionnelle?


20

Disons que nous avons deux facteurs (A et B), chacun avec deux niveaux (A1, A2 et B1, B2) et une variable de réponse (y).

Lorsque vous effectuez une ANOVA bidirectionnelle du type:

y~A+B+A*B

Nous testons trois hypothèses nulles:

  1. Il n'y a pas de différence dans les moyennes du facteur A
  2. Il n'y a pas de différence dans les moyennes du facteur B
  3. Il n'y a pas d'interaction entre les facteurs A et B

Lorsqu'elles sont écrites, les deux premières hypothèses sont faciles à formuler (pour 1 c'est )H0:μA1=μA2

Mais comment formuler l'hypothèse 3?

edit : et comment serait-il formulé pour le cas de plus de deux niveaux?

Merci.


3
Je n'ai pas la réputation de me permettre d'éditer, mais je pense que tu veux (ou μ A 1 si vous voulez un double indice) [oups, il a automatiquement tex-ifié que : ou]H0=μA1=μA2μA1H_0 = \mu_{A1}=\mu_{A2}\mu_{A_1}
Ben Bolker

1
Oups, vous n'avez pas vu que vous utilisez des majuscules pour désigner le nom du facteur et ses niveaux - corrigez-le (en suivant la notation @Ben).
chl

Réponses:


18

Je pense qu'il est important de séparer clairement l'hypothèse et son test correspondant. Pour ce qui suit, je suppose un plan CRF- équilibré entre les sujets (tailles de cellules égales, notation de Kirk: plan factoriel complètement aléatoire).pq

est l'observation i dans le traitement j du facteur A et le traitement k du facteur B avec 1 i n , 1 j pYijkijAkB1in1jp et . Le modèle est Y i j k = μ j k + ϵ i ( j k ) ,1kqYijk=μjk+ϵi(jk),ϵi(jk)N(0,σϵ2)

Conception: B 1 ... B k ... B q A 1 μ 11 ... μ 1 k ... μ 1 q μ 1. ... ... ... ... ... ... ... A j μ j 1 ... μ j k ... μ j q μ j . ... ... ... ... ... ... ... A p μ p 1 ... μ B1BkBq A1μ11μ1kμ1qμ1.Ajμj1μjkμjqμj.Apμp1μpkμpqμp. μ.1μ.kμ.qμ

est la valeur attendue dans la cellule j k , ϵ i ( j k ) est l'erreur associée à la mesure de la personne i dans cette cellule. Lanotation ( ) indique que les indices j k sont fixes pour toute personne i donnéecar cette personne est observée dans une seule condition. Quelques définitions des effets:μjkjkϵi(jk)i()jki

(valeur moyenne attendue pour le traitementjdu facteurA)μj.=1qk=1qμjkjA

(valeur moyenne attendue pour le traitementkdu facteurB)μ.k=1pj=1pμjkkB

(effet du traitement j du facteur A ,p j = 1 α j = 0 )αj=μj.μjAj=1pαj=0

(effet du traitement k du facteur B ,q k = 1 β k = 0 )βk=μ.kμkBk=1qβk=0

(effet d'interaction pour la combinaison du traitement j du facteur A avec le traitement k du facteur B ,p j = 1 ( α β ) j k =(αβ)jk=μjk(μ+αj+βk)=μjkμj.μ.k+μ
jAkBj=1p(αβ)jk=0k=1q(αβ)jk=0)

αj(k)=μjkμ.k
(conditional main effect for treatment j of factor A within fixed treatment k of factor B, j=1pαj(k)=01qk=1qαj(k)=αjj,k)

βk(j)=μjkμj.
(conditional main effect for treatment k of factor B within fixed treatment j of factor A, k=1qβk(j)=01pj=1pβk(j)=βkj,k)

With these definitions, the model can also be written as: Yijk=μ+αj+βk+(αβ)jk+ϵi(jk)

This allows us to express the null hypothesis of no interaction in several equivalent ways:

  1. H0I:jk(αβ)jk2=0
    (all individual interaction terms are 0, such that μjk=μ+αj+βkj,k. This means that treatment effects of both factors - as defined above - are additive everywhere.)

  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: In a diagramm which shows the expected values μjk with the levels of factor A on the x-axis and the levels of factor B drawn as separate lines, the q different lines are parallel.


1
A really impressive answer Caracal - thank you.
Tal Galili

9

An interaction tells us that the levels of factor A have different effects based on what level of factor B you're applying. So we can test this through a linear contrast. Let C = (A1B1 - A1B2) - (A2B1 - A2B2) where A1B1 stands for the mean of the group that received A1 and B1 and so on. So here we're looking at A1B1 - A1B2 which is the effect that factor B is having when we're applying A1. If there is no interaction this should be the same as the effect B is having when we apply A2: A2B1 - A2B2. If those are the same then their difference should be 0 so we could use the tests:

H0:C=0vs.HA:C0.


1
Thanks Dason, that helped. Also, after reading your reply, it suddenly became clear to me that I am not fully sure how this generalizes in case we are having more factors. Could you advise? Thanks again. Tal
Tal Galili

2
You can test multiple contrasts simultaneously. So for example if A had three levels and B had 2 we could use the two contrasts: C1 = (A1B1 - A2B1) - (A2B1 - A2B2) and C2 = (A2B1 - A2B2) - (A3B1 - A3B2) and use a 2 degree of freedom test to simultaneously test if C1 = C2 = 0. It's also interesting to note that C2 could equally have been (A1B1 - A1B2) - (A3B1 - A3B2) and we would come up with the same thing.
Dason

Hi @Dason: you seem to have multiple accounts. Could you please complete the form at stats.stackexchange.com/contact and request that they be merged? That will simplify your use of this site (and give you the combined net reputation of both accounts).
whuber
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.