我使用optim()进行ML估计,我想知道是否有控制选项来限制功能评估的数量?否则就会继续前进。
我目前的代码如下:
fit <- mle(LL, start = sv, method = "BFGS", control=list(maxit=10, reltol=1))
谢谢!
答案 0 :(得分:0)
您可以在stop()
功能中添加全局计数器变量和browser()
或LL
来电:
maxEvals<-0
LL<-function(...) {
maxEvals<-maxEvals+1
if (maxEvals>100)
stop('Maximum evaluations reached')
}
但即使ML的单次迭代花费太长时间,您也可能应该检查您的功能和数据。