So I regularly make use of lme4 and lmerTest to fit Linear Mixed Models. I then do the post ANOVA contrasts via the emmeans (Estimated Marginal Means) package.
However, I understand that for Generalized Linear Mixed Models, the ANOVA (or analysis of deviance) table does not provide any p-values at all. Is it still valid to perform the contrasts without this?
One method I read about is Cluster Bootstrapping the GLMM. I found a nice package called ClusterBoot which seems to perform the cluster bootstrap. It was released recently: https://link.springer.com/article/10.3758/s13428-019-01252-y
I am able to get bootstrap coefficients, SEs, and percentile/parametric intervals for all the coefficients with this. However, this is difficult to interpret since I don't know how to use all of this to get ANOVA contrasts for my pairwise comparisons from this. Emmeans interfaces nicely with lme4 and provides this but can't be used with this.
How would I go about getting contrasts knowing just the Bootstrap Coefficients/SEs? Along with corresponding p-values.