我正在尝试使用lme4的开发版本和本教程对混合效果模型进行功效分析.我在教程中注意到lme4引发了收敛错误:
## Warning: Model failed to converge with max|grad| = 0.00187101 (tol =
## 0.001)
Run Code Online (Sandbox Code Playgroud)
当我运行数据集的代码时出现相同的警告,其中:
## Warning message: In checkConv(attr(opt, "derivs"), opt$par, checkCtrl =
control$checkConv, :
Model failed to converge with max|grad| = 0.774131 (tol = 0.001)
Run Code Online (Sandbox Code Playgroud)
来自此更新版本的常规glmer调用的估计值与我使用更新的CRAN版本时略有不同(在这种情况下没有警告).知道为什么会这样吗?
编辑
我试图指定的模型是:
glmer(resp ~ months.c * similarity * percSem + (similarity | subj), family = binomial, data = myData)
Run Code Online (Sandbox Code Playgroud)
我拥有的数据集有一个主体间(年龄,居中)和两个主体内变量(相似性:2个水平,percSem:3个水平)预测二元结果(虚假记忆/猜测).另外,每个受试者内细胞具有3个重复测量.因此,对于每个个体总共存在2×3×3 = 18个二元响应,总共38个参与者.
structure(list(subj = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, …Run Code Online (Sandbox Code Playgroud) 我一直在努力将缩放和居中的模型系数从glmer模型转换回未经中断和未缩放的值.
我在lme4(v1.1.7)包中使用GLMM分析了一个数据集.它涉及声学接收器的最大检测范围的计算和环境变量的影响.
样本数据:
dd <- structure(list(SUR.ID = c(10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, …Run Code Online (Sandbox Code Playgroud)