[英]Combining boxplots in R
我有一个关于箱线图的绘图问题(使用基本图形)。
我有几个数据阵列,希望将它们转换成箱形图并进行比较。 数组反映了不同的实验,我想展示的是基本结果和实验的百分比差异(在一个图上!)。 即在第一个y轴上的基本结果,在第二个y轴上的%diff:
base <- array(runif(12*24*3), dim=c(12,24,3))
exp1 <- array(runif(12*24*3), dim=c(12,24,3))
exp2 <- array(runif(12*24*3), dim=c(12,24,3))
exp3 <- array(runif(12*24*3), dim=c(12,24,3))
exp4 <- array(runif(12*24*3), dim=c(12,24,3))
# calc p.diff
p.diff <- function(mod,base) {
100.0*((mod-base)/base) }
a <- p.diff(exp1,base)
b <- p.diff(exp2,base)
c <- p.diff(exp3,base)
# combine the % diff arrays
exps <- list(a,b,c)
# plot the results
boxplot(base, xlim=c(1,4), col="gray", xaxt="n", ylab="Base values", outline=FALSE)
axis(side=1, 1:4, labels=c("base","% exp1","% exp2","% exp3") )
par(new=TRUE)
boxplot(exps, col="red", ylim=c(-200,200), outline=FALSE, axes=FALSE)
axis(4)
grid()
这几乎可以用,但是我无法正确定位不同箱形图的位置(如果运行我的示例,您将明白我的意思)。 那么,有没有更好的方法来控制箱形图的放置? 还是产生类似类型人物的更好方法?
编辑(1):您需要定义X轴的刚性序列。 这样地块就不会重叠。 只是尝试玩它。
我认为X轴的标签位置不正确? 我不知道一种更优雅的方式,但这是一个解决方案:
# plot the results
boxplot(base, xlim=c(1,4), col="gray", xaxt="n", ylab="Base values", outline=FALSE)
axis(side=1,1,labels=('base'))
par(new=TRUE)
boxplot(exps, col="red", ylim=c(-200,200), outline=FALSE, axes=FALSE)
axis(4)
axis(side=1,1:3,labels=c("% exp1","% exp2","% exp3"))
grid()
所以我在创建boxplot
之后添加了每个标签。 首先绘制base
并标记,然后绘制exps
并标记。 它能解决您的问题吗?
编辑:为了更清楚,您要添加一个具有3个值的新图,这就是为什么axis(side=1,1:3,labels=c("% exp1","% exp2","% exp3"))
是1到3 ...
编辑(2):
为什么不在图表中使用多行并尝试绘制2个图形? 这是您的数据的示例:
#divide your plottin area into 2 columns with one row.
par(mfrow = c(1, 2))
# plot the results
boxplot(base, col="gray", xaxt="n", ylab="Base values", outline=FALSE,axes=FALSE)
axis(2)
axis(side=1,1,labels=('base'))
segments(0,0,1,0)
boxplot(exps,col="red", xaxt="n", ylim=c(-200,200), outline=FALSE, axes=FALSE)
axis(4)
axis(side=1,at=(1:3),labels=c("% exp1","% exp2","% exp3"))
您可以从这里获得有关它的更多信息
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.