使用dplyr包生成绘图文件失败

Question

我使用R中的dplyr包来制作地块，并且很难用pdf格式保存地块。

我有一个有几组数据的数据框架，我想要绘制每组数据并以pdf格式保存这些图。我试图在dplyr包中使用do和group_by，但是没有生成pdf文件。没有错误信息。

rm(list = ls())

library(dplyr)

df=data.frame(x = c(rep('a',10),rep('b',10)),y = c(1:20))

plt = function(df)
{
  pdf('plt1.pdf')

  plot(df$x,df$y)

  dev.off()

  return(data.frame())
}

plt(df)

df %>% group_by(x) %>% do(data.frame(plot(.$y)))

我也尝试了/，类似的代码也不起作用。

有人能帮我吗？

Answer 1

您需要在plt中调用do函数；如果要为每组数据生成一个单独的pdf文件，则必须根据数据组相应地修改文件名，否则文件将被最后一组数据中的一个情节覆盖：

plt = function(df) {
    pdf(sprintf('plt_%s.pdf', df$x[1]))
    plot(df$x,df$y)
    dev.off()
    return(data.frame())
}

df %>% group_by(x) %>% do(plt(.))

此外，如果您只想拥有按x分组的y的单个方框图，则可以使用一个boxplot。

boxplot(y ~ x, data = df)

Answer 2

以下是解决您的问题的解决方案：

rm(list = ls())
library(dplyr)
df=data.frame(x = c(rep('a',10),rep('b',10)),y = c(1:20))

plt = function(k, dtset) {
 pdf(paste('plt',k,'.pdf',sep=""))
 plot(dtset$y)
 dev.off()
}
plt(1, df)

grps <- unique(df$x)
for (k in seq_along(grps)) {
 plt(k, subset(df, df$x==grps[k]))
}

Answer 3

如果您想将所有的地块放在一个PDF中，并且不介意使用ggplot，您可以这样做：

library(dplyr)
library(ggplot2)
library(gridExtra)

# i'm using a different test data set to make scatter plots
set.seed(1)
df <- data.frame(id = rep(c("a", "b"), each = 10),
                 var1 = rnorm(20),
                 var2 = rnorm(20),
                 stringsAsFactors = FALSE)

# make a list of plots, one for each group (id)
plotlist <- lapply(unique(df$id), function(x)

  df %>%
    filter(id == x) %>%
    qplot(x = var1, y = var2, data = ., main = x)

)

# use marrangeGrob from gridExtra to get an arranged version of that list
mg = marrangeGrob(plotlist, nrow = 2, ncol = 1)

# now use ggsave to write a pdf with the plots arranged as specified in the
# previous step to the working directory; see documentation for ways to tinker
# with layout, file type, etc.
ggsave("myplots.pdf", mg)