Bonjour, et merci de l'aide que vous apportez aux utilisateurs !
Je suis en train d'écrire un script pour analyser des données provenant d'un appareillage. Je pars d'une DF qui comporte une colonne ('GROUP') qui encode le groupe expérimental (3 groupes au total), une colonne ('Time') qui encode le moment de l'analyse (12 moments différents) et une colonne ('OCR') qui encode la donnée mesurée par l'appareil dans ces 3 groupes et à ces 12 moments, chaque groupe contient 3 échantillons indépendants, soit 108 observations au total. Je représente ces données sous la forme de séries temporelles [3 courbes sur un graphique OCR = f(temps)]... c'est très beau !
Je cherche à comparer les moyennes de 2 (hors des 3) groupes à chaque temps... et bien bizarrement je n'y arrive pas ! :)
Voilà comment j'ai procédé :
TWANOVA <- aov(OCR ~ Group * Time, data = RATE_data)
TukeyHSD(TWANOVA, "Group:Time")
et j'ai donc toutes les comparaisons par paires, comment puis-je limiter aux 12 paires qui m'intéressent :
CONTROL:1-EXPERIMENTAL:1
CONTROL:2-EXPERIMENTAL:2
.../...
CONTROL:12-EXPERIMENTAL:12
Merci pour votre aide !