创建一个计算正/负ggplot的geom直方图

Question

我无法在正y轴上绘制正log2_ratio计数，在负y轴上绘制负log2_ratio计数。

本质上，我希望正数在x轴之上，而负数在x轴之下。

这是数据框和代码：

   chrom      chr_start   chr_stop     num_positions normal_depth tumor_depth log2_ratio gc_content sample
   324202     1 156249804 156249858            55         12.3         4.7     -1.399       34.5     10
   324203     1 156250463 156250473            11         10.0         4.6     -1.109       27.3     10
   324204     1 156250664 156250705            42         12.0         7.4     -0.704       19.0     10
   324205     1 156250816 156250847            32         11.7         4.6     -1.343       40.6     10
   324206     1 156251108 156251132            25         10.6         3.6     -1.569       60.0     10
   324207     1 156251411 156251464            54         12.3         6.8     -0.863       46.3     10

newHist = ggplot(resultsPileup1COMBINED[resultsPileup1COMBINED$sample <= 25,],
          aes(x=sample)) +
          geom_histogram(fill="blue" , bindwidth = 1) +
          geom_histogram(data=resultsPileup1COMBINED[resultsPileup1COMBINED$sample > 25,],
            fill="gray50" , binwidth = 1) +
            scale_x_continuous(breaks = seq(from = 1, to = 50, by = 3))

这是当前的图表：

Answer 1

如果您要求全新的图表，请尝试：

-..count..

在这个情节中，我们有＆lt;底部为0 log2_ratios，x值为负，与上面对齐，使用 - ..count ..方法

编辑：要求稍微不同的图表，留下上面的后代。

要绘制每个bin的+/-值的数量，我们再次使用ggplot() + geom_histogram(data = resultsPileup1COMBINED[resultsPileup1COMBINED$log2_ratio < 0, ], aes(x = sample, y = -..count..)) + geom_histogram(data = resultsPileup1COMBINED[resultsPileup1COMBINED$log2_ratio > 0, ], aes(x = sample, y = ..count..)) 技巧绘制出来：

scale_x_discrete(breaks = seq(from = 1, to = 50, by = 3))

再一次，休息和颜色取决于你。

要使其与原始图表类似，请确保包括：

close()