全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
4389 3
2014-09-24
有一个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

求问大神,各位有会的吗?

急急急!
二维码

扫码加我 拉你入群

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

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

全部回复
2014-9-24 09:03:19
说得有理哦       
              
                          
                  
                   
                                 
        
二维码

扫码加我 拉你入群

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

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

2015-1-8 20:42:39
拟合伽马分布可以这么来做——

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")

二维码

扫码加我 拉你入群

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

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

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

x.gam
谢谢!
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群