Ручная установка цветов для графика Воронио при раскраске по фактору с использованием ggvoronoi

Question

утром, днем ​​или вечером.

# Reproducible data
df <- quakes[1:20, 1:2]
df$years <-  as.factor(rep(c("2000","2020"), each=10))
df$cluster <- as.factor(c("1","1","1","1","1","1","2","2","2","2",
                          "2","2","2","2","2","3","3","3","3","3"))

Я использую данные GPS для создания графиков вороной и раскраски их по коэффициенту (вывод кластеризации k-средних). Мне нужно создать довольно много графиков, поэтому я запускаю их в al oop следующим образом:

years <- levels(df$years)

library(dplyr)
library(ggplot2)
library(ggvoronoi)

for(i in years){
  #
  single_year <- df %>% 
    filter(years == i)
  #
  #
  plot <- ggplot(single_year,
                 aes(x=lat,
                     y=long)) +
    #
    geom_voronoi(aes(fill=(cluster))) +
    #
    stat_voronoi(geom="path" )+
    #
    geom_point() +
    #
    labs(title = paste(i))
  #
  #
  ggsave(paste0(i,".jpeg"), plot = last_plot(), # Watch out for the SAVE!!!
         device = 'jpeg')
  #
}

Это дает мне следующие (удивительные) графики:

Эта проблема в цвете. Я хочу иметь согласованность между сюжетами. Например, для любого графика кластер 2 будет синим, кластер 3 = красным и т. Д. c.

Меня смущает, какой из множества вариантов цвета ggplot использовать здесь для обеспечения согласованности. Большое спасибо!!

Answer 1

Вы можете определить вектор для присвоения цвета каждому значению переменной «cluster», а затем передать их в аргумент values = функции scale_fill_manual, как показано ниже:

library(ggplot2)
library(ggvoronoi)
library(dplyr)
for(i in df$years){
  #
  col  = c("1" = "green", "2" = "blue", "3" = "red")
  single_year <- df %>% 
    filter(years == i)
  #
  #
  plot <- ggplot(single_year,
                 aes(x=lat,
                     y=long)) +
    #
    geom_voronoi(aes(fill = cluster)) +
    #
    stat_voronoi(geom="path" )+
    #
    geom_point() +
    #
    labs(title = paste(i))+
    scale_fill_manual(values = col)
  #
  #
  ggsave(paste0(i,".jpeg"), plot = last_plot(), # Watch out for the SAVE!!!
         device = 'jpeg')
  #
}

и

Отвечает ли он на ваш вопрос?