如何用R求伽马分布的极大似然估计 - 经管之家

› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛

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

5641

5

收藏 2016-04-19

参考帖子上的自己编的如下，总会出现问题，还有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)

二维码

扫码加我拉你入群

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

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

全部回复

2016-4-19 21:43:51

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)

二维码

扫码加我拉你入群

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

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

2016-4-19 21:44:33

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)

二维码

扫码加我拉你入群

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

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

2016-4-25 09:52:17

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

二维码

扫码加我拉你入群

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

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

2016-5-3 19:42:45

路过，帮顶

二维码

扫码加我拉你入群

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

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

2016-5-3 21:26:25

路过，帮顶！

二维码

扫码加我拉你入群

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

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

相关推荐

栏目导航

热门文章

推荐文章

扫码加好友，拉您进群

各岗位、行业、专业交流群