J'ai des problèmes à trouver une solution concernant la façon d'exécuter un test post-hoc (Tukey HSD) après une ANOVA à mesures répétées à 2 facteurs (tous deux intra-sujets) en R. Pour l'ANOVA, j'ai utilisé la fonction aov:
summary(aov(dv ~ x1 * x2 + Error(subject/(x1*x2)), data=df1))
Après avoir lu les réponses à d'autres questions, j'ai compris que je devrais d'abord relancer l'ANOVA en utilisant une autre fonction (par exemple, lme). C'est ce que j'ai trouvé.
Lme.mod <- lme(dv ~ x1*x2, random=list(subject=pdBlocked(list(~1, pdIdent(~x1-1), pdIdent(~x2-1)))), data=df1)
anova(Lme.mod)
Les deux effets principaux étaient significatifs, mais il n'y avait aucun effet d'interaction. Ensuite, j'ai utilisé ces fonctions pour les comparaisons post-hoc:
summary(glht(Lme.mod, linfct=mcp(x1="Tukey")))
summary(glht(Lme.mod, linfct=mcp(x2="Tukey")))
Cependant, il y avait quelques problèmes:
Tout d'abord, le fichier d'aide R indique que «La fonction mcp doit être utilisée avec précaution lors de la définition des paramètres d'intérêt dans les modèles ANOVA ou ANCOVA bidirectionnels (...) multcomp version 1.0-0 et supérieure génère des comparaisons pour les effets principaux uniquement, en ignorant les covariables et les interactions (les anciennes versions font automatiquement la moyenne des termes d'interaction). Un avertissement est donné. " Et bien sûr, j'ai reçu le message d'avertissement suivant:
Warning message:
In mcp2matrix(model, linfct = linfct) :
covariate interactions found -- default contrast might be inappropriate
Une autre chose déconcertante était que, bien que les deux effets principaux soient significatifs, il n'y avait pas de différences significatives dans les comparaisons post-hoc pour l'un des facteurs (x1). Je n'ai jamais rencontré ça auparavant. Les scripts / analyses sont-ils corrects / appropriés ou y a-t-il quelque chose qui me manque? Toute aide serait appréciée!
la source
Réponses:
Aurait
être ce que vous recherchez, c'est-à-dire faire des tests posthoc parmi toutes les combinaisons de niveaux de mesure des deux facteurs x1 et x2? (J'ai également imposé une symétrie composée, pour que le résultat lme corresponde à celui des mesures répétées appelées aov)
la source
Test de multicomparaison de Tukey
Installer le package multcomp install.packages ("multcomp")
Rendre multcomp disponible pour une bibliothèque d'utilisation ("multcomp")
Vérifier qu'il fonctionne - Explique quels packages sont actuellement ouverts dans R search ()
Utilisez ensuite la fonction glht ()
la source