[R] contrasts for lmer model

Kay Cichini Kay.Cichini at uibk.ac.at
Wed May 19 14:29:36 CEST 2010


.. i found this helpful thread:

http://r.789695.n4.nabble.com/Please-Post-Planned-Contrasts-Example-in-lme-nlme-td1598717.html#a1599031

m0<-lm(rich ~ gap*stage,data=richness)

cc<-contrast(m0,
a = list(stage="A",gap="0"),
b = list(stage="A",gap="1"))

m2<-lmer(rich~gap*stage+(1|site),family=poisson,data=richness)

# plug in matrix from contrast(m0)
summary(glht(m2, linfct = cc$X)) 

ok, this works fine - for each stage seperatly..

but what if i needed all 4 gap effects within each of the stages?
i think that then it may be needed to test "gap" within each stage
simultanously so that the p-values will be adjusted accordingly, isn't it?

is this possible? - if so, how would this be done?
anyone who can help with this..

best,
kay



-----
------------------------
Kay Cichini
Postgraduate student
Institute of Botany
Univ. of Innsbruck
------------------------

-- 
View this message in context: http://r.789695.n4.nabble.com/contrasts-for-lmer-model-tp2222682p2222743.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list