我使用下面的代码得到了这个图
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"))
在我的原始数据中，我有 30 个级别id。如果我在一页上绘制所有方面，它将变得如此凌乱且难以阅读。因此，我希望有一个类似于上面的情节的情节，但每页上只有 4 个方面，所以它是可读的。
我使用ggplus 包来做到这一点。
#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 手册（第 5 页）和这个问题中获得的。
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）：
无效的'nrow'值（太大或NA）此外：警告消息：1：删除了包含缺失值的6681行（geom_path）。2：删除了包含缺失值的 6681 行 (geom_path)。3：删除了包含缺失值的 6681 行 (geom_path)。4：在 min(layout$ROW) 中：min 没有非缺失参数；返回 Inf 5：在 max(layout$COL) 中：max 没有非缺失参数；返回 -Inf 6：在 max(layout$ROW) 中：max 没有非缺失参数；返回 -Inf 7: In matrix(list(zeroGrob()), nrow = nrow, ncol = ncol) : 由强制引入到整数范围的 NA
有什么建议可以解决这些错误吗？
