全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2330 4
2013-02-03
请问一下,怎么样把卡方分布的图加到直方分布图上面啊? 我要目测该分布是否是自由度为2的卡方分布。 谢谢大家
二维码

扫码加我 拉你入群

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

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

全部回复
2013-2-3 07:35:13
最笨的方法就是用lines(x, dchisq(x, 2))了,x取得足够密即可。
二维码

扫码加我 拉你入群

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

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

2013-2-3 08:14:30
hugebear 发表于 2013-2-3 07:35
最笨的方法就是用lines(x, dchisq(x, 2))了,x取得足够密即可。
你好。 我的问题是,我用hist()做出了一个直方图,然后我要在直方图上加一条自由度为2的卡方分布线,请问如何实现。
二维码

扫码加我 拉你入群

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

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

2013-2-3 13:24:33
tiamojuan 发表于 2013-2-3 08:14
你好。 我的问题是,我用hist()做出了一个直方图,然后我要在直方图上加一条自由度为2的卡方分布线,请问 ...
。。。我不是已经在上一帖中告诉你了么。。。
好吧,我再说详细点。
当你用hist()画完直方图后。我假设你横坐标的范围是c(1, 50)。然后:
x <- seq(1, 50, length = 1000)  #取1000个点,用以画卡方分布的密度曲线
lines(x, dchisq(x, df = 2)) #lines绘出这条曲线,如有疑问请使用帮助文件
二维码

扫码加我 拉你入群

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

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

2013-2-3 13:29:17
参:http://www.r-bloggers.com/multiple-y-axis-in-a-r-plot/
par(mar=c(5, 6, 4, 4) + 0.1)
x<-rnorm(100)
hist(x,xlab="",ylab="freq.",main="")
par(new=T)
x<-sort(x)
y<-dchisq(x,2)
plot(x,y,"l",axes=F,
xlab="x",ylab="",
lty=2)
axis(4,ylim=c(0,max(y)))
mtext(side=4,text="chi-sq. density",
line=2)

fig.jpeg
小advice:主贴文字不要用晃眼的颜色

二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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