[R] contrasts for lmer model
Kay Cichini
Kay.Cichini at uibk.ac.at
Wed May 19 15:23:33 CEST 2010
plugged it together myself -
obviously this should work:
m0<-lm(rich ~ gap*stage,data=richness)
cc<-contrast(m0,
list(stage = levels(stage), gap = "1"),
list(stage = levels(stage), gap = "0"))
# plug in matrix from contrast(m0)
summary(glht(m3, linfct = cc$X))
Simultaneous Tests for General Linear Hypotheses
Fit: glmer(formula = rich ~ gap * stage + (1 | site), family = poisson,
REML = 1)
Linear Hypotheses:
Estimate Std. Error z value Pr(>|z|)
1 == 0 0.04820 0.15594 0.309 0.9965
2 == 0 -0.44679 0.10204 -4.379 4.78e-05 ***
3 == 0 -0.52847 0.09133 -5.786 2.88e-08 ***
4 == 0 -0.44396 0.15549 -2.855 0.0171 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Adjusted p values reported -- single-step method)
-----
------------------------
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-tp2222682p2222820.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list