rsolnp等式约束不满足

Question

为什么rsolnp不满足我的等式约束？

require(Rsolnp)
f=function(p){sum(p*log(p))}
g=function(p){sum(p)-1}
p0=runif(10)
p0=p0/sum(p0)
res=solnp(p0,f,g)
print(g(res$par))

返回

[1] 2.678794

Answer 1

您需要定义eqB的{​​{1}}参数：

solnp