R Пользовательская цветовая шкала для печати

Question

Я пытаюсь создать пользовательский масштаб, который будет лучше отображать имеющиеся у меня данные.Для первого графика данные были относительно однородными и находились в диапазоне от 0 до 3. Однако для моего второго набора данных большинство данных меньше 100, а примерно 90% данных меньше 1000. (графики внизу)

Столбец со значениями, которые я пытаюсь построить, называется «data1», а значения представляют собой целые числа в диапазоне от 0 до 10000.

Я пытался применить код на https://gis.stackexchange.com/questions/178597/how-to-create-a-continuous-scale-with-distinct-custom-color-and-value-breaks-wit, вставив

map.df$brks <- cut(map.df$data1, 
                 breaks=c(0, 1, 10, 25, 100, 250, 1000, 10000), 
                 labels=c("0", "1 - 10", "11 - 25", "26 - 100", 
                          "101 - 250", "251 - 1000", "1001 - 10000"))

, где находится закомментированная строка, и заменив аргумент заполнения на brks в вызове функции ggplot, но я получаю следующую ошибку: Ошибка: дискретное значение передается в непрерывном масштабе.

stateNames <- c("alabama", "alabama", "arizona", "arkansas", "california", 
                "colorado", "connecticut", "florida", "georgia")
countyNames <- c("autauga", "blount", "apache", "bradley", "orange", 
                 "boulder", "new haven", "flagler", "turner")
dataCol <- c(0, 5, 15, 50, 150, 1000, 5000, 249, 30)

library(ggplot2)
map.county <- map_data('county')
counties   <- unique(map.county[,5:6])
data_map <- merge(counties, data.frame(region=stateNames, 
                                     subregion=countyNames, 
                                     data1= dataCol),
                  by=c("region","subregion"), all.x=T, all.y=F
)

data_map$data1[which(is.na(data_map$data1))] <- fillValue
for(i in 1:length(data_map$data1)) {
  if(is.na(data_map$data1[i])) {
    data_map$data1[i] <- 0
  }
}
library(data.table)
map.county <- data.table(map_data('county'))
setkey(map.county,region,subregion)
data_map <- data.table(data_map)
setkey(data_map,region,subregion)
map.df <- map.county[data_map]

############################################################# 

ggplot(map.df, aes(x=long, y=lat, group=group, fill=data1)) + 
  scale_fill_gradientn("",colours=brewer.pal(9,"YlOrRd"))+
  geom_polygon(colour = borderColor, size = 0.001) +
  coord_map()  +
  expand_limits(x = map.county$long, y = map.county$lat) +
  scale_x_continuous(breaks = NULL) + 
  scale_y_continuous(breaks = NULL) +
  labs(x = "", y = "") +
  labs(fill=dataName)  +
  theme(legend.position = "bottom", panel.background = element_blank())

У кого-нибудь есть советы, как мне создать собственный масштаб?Заранее спасибо!

Answer 1

Вы почти у цели.Пояснения в аннотированном коде ниже:

# add include.lowest = TRUE, otherwise the 0 values will become NAs
map.df$brks <- cut(map.df$data1, 
                   breaks = c(0, 1, 10, 25, 100, 250, 1000, 10000), 
                   labels = c("0", "1 - 10", "11 - 25", "26 - 100", 
                            "101 - 250", "251 - 1000", "1001 - 10000"),
                   include.lowest = TRUE)

# use scale_fill_brewer rather than scale_fill_gradientn; it expects
# a categorical variable (e.g. brks) while scale_fill_gradientn expects
# a continuous variable (e.g. data1).
# also, I suggest placing the legend at the side, as a vertical stack of 
# boxes, which resembles a color bar.
# theme_void can be used to remove all backgrounds, axis labels, etc.
ggplot(map.df, aes(x=long, y=lat, group=group, fill=brks)) + 
  geom_polygon(colour = "grey", size = 0.001) +
  coord_map()  +
  expand_limits(x = map.county$long, y = map.county$lat) +
  scale_fill_brewer(palette = "YlOrRd") +
  theme_void() +
  theme(legend.position = "right")