在R中为时间序列绘图编写循环的技术求助

原代码的核心问题 #

• 循环变量引用错误：'i'是固定字符串，不是循环变量，无法匹配当前迭代的国家
• 缺少按日期统计对象数量的关键步骤：原代码未对每个国家的日期分组计数，直接调用n会报错
• 绘图后未执行打印操作：循环中生成的ggplot对象必须通过print()才能显示
• 标题未动态替换：标题里的i是静态字符串，无法对应当前国家名称

修正后的完整代码 #

# 加载所需依赖包
library(dplyr)
library(tidyr)
library(ggplot2)
library(hrbrthemes) # theme_ipsum依赖此包

# 第一步：全局预处理数据，按国家+日期统计对象数量
df_summary <- df_period5 %>%
  mutate(period5 = as.Date(period5)) %>% # 统一转换日期格式
  group_by(Country, period5) %>%
  summarise(n = n(), .groups = "drop") # 统计每个国家-日期的对象总数

# 第二步：循环处理每个国家并绘图
for (country in unique(df_summary$Country)) {
  # 筛选当前国家的数据
  country_data <- df_summary %>%
    filter(Country == country)
  
  # 补全缺失年份，直接填充NA为0
  country_data_complete <- country_data %>%
    complete(period5 = seq.Date(min(period5), max(period5), by = "year"),
             fill = list(n = 0))
  
  # 生成时间序列图
  plot_obj <- ggplot(country_data_complete, aes(x = period5, y = n)) +
    geom_area(fill = "blue", alpha = 0.5) +
    geom_line(color = "blue") +
    ylab("Number of Objects") +
    theme_ipsum() +
    labs(title = paste("Number of Acquired Objects Over Time in", country))
  
  # 打印图形（循环中必须执行此步骤才会显示图）
  print(plot_obj)
}

关键改进说明 #

• 提前全局预处理：将日期转换、计数统计放在循环外完成，避免重复操作，提升运行效率
• 正确引用循环变量：用Country == country（变量）替代原代码的Country == 'i'（静态字符串），匹配当前迭代的国家
• 简化NA填充：利用complete()的fill参数直接填充缺失值为0，比ifelse更简洁高效
• 动态生成标题：通过paste()拼接国家名称，实现每个图的标题与对应国家匹配
• 强制打印图形：循环中生成的ggplot对象必须调用print()才能在控制台或RStudio中显示

内容的提问来源于stack exchange，提问作者Aby Ramata Almamy Mbaye