An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20120320/45ddbe3a/attachment-0001.pl>
Anova/deviance table in glmmADMB
3 messages · Thomas Merkling, John Fox, Ben Bolker
Dear Thomas, The Anova() function in the car package doesn't have a specific method for objects produced by glmmADMB() (with which I'm unfamiliar). If there are coef() and vcov() methods for those objects, then the default Anova() method should produce Wald tests. Best, John ------------------------------------------------ John Fox Sen. William McMaster Prof. of Social Statistics Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox/ On Tue, 20 Mar 2012 19:32:11 +0100
Thomas Merkling <thomas.merkling at cict.fr> 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. - I also tried to compare nested models via LRT using anova(model1,model2) for glmmADMB models but I didn't get any p-values, it produced NaN, although it apparently calculate the differences in likelihood. Where can this come from ? - I tried to use the coefplot2 library (downloaded from http://www.math.mcmaster.ca/bolker/R/src/contrib) as showed on the glmmADMB home page ("getting started with glmmADMB"), but I couldn't use it and got this error message : unable to find an inherited method for function "coefplot2", for signature "glmmadmb". Is it normal ? It worked fine with glmer output. Thanks for your help ! Thomas -- ****NEW ADDRESS AND PHONE NUMBER **** Thomas Merkling, Doctorant (PhD Student) Web Page <http://www.edb.ups-tlse.fr/Merkling-Thomas.html> Laboratoire "Evolution et Diversit? Biologique" -EDB UMR 5174 - b?t 4R1 - bureau 33 RDC Universit? Paul Sabatier Toulouse 3 118, route de Narbonne 31062 TOULOUSE Cedex O9, FRANCE T?l: 33 5-61-55-67-56 Fax: 33 5-61-55-73-27 [[alternative HTML version deleted]]
-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1
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.
Hmmm. 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) The newest version of glmmADMB is 0.7.5.1 -- it may take a little while for r-forge to catch up, and I was having a few dependency issues. If you don't see it there in ~ 24 hours, drop me a line.
- I also tried to compare nested models via LRT using anova(model1,model2) for glmmADMB models but I didn't get any p-values, it produced NaN, although it apparently calculate the differences in likelihood. Where can this come from ?
This ought to have worked: perhaps you had the models in the wrong order?
- I tried to use the coefplot2 library (downloaded from http://www.math.mcmaster.ca/bolker/R/src/contrib) as showed on the glmmADMB home page ("getting started with glmmADMB"), but I couldn't use it and got this error message : unable to find an inherited method for function "coefplot2", for signature "glmmadmb". Is it normal ? It worked fine with glmer output.
Can you try getting it from r-forge? Again, if the problem persists let me know. Ben
Thanks for your help ! Thomas
-----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.10 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/ iQEcBAEBAgAGBQJPaPoCAAoJED2whTVMEyK9qAYIAJZKoj3hh8fdTR7U22I64HhM SmLeTSCmauRzlQXbgSeKx720UGsJtGUzu9Sw9AHiTKqpGrujOj48balEf1MF3Crk vkTsdDM7Pp5y36xsqw6Ps122DnDQ30ctYh6IDD0+XJX7BebEK5P9or4xoy1lvbOj L6STkB3wmCYKtJMhFJkvvmCu/+S4dSzOb148edSuO0WFwIC4+Gax7UtHhoj1vgkm 2Gfz3uHre+1aEh0CantGyFgAYBxEIA0F/5AloWQLCRCtqSIi8eNCS0HBt9hmzxLW 455R8PXJyVmtrRn2fWO8Fa5/Two56WMeQwYV8DenBtydM+bGe0NHTEiirHtvf/g= =dVAk -----END PGP SIGNATURE-----