The glmmTMB package has a function Anova.glmmTMB that reports the anova table with chi-square tests. Is this not sufficient for your purposes?

Example:

library(glmmTMB)

m2

Anova.glmmTMB(m1)

does not work.

my code is:

m1

Anova(glmmTMB(m1))

Error in I.p[c(subs.relatives, subs.term), , drop = FALSE] :

subscript out of bounds

In addition: Warning messages:

1: In is.na(coef(mod)) :

is.na() applied to non-(list or vector) of type 'NULL'

2: In model.matrix.default(model) :

the response appeared on the right-hand side and was dropped

3: In model.matrix.default(model) :

problem with term 11 in model.matrix: no columns are assigned

Does anyone know how to fix any of these errors?

Does the example in my response, m2 with the Salamanders data, work?

An update,

I updated ALL my packages, and it was worth it. It now works.

Thank you Sal Mangiafico

m1

After updating all my packages, it still does not work for me. Anyone can help me?

I'd ask the same question, can you get the Salamanders example above to work?

Am probably too late by now, glmmTMB:::Anova.glmmTMB(model, ...) works.

https://rdrr.io/cran/glmmTMB/man/downstream_methods.html

Just using Anova() from library(car) works for me.

mod

glmmTMB:::Anova.glmmTMB(m1)

@ Laurent De Vriendt you are a life saver, my anova for glmmtmb wasn;t working until I saw your comment with the interactions type 3

You can get the chi-squared result using the following:

glmmTMB:::Anova.glmmTMB(model, type = c("II", "III", 2, 3), test.statistic = c("Chisq", "F"), component = "cond", vcov. = vcov(model))

This will give you both the chi-squared result and the f-stat.

You can also use the "car" package as follows:

library(car)

Anova(model)

There you go!

Best wishes

Christoph