I would like to see the joint effects of an interaction term. I tried joint_tests() but it gives f ratios. However, I would prefer Incidence Rate Ratios (IRR) because the rest of my results section uses those. My model is a GLMM with Poisson distribution. Any ideas?
I also tried the following but the outcome remains the same:
EMM <- emmeans(model1, ~ A | B) PRS <- contrast(EMM, "revpairwise", interaction = "consec") summary(PRS) summary(PRS, type = "response", infer = c(TRUE, TRUE))
Thank you in advance!