I am using glmmTMB to analyse count data that suffers from both overdispersion and zero inflation, thus I am using a negative binomial GLMM with ziformula=~1
When I tried to check the residuals for model validation I got:
Error in residuals.glmmTMB(Mod_5, type = "pearson") : pearson residuals are not >implemented for models with zero-inflation or variable dispersion
Is there a way to obtain Pearson residuals by hand in R?
Thank you, Livia