GLMM 和两个斜坡

机器算法验证 r 物流 混合模式 咕噜咕噜 lme4-nlme
2022-03-22 01:35:42

我的结果变量是二项式的，我有 11 个自变量和一个时间变量。时间变量有不同的斜率，所以我将它固定为time-before和time-after。我使用了lme4包（glmer功能）。我有一个随机截距和两个随机斜率。我创建了这样的模型：

m3.glmm <- glmer(y ~ timebefore + timeafter + x1 + x2 +...+ x11 +     
(1+timebefore+timeafter|id),
             data = data, family = binomial (link="logit"), nAGQ=3)

当我使用这个模型时，我遇到了这个错误：

Error in updateGlmerDevfun(devfun, glmod$reTrms, nAGQ = nAGQ) : 
  nAGQ > 1 is only available for models with a single, scalar random-effects term

任何人都对如何拟合（或编码）这个模型有一个简单的解释？

1个回答

问题是您为对数似然的自适应高斯正交 (AGQ) 近似选择的横坐标或节点数，由 指定nAGQ。默认值为 1（相当于拉普拉斯近似）。

该glmer函数的详细信息部分（lme4 帮助页面中的第 29 页）指出：

GLMM 最可靠的近似是自适应 Gauss-Hermite 正交，目前仅用于具有单个标量随机效应的模型。

将 AGQ 近似限制为单个标量随机效应不是 AGQ 的限制，但似乎是lme4包作者做出的决定，正如Douglas Bates 在 2011 年在此指出的那样（以下引用的相关文章）：

似乎可以通过在 glmer 中加入自适应 Gauss-Hermite 方法来解决这个问题……在 glmer 的版本中有这样的方法，但仅适用于非常特定的模型。我们将添加它，但现在我们正专注于开发中的其他问题。

所以，为了让你的代码执行，我相信设置nAGQ为 1 会起作用。

其它你可能感兴趣的问题
上一篇没有距离矩阵的聚类 下一篇Calinski-Harabasz 指数背后的直觉