R：nlme fitで冗長な出力から情報を抽出する

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) 

はあなたのコードを実行し、最も深いフレームを選択します。エラーをスローする関数の環境で終わるはずです。 nlme - capture.outputは（（FM2 <を試してみてください - あなたはアウト< `使用して文字列をつかむことができ

LL <<- nlmeFit$logLik