I am a newbie here and my question is whether I should use a parametric or non-parametric post-hoc test based on the results from a generalized linear model; and if non-parametric is appropriate, how to conduct it.
Because this is about concepts, I'm not posting any data (unless it's necessary). I fitted generalized linear mixed-effect models with negative binomial function in Rstudio with lme4 package (glmer.nb model).
I understand that the GLM is for non-parametric data, but if I want to run a follow-up post-hoc test based on the model results, do I use parametric or non-parametric test?
I am able to conduct Tukey-adjusted parametric test of the model results using emmeans package (emmeans command) and multcomp package (cld command) to analyse significant differences between treatments. I.e. emmeans(modelX, "treatments").
However, if this is not appropriate and non-parametric test should be conducted instead, can someone enlighten me on how to do it using the model results and/or accounting for the random effects?
Thank you in advance!