[R-sig-ME] lmer/glmer standard error interpretation and visualization

Colin Wahl biowahl at gmail.com
Wed Mar 7 21:20:32 CET 2012


Hello,
I am in the process of finalizing figures for my thesis on stream
invertebrate distributions among watershed and riparian types. See
below for additional information on the design. I'm having difficulty
including standard errors from the lmer modeling as error bars in the
figures. Here is the table I've created from the lmer output:
estimates of %EPT and St Error are back transformed from logits and
converted from fractions to percents. Estimates are also absolute (not
relative to the intercept).

 Watershed     Effect                  Estimate  St. Error  z score  p value

 Forested      Intercept: F vs. 0       28.23      59.6     -2.346     0.019*
                    Riparian:  F vs. NF   16.017     62.3     -1.436     0.151
 Cultivated    Watershed: C vs. F    1.351     65.3     -5.297    <0.000*
                    Riparian:  F vs. NF    1.555      69.2      1.071    0.284
 Developed   Watershed: D vs. F    0.175      66.8     -7.714    <0.000*
                    Riparian:  F vs. NF    0.292      70.9      1.391     0.164
 Grassland    Watershed: G vs. F    28.94     66.6       0.05      0.960
                    Riparian:  F vs. NF     1.967     70.7     -2.595     0.009*

The st. errors are huge. I initially used standard error calculations
in excel for error bars (stdev(x)/sqrt(n(x))), which look very
reasonable, and are reflective of significant differences.

Does anyone have any advice to offer for visualizing these glmer
results? Should I use the huge model St. Errors? My inclination is
yes, because they are used to calculated significant differences, but
28 + or - 59.6 with a significant p value seems ridiculous.

Thank you,
Colin Wahl
M.S. Candidate
Dept. of Biology
Western Washington University




More information about the R-sig-mixed-models mailing list