将误差条添加到条形图中

时间:2015-10-07 06:13:55

标签: r vector plot bar-chart

我有两个向量。我想制作第一个矢量的条形图(简单,正确)。扭曲是第二个向量的每个元素都是第一个向量的每个元素的标准偏差(它本身是4个其他值的平均值)。我怎么能这样做?

有问题的载体:

-4.6521175 0.145839723
 1.1744100 0.342278694
-0.2581400 0.003776341
-0.3452675 0.073241199
-2.3823650 0.095008502
 0.5625125 0.021627196

即,如何将第二列向量的元素作为误差条添加到第一列向量中的相应元素?

enter image description here

注意:在你问之前,是的,我在网站上进行了大量搜索,并进行了大量的谷歌搜索,但我的问题更具体,即我发现的并不匹配我需要什么。

2 个答案:

答案 0 :(得分:10)

我个人最喜欢arrows()这种图形:

df <- data.frame(bar = c(-4.6521175, 1.1744100, -0.2581400,  -0.3452675, -2.3823650, 0.5625125),
error = c(0.145839723, 0.342278694, 0.003776341, 0.073241199, 0.095008502, 0.021627196))

foo <- barplot(df$bar,ylim=c(-6,2),border=NA)
arrows(x0=foo,y0=df$bar+df$error,y1=df$bar-df$error,angle=90,code=3,length=0.1)

enter image description here

两个细节:

    {li>

    border=NA barplot()删除条形周围的边框,因此您可以在第三个条形图块周围看到错误的胡须。由于第三个误差很小,所以晶须几乎位于条形边缘的顶部。

  1. 我使用length中的arrows()参数来减少水平胡须的宽度,如果我们有更多的条形,这尤其重要。默认值为length=0.25

  2. 但请注意"dynamite plots" have major disadvantages。您写道,您的数据仅来自每个条形的四个原始点。在这种情况下,只需肯定就可以更好地绘制原始数据的(抖动)点图。

答案 1 :(得分:8)

geom_bar的{​​{1}}和geom_errorbar的实施:

ggplot2

这导致:

enter image description here

如果要删除x轴上的数字,可以添加:

library(ggplot2)
ggplot(df, aes(x=row.names(df), y=V1)) +
  geom_bar(stat="identity", fill="grey") +
  geom_errorbar(aes(ymin = V1 - V2, ymax = V1 + V2), width=0.6) +
  theme_classic() 

到你的ggplot代码。

使用过的数据:

  theme(axis.title.x=element_blank(),
        axis.text.x=element_blank(),
        axis.ticks.x=element_blank())

作为对your comment的回应,当您想要绘制如此大量的柱子时,有两种可能的解决方案:

1:仅包含一系列轴标签:

df <- read.table(text="-4.6521175 0.145839723
 1.1744100 0.342278694
-0.2581400 0.003776341
-0.3452675 0.073241199
-2.3823650 0.095008502
 0.5625125 0.021627196", header=FALSE)

这给出了:

enter image description here

可以看出,在情节中填充如此多的条形图是不理想的。因此,请参阅备选方案2。

2:创建一个可用于创建构面的分组变量:

ggplot(df2, aes(x=as.numeric(row.names(df2)), y=V1)) +
  geom_bar(stat="identity", fill="grey", width=0.7) +
  geom_errorbar(aes(ymin = V1 - V2, ymax = V1 + V2), width=0.5) +
  scale_x_continuous(breaks=c(1,seq(10,200,10)), expand=c(0,0)) +
  theme_classic() +
  theme(axis.text.x=element_text(size = 6, angle = 90, vjust = 0.5))

这给出了:

enter image description here

最后两个例子的使用数据:

df2$id <- rep(letters[1:20], each=10)

ggplot(df2, aes(x=as.numeric(row.names(df2)), y=V1)) +
  geom_bar(stat="identity", fill="grey", width=0.7) +
  geom_errorbar(aes(ymin = V1 - V2, ymax = V1 + V2), width=0.5) +
  scale_x_continuous(breaks=as.numeric(row.names(df2))) +
  facet_wrap(~ id, scales = "free_x") +
  theme_bw() +
  theme(axis.text.x=element_text(angle = 90, vjust = 0.5))