Modele effets imbriqués

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Jérémy Jachacz
Messages : 153
Enregistré le : 15 Avr 2014, 12:56

Modele effets imbriqués

Messagepar Jérémy Jachacz » 30 Jan 2017, 11:07

Bonjour,

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
Statisticien (69)

Maxime Hervé
Messages : 427
Enregistré le : 03 Mar 2010, 14:21
Contact :

Re: Modele effets imbriqués

Messagepar Maxime Hervé » 30 Jan 2017, 15:00

Bonjour,
en utilisant le package lsmeans ça donnerait :

Code : Tout sélectionner

lsm <- lsmeans(model_AC,~Operator|Lab)
contrast(lsm,"pairwise")

La syntaxe Operator|Lab précisant que vous souhaitez faire des comparaisons multiples entre niveaux d'Operator au sein de chaque niveau de Lab.
Maxime

Jérémy Jachacz
Messages : 153
Enregistré le : 15 Avr 2014, 12:56

Re: Modele effets imbriqués [Résolu]

Messagepar Jérémy Jachacz » 30 Jan 2017, 16:06

Bonjour Maxime,

merci, ça fonctionne parfaitement

Jérémy
Statisticien (69)


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité