4

我想知道是否有人知道如何将 lme() 中的优化器精确设置为 opt='optim' 更改参数估计。

如本例所示:

ctrl <- lmeControl(opt='optim');
flow.lme <- lme(rate ~ nozzle, error= nozzle|operator, control=ctrl, data=Flow)

此处提出并回答了一个相关问题(https://stats.stackexchange.com/questions/40647/lme-error-iteration-limit-reached),但我没有声誉点来评论它。:)

4

1 回答 1

2

来自?lmeControl

opt:要使用的优化器,'"nlminb"'(默认)或'"optim"'。

optimMethod: character - 与 'optim' 优化器一起使用的优化方法。默认值为“BFGS”。另一种选择是“L-BFGS-B”。

作为估计过程的一部分,lme必须使用非线性优化函数来估计方差-协方差参数。 nlminb()并且optim()是 R 中两个主要的内置优化器:虽然nlminb使用单一的底层算法,optim但可以选择算法。

很难先验地知道哪种非线性优化函数对特定数据集最有效。

于 2014-10-14T03:19:06.227 回答