Question

我正在尝试使用以下代码优化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])

Answer 1

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

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

您可以使用

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

使用Julia Optim进行优化时出错

1 个答案: