R语言中mle如何修改？急！谢谢！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

LL<-function(params,data)
  {
theta<-params[1]
a<-params[2]
b<-params[3]
ll<-log(1-theta+theta/sum(prod(a+b/i)))+n*log(1-theta+theta/sum(prod(a+b/i))-1+theta)+sum((n-i+1)*(a+b/i))
#log-likelihood函数
return(-ll)

  }

geyser.res<-nlminb(c(0.5,0.5,0.2),LL=0,data=1:100,lower=c(-1,-1,0),upper=c(1,1,1))
Error in objective(.par, ...) : 缺少参数"objective",也没有缺省值
geyser.res$par

这个地方要如何修改呢？？？？谢谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MLE R语言

Sunzhuowei789 发表于 2018-4-1 17:02
LL

你这没有函数啊

Sunzhuowei789 发表于 2018-4-1 17:02
LL

LL的地方不用赋值吧

中国梦丶 发表于 2018-4-1 17:26
你这没有函数啊

请问如果函数是r中没有的 就是自己定义的那种该怎么办呢
密度函数是：p(i)=p(i-1)(a+b/i)

Sunzhuowei789 发表于 2018-4-1 18:30
请问如果函数是r中没有的 就是自己定义的那种该怎么办呢
密度函数是：p(i)=p(i-1)(a+b/i)

> LL<-function(params,data)
+   {
+ theta<-params[1]
+ a<-params[2]
+ b<-params[3]
+ n<-1000
+ ll<-log(1-theta+theta/sum(prod(a+b/i)))+n*log(1-theta+theta/sum(prod(a+b/i))-1+theta)+sum((n-i+1)*(a+b/i))
+ #log-likelihood函数
+ return(-ll)
+
+   }
>
> geyser.res<-nlminb(c(0.5,0.5,0.2),LL,data=1:100,lower=c(-1,-1,0),upper=c(1,1,1))
There were 50 or more warnings (use warnings() to see the first 50)
> geyser.res$par
[1]  0.758887894 -0.000199301  0.199500293

这个该怎么修改呢？

Sunzhuowei789 发表于 2018-4-1 18:33
&gt; LL

没懂你意思

Sunzhuowei789 发表于 2018-4-1 17:02
LL

函数那里就是你自己写的似然函数啊