ggpubr::ggarrange处理未知数量子图列表失效问题求助

需要通过双层循环，借助ggpubr包的ggarrange函数生成绘图：外层按places分组，每个分组内按dates生成对应散点图，再将同一places下的子图排列。但由于每个places对应的dates数量不固定，尝试用Figure1[[1]]:Figure1[[length(Figure1)]]的语法传递子图列表时触发报错：

Error in Figure1[[1]]:Figure1[[length(Figure1)]] : NA/NaN argument
In addition: Warning messages:
1: In Figure1[[1]]:Figure1[[length(Figure1)]] :
  numerical expression has 9 elements: only the first used
2: In Figure1[[1]]:Figure1[[length(Figure1)]] :
  numerical expression has 9 elements: only the first used

1. 列表元素提取语法错误：:是用于生成数值序列的运算符，不能用来提取列表的连续元素。Figure1[[1]]:Figure1[[length(Figure1)]]会把两个ggplot对象当成数值处理，直接引发报错。
2. 数据筛选变量名错误：内层循环中筛选日期时误用了my.ds.sub$date，实际数据列名为dates。
3. 未定义变量：ggsave中的all_stands[g]是未定义的变量，与外层循环的变量不匹配。
4. ggplot语法缺失：theme_bw()前缺少+号，导致绘图语法不完整。

关键修改点 #

• 直接将子图列表Figure1通过plotlist参数传给ggarrange，自动适配不同长度的子图集合。
• 修正数据筛选和变量名错误，补全ggplot语法。
• 用lapply批量统一修改子图样式，避免重复代码。

完整修正代码 #

# 加载依赖包
library(ggplot2)
library(ggpubr)

# 生成模拟数据
places <-c("A","A","A","B","B","B","B","A","A","A","B","B","B","B","A","A","A","B","B","B","B","A","A","A","B","B","B","B")
dates <-c("2019-01-01","2019-02-01","2019-03-01","2019-01-01","2019-02-01","2019-03-01","2019-04-01",
          "2019-01-01","2019-02-01","2019-03-01","2019-01-01","2019-02-01","2019-03-01","2019-04-01",
          "2019-01-01","2019-02-01","2019-03-01","2019-01-01","2019-02-01","2019-03-01","2019-04-01",
          "2019-01-01","2019-02-01","2019-03-01","2019-01-01","2019-02-01","2019-03-01","2019-04-01")
areas <-rnorm(56,mean=100,sd=10)
my.ds <- data.frame(places,dates,areas)

# 按places分组生成并排列子图
all_places <- unique(my.ds$places)

for (a in 1:length(all_places)) {
  Figure1 <- list()
  my.ds.sub <- my.ds[my.ds$places == all_places[a],]
  unique_dates <- unique(my.ds.sub$dates)
  
  # 内层循环生成每个日期的散点图
  for (d in 1:length(unique_dates)) {
    my.ds.sub.date <- my.ds.sub[my.ds.sub$dates == unique_dates[d],]
    Figure1[[d]] <- ggplot(data = my.ds.sub.date, aes(x = dates, y = areas)) +
      geom_point() +
      theme_bw()
  }
  
  # 统一移除所有子图的x轴文本
  Figure1 <- lapply(Figure1, function(p) p + rremove("x.text"))
  
  # 排列子图
  combined_plot <- ggarrange(plotlist = Figure1,
                             font.label = list(size = 18, color = "black", face = "plain", family = "serif"),
                             widths = rep(1, length(Figure1)),
                             heights = rep(1, length(Figure1)),
                             ncol = ceiling(sqrt(length(Figure1))), 
                             nrow = ceiling(sqrt(length(Figure1))),
                             common.legend = TRUE, 
                             legend = "bottom")
  
  # 保存图片
  ggsave(paste0(all_places[a], ".svg"), plot = combined_plot, dpi = 600, width = 15, height = 15, units = "in")
}

代码说明 #

• 通过plotlist = Figure1直接传递子图列表，无需手动指定元素范围，完美适配不同数量的子图。
• 用lapply批量修改子图样式，统一移除x轴文本，简化重复操作。
• 修正了所有语法和变量错误，确保代码可正常运行。

内容的提问来源于stack exchange，提问作者Leprechault