[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