ggplot，在具有不同行的多个页面中绘制

Question

我想使用 facet_wrap 制作图表并将其绘制在 pdf 文件的不同页面中。我读过儿子很多选择，这很有效： R + ggplot: plotting over multiple pages 但只有当你在每一页中有相同的行时。

我有这个演示数据来尝试解释我的案例：

A <- data.frame(TIME = rep(c(0, 5, 10, 15, 30, 45, 60), 5))
A$C <- (1 - exp(-0.2*A$TIME))
A$ID <- rep(1:5, each = 7)
A$R <- rnorm(35, mean = 1, sd = 0.01)
A$C2 <- A$C*A$R
Pages <- 5
A2 <- A[c(1,4:8,10:22,24:35),]

因此，我的 ID 具有不同的观察次数。我试图用每个 ID 中的观察次数制作一个向量（我想要每页一个 ID），但它不起作用。

nrws <- ddply(A2, .(ID), "nrow")
nsamp <- nrws[,2]

pdf("Test.pdf")
for (i in seq(Pages))
  {
    slice = seq(((i-1)*nsamp[i]),(i*nsamp[i]))
    slice2 = slice[!(slice > nrow(A2))]
    A3 = A2[slice2,]
    p1 <- ggplot(A3, aes(x = TIME, y = C2)) +  
      geom_line(size = 0.5) +
      geom_point(size = 1) +
      facet_wrap(~ID)
print(p1)
  }
dev.off()

你能帮帮我吗？

提前致谢，纳乔

Answer 1

我认为您在计算 "slices" 时想得太多了。也许你想要这个？

不太确定。如果您只希望每页有一个 ID，则不需要 facet_wrap，并且您可能需要明确设置比例以使其在每一页之间保持相同。

library(plyr)

A <- data.frame(TIME = rep(c(0, 5, 10, 15, 30, 45, 60), 5))
A$C <- (1 - exp(-0.2*A$TIME))
A$ID <- rep(1:5, each = 7)
A$R <- rnorm(35, mean = 1, sd = 0.01)
A$C2 <- A$C*A$R
Pages <- 5
A2 <- A[c(1,4:8,10:22,24:35),]

nrws <- ddply(A2, .(ID), "nrow")
nsamp <- nrws[,2]

pdf("Test.pdf")
for (i in seq(Pages))
{
#  slice = seq(((i-1)*nsamp[i]),(i*nsamp[i]))
#  slice2 = slice[!(slice > nrow(A2))]
#  A3 = A2[slice2,]

  A3 = A2[A2$ID==i,]
  p1 <- ggplot(A3, aes(x = TIME, y = C2)) +  
    geom_line(size = 0.5) +
    geom_point(size = 1) +
    facet_wrap(~ID)
  print(p1)
}
dev.off()

ggplot，在具有不同行的多个页面中绘制

ggplot, plotting in multiple pages with different rows

pdf

loops

r

ggplot2

facet-wrap