Basically, I haven't gotten around to adding a contrast.lmer methods
yet. It's only my list, but in the short term, you are welcome to
contribute code to the package.
There have also been discussions on this list about the
appropriateness of such contrats with nonlinear models, but I need to
imple
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
.. 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=pois
hello,
i found it most convenient to use package contrast for planned comparisons
on mixed models.
for instance i have a model with 2 fixed factors, one with 4 levels (stage)
and one with 2 levels (gap) and a nested random factor (site) and i tested
gap within level A of factor stage:
library(co
4 matches
Mail list logo
