将图添加到列表ggplot2

时间:2018-12-03 17:31:38

标签: r ggplot2 ggpubr

我正在尝试将ggplot2地块添加到列表中,以便可以使用ggpubr的ggarrange函数将它们组织在多个页面上。

由于有数百个图,因此我正在使用一个函数来生成和保存图,但是我无法获得将图返回到环境或将名称写入列表的功能。

我敢肯定,这很简单,我很想念但找不到它。

我正在使用的绘图功能是:

histFacet.plot <- function(x, results, info, ...) {

  md<- names(x) %in% c("rn","Taxa","year","rep","block","column",
                       "range", "entity_id")
  traits <- names(x[ , !md])
  for (i in traits) {
    i <-ggplot(data = x, aes_string(x = i)) + 
      geom_histogram(colour="black", fill="white") + 
      #facet_grid(x$year ~ .) +
      theme_bw() +
      xlab(paste0(i)) +
      ylab("Frequency") +
      theme(panel.grid.major = element_blank()) +
      theme(panel.grid.minor = element_blank()) +
      theme(axis.text = element_text(size = 15)) +
      theme(axis.title = element_text(size = 15)) +
      theme(strip.text = element_text(size = 15)) 
    #ggsave(paste0(i,"_",info,".pdf"),path=paste(results, sep=''))
    plotList<- list(plotList, list(i))
    print(i)

  }
  return(i)
}

histFacet.plot(pd,'~/Dropbox/Research_Poland_Lab/AM Panel/Figures/Hist/',
               "_raw_2018")

1 个答案:

答案 0 :(得分:0)

您的大问题是您return(i)而不是return(plotList)。在迭代器中使用ifor循环内重新分配i很奇怪。尤其是当情节使用i作为字符串时,我会尝试这样做:

histFacet.plot <- function(x, results, info, ...) {

  md <- names(x) %in% c("rn","Taxa","year","rep","block","column",
                       "range", "entity_id")
  traits <- names(x[ , !md])
    plotList = list()
    for (i in traits) {
        thisPlot <- ggplot(data = x, aes_string(x = i)) + 
          geom_histogram(colour="black", fill="white") + 
          #facet_grid(x$year ~ .) +
          theme_bw() +
          xlab(i) +
          ylab("Frequency") +
          theme(panel.grid.major = element_blank()) +
          theme(panel.grid.minor = element_blank()) +
          theme(axis.text = element_text(size = 15)) +
          theme(axis.title = element_text(size = 15)) +
          theme(strip.text = element_text(size = 15)) 

        plotList[[i]] = thisPlot
        print(i)

      }
      return(plotList)
}

由于您没有共享示例数据,因此未经测试,当然。如果您提供了一个小的,可重现的示例数据集,则很乐意进行测试/调试。

我不确定是否要print(i)将绘图打印到图形设备上(如果是,请更改为print(thisPlot))还是将当前特征打印到控制台上以更新进度循环(如果是这样,请更改为message(i)以使其易于禁用)。

其他一些注意事项:如果您是多方面的,请使用year ~ .作为公式,而不要使用x$year ~ .。如果i已经成为书迷,则paste0(i)i(在您的xlab中)相同。