我在尝试最大化R中的对数似然函数时遇到问题。
我的函数有5个参数,其中我想约束3为正,而留下其他2(应该非常小,接近0,但正或负)作为无约束。
我怎样才能在R中这样做?到目前为止,我简单地在R中运行约束最大化,要求前3个参数大于0,并要求其他2大于非常小的负数(例如,-300)。 R中的代码如下所示:
optcinstMO<-constrOptim(c(1,1,1,1,1),f=neg.loglik2dLO, data1=x, data2=z, data3=y, method="Nelder-MEad",ui=rbind(c(1,0,0,0,0),c(0,1,0,0,0),c(0,0,1,0,0),c(0,0,0,1,0),c(0,0,0,0,1))ci=c(0,0,0,-300,-300))
其中 neg.loglik2dLO 是我正在最小化的负对数似然函数。在R中还有其他方法可以解决这个“混合”最大化问题吗?
此外,任何人都可以告诉我,通过使用我的“自制”方法,我向我的最大化介绍了什么样的错误(如果有的话)?