问题描述
我正在拟合具有二项分布(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
的模型的诊断图:
不存在离散或零通胀问题,但分位数偏差很大。我可以在残差中看到一个轻微的模式,但我不知道我应该如何解释这个或我应该采取哪些步骤来改进模型拟合。
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)