问lmer的替代优化算法

EN

您不能使用newuoa或uobyqa，因为它们都不允许对参数进行约束。来自?lmerControl (重点添加)

任何允许框约束的最小化函数都可以使用，只要它 (1)接受输入参数fn(要优化的函数)、“par”(起始参数值)、“下限”和“控制参数”(从“控制”参数传递过来的控制参数)和 (2)返回一个包含(至少)元素‘par’(最佳拟合参数)、‘fval’(最佳拟合函数值)、‘conv’(收敛码等于零的成功收敛)和(可选)‘message’(信息消息或收敛失败解释)的列表。

"bobyqa“开头的b表示”绑定“(如在约束中)，我假设其他算法中的u类似地表示”无约束“。您可以查看此文件，以便某些机器(重新)可以将同一模型与一组不同的优化器相匹配：

allFit <- system.file("utils", "allFit.R", package="lme4")
file.show(allFit)

我目前所知道的允许框约束的所有优化器列表不需要指定显式梯度函数( optimx包中的大多数绑定约束优化器都需要)，如上面的文件所示

• BOBYQA (minqa和nloptr包实现)
• Nelder (lme4、nloptr和dfoptim包实现)
• 来自R基的nlminb (来自贝尔实验室端口库)
• L-BFGS-B来自R基地，途经optimx (Broyden Goldfarb-Shanno，经纳什)

除了这些内置于allFit.R的功能之外，您还可以使用来自nloptr的COBYLA或子丛优化器:参见?nloptwrap。在subplex包中还有另一个子丛的实现:可能还有其他几个我错过了。