[R] GLMM with poisson distribution (lme4)
Jan Michael Kaesler
jan_kaesler at gmx.de
Tue Mar 24 23:20:09 CET 2009
Dear R-Users,
I have a question to the GLMM via the lme4 package.
I have 160 nest-boxes which are placed at 8 different localities.
"Count" is the number of animals which were found inside the boxes during the observation time.
The independent variables are factors which are supposed to influence
the occurence of the animal.
1. Is the formula right like this?
Here are my outcomes of running the glmer function:
m1<-glmer(Count~A1+A2+A3+A4+A5+(1|locality), family=poisson(log))
summary(m1)
Generalized linear mixed model fit by the Laplace approximation
Formula: Count ~ A1 + A2 + A3 + A5 + A5 + (1 | locality)
AIC BIC logLik deviance
477.8 499.4 -231.9 463.8
Random effects:
Groups Name Variance Std.Dev.
locality (Intercept) 0.15181 0.38963
Number of obs: 160, groups: locality, 8
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -4.410558 1.854530 -2.378 0.017394 *
A1 -0.003627 0.002199 -1.649 0.099089 .
A2 0.037399 0.018599 2.011 0.044350 *
A3 0.041520 0.018579 2.235 0.025432 *
A4 0.115765 0.034663 3.340 0.000839 ***
A5 0.532708 0.073443 7.253 4.06e-13 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Correlation of Fixed Effects:
(Intr) S2 K2 SP2 SUMSP2
A1 0.005
A2 -0.989 -0.030
A3 -0.984 -0.123 0.992
A4 0.001 -0.076 -0.068 -0.036
A5 -0.037 0.123 -0.058 -0.060 -0.041
2. Does it make sense to manually cut out these values which are not significant and to it compare the AIC's of the different models?
There are still more variables which I didn't took inside to keep the example little bit shorter.
Because of an advise I wanted to build up the model by stepwise regression.
I actually wanted to use the functions drop1/add1 which are not possible for this object.
Thanks a lot for every help!
Kind regards,
Jan Kaesler
--
Aufgepasst: Sind Ihre Daten beim Online-Banking auch optimal geschützt?
Jetzt absichern: https://homebanking.gmx.net/?mc=mail@footer.hb
More information about the R-help
mailing list