[R-sig-ME] Calculating SE from GLMM results (using glmer{lme4})
Raldo Kruger
raldo.kruger at gmail.com
Thu Sep 24 08:05:32 CEST 2009
Dear R users,
Please excuse the basic questions, but I’m new to GLMMs and R!
I’m analyzing an experiment where Seedling numbers in plots where seed
has been sown on restoration sites is the response variable. I’m most
interested in determining whether the Nutrients (N) and water
absorbing polymer Gel (Ge) additions to the soil substrate contribute
positively to the survival of the seedlings, over a 3 year time period
(for simplicity I'm just using 3 time periods, each in the same season
for the 3
successive years).
Fixed factors: Nutrients (0 and 1), Gel (0 and 1)
Random factors: Site (4 non replicate sites), Year (3 time periods)
Response variable: Seedling numbers (counts) / 0.25m2 plot
The results are as follows:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 4.52982 0.24486 18.5 <2.00E-16 ***
N -0.07922 0.08415 -0.94 0.346489
Ge 0.20766 0.08428 2.46 0.013744 *
Year -1.62937 0.04672 -34.88 <2.00E-16 ***
N:Ge -0.44213 0.11898 -3.72 0.000202 ***
N:Year 0.11705 0.06322 1.85 0.064125 .
Ge:Year -0.04861 0.0645 -0.75 0.451132
N:Ge:Year 0.11458 0.08917 1.28 0.198821
1) So as I understand (from previous correspondence with R-users) the
number of seedlings in the control plots in year 0 is
exp(4.53) = 92.7. Is the standard error calculated with
0.24486 (i.e. 92.7*0.24), or with 92.7*exp(0.24).
2) And for the N:Ge treatment, the effect is exp(-0.08+0.21-0.44)
=0.73 (I.e. a 27% reduction compared to the control), right? So
is the SE for the N:Ge effect calculated as the sum of the
SE’s too, i.e. 0.08+0.08+0.12, or is it just 0.12?
3) Lastly, is it possible to fit two response variables in one GLMM?
E.g. seedling numbers and height.
Many thanks,
Raldo Krüger
Msc student
University of Cape Town
--
Raldo
More information about the R-sig-mixed-models
mailing list