18

我已经四处寻找解决方案，但似乎大多数处理单独生成的绘图被组合成 PDF 格式，而不是将使用刻面生成的绘图分离到 PDF 的单独页面上。

示例数据

在上述数据中使用以下代码，选择生成列表中的所有项目：

Ex<-read.csv("StackOverflowEx (3).csv")
library(ggplot2)
library(reshape2)
vars <- select.list(names(Ex),multiple=TRUE,graphics=TRUE)
Cases<-subset(Ex,select=vars)

gg<-melt(Cases,id=c("Item","Total","Admin"))
print(ggplot(gg, aes(x=Total,y=Admin))+
  geom_point(colour="dark green",size=1.5)+
  geom_point(aes(y=value,color=variable))+
  geom_smooth(aes(y=value,fill=variable),
              method=loess,size=1,linetype=1,se=T)+
  facet_wrap(~variable,ncol=2,nrow=1000)+
  ylim(0,1)+
  labs(x="Expected",y="Admin",title=vars))

...应该生成所有 8 (AH) Case 的 facet wrap。然而，生成这个往往会限制图表并降低它们的可读性，实际上，我打算在 500 多个案例中使用它（它只返回标有列名且没有可读图表的条形图）。

是否可以在转换为 PDF 时指定分面中显示在单个页面上的图表数量，而不是将所有图表压缩到单个页面中？例如，使用上述数据，在包含所有 8 个案例的单独页面上生成两个 2x2 图（即第 1 页上的案例 AD，第 2 页上的案例 EH）。

我可以通过突出显示 4 个案例 +“项目”、“总计”和“管理员”并为接下来的 4 个案例重复并组合生成的 PDF 来做到这一点。然而，在实践中超过 500 个案例，这意味着超过 100 次迭代，很可能会出现人为错误。一些自动化过程的帮助会很棒。

4

1 回答 1

8

因为还没有答案。这是一个：

library(ggplot2)
library(reshape2)

Ex <- read.csv("C:/Users/Thomas/Desktop/StackOverflowEx (3).csv")
gg <- melt(Ex,  id = c("Item", "Total", "Admin"))

# put in number of plots here
noPlots <- 4
# save all variables in a seperate vector to select in for-loop
allVars <- unique(gg$variable)
noVars <- length(allVars)

# indices for plotting variables
plotSequence <- c(seq(0, noVars-1, by = noPlots), noVars)

# pdf("plotpath.pdf") # uncomment to save the resulting plots in a pdf file
# loop over the variables to plot
for(ii in 2:length(plotSequence)){
  # select start and end of variables to plot
  start <- plotSequence[ii-1] + 1
  end <- plotSequence[ii]

  # subset the variables and save new temporary data.frame
  tmp <- subset(gg, variable %in% allVars[start:end])
  cat(unique(tmp$variable), "\n")

  # generate plot
  p <- ggplot(tmp,  aes(x = Total, y = Admin))+
    geom_point(colour = "dark green", size = 1.5)+
    geom_point(aes(y = value, color = variable))+
    geom_smooth(aes(y = value, fill = variable), 
                method = loess, size = 1, linetype = 1, se = T)+
    facet_wrap(~variable, ncol = 2, nrow = 1000)+
    ylim(0, 1)+
    labs(x = "Expected",  y = "Admin")
  print(p)
}
# dev.off()
于 2015-04-27T20:20:22.653 回答