Je travaille sur un modèle à effets imbriqués (dont je ne connais pas grand chose) :
J'ai 4 labos et 2 opérateurs par labo. Parmis ces 4 labos il y en a un qui est " la référence".
Avec le modèle qui suis, j'obtiens cette table d'ANOVA :
Code : Tout sélectionner
lm(CFU.AD~Lab+Operator%in%Lab,data=don_AD)
anova(model_AD)
Analysis of Variance Table
Response: CFU.AD
Df Sum Sq Mean Sq F value Pr(>F)
Lab 3 1.3639 0.45465 14.3729 3.573e-08 ***
Lab:Operator 4 1.2562 0.31405 9.9283 4.569e-07 ***
Residuals 134 4.2387 0.03163
j'ai pu faire les comparaisons multiples associées au labos avec :
Code : Tout sélectionner
summary(glht(model_AC, linfct=mcp(Lab="Dunnett")))
Mais comment obtenir ces comparaison pour l'effet imbriqué ?
Cordialement
Jérémy