Question

我使用下面的代码获得了这个情节

library(tidyverse)
ggplot(df2, aes(x =Date, y = Sales, color = id))+ 
  geom_line(size = 0.01, alpha = 0.3)+
  +facet_wrap(~id)
  scale_x_date(breaks = seq(as.Date("2001-01-01"),
                            as.Date("2007-01-01"), by="1 year"),
               labels = date_format("%Y"))

在我的原始数据中，我有id的30个级别。如果我在一个页面上绘制所有方面，它将是如此混乱和难以阅读。因此，我希望有一个类似于上图的情节，但每个页面只有4个面，因此它是可读的。

我用ggplus package这样做了。

#devtools::install_github("guiastrennec/ggplus")
library(ggplus)
p <- ggplot(df2, aes(x =Date, y = Sales))+ 
  geom_line(size = 0.01, alpha = 0.3)+
  scale_x_date(breaks = seq(as.Date("2001-01-01"),
                            as.Date("2007-01-01"), by="1 year"),
               labels = date_format("%Y"))

facet_multiple(plot = p, 
               facets = "id", 
               ncol = 2, 
               nrow = 2)

如果没有NA，它工作正常，我得到两页的情节。但是，由于我的数据有NA，我收到了这个错误

if（zero_range（range））{：缺失值为TRUE / FALSE时出错需要

我将非常感谢您解决此错误的任何建议。

数据

Date <- as.Date(c(seq(as.Date("2001-01-03"), as.Date("2006-10-17"), by = 1), 
                  seq(as.Date("2001-05-10"), as.Date("2006-12-17"), by = 1),  
                  seq(as.Date("2001-04-12"), as.Date("2006-11-17"), by = 1),
                  seq(as.Date("2001-03-12"), as.Date("2006-10-12"), by = 1),
                  seq(as.Date("2001-01-12"), as.Date("2006-11-03"), by = 1),
                  seq(as.Date("2001-03-11"), as.Date("2006-10-21"), by = 1),
                  seq(as.Date("2001-02-14"), as.Date("2006-11-25"), by = 1),
                  seq(as.Date("2001-04-22"), as.Date("2006-12-27"), by = 1)))
id  <- c(rep("AAA", 2114), rep("BBB", 2048), rep("CCC", 2046), rep ("DDD", 2041), rep ("EEE", 2122), rep ("FFF", 2051), rep ("GGG", 2111) , rep ("HHH", 2076))
Sales <- c(sample(10:20, 2114, replace = T), sample(50:60, 2048, replace = T), sample(80:90, 2046, replace = T), sample(80:90, 2041, replace = T),
           sample(12:70, 2122, replace = T), sample(30:90, 2051, replace = T), sample(15:70, 2111, replace = T), sample(100:120, 2076, replace = T))
df <- data.frame(id , Date, Sales)

df1 <- df[c(1:50, 1050:2000, 3000:3600, 4000:4350, 6000:6400, 9000:9700, 10456:11254, 12000:12850, 14500:15200),] %>% 
  tidyr::spread(id, Sales)

df2 <- data.frame(Date = seq(as.Date("2001-01-01"), as.Date("2006-12-31"), by = 1)) %>% 
  dplyr::left_join(., df1, by ="Date") %>% 
  tidyr::gather("id", "Sales", 2:9)

更新

我尝试使用ggforce::facet_wrap_paginate在许多页面上绘制构面。

以下部分代码是从ggforce manual（第5页）和this question获得的。

n_pages_needed_df2 <- ceiling(
  length(levels(df2$id)) * length(levels(df2$id)) / 4
)


for (i in seq_len(n_pages_needed_df2)) {
  pl<-ggplot(df2, aes(x =Date, y = Sales , color = id))+ 
    geom_line(size = 0.2, alpha = 0.3)+
    scale_x_date(breaks = seq(as.Date("2001-01-01"),
                              as.Date("2007-01-01"), by="1 year"),
                 labels = date_format("%Y"))+
    facet_wrap_paginate(~id, ncol = 2, nrow = 2, page = i,
                        strip.position="top", scales="free_y")
  ggsave(paste("fig-", i, ".png", sep=""), height = 5.8, width = 10.04 , dpi = 600)
}

它工作正常，我得到了以下两个.png图（我更喜欢将其保存为.png，因为将我更容易/更快地导入到Microsoft Office中）。

第一个情节

第二个情节

但是，它导致了第三个空图和这些错误

矩阵错误（list（zeroGrob（）），nrow = nrow，ncol = ncol）：
无效＆＃39; nrow＆＃39;值（太大或NA）另外：警告消息： 1：删除了包含缺失值的6681行（geom_path）。 2：删除了包含缺失值（geom_path）的6681行。 3：删除了 6681行包含缺失值（geom_path）。 4进 min（layout $ ROW）：min的非缺失参数;返回Inf 5：In max（layout $ COL）：max的非缺失参数;返回--Inf 6：在max（layout $ ROW）中：max没有非缺失参数;返回-Inf 7：在矩阵中（list（zeroGrob（）），nrow = nrow，ncol = ncol）：NA 通过强制引入整数范围

我是如何修复这些错误的？

Answer 1

删除NAs效果很好：

library(ggplus)
df2 <- df2[!is.na(df2$Sales),]
pdf("C:\\1\\test.pdf", 7, 5)
p <- ggplot(df2, aes(x =Date, y = Sales)) +
  geom_line(aes(colour=id),size = 0.01)+
  scale_x_date(breaks = seq(as.Date("2001-01-01"),
                            as.Date("2007-01-01"), by="1 year"),
               labels = date_format("%Y"))
facet_multiple(plot = p, facets = 'id', ncol = 2, nrow = 2)
dev.off()

PS：另一种选择是使用gridExtra包。

使用ggforce / ggplus

1 个答案: