全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
8420 9
2016-01-12
复制代码

#Step1:生成10万个服从gamma分布的数,shape=2,rate=0.2
复制代码
QQ截图20160112150438.jpg

#Step2:运用MASS包中的fitdistr()求出参数shape、rate
复制代码
#Step3:将actual graph和fit graph画在一张图上
复制代码
QQ截图20160112151154.jpg
#问题就出现了,参数值相当完美,画出的图为什么重合度那么差,所以,我的问题就是:
1.拟合曲线是这样画吗,我感觉应该不是,还希望能请教下大家;

2.对于任何一个模型,求出参数之后,还请以这个例子为例,帮我大概解答下怎么评估所建的模型的拟合好坏,涉及到哪些常用的检验

最后,感觉这个问题应该不难,自己谷歌了一两天也没解决,总结一下,应该是基础不够扎实,对于这类问题还没有建立系统的思想来一步步深入,所以在边科研边补充统计学的知识,希望大家能多指教一些,解决这个问题固然重要,还是希望从问题中学到更多衍生的知识,在此先谢谢各位了

二维码

扫码加我 拉你入群

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

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

全部回复
2016-1-13 18:21:45
你的也太复杂了,不对呀,拟合的数据直接用核密度估计曲线啊,应该是你求拟合数据的频率有问题
x <- rgamma(100000,shape = 2,rate = 0.2)
hist(x,freq=F,xlim=c(0,60),ylim=c(0,0.1))
lines(density(x),col="red")
curve(dgamma(x,2,rate=0.2),add=T)
二维码

扫码加我 拉你入群

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

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

2016-1-13 18:23:58
你试试这个吧,检验的话,用分布拟合检验,有点复杂,简单的就是ks检验
二维码

扫码加我 拉你入群

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

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

2016-1-13 20:34:39
sunshine668890 发表于 2016-1-13 18:21
你的也太复杂了,不对呀,拟合的数据直接用核密度估计曲线啊,应该是你求拟合数据的频率有问题
x
嗯,自己对于密度的求法确实存在问题,应该就是个错误,学习了,谢谢你的解答
二维码

扫码加我 拉你入群

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

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

2016-1-13 20:35:34
sunshine668890 发表于 2016-1-13 18:23
你试试这个吧,检验的话,用分布拟合检验,有点复杂,简单的就是ks检验
嗯,明白了,ks.test(),已经慢慢摸索到了,谢谢你
二维码

扫码加我 拉你入群

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

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

2016-1-14 10:57:06
kizgel 发表于 2016-1-13 20:35
嗯,明白了,ks.test(),已经慢慢摸索到了,谢谢你
不客气,我也是刚学完R这门课
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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