方框和晶须图显示以下信息:最大值、最小、平均值、第75百分位数、第25百分位数。如果我有这些信息,我能画出相应的B&W图吗?
我有一个名为TP.df的数据框架:
pb1 ag1 pb2 ag2 pb3 ag3
Nb 498 498 85 85 68 68
Min 0 0 0 0 0 0
Max 1.72 461 2.641 260.8 0.3 144
Mean 0.06 19.2 0.15 35.35 0.02 9.11
75_p 0.06 20 0.08 33 0.02 8
25_p 0.01 10 0 14 0.01 4
文件:
,pb1,ag1,pb2,ag2,pb3,ag3
Nb,498,498,85,85,68,68
Min,0,0,0,0,0,0
Max,1.72,461,2.641,260.8,0.3,144
Mean,0.06,19.2,0.15,35.35,0.02,9.11
75_p,0.06,20,0.08,33,0.02,8
25_p,0.01,10,0,14,0.01,4
我怎样才能得到相应的盒子和晶须图:
pb1
,ag1
,pb2
,ag2
,pb3
,ag3
0
到max(TP.df[Max,])
发布于 2013-08-02 07:38:15
使用bxp
从一个五个数字的摘要中绘制一个盒子图.
这能做你想做的事吗?
testdata=data.frame(R1=c(0,5,3,2,4),R2=c(1,7,3,2.8,6))
o=c(1,5,3,4,2) # the rows in increasing order
bxp.data=list(stats=data.matrix(testdata[o,]),n=rep(1,ncol(testdata)))
# the n=... parameter doesn't affect the plot, but it still needs to be there
bxp(bxp.data)
您可以使用标准的图形参数(参见help(par)
)来标记轴并使其看起来漂亮。
(注:排号为1,因为问题是在我发布这个答案后编辑的。)我将保留原样的答案,以便与前两个注释相匹配。)
发布于 2013-08-02 07:25:59
这可能是一种解决办法:
#create a dummy boxplot that you can modify the data easily
z<- boxplot(1:10)
#look at the outbut an assign yout data to stats
z$stats<- your_data
#use bxp to plot, via add you can combine all three
bxp(z)
或者从一开始就使用bxp函数。另见:http://r.789695.n4.nabble.com/Box-plot-with-5th-and-95th-percentiles-instead-of-1-5-IQR-problems-implementing-an-existing-solution-td3456123.html
本
发布于 2019-07-11 23:36:07
使用ggplot
使用手动盒绘图,您可以操作绘图美学并添加其他参数:
ggplot()+geom_boxplot(aes(x=1, y = 2:5, lower = 3, upper = 4, middle = 3.5, ymin=2, ymax=5))
https://stackoverflow.com/questions/18010347
复制相似问题