如何使用ggplot2绘制重叠范围

Question

我试图理解ggplot2。特别是，我试图找出是否有更好（更优雅，更简单）的方法来创建Bioconductor IRanges包装插图（找到here中找到的情节，如图所示第12页，代码在第11页）。

在插图中，使用以下代码生成绘图：

plotRanges <- function(x, xlim = x, main = deparse(substitute(x)),
+ col = "black", sep = 0.5, ...) +{
+ height <- 1
+   if (is(xlim, "Ranges"))
+     xlim <- c(min(start(xlim)), max(end(xlim)))
+   bins <- disjointBins(IRanges(start(x), end(x) + 1))
+ plot.new()
+   plot.window(xlim, c(0, max(bins)*(height + sep)))
+   ybottom <- bins * (sep + height) - height
+   rect(start(x)-0.5, ybottom, end(x)+0.5, ybottom + height, col = col, ...)
+   title(main)
+ axis(1) +}

ir <- IRanges(c(1, 8, 14, 15, 19, 34, 40),
+   width = c(12, 6, 6, 15, 6, 2, 7))

plotRanges(ir)

堆叠条形图是通过绘制矩形创建的，而且必须计算每个矩形的角点，高度和宽度，这让我觉得不是很优雅，ggplot2有更优雅的方式来做这个？我知道'优雅'不是一个非常精确的描述，但我希望你明白我的意思（如果不是我会尝试更好地解释）。

Answer 1

这是使用ggplot2生成类似情节的方法。我使用IRanges的示例数据。

library(IRanges)   

# example data
ir <- IRanges(c(1, 8, 14, 15, 19, 34, 40),
              width = c(12, 6, 6, 15, 6, 2, 7))
# IRanges of length 7
#     start end width
# [1]     1  12    12
# [2]     8  13     6
# [3]    14  19     6
# [4]    15  29    15
# [5]    19  24     6
# [6]    34  35     2
# [7]    40  46     7    

bins <- disjointBins(IRanges(start(ir), end(ir) + 1))
# [1] 1 2 1 2 3 1 1

dat <- cbind(as.data.frame(ir), bin = bins)

library(ggplot2)
ggplot(dat) + 
         geom_rect(aes(xmin = start, xmax = end,
                       ymin = bin, ymax = bin + 0.9)) +
  theme_bw()