【急急急】R 语言根据一个table 模拟 Gamma 分布

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

有一个table，如下：

"table"|"lower"|"upper"|"frequency"|"average"
1|0|2.5|41|1.389
2|2.5|7.5|48|4.661
3|7.5|12.5|24|9.991
4|12.5|17.5|18|15.482
5|17.5|22.5|15|20.232
6|22.5|32.5|14|26.616
7|32.5|47.5|16|40.278
8|47.5|67.5|12|56.414
9|67.5|87.5|6|74.985
10|87.5|125|11|106.851
11|125|225|5|184.735
12|225|300|4|264.025
13|300|NA|3|300


其中，| 代表 分隔符号。
请据此模拟一个gamma 分布出来。（用最大似然估计）

老师提示： multinomial

求问大神，各位有会的吗？

急急急！

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：gamma Table ABLE GAM ABL frequency average

说得有理哦       
              
                          
                  
                   
                                 
        

拟合伽马分布可以这么来做——

x.gam<-rgamma(200,rate=0.5,shape=3.5) library(stats4) ## loading package
stats4 ll<-function(lambda,alfa) {n<-200 x<-x.gam -n*alfa*log(lambda)+n*log(gamma(alfa))-(alfa- 1)*sum(log(x))+lambda*sum(x)} ## -log-likelihood function
est<-mle(minuslog=ll, start=list(lambda=2,alfa=1))
summary(est)
  
对于常用的分布形式，也可以直接用MASS包里的fitdistr函数直接拟合 如
library(MASS) ## loading package MASS
fitdistr(x.gam,"gamma")

DM小菜鸟 发表于 2015-1-8 20:42
拟合伽马分布可以这么来做——

x.gam

谢谢！