[ADMB Users] Anova table for glmm.admb objects

[Paula Rosewarne]

Dear All,

I came across this posting (below) on a different list- how to get an Anova table from glmm.admb model objects- but I cannot get the car::Anova function to work for my model, or for the owl example I worked through from the glmmadmb package help notes, so I am guessing my code is wrong

When I tried it with the owl example I get the following message, the same as when I try for my model:

> car::Anova(fit_zipoiss)
Error in data.frame(df, teststat, p) :
  arguments imply differing number of rows: 5, 6
In addition: Warning message:
In Ops.factor(1, Nest) : | not meaningful for factors



(I am using glmm.admb v0.7.2.4)



Please could you advise, many thanks,

Paula



>From the r-sig-mixed-models list:
Le 20/03/2012 22:43, Ben Bolker a écrit : On 12-03-20 02:32 PM,
> Thomas Merkling wrote:
>>>> Dear Ben and other list members,
>>>>
>>>> - Is there any way to produce a Anova/deviance table for a
>>>> model fitted with glmmADMB ? I used the Anova() function from
>>>> the car library for glmer models, but it does not seem to
>>>> work with glmmadmb (I'm using glmmADMB 0.7) and I would like
>>>> only one p-value for each term and interaction and NOT one
>>>> p-value for each level of the interaction.
Ben's reply:
> car::Anova() should work now -- I had to add a model.frame() and a
> df.residual() method for glmmadmb objects.  (The df.residual
> number may be a little dodgy -- I'm not sure I counted the
> parameters right -- but I don't think it's actually used for much
> by default, cause you get Wald chi-square tests)




Paula Rosewarne, PhD researcher,
Faculty of Biological Sciences
Manton Building 8.17
Clarendon Way
University of Leeds
LS2 9JT
UK

Email: bspjr at leeds.ac.uk<mailto:bspjr at leeds.ac.uk>