Skip to content
Prev 4835 / 7420 Next

How do I interpret linear mixed model contrast estimates from multcomp::glht()?

So, what do the rows correspond to in the summary (e.g. "1 == 0")? I was
thinking the answer was buried *cc*, but I can't figure it out. Consider
this modified example I stole from here
<https://stat.ethz.ch/pipermail/r-sig-mixed-models/2009q4/003061.html>...
(1|schoolid), egsingle)
= levels(egsingle$female)), +
                                                b = list(year = c(.5, 1.5,
2.5), size = 800, female = levels(egsingle$female)))
Simultaneous Tests for General Linear Hypotheses

Fit: lme4::lmer(formula = math ~ year * size + female + (1 | childid) +
    (1 | schoolid), data = egsingle)

Linear Hypotheses:
              Estimate   Std. Error   z value   Pr(>|z|)
1 == 0  0.12774    0.08020     1.593     0.1272
2 == 0  0.15322    0.08066     1.900    0.0669 .
3 == 0  0.17870    0.08178     2.185    0.0341 *
4 == 0  0.12774    0.08020     1.593    0.1273
5 == 0  0.15322    0.08066     1.900    0.0669 .
6 == 0  0.17870    0.08178     2.185    0.0342 *
---
Signif. codes:  0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1
(Adjusted p values reported -- single-step method)

Ultimately I would like to create a dataframe so I can plot the contrasts,
something like this...
x[["test"]]$coefficients, Std.Error = x[["test"]]$sigma)

Thanks for the help!
Cheers,
MVS
=====
Matthew Van Scoyoc

<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=mvanscoyoc at aggiemail.usu.edu>
https://sites.google.com/site/scoyoc/
=====
Think SNOW!