Question

library(nlme)
fm2 <- nlme(height ~ SSasymp(age, Asym, R0, lrc),
            data = Loblolly,
            fixed = Asym + R0 + lrc ~ 1,
            random = Asym ~ 1,
            start = c(Asym = 103, R0 = -10, lrc = -8), verbose = TRUE)

**Iteration 1
LME step: Loglik: -115.9653, nlminb iterations: 1
reStruct  parameters:
     Seed 
-6.000592 
Error in nlme.formula(height ~ SSasymp(age, Asym, R0, lrc), data = Loblolly,  : 
  step halving factor reduced below minimum in PNLS step

在这个例子中，我稍微用起始值进行了融合，因此算法没有收敛。但是，我有兴趣从详细输出中提取对数似然值-115.9653。有没有办法做到这一点？

Answer 1

我只有一个手动解决方案：

设置浏览错误：

options(error = recover)

运行您的代码并选择最深的框架。您应该最终在抛出错误的函数的环境中。您可以通过以下方式访问可能性：

nlmeFit$logLik

并使用以下命令将其分配给全局环境

LL <<- nlmeFit$logLik

R：从nlme拟合的详细输出中提取信息

1 个答案: