LME（）错误-达到迭代限制

在指定交叉混合效果模型时，我尝试包括交互。但是，我收到以下错误消息：

Error in lme.formula(rate ~ nozzle, random = ~nozzle | operator, data = Flow) : 
nlminb problem, convergence error code = 1
message = iteration limit reached without convergence (10)

该模型具有以下特征：1. 3种喷嘴类型（固定效果）2. 5个操作员，每个操作员对3种喷嘴类型的燃油流量进行3次重复测量。

我被要求在模型中包括喷嘴类型和操作员之间的相互作用。这是我的模型代码：

flow.lme <- lme(rate ~ nozzle, error= nozzle|operator, data=Flow)

为什么我会收到此错误消息？

我还没有听说过error参数，lme也没有在文档中看到它。您确定这不是错字吗？但是，要回答您提出的问题：

尝试 ?lmeControl

设置maxIter，msMaxIter，niterEM，和/或msMaxEval参数高于默认值，可以解决这个问题。捕获来自的输出lmeControl，然后将该对象传递给的control参数lme。

要么...

新的默认优化器lme使用的是片状的。当我将其改回旧的优化器时，这类问题有一半的时间为我解决了。您可以通过将opt参数设置为lmeControl来实现'optim'。

因此，将其放在一起：

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

首先，这是一个方差分析模型，而不是混合模型。

其次，在我看来，您的模型未被识别。在方程式中，您有

最后一项对您拥有的15个观测值有15个单独的值。没有任何自由度可以在模型中获得任何其他条件。包括互动是一个不好的建议。您必须丢弃它们。即使将它们包括为交叉效果也无济于事，因为它们将与固定效果完全共线，并且无法估计。具有15个观测值的最大似然或REML模型没有意义；最大似然理论的渐近结果根本行不通：这是您要在犁过的田野上行驶的法拉利。