[ADMB Users] Nested random effects syntax in glmmADMB 0.6.2 alpha
frocga at gmail.com
Thu Aug 4 21:52:57 PDT 2011
Today I managed to get the alpha version of glmmADMB 0.6.2 loaded and
running. I am running negative binomial model with nested random effects.
I decided to start with a reduced model of my data, so I began with the
following syntax and R would just hang and I get the pinwheel of death (Mac
I can't find the exact mailing list message, but somewhere I believe I read
that +(1|A)+(1|B) is the same as the above syntax in lme4 as long as factor
B is coded with unique individual labels.
I then tried the following model syntax and it worked just fine:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.841 0.234 -3.59 0.00034 ***
treaty -0.702 0.311 -2.26 0.02413 *
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Number of observations: total=160, =32, =160
Random effect variance(s):
Negative binomial dispersion parameter: 403.43 (std. err.: 2.1584)
My data structure:
A =(randomly selected plot) B =(randomly selected sub-plot within A)
X =(disturbed=Y or N) Y =(total count of individuals present)
Is glmmADMB treating B as being nested in A like it does in lme4? To the
best of your knowledge do the results look like they should?
I'll have more questions later but gotta get this sorted first.
Thanks for the help,
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Users