疏通全局模型无法收敛

Question

我正在尝试运行在疏浚 (MuMIn) 框架中使用多个变量组合的 gamm 模型，并使用基于截止 TRUE/FALSE 的相关矩阵作为子集。

问题是，我的完整模型非常复杂，具有随机效应，>20 预测变量包括 3 个固定效应，因此不会收敛。我的预测变量是高度相关的，如前所述，我包含了一个 TRUE/FALSE 相关矩阵作为子集。由于相关性截止值较低，我预计组合很少，对完整模型并不感兴趣。

但是，我构建的包含 20 多个预测变量（平滑项和固定效应、随机效应）的 gamm（或 uGamm）完整模型并未收敛。这实质上意味着 dredge 不会运行，因为完整模型不会收敛，即使我对完整模型本身不感兴趣，而是对单个模型和一些组合感兴趣。

我似乎也无法使用 fm1 <- gamm(y ~ (.), data = Cement) 选项，因为这会产生以下错误：

Error in s(.) : s(.) not yet supported

有没有办法绕过不会收敛的完整模型并且不管我的变量组合如何？ and/Or 有没有办法避免 (.) 的上述错误？

有什么建议吗？我的数据太大，无法包含在这里。

Answer 1

您可以通过在更简单的模型中替换公式来欺骗 dredge，然后将该模型用作 "global model"。将完整的模型公式分配给 <gam.model> $ gam $ formula。例如：

# fit some simple model
fm <- uGamm(y ~ s(x1), ...) # add other arguments as in the proper full model
fm$gam$formula <- y ~ x0 + s(x1) + s(x2) + ... + s(xn)
dredge(fm)

疏通全局模型无法收敛

Dredge with the global model failing to converge

r

mixed-models