[R-sig-ME] Anova (type III-tests) table based on LRT for glmmTMB models

Sorkin, John j@orkin @ending from @om@um@ryl@nd@edu
Mon Jul 30 17:58:08 CEST 2018


Although not a perfect answer to your question, and I am not certain it will work with glmmTMB, I suggest you look at the drop1 function. It may give you an answer that you can use.


John David Sorkin M.D., Ph.D.
Professor of Medicine
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division of Gerontology and Geriatric Medicine
Baltimore VA Medical Center
10 North Greene Street
Baltimore, MD 21201-1524
(Phone) 410-605-7119
(Fax) 410-605-7913 (Please call phone number above prior to faxing)

From: R-sig-mixed-models <r-sig-mixed-models-bounces using r-project.org> on behalf of Guillaume Adeux <guillaumesimon.a2 using gmail.com>
Sent: Monday, July 30, 2018 4:59 AM
To: r-sig-mixed-models using r-project.org
Subject: [R-sig-ME] Anova (type III-tests) table based on LRT for glmmTMB models

CAUTION: This message originated from a non UMB, UMSOM, FPI, or UMMS email system. Whether the sender is known or not known, hover over any links before clicking and use caution opening attachments.

Hello everyone,

I'm looking for a method/function in order to produce an Anova table based
on Likelihood Ratio Tests (LRT) for a glmmTMB model (R software). In my
case it is with a beta distribution and log link. My response is a ratio
(%) in a repeated measures design.


   the function Anova() from the {car} package doesn't not run on single
   models (i.e. Anova(mod)). It only allows comparison of two models (i.e.

   for glmer models, I was used to using the mixed() function from the
   {afex} packages which produced Anova tables (type III tests) based on LRT
   (or parametric bootstrap) for glmms.

Could anyone shed their on light on a function like mixed() which would run
on glmmTMB objects or on a procedure to do this by hand?

I suppose if only one fixed predictor was present in the model, this would
be simple by comparing it to a null model but my model contains an
interaction. Hence, I am incable of comparing a model A+B+B:C with a model
containing A+B:C.

Thanks for your interest.

Guillaume ADEUX

        [[alternative HTML version deleted]]

R-sig-mixed-models using r-project.org mailing list

	[[alternative HTML version deleted]]

More information about the R-sig-mixed-models mailing list