如何用R求伽马分布的极大似然估计

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

参考帖子上的自己编的如下，总会出现问题，还有optim（c（2,1）。。。。）这个表示什么，新手求大神教！！
dat<-c(90,10,60,186,61,49,14,24,56,20,79,84,44,59,29,118,25,156,310,76,26,44,23,62,130,208,70,101,208)
f<-function(theta,dat){
sd<-theta[1]
beta<-theta[2]
n<-length(dat)
return(-n*log(beta^sd*gamma(sd))+(sd-1)*sum(log(dat))-beta^(-1)*sum(dat))
}
optim(c(2,1),f,dat=dat)

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：极大似然估计 极大似然 伽马分布 似然估计 如何用 如何

dat=c(90,10,60,186,61,49,14,24,56,20,79,84,44,59,29,118,25,156,310,76,26,44,23,62,130,208,70,101,208)
f<-function(theta){
sd<-theta[1]
beta<-theta[2]
n<-length(dat)
return(-n*log(beta^sd*gamma(sd))+(sd-1)*sum(log(dat))-beta^(-1)*sum(dat))
}

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

dat=c(90,10,60,186,61,49,14,24,56,20,79,84,44,59,29,118,25,156,310,76,26,44,23,62,130,208,70,101,208)
f<-function(theta){
sd<-theta[1]
beta<-theta[2]
n<-length(dat)
return(-n*log(beta^sd*gamma(sd))+(sd-1)*sum(log(dat))-beta^(-1)*sum(dat))
}

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

zhou1_20 发表于 2016-4-19 21:44
dat=c(90,10,60,186,61,49,14,24,56,20,79,84,44,59,29,118,25,156,310,76,26,44,23,62,130,208,70,101,208 ...

算出来的结果还是跟文献中不一样，而且conv是1

路过，帮顶

路过，帮顶！