在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。