我想知道是否有人知道如何精确地将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)，但是我没有对此发表评论的声誉点。 :)
最佳答案
从?lmeControl:
作为估计过程的一部分，lme必须使用非线性优化函数来估计方差-协方差参数。 nlminb()和optim()是R中两个主要的内置优化器:nlminb使用单个基础算法，而optim提供了多种算法。
很难先验地知道哪种非线性优化函数对特定数据集效果最佳。