R求极大似然估计

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

正在试图入门，求大神指导，怎么在R中求GAMMA函数的极大似然估计。
我写了这个之后提示错误：
ilprec<-scan("E:\\作业\\非参\\illinois601234.txt")
Read 227 items
> a=sum(ilprec)
> a
[1] 50.937
> b=prod(ilprec)
> b
[1] 5.974125e-293
> f<-function(params,ilprec){
+ lambda=params[1]
+ alpha=params[2]
+ n=227
+ L<-n*alpha*log(lambda)-n*log(Gamma(alpha))+(alpha-1)*log(b)-lambda*a
+ return(-L)}
> result=optim(c(1,2),f,ilprec)
提示说：Error in Gamma(alpha) :
  link "alpha" not available for gamma family; available links are ‘inverse’, ‘log’, ‘identity’
求各位指导~
小的感激不尽。。

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：极大似然估计 极大似然 似然估计 gamma prod

路过

dat <- rgamma(100, 3, 2)
f <- function(x, dat) {
    a <- x[1]
    b <- x[2]
    return(-sum(a*log(b)-log(gamma(a))+(a-1)*log(dat)-b*dat))
}

optim(c(2, 1), f, dat=dat)