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

Question

我正在拟合具有二项分布 (lme4::glmer和glmmTMB::glmmTMB) (n=341) 的 GLMM。我根据这篇文章中的建议使用了这两个包：解释 Gamma GLMM 的 DHARMa 残差。这是我安装的模型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), na.action = na.omit, data=Q1b.df, family = "binomial", control = glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 100000)))

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

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