在R中使用Optim命令时出错

Question

logl = function(theta, x, y)
{
b0 <- theta[1]
b1 <- theta[2]
sig.hat <- theta[3] 
b2 <- theta[4]
n = length(x)
logl <- .5*n*log(2*pi)-.5*n*(sig.hat) - (1/(2*(exp(sig.hat))))*sum((y - (b0+b1*x+b2*(x^2))^2))
}
theta.start = c(.015,.02,3,1)
optim(theta.start, fn=logl, x=x, y=y, control=list(fnscale=-1), hessian = T)

这给了我一个错误：

Error in optim(theta.start, fn = logl, x = x, y = y, control = list(fnscale = -1),  : 
  non-finite finite-difference value [1]

在我不小心让R尝试记录负数之前我得到了这个错误，但现在我不知道是什么问题？

Answer 1

设置fnscale = -1时，它会成为最大化问题。快速浏览目标函数似乎表明您的问题是无限制的，即您可以继续增加b0以增加logl。