使用 Julia Optim 进行优化时出错答案

【问题标题】：Error while doing optimization using Julia Optim使用 Julia Optim 进行优化时出错
【发布时间】：2014-12-30 09:27:39
【问题描述】：

我正在尝试使用以下代码优化 Julia 中的一个函数，但得到了 LoadError。谁能告诉我我的代码有什么问题？

data = readdlm("RGPSLessData.csv", ';')
using Optim
function f(x::Vector,n11,E)
     return sum(-log((x[5] * dnbinom(n11, x[1], x[2]/(x[2]+E)) + (1-x[5]) * dnbinom(n11, x[3], x[4]/(x[4]+E)))))
end
res = optimize(x -> f(x,data[:,1],data[:,2]),[0.2,0.06,1.4,1.8,0.1])

【问题讨论】：

看看newton.jl test from the Optim.jl package。
@rickhg12hs，感谢您的回复。我对这一切都很陌生。你能告诉我我必须做的所有改变吗？提前致谢
鉴于newton 方法示例与您的示例代码非常不同，我不清楚您想做什么。
相关问题：stackoverflow.com/questions/26158142/…
@Vincent Zoonekynd 感谢您的回复，我已经更新了我的代码.. 我已经在我的问题中更新了它，但它仍然显示 loadError 请帮帮我..

标签： optimization julia

【解决方案1】：

dnbinom 不是 Julia 中的函数（在 cmets 中它似乎来自 R）。概率分布由Distributions.jl 包提供。而不是

dnbinom(n11, x[1], x[2]/(x[2]+E))

你可以使用

pdf(NegativeBinomial(x[1], x[2]/(x[2]+E)), n11)

【讨论】：