我使用 R 中的 glmer 命令(来自 lme4 包)估计了随机系数风险模型。该命令如下所示。
(logit.full <-
glmer(event ~
+ V12 * I(V1 - 2)
+ V13
+ V9 * I(V5 - 2)
+ V11
+ V10
+ V2
+ V3
+ V4
+ V6 + V7 + V8
+ (1 + V6 + V7 + V8 | V14),
family=binomial("logit"), data=dataset,
verbose=TRUE, control=list(maxIter=400)))
该模型在过去三个月中一直运行良好。现在,将软件包更新到 1.0-4 后,“control”命令似乎有问题,我收到以下错误消息:
Warning in glmer(event ~ a1+a2+a3 :
Use control=glmerControl(..) instead of passing a list of class “list”
Error in function (optimizer = c("bobyqa", "Nelder_Mead"), restart_edge = FALSE, :
unused argument(s) (maxIter = 400)
有谁知道如何解决这个问题?
From ?glmerControl
:
optCtrl:要传递给的附加参数的“列表”
非线性优化器(参见“Nelder_Mead”、“bobyqa”)。在
特别是,“Nelder_Mead”和“bobyqa”都使用“maxfun”来
指定他们将进行的函数评估的最大数量
放弃之前先尝试 - 与“optim”相反
“optimx”包装的优化器,使用“maxit”。
诚然,这只是相当长且复杂的帮助页面的一小部分。
so control=glmerControl(optCtrl=list(maxfun=...))
应该这样做。
我可以看到这可能是一个常见问题解答,因此我们可以添加一些特殊代码来检测这种情况(和/或在文档中添加更突出的注释)。
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)