R中合并图的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请问一下，怎么样把卡方分布的图加到直方分布图上面啊？ 我要目测该分布是否是自由度为2的卡方分布。 谢谢大家

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：卡方分布 自由度 怎么样 分布图 分布图 自由度

最笨的方法就是用lines(x, dchisq(x, 2))了，x取得足够密即可。

hugebear 发表于 2013-2-3 07:35
最笨的方法就是用lines(x, dchisq(x, 2))了，x取得足够密即可。

你好。 我的问题是，我用hist()做出了一个直方图，然后我要在直方图上加一条自由度为2的卡方分布线，请问如何实现。

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绘出这条曲线，如有疑问请使用帮助文件

参：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)


小advice：主贴文字不要用晃眼的颜色