Discussion :

[niaboc]

Bonjour,

je ne comprends les résultats de la proc ANOVA :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
data exp;
   input A $ B $ Y @@;
   datalines;
A1 B1 12 A1 B1 14     A1 B2 11 A1 B2 9
A2 B1 20 A2 B1 18     A2 B2 17
;
proc anova data=exp;
   class A B;
   model Y=A B;
run;
quit;

Même si on a un plan déséquilibré??

Pourquoi le test global "F" ne donne pas le même résultat qu'avec une proc GLM?

Et pourquoi on a infinity ici??

Merci

[BayHaym]

Bonjour

Selon la doc SAS :

Use PROC ANOVA for the analysis of balanced data only

Source : http://support.sas.com/documentation...va_sect001.htm

Ton exemple étant une unbalanced ANOVA avec deux facteurs, on retrouve le warning à l'exécution de ton code. Le RMSE calculé est nul, la stat du test de Fisher faisant intervenir le RMSE au dénominateur, tu te retrouves avec une stat égale à l'infinie (=> division par 0). Il faut donc se tourner vers la PROC GLM.

La différence observée entre PROC ANOVA et PROC GLM doit venir de l'algo de calcul qui diffère entre les deux. Visiblement ANOVA n'utilise que les diagonales des matrices alors que GLM non.
Sources : http://support.sas.com/documentation...va_sect021.htm
http://support.sas.com/documentation...lm_sect040.htm

[niaboc]

Merci pour ta réponse.

Si j'ai bien compris, la proc ANOVA ne fait pas vraiment une estimation par les moindres carrés ; c'est simplement des différences de moyennes qui permettent d'estimer les paramètres... et c'est pourquoi ça ne marche pas pour les plans déséquilibrés, car on voit clairement que des différences de moyenne ne peuvent pas estimer les paramètres dans un tel cas.

Ai-je bien compris?

[BayHaym]

Honnêtement je suis pas assez à l'aise avec les formules pour confirmer toute ta phrase, mais je suis d'accord avec le début.

Dans les moindres carrés classiques, le vecteur de paramètres est calculé via inversion de la matrice X'X, alors que la PROC ANOVA ne travaille que sur la diagonale de X'X. Ce qui revient au même uniquement en cas de plan équilibré (+ autres exceptions listées dans la doc SAS).

[niaboc]

On est d'accord

Merci pour ta réponse et pour les sources