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

[Mollie Brooks]

Hi Guillaume,

I’m not very experienced with Anova, but I know the development version of glmmTMB supports car::Anova. Ben recently added this functionality.

To install the development version try devtools::install_github("glmmTMB/glmmTMB/glmmTMB")
or see further instructions here https://github.com/glmmTMB/glmmTMB

After installing, check out vignette("model_evaluation")

cheers,
Mollie

> On 30Jul 2018, at 10:59, Guillaume Adeux <guillaumesimon.a2 using gmail.com> wrote:
> 
> 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 using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models


	[[alternative HTML version deleted]]