求助：使用极大似然估计，如何设定估计参数的上下界？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

论坛里面的各位高手，菜鸟请教个问题：
用极大似然估计，如何设定估计参数的上下界？  下面是写的程序，我想设定要估计的参数param为[-1,1],可是R里面的极大似然估计包里面的函数好像都没有设置参数的上下界！在这里请赐教啦！
f=read.table("f:/data1.txt",header=TRUE)   #导入数据
cdf1=ecdf(f$x)
cdf2=ecdf(f$y)
u1=cdf1(f$x)
u2=cdf2(f$y)  #得出经验分布函数值
n=length(u1)
  #对u1,u2的数据进行处理，防止出现求逆的时候出现无穷大，这里把1约等于0.99
     for(i in 1:n)
      { if(u1[i]==1)  u1[i]=0.99
        if(u2[i]==1)  u2[i]=0.99
       }      
y1=qnorm(u1)
y2=qnorm(u2)
# 进行极大似然估计
library('maxLik')
loglik=function(param){
ll=-n*log(1-param^2)/2-(1/(2*(1-param^2)))*(sum(y1^2)-2*param*sum(y1*y2)+sum(y2^2))
  ll
}
res=maxLik(loglik,start=0,method="BFGS")

运行后出现：
Warning messages:
1: In log(1 - param^2) : 产生了NaNs
2: In log(1 - param^2) : 产生了NaNs
3: In log(1 - param^2) : 产生了NaNs
4: In log(1 - param^2) : 产生了NaNs
5: In log(1 - param^2) : 产生了NaNs
6: In log(1 - param^2) : 产生了NaNs
7: In log(1 - param^2) : 产生了NaNs
8: In log(1 - param^2) : 产生了NaNs
9: In log(1 - param^2) : 产生了NaNs
所以这里肯定是参数的界限出了问题！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：极大似然估计 似然估计 极大似然 messages function 求助 参数 极大似然估计 下界

难道没有人知道么？看来得用其他函数啦！

同问这个问题

让我想想......

楼主，这个问题解决了吗，同求！

qinglingxuanfen 发表于 2014-11-24 11:12
楼主，这个问题解决了吗，同求！

亲，你这个问题解决了吗？我也想问！！

用optim函数可以