ggplot2:在同一页面中绘制多个直方图,但是一个坐标反转

时间:2013-06-19 09:45:53

标签: r ggplot2 histogram

我想从相同的数据中绘制两个或更多个直方图,但是反转一个组的坐标,如下所示:

enter image description here

(注意:这个数字取自Grossman et al 2011, A composite of Multiple Signals Distinguishes Causal Variants in Regions of Positive Selection

例如,让我们从ggplot2中获取钻石数据集:

> library(ggplot2)
> head(diamonds)
  carat       cut color clarity depth table price    x    y    z
1  0.23     Ideal     E     SI2  61.5    55   326 3.95 3.98 2.43
2  0.21   Premium     E     SI1  59.8    61   326 3.89 3.84 2.31
3  0.23      Good     E     VS1  56.9    65   327 4.05 4.07 2.31
4  0.29   Premium     I     VS2  62.4    58   334 4.20 4.23 2.63
5  0.31      Good     J     SI2  63.3    58   335 4.34 4.35 2.75
6  0.24 Very Good     J    VVS2  62.8    57   336 3.94 3.96 2.48

我尝试过的一种方法是使用stat_bin计算直方图而不绘制它,更改它返回的计数列的值。

> my_sb <-stat_bin(data=diamonds, mapping=aes(x=x)

stat_bin的文档说这个函数应该返回一个等于映射的data.frame,但是添加了四个新列(count,density,ncount,ndensity)。但是,我无法在任何地方找到这些列:

# I supposed that this should contain a count, density columns, but it does not.
> print(head(my_sb$data))  
  carat       cut color clarity depth table price    x    y    z
1  0.23     Ideal     E     SI2  61.5    55   326 3.95 3.98 2.43
2  0.21   Premium     E     SI1  59.8    61   326 3.89 3.84 2.31
3  0.23      Good     E     VS1  56.9    65   327 4.05 4.07 2.31
4  0.29   Premium     I     VS2  62.4    58   334 4.20 4.23 2.63
5  0.31      Good     J     SI2  63.3    58   335 4.34 4.35 2.75
6  0.24 Very Good     J    VVS2  62.8    57   336 3.94 3.96 2.48

另一种可能的方法是使用scale_y_reverse(),但我不知道如何将它应用于单个数据集。

我能想到的第三种方法是使用viewPorts,但我不太清楚如何实现它。

1 个答案:

答案 0 :(得分:4)

喜欢这样:

ggplot(data = diamonds) + 
    geom_histogram(aes(x = x,y = ..count..)) + 
    geom_histogram(aes(x = x,y = -..count..))

仅供参考 - 我不记得过去我是怎么做到的,所以我用Google搜索了“ggplot2倒排直方图”并点击了first hit,一个StackOverflow问题。

我不确定stat_bin返回的proto对象究竟是如何构造的,但新变量就在那里。这种方式的工作方式是geom_histogram本身调用stat_bin来执行分箱,因此它可以访问计算变量,我们可以将其映射到y变量。