Разделение facet_wraps на несколько выходов .png в ggplot2 - PullRequest
Новая
       25

Разделение facet_wraps на несколько выходов .png в ggplot2

0 голосов
/ 12 мая 2018

Я пытаюсь создать несколько графиков фасетов на группу в данных. Проблема, с которой я сталкиваюсь, заключается в том, что в моих реальных данных иногда число граней на группу превышает 100, и ggplot2 пытается поместить все аспекты в один график. Это не хорошо для визуальности.

По сути, я хочу добиться того, чтобы при превышении определенного числа окон фасетов на странице (скажем, рисовать 10 окон на странице) я хотел бы продолжать отображать их на следующих графиках для каждой группы.

Минимальный пример того, как это охватывает только грани для каждой группы. 

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()
}

Я нашел несколько соответствующих сообщений, таких как

Несколько графиков на нескольких страницах с использованием ggplot

ggplot2: графики на нескольких страницах

но они используют вывод в формате .pdf, и этот пример является наиболее актуальным, но он работает с выводом .pdf

@ гибкое решение с выходной версией .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()

1 Ответ

0 голосов
/ 14 мая 2018

Вот пример, который использует 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() извлекаете эти данные и сохраняете файл .png с помощью ggsave().

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
  )

}

Example output

Я упростил ваш сюжет, чтобы сделать код более компактным.

...