使用 dHARMA 进行二项式 GLMM 的残差诊断

如何解决使用 dHARMA 进行二项式 GLMM 的残差诊断

我正在拟合具有二项分布（lme4::glmer 和 glmmTMB::glmmTMB）（n=341）的 GLMM。我根据这篇文章中的建议使用了这两个软件包：Interpretation of DHARMa residuals for Gamma GLMM。这是我安装 glmmTMB 的模型：

mod.ses1 <- glmmTMB(blast_stop ~
                  log(fshg_hhinc7days4.1e) +
                  ( 1 | village.1.01c),na.action = na.omit,data=Q1b.df,family = "binomial"(link = "probit")) 

还有 lme4：

mod.ses1A <- glmer(blast_stop ~
                 log(fshg_hhinc7days4.1e) +
                 ( 1 | village.1.01c),family = "binomial",control = glmerControl(optimizer = "bobyqa",optCtrl = list(maxfun = 100000)))

具有这两个软件包的模型产生了类似的问题。以下是装有 lme4 的模型的诊断图：

不存在离散或零通胀问题，但分位数偏差很大。我可以在残差中看到一个轻微的模式，但我不知道我应该如何解释这个或我应该采取哪些步骤来改进模型拟合。