这里有不等式约束条件哈，用Rsolnp包，里面的solnp可以加不等式——

solnp(pars, fun, eqfun = NULL, eqB = NULL, ineqfun = NULL, ineqLB = NULL,
ineqUB = NULL, LB = NULL, UB = NULL, control = list(), ...)
   
这里面的——
ineqLB = NULL,这是lower那个点
ineqUB = NULL,这个是upper那个点

我也碰到过这样得情况，我个人觉得可以改变一下初值，是因为程序从这个程序出发在做最优化的时候，正好碰上某个值使得你的目标函数无法求解，比如说有可能使得分母等于0了，所以你可以试试调试初值。