Les posteurs les plus actifs de la semaine
Eric Wajnberg
 
zezima
 
Nik
 
margotte185
 
c@ssoulet
 
Coco
 
Ayana
 


J'ai l'impression que l'option "Satterthwaite" ne marche pas

Voir le sujet précédent Voir le sujet suivant Aller en bas

J'ai l'impression que l'option "Satterthwaite" ne marche pas

Message par zezima le Lun 19 Juin 2017 - 13:13

Bonjour,

J'essaie de faire une correction de Satterthwaite dans une ANOVA afin de prendre en compte l'hétérogénéité des variances et d'évaluer les différences de résultats avec une ANOVA classique (homoscédastique).

Voici les 2 tests que j'ai réalisé :

Code:
/*data*/
data data; input trt $ y;
datalines;
A   21.125
A   20.575
A   40.525
A   26.15
A   26.35
A   44.275
A   49.725
A   63.6
A   49.35
B   68.875
B   58.025
B   58.075
B   68.025
B   52.3
B   15.2
B   14.6
B   18.7
B   20.2
;

/*ANOVA*/
proc mixed data=data;
class trt;
model y = trt;
ODS OUTPUT Tests3=Tests3;
run;

/*Satterthwaite*/
proc mixed data=data;
class trt;
model y = trt / ddfm=SATTERTHWAITE;
ODS OUTPUT Tests3=Tests3;
run;

Quelqu'un saurait-il comment appliquer cette option, les résultats sont systématiquement similaires au niveau de la p-valeur.

Merci d'avance
avatar
zezima

Nombre de messages : 747
Date d'inscription : 26/02/2013

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: J'ai l'impression que l'option "Satterthwaite" ne marche pas

Message par Eric Wajnberg le Lun 19 Juin 2017 - 15:46

zezima a écrit:Quelqu'un saurait-il comment appliquer cette option, les résultats sont systématiquement similaires au niveau de la p-valeur.
Formellement, l'approximation de Satterthwaite ne modifie pas (forcément) la p-value, mais les df. Je commencerais pas vérifier ceci. Par ailleurs, il y modification que si les variances sont différentes. est-ce le cas sur ce jeux de données (je n'ai pas vérifié) ?

HTH, Eric.
avatar
Eric Wajnberg

Nombre de messages : 694
Date d'inscription : 14/09/2012

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: J'ai l'impression que l'option "Satterthwaite" ne marche pas

Message par Ayana le Mar 20 Juin 2017 - 6:07

Bonjour,
Eric a raison, la correction de Satterthwaite ne modifie que le nombre de degres de liberte.
Cependant, dans le cas de modeles mixtes, si la variance de l'effet aleatoire est nulle ou extremement proche de 0, aucune correction n'est appliquee dans R. Donc je suppose qu'un principe similaire est applique ici si les variances sont egales, afin de ne pas perdre en puissance
Ayana
avatar
Ayana

Nombre de messages : 458
Localisation : Londres
Date d'inscription : 18/08/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: J'ai l'impression que l'option "Satterthwaite" ne marche pas

Message par Ayana le Mar 20 Juin 2017 - 6:18

Essaie l'option method=type 3 en option de ta proc mixed afin d'avoir la table d'ANOVA pour comparer les valeurs. Et surtout, dans SAS, es-tu sur que ddfm fasse quelque chose lorsque la commande random est vide? Car initialement, proc mixed a ete developpee pour des modeles a effets aleatoires, et l'option ddfm pour corriger les DoF pour eviter une inflation de l'erreur de type I observee qand le nombre de "clusters" est faible.
Ayana


Dernière édition par Ayana le Mar 20 Juin 2017 - 6:25, édité 1 fois (Raison : q)
avatar
Ayana

Nombre de messages : 458
Localisation : Londres
Date d'inscription : 18/08/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: J'ai l'impression que l'option "Satterthwaite" ne marche pas

Message par zezima le Mar 20 Juin 2017 - 9:30

Merci pour vos réponses, qui expliquent pourquoi j'obtenais des résultats similaires dans les deux cas de figures.

Cela m'a mis la puce à l'oreille et en fouillant plus profondément, le Satterthwaite ne change en effet que le ddlc'est la correction de Welch qui prend en compte l'hétérogénéité des variances.

Une façon d'obtenir des résultats différents avec le Satterthwaite est :
Code:

proc mixed data=data;
class trt;
model y = trt / ddfm=SATTERTHWAITE;
ODS OUTPUT Tests3=Tests3;
repeated / GRP=trt;
lsmeans trt;
run;

Après, au niveau de quel correction prendre entre le Welch et le Satterthwaite en cas d'hétéscédasticité, c'est une autre question.
avatar
zezima

Nombre de messages : 747
Date d'inscription : 26/02/2013

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: J'ai l'impression que l'option "Satterthwaite" ne marche pas

Message par Contenu sponsorisé


Contenu sponsorisé


Revenir en haut Aller en bas

Voir le sujet précédent Voir le sujet suivant Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum