r - lme 函数中的 R 错误：外部函数调用中的 NA/NaN/Inf (arg 3)

Question

我在 nlme 包中运行线性混合模型。

control <- lmeControl(maxIter=100,opt = c("optim"))
lme(response ~ 0+factorA+covariate,random=~1|factorB,
               weights=varIdent(form= ~1|factorA),control=control),

并且，它出现以下错误。

Error in logLik.reStruct(object, conLin) :NA/NaN/Inf in foreign function call (arg 3)

和收敛误差一样吗？还是其他人？

谢谢

Answer 1

这不是一个完整的答案，但在评论中的格式会很糟糕。（如果问题得到解决，稍后将删除或编辑。）

到目前为止，我无法复制。这是一个可重现的示例：

set.seed(101)
d <- data.frame(covariate=runif(500),response=rnorm(500),
                factorA=factor(sample(1:5,size=500,replace=TRUE)),
                factorB=factor(sample(1:5,size=500,replace=TRUE)))
library("nlme")
control <- lmeControl(maxIter=100,opt = c("optim"))
m1 <- lme(response ~ 0+factorA+covariate,random=~1|factorB,
          weights=varIdent(form= ~1|factorA),control=control,
          data=d)

由于这有效，因此很难走得更远。

最好能提供一个可重现的例子。除此之外，您可以通过设置options(error=recover)或来调试/帮助调试debug(nlme:::logLik.reStruct)。进入浏览器后，您可以尝试以下命令

ls()
str(object)
str(conLin)
dput(object)
dput(conLin)

并酌情提供结果。