Anova (type III-tests) table based on LRT for glmmTMB models
Guillaume, 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 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 GRECC (BT/18/GR) 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 at r-project.org> on behalf of Guillaume Adeux <guillaumesimon.a2 at gmail.com>
Sent: Monday, July 30, 2018 4:59 AM
To: r-sig-mixed-models at r-project.org
Subject: [R-sig-ME] Anova (type III-tests) table based on LRT for glmmTMB models
Sent: Monday, July 30, 2018 4:59 AM
To: r-sig-mixed-models at 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.
Anova(mod,mod1)).
-
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 at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models