缺少优化函数参数答案

【问题标题】：optim function argument missing缺少优化函数参数
【发布时间】：2013-12-27 01:48:34
【问题描述】：

这是我的代码。 kum.loglik 函数返回负对数似然并接受两个参数 a 和 b。我需要找到使用 optim 函数最小化此函数的 a 和 b。（n1,n2,n3 是预先指定的并传递给优化函数。

kum.loglik = function(a, b, n1, n2, n3) {
  loglik = n1*log(b*beta(1+2/a,b)) + n2 * log(b*beta(1+2/a,b)-2*b*beta(1+1/a,b)+1) +
    n3 * log(b*beta(1+1/a,b)-b*beta(1+2/a,b))
  return(-loglik)
}
optim(par=c(1,1), kum.loglik, method="L-BFGS-B",
      n1=n1, n2=n2, n3=n3,
      control=list(ndeps=c(5e-4,5e-4)))

这段代码应该可以正常工作，但它给出了错误消息

Error in b * beta(1 + 2/a, b) : 'b' is missing

这段代码有什么问题？

【问题讨论】：

标签： r optimization

【解决方案1】：

问题是（直接来自优化帮助）：

fn: A function to be minimized (or maximized), with first
argument the vector of parameters over which minimization is
to take place.

您的kum.loglik 函数需要采用一个向量v，您可以从中提取参数，例如：

kum.loglik=function(v) { a = v[1]; b = v[2]; ...}

【讨论】：

【解决方案2】：

我总是使用以下，它会给你最好的结果

p0 <- c(a,b) #example vector of starting values
m <- optim(p0, loglik, method="BFGS", control=list(fnscale=-1, trace=10),
hessian=TRUE, x=data.frame)
#for table of results
rbind(m$par, sqrt(diag(solve(-m$hessian))))

【讨论】：