将参考线添加到R中带有ggplot的条形图中

Question

这是一个最小的示例，显示了我要绘制的图。

数据看起来像这样：

plot1 = data.frame(
    Factor1 = as.factor(rep('A', 4)),
    Factor2 = as.factor(rep(c('C', 'D'), 2)),
    Factor3 = as.factor(c( rep('E', 2), rep('F', 2))),
    Y = c(0.225490, 0.121958, 0.218182, 0.269789)
)

plot2 = data.frame(
    Factor1 = as.factor(rep('B', 4)),
    Factor2 = as.factor(rep(c('C', 'D'), 2)),
    Factor3 = as.factor(c( rep('E', 2), rep('F', 2))),
    Y = c(-0.058585, -0.031686, 0.013141, 0.016249)
)

虽然绘图的基本代码如下：

require(ggplot2)
require(grid)

p1 <- ggplot(data=plot1, aes(x=Factor2, y=Y, fill=factor(Factor3))) +
    ggtitle('Type: A') +
    coord_cartesian(ylim = c(-0.10, 0.30)) +
    geom_bar(position=position_dodge(.9), width=0.5, stat='identity') +
    scale_x_discrete(name='Regime',
        labels=c('C', 'D')) +
    scale_y_continuous('Activations') +
    scale_fill_brewer(palette='Dark2', name='Background:',
        breaks=c('E','F'),
        labels=c('E','F')) +
    theme(axis.text=element_text(size=11),
        axis.title.x=element_text(size=13, vjust=-0.75),
        axis.title.y=element_text(size=13, vjust=0.75),
        legend.text=element_blank(),
        legend.title=element_blank(),
        legend.position='none',
        plot.title=element_text(hjust=0.5))

p2 <- ggplot(data=plot2, aes(x=Factor2, y=Y, fill=factor(Factor3))) +
    ggtitle('Type: B') +
    coord_cartesian(ylim = c(-0.10, 0.30)) +
    geom_bar(position=position_dodge(.9), width=0.5, stat='identity') +
    scale_x_discrete(name='Regime',
        labels=c('C', 'D')) +
    scale_y_continuous('Activations') +
    scale_fill_brewer(palette='Dark2', name='Background:',
        breaks=c('E','F'),
        labels=c('E','F')) +
    theme(axis.text=element_text(size=11),
        axis.title.x=element_text(size=13, vjust=-0.75),
        axis.title.y=element_blank(),
        legend.text=element_text(size=11),
        legend.title=element_text(size=13),
        plot.title=element_text(hjust=0.5))

pushViewport(viewport(
    layout=grid.layout(1, 2, heights=unit(4, 'null'),
        widths=unit(c(1,1.17), 'null'))))
print(p1, vp=viewport(layout.pos.row=1, layout.pos.col=1))
print(p2, vp=viewport(layout.pos.row=1, layout.pos.col=2))

这个数字看起来像这样：

但是，我需要这样的东西：

黑色粗线是参考值。它们是恒定的，图中表示该“参考情况”。但是，在其他需要制作条形图的图中，该图将发生变化，但参考值应保持不变，以使比较简单明了。我知道我应该使用geom_segment()，但是在我尝试进行此工作时这些行只是缺失了小节。

有任何帮助/建议吗？谢谢！

Answer 1

我能够使用geom_errorbarh执行此操作。例如，对于第二个数字：

p1 + 
  geom_errorbarh(
    aes(xmin = as.numeric(Factor2)-.2,xmax = as.numeric(Factor2)+.2), #+/-.2 for width
    position = position_dodge(0.9), size = 2, height = 0
  )

输出：而且，如果我了解您描述的其他图，则可以在其中指定参考数据，例如data = plot1

Answer 2

如果您的参考不会改变，您可以创建另一个数据集并将其合并到要绘制的数据集中。

[这里，我首先添加plot1和plot2。然后，我创建一个新的数据集作为参考数据集。

library(dplyr)
new_df = rbind(plot1, plot2)
ref_plot = new_df
ref_plot <- ref_plot %>% rename(Ref_value = Y)

然后，您现在具有要绘制的数据集new_df和包含每个条件的参考值的ref_plot。

[而不是使用grid并创建两个不同的图，然后将其合并，我更喜欢使用facet_wrap将所有图放在同一图形上。这更加方便，不需要写两次相同的东西。

正如@AHart在我面前提到的几分钟，您可以使用geom_errorbar在绘图上定义参考值。区别在于我更喜欢​​使用geom_errorbar而不是geom_errobarh。

这里是剧情：

library(ggplot2)
new_df %>% left_join(ref_plot) %>% 
  ggplot(aes(x = Factor2, y = Y, fill = Factor3))+
  geom_bar(stat = "identity", position = position_dodge())+
  geom_errorbar(aes(ymin = Ref_value-0.00001, ymax = Ref_value+0.0001, group = Factor3), position = position_dodge(.9),width = 0.2)+
  facet_wrap(.~Factor1, labeller = labeller(Factor1 = c(A = "Type A", B = "Type B"))) +
  scale_x_discrete(name='Regime',
                   labels=c('C', 'D')) +
  scale_fill_brewer(palette='Dark2', name='Background:',
                    breaks=c('E','F'),
                    labels=c('E','F')) +
  theme(axis.text=element_text(size=11),
        axis.title.x=element_text(size=13, vjust=-0.75),
        axis.title.y=element_blank(),
        legend.text=element_text(size=11),
        legend.title=element_text(size=13),
        plot.title=element_text(hjust=0.5))