在ggplot2循环中基于分类变量Colonizers修改线条颜色

原代码核心问题 #

• 颜色映射错误：将Colonizers列用单引号包裹，导致ggplot把它当成固定字符串，而非数据列引用
• 重复定义美学映射：主ggplot()中已声明aes，geom_line()无需重复指定x、y、colour
• 日期轴处理错误：period5是Date类型，却使用scale_x_continuous()，且日期序列生成方式错误
• 刻度设置不合理：y轴刻度seq(0,5000,10)会生成500个刻度，导致图表拥挤

修正后的代码 #

library(tidyverse)

# 循环生成每个国家的折线图
for (i in country_count$Country) {
  plot_data <- df_period5 %>%
    filter(Country == i) %>%
    complete(period5 = seq.Date(min(period5), max(period5), by = "year")) %>%
    mutate(n = ifelse(is.na(n), 0, n))
  
  print(
    ggplot(plot_data, aes(x = period5, y = n, colour = Colonizers)) +
      geom_line(linewidth = 1) +  # 无需重复aes，直接继承主映射
      ylab("Number of Objects") +
      ggtitle(str_to_title(i)) +  # 标题首字母大写更美观
      theme(text = element_text(size = 16, family = "serif")) +
      scale_x_date(
        breaks = seq.Date(as.Date("1900-01-01"), as.Date("2019-01-01"), by = "10 years"),
        date_labels = "%Y"  # 只显示年份，让轴标签更清晰
      ) +
      scale_y_continuous(breaks = seq(0, 5000, 500))  # 调整y轴刻度间隔，避免拥挤
  )
}

关键修复说明 #

1. 正确映射颜色变量：移除Colonizers的单引号，让ggplot识别为数据框中的分类列，自动按殖民者分组着色
2. 简化美学映射：geom_line()直接继承ggplot()中的aes设置，避免冗余代码
3. 日期轴正确配置：使用scale_x_date()处理日期类型轴，用seq.Date()生成正确的日期序列，搭配date_labels优化标签显示
4. 优化刻度设置：调大y轴刻度间隔，避免因刻度过多导致的图表混乱

测试数据结构 #

df_period5 <- structure(list(period5 = structure(c(-11688, -11688, -7670, -14245, 
-7305, -13149, 7305, -2557, -10592, -14610), tzone = "America/Los_Angeles", class = "Date"), 
    Country = c("cameroun", "tchad", "tchad", "mali", "tchad", 
    "cameroun", "madagascar", "madagascar", "mali", "mali"), 
    n = c(4135L, 2725L, 2479L, 1803L, 1674L, 1264L, 1088L, 1032L, 
    863L, 761L), Independence.x = c(1960, 1960, 1960, 1960, 1960, 
    1960, 1960, 1960, 1960, 1960), Colonizers.x = c("France", 
    "France", "France", "France", "France", "France", "France", 
    "France", "France", "France"), Independence.y = c(1960, 1960, 
    1960, 1960, 1960, 1960, 1960, 1960, 1960, 1960), Colonizers.y = c("France", 
    "France", "France", "France", "France", "France", "France", 
    "France", "France", "France"), Independence.x.x = c(1960, 
    1960, 1960, 1960, 1960, 1960, 1960, 1960, 1960, 1960), Colonizers.x.x = c("France", 
    "France", "France", "France", "France", "France", "France", 
    "France", "France", "France"), Independence.y.y = c(1960, 
    1960, 1960, 1960, 1960, 1960, 1960, 1960, 1960, 1960), Colonizers.y.y = c("France", 
    "France", "France", "France", "France", "France", "France", 
    "France", "France", "France"), Independence = c(1960, 1960, 
    1960, 1960, 1960, 1960, 1960, 1960, 1960, 1960), Colonizers = c("France", 
    "France", "France", "France", "France", "France", "France", 
    "France", "France", "France")), row.names = c(NA, -10L), class = c("tbl_df", 
"tbl", "data.frame"))

# 假设country_count是包含唯一国家的数据集
country_count <- df_period5 %>% distinct(Country)

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