0

我正在尝试gridGraphics::grid.echo根据此链接使用,因此我可以将Gviz plotTracks绘图与ggplotusing cowplot's结合起来plot_grid

以下Gviz小插图这个链接,这就是我所做的:

require(Gviz)
data(geneModels)
gtrack <- GenomeAxisTrack()
itrack <- IdeogramTrack(genome = "hg19", chromosome = as.character(geneModels$chromosome[1]))
grtrack <- GeneRegionTrack(geneModels, genome = "hg19",chromosome = as.character(geneModels$chromosome[1]), name = "Gene Model")

require(gridGraphics)

gwrap_plot <- function(x) {
  tree <- grid::grid.grabExpr(gridGraphics::grid.echo(x))
  u <- grid::unit(1, 'null')
  gtable::gtable_col(NULL, list(tree), u, u)
}

graphics.off()
plotTracks(list(itrack, gtrack, grtrack))
track.plot <- recordPlot()
gwrap_plot(track.plot)

我得到这个错误:

 Error in switch(x[[2]][[1]]$name, C_abline = C_abline(x[[2]]), C_plot_new = C_plot_new(x[[2]]),  : 
  EXPR must be a length 1 vector 

知道gridGraphics's grid.echowith Gviz' 的plotTracks情节有什么问题吗?

4

1 回答 1

1

grid.echo用于基本图形,Gviz 似乎使用grid图形,

p1 = grid::grid.grabExpr(plotTracks(list(itrack, gtrack, grtrack), add = TRUE))
p2 = ggplot2::qplot(1:10, 1:10)

gridExtra::grid.arrange(p1, p2, ncol=2)

在此处输入图像描述

于 2017-01-02T08:24:05.523 回答