1

我有一个大数据框，我正在使用 ggplot。

ggplot()+
  geom_line(data=DATA,aes(logl,PercPos,group=name),col="blue")+
  geom_line(data=DATA,aes(logl,PercNeg,group=name),col="green")+ 
  geom_histogram(data=DATA, aes(x=logl,alpha=.5,group=name),bindwidth=0.1,density=T)+
  facet_wrap(~ name,as.table=T)+
  theme_bw()+
  xlim(min(DATA$logl),max(DATA$logl))

它生成 180 个不同的情节。我希望将结果绘制在不同的文件中，例如每个文件中不超过 30 个绘图。

有没有一种简单的自动方法可以做到这一点？

感谢帮助。

4

1 回答 1

2

这是一种方法：

library(ggplot2)
library(plyr)

# dummy data with 181 different levels
N <- 181
d <- data.frame(x=rnorm(10*N), y=rnorm(10*N), f=gl(N, 10))

# group levels to fit ppp plots per page
split_df_ppp <- function(d, f, id.var = ".page_split", ppp=30){

  stopifnot(is.factor(f))
  n <- length(levels(f))
  pages <- n%/%ppp + as.logical(n%%ppp)
  groups <- split(levels(f), gl(pages, ppp, n))

  d[[id.var]] <- f
  levels(d[[id.var]]) <- groups

  invisible(d)

}

d2 <- split_df_ppp(d, d$f)
str(d2)

# ggsave() could go in this function, if multiple files are preferred
plot_one_page <- function(.d){
  qplot(x, y, data=.d) +facet_wrap(~f)
}

# here saving all plots as separate pages in one file
pdf("multipage.pdf")
# loop through pages
 d_ply(d2, ".page_split", plot_one_page, .print=TRUE)
dev.off()

您也可以使用marrangeGrob来自 gridExtra，在这种情况下不需要刻面：

lp = dlply(d, "f", function(d) qplot(x, y, data=d))
g = do.call(gridExtra::marrangeGrob, c(lp, ncol=5, nrow=6))

ggsave("multipage.pdf", g)
于 2013-05-29T12:39:33.150 回答