是否有使用 MuMin 的疏浚功能解决 glmmTMB 模型中的非收敛问题的解决方法?
Is there a workaround for non-convergence in glmmTMB models using the dredge function of MuMin?
我已经尝试使用 MuMin-Package 的 dredge-function 来构建负二项式广义线性混合模型,该模型装有 glmmTMB 包。
因为我的完整模型无法收敛,我已经尝试了此处描述的解决方法:
但是当我使用简化模型并重写 model$call$function 中的函数时,dredge 会忽略此更改并使用简化模型而不是完整模型。
对于 glmmTMB 包的功能是否有其他解决方法?
下面是一些示例代码:
# 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"
,但对于您的模型,只使用第一个。
我已经尝试使用 MuMin-Package 的 dredge-function 来构建负二项式广义线性混合模型,该模型装有 glmmTMB 包。
因为我的完整模型无法收敛,我已经尝试了此处描述的解决方法:
但是当我使用简化模型并重写 model$call$function 中的函数时,dredge 会忽略此更改并使用简化模型而不是完整模型。
对于 glmmTMB 包的功能是否有其他解决方法?
下面是一些示例代码:
# 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"
,但对于您的模型,只使用第一个。