已有变量的置信区间数据，用什么绘图函数画置信区间阴影图？

5226

收藏 2016-09-14

0.1_pre	0.5_pre	0.9_pre
10.03153	10.377	10.72248
10.00079	10.35178	10.70288
10.04735	10.39759	10.74797
10.03385	10.3793	10.72478
10.0155	10.36892	10.72225
10.05432	10.40815	10.76136
10.05113	10.39674	10.74235
10.08417	10.43517	10.78625
10.04867	10.39353	10.73839
10.03898	10.3868	10.73451
10.06012	10.40519	10.75047
10.04588	10.39396	10.74205
10.00231	10.35539	10.70855
10.00163	10.35535	10.70915
10.00489	10.35688	10.70893
10.03977	10.38506	10.73037
10.0033	10.35598	10.70871
10.01921	10.36751	10.71587
10.05194	10.39564	10.73928
10.06064	10.40316	10.74564
10.09896	10.45761	10.81636
10.06952	10.41527	10.76105
10.07373	10.44153	10.80911
10.11984	10.50313	10.88585
10.08135	10.42877	10.77609
10.04901	10.39541	10.7418
10.07467	10.42503	10.77525
10.10553	10.47747	10.84912
10.08382	10.4503	10.816
10.08841	10.4423	10.79611
10.06373	10.40829	10.75289
10.08229	10.42926	10.77631
10.06692	10.41033	10.75372
10.06184	10.41021	10.75864
10.06074	10.41524	10.77027

数据如上所示，第一列是0.1置信区间，第二列是均值，第三列是0.9置信区间，求问用什么绘图函数可以把这三列数据绘制在一副图中，并且用阴影来表示置信区间。
期待众位大神解答。

扫码加我拉你入群

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

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

全部回复

luzifer_li

2016-9-14 17:22:33

试试这个吧，自己写了个函数。只弄了前五行。

#read.csv('1.csv',head=T,stringsAsFactors=F)->inter1
inter1是你的data.frame，三列

library(grid)
grid.newpage()
pushViewport(viewport(h=0.8,w=0.8))
pushViewport(dataViewport(xData=c(0,1),yData=range(as.matrix(inter1)),extension=c(0,0.08)))
ppoints(nrow(inter1),1/2)->interx
inter.wid<-0.75/nrow(inter1)
grid.rect(x=interx, y=inter1[,1],width=inter.wid,just='bottom',height=inter1[,3]-inter1[,1],gp=gpar(fill='gray',col='gray'),default='native')
grid.segments(x0=interx-inter.wid*0.5, x1=interx+0.5*inter.wid,y0=inter1[,2],y1=inter1[,2],gp=gpar(col='red',lwd=1.5),default='native')
grid.rect()
grid.yaxis()

附件列表

try.png

原图尺寸 3.42 KB

扫码加我拉你入群

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

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

慕清雪

2016-9-14 17:40:54

luzifer_li 发表于 2016-9-14 17:22
试试这个吧，自己写了个函数。只弄了前五行。

#read.csv('1.csv',head=T,stringsAsFactors=F)->inter1

谢谢大神，图中您使用的是条形图，请问如果是使用折线图，应该修改什么参数呢？

扫码加我拉你入群

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

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

luzifer_li

2016-9-18 09:37:07

<code>
grid.newpage()
pushViewport(viewport(h=0.8,w=0.8))
pushViewport(dataViewport(xData=c(0,1),yData=range(as.matrix(inter1)),extension=c(0,0.08)))
ppoints(nrow(inter1),1/2)->interx
inter.wid<-0.75/nrow(inter1)
grid.rect()
grid.rect(x=interx, y=inter1[,1],width=inter.wid,just='bottom',height=inter1[,3]-inter1[,1],gp=gpar(fill='gray',col='gray'),default='native')
#grid.segments(x0=interx-inter.wid*0.5, x1=interx+0.5*inter.wid,y0=inter1[,2],y1=inter1[,2],gp=gpar(col='red',lwd=1.5),default='native')
grid.points(interx,inter1[,2],pch=20,default='native',gp=gpar(col='red'))
grid.lines(interx,inter1[,2],default='native',gp=gpar(col='red',lwd=2))
grid.yaxis()
</code>

附件列表

try.png

原图尺寸 7.21 KB

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群