[英]Is there a workaround for non-convergence in glmmTMB models using the dredge function of MuMin?
我嘗試將 MuMin-Package 的疏通功能用於配備 package glmmTMB 的負二項式廣義線性混合 model。
因為我的完整模型未能收斂,所以我嘗試了此處描述的解決方法: Dredge with the global model failed to blend
But when I use a simplified model and rewrite the function in model$call$function, dredge ignores this change and uses the simplified model instead of the full model.
glmmTMB-package 的功能是否有另一種解決方法?
下面是一些示例代碼:
# The full_model does not converge
full_model <- glmmTMB(y ~ x1 * x2 * (x3 + x4 + x5 + x6) + (1|RE1/RE1.1/RE1.2) + (1|RE2), data = df, family = "nbinom2")
# The simple_model does converge
simple_model <- glmmTMB(y ~ x1 + x2 + x3 + x4 + x5 + x6 + (1|RE1/RE1.1/RE1.2) + (1|RE2), data = df, family = "nbinom2")
# Change formula in the model
simple_model$call$formula <- y ~ x1 * x2 * (x3 + x4 + x5 + x6) + (1|RE1/RE1.1/RE1.2) + (1|RE2)
# use dredge, but this ignores the changed formula
dredge(simple_model)
謝謝!
您必須替換simple_model $ modelInfo $ allForm
的元素。 這是三個公式: "formula"
、 "ziformula"
和"dispformula"
,但在您的 model 的情況下,僅使用第一個。
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.