在 ggplot2 中拆分 facet_wraps 多个 .png 输出

Splitting facet_wraps over multiple .png outputs in ggplot2

我正在尝试为数据中的每组生成多个方面图。我面临的问题是,在我的真实数据中,有时每组的方面超过 100,并且 ggplot2 试图将所有方面放入一个图中。这不利于视觉效果。

基本上,我想要实现的是,当每页超过一定数量的 windows 时(假设每页绘制 10 window),我想继续绘制它们在每组的下一张图中。

制作此内容的最小示例仅涵盖每组的绘图构面。

x <- rnorm(25,0,5)
y <- rnorm(25,5,10)

gr=rep(c(1,2),c(20,5))

set <- c(rep(1:20,1),1:5)

test_data <- data.frame(x,y,set,gr)

utils::View(test_data)

对于这个最小的例子,gr 1 的集数是 20。但我想将它分开,也就是说每页 10 个面。

library(ggplot2)
library(gtable)
library(grid)
library(gridExtra)

plot_list=list()

for (i in unique(test_data$gr)){


plot_sets<- ggplot(data = test_data[test_data$gr==i,], aes(y = y, x=x,fill=factor(set))) +


geom_point(alpha = 0.8,size=8,shape=21,col="black",aes(group=set)) +


  facet_wrap(~set)+


  theme_bw()+
  theme(legend.position = "top",plot.title=element_text(face="bold",colour ="Blue",size=25),
        axis.title.x = element_text(face="bold",size=24), 
        axis.title.y = element_text(face="bold",size=24),
        #axis.text.x = element_blank(),
        axis.text.x = element_text(size = 22 , colour = "black",angle=0,vjust = 0.0, hjust = 0),
        axis.text.y = element_text(size = 22 , colour = "black"),
        strip.text = element_text(size=22, lineheight=1.0,face="bold"),
        strip.background = element_rect(fill="#FFFFCC", colour="black", size=0.1),
        legend.title = element_text(colour="black", size=20,face="bold"),
        legend.text = element_text(colour="black", size=26),
        panel.spacing = unit(0.2, "lines"))+
  labs(fill='set',col="set",x = "x",y="y",title=paste(i,j,"my_Sets",sep="_"))+
  guides(fill = guide_legend(override.aes = list(alpha=1,nrow=1,size=5)))  

  plot_list[[i]]=plot_sets

}

以 .png 格式保存图表

for (i in 1:length(unique(test_data$gr))) {

  file_name = paste(unique(test_data$gr), "my_sets.png", sep="_")

  png(paste("C:/test","/",file_name[i],sep=""), width =6400, height =4400, pointsize = 16, units = 'px', res = 200)

  grid.draw(plot_list[[i]])
  dev.off()
}

我找到了一些相关的帖子,例如

Multiple graphs over multiple pages using ggplot

但他们使用 .pdf 输出,这个例子是最相关的,但它适用于 .pdf 输出

@agstudy 解决方案与 .pdf 输出版本

tab <- data.frame(Date = rep(seq(as.Date('2013-03-01'),
                       by = 'day', length.out = ii), nn),
                   value = rep(runif(nn, 100, 200)))

tab$Station <- rep(1:nn, each = ii)

*library(gridExtra)
library(ggplot2)
pdf('test.pdf', width=21, height=27)
i = 1
plot = list() 
for (n in unique(tab$Station)){
  ### process data for plotting here ####
  plot[[i]] = ggplot(tab[tab$Station==n,], aes(x=Date,y = value)) +

    geom_line() +
    facet_wrap(~ Station, ncol = 5)

  if (i %% 20 == 0) { ## print 8 plots on a page
    print (do.call(grid.arrange,  plot))
    plot = list() # reset plot 
    i = 0 # reset index
  }
  i = i + 1
}

if (length(plot) != 0) { 
  print (do.call(grid.arrange,  plot))
}
dev.off()

这是一个使用 cut()ggsave() 的示例。

假设您有 80 套想要刻面包裹。

library(tidyverse)

param_n <- 800

data <- data_frame(
  x = rnorm(param_n, 0, 5),
  y = rnorm(param_n, 5, 10),
  group = factor(sample(c(1, 0), param_n, TRUE)),
  set = sample(1:(param_n / 10), param_n, TRUE)
) 

一个情节太多了,所以你用 cut() 把它分成十页。

data_cut <- 
  data %>%
  mutate(page = cut(set, 10))           # split into 10 groups

然后您只需遍历每一页,filter() 输出数据,并用 ggsave() 保存一个 .png 文件。

for (current_page in levels(data_cut$page)) {

  fig <- 
    data_cut %>%
    filter(page == current_page) %>%
    ggplot(
      aes(y = y, x = x,
          colour = group,
          group = set)) +
    facet_wrap(~ set) +
    geom_point() +
    labs(title = paste("Sets", current_page)) +
    theme_bw() +
    theme(legend.position = "bottom")

  ggsave(
    paste0("figures_sets_", current_page, ".png"),
    plot = fig
  )

}

我简化了你的情节使代码更紧凑。