挖掘全球模型未能收敛

时间:2015-02-26 14:26:22

标签: r mixed-models

我正在尝试在dredge(MuMIn)框架中运行具有多个变量组合的gamm模型,其中基于截断的TRUE / FALSE相关矩阵作为子集。

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

然而,我构建的gamm(或uGamm)完整模型,包括20多个预测变量(平滑项和固定效果,随机效应)不会收敛。这实际上意味着挖泥机不会运行,因为完整模型不会收敛,即使我对完整模型本身不感兴趣,而是在单个模型和一些组合中。

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

  

s(。)中的错误:s(。)尚不支持

有没有办法绕过不会收敛的完整模型并且无论如何都有我的变量组合?和/或有没有办法避免(.)的上述错误?

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

1 个答案:

答案 0 :(得分:2)

您可以通过在更简单的模型中替换公式来欺骗dredge,然后将该模型用作“全局模型”。将完整的模型公式分配给<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)