问题描述
我正在尝试绘制纽约州的县级地图。我想根据每个县的工会化程度为其着色。我需要地图和图例具有四种离散的红色，而不是红色渐变。我需要图例以不重叠的标签/范围（例如0-25、26-50、51-75、76-100）显示这四种不同的颜色。
这是我的数据：
fips unionized
1 36001 33.33333
2 36005 86.11111
3 36007 0.00000
4 36017 0.00000
5 36021 0.00000
6 36027 66.66667
7 36029 40.00000
8 36035 50.00000
9 36039 0.00000
10 36047 82.85714
11 36051 0.00000
12 36053 100.00000
13 36055 30.76923
14 36057 0.00000
15 36059 84.37500
16 36061 81.81818
17 36063 60.00000
18 36065 50.00000
19 36067 71.42857
20 36069 0.00000
21 36071 55.55556
22 36073 0.00000
23 36079 100.00000
24 36081 92.15686
25 36083 50.00000
26 36085 100.00000
27 36087 87.50000
28 36101 0.00000
29 36103 63.88889
30 36105 0.00000
31 36107 0.00000
32 36111 50.00000
33 36113 50.00000
34 36115 100.00000
35 36117 0.00000
36 36119 73.33333
37 36121 0.00000
38 36123 0.00000
我已经成功地用渐变色制作了地图，但无法弄清楚如何在地图和图例中制作离散的颜色。
这是我的代码：
library(usmap)
library(ggplot2)
plot_usmap(regions = "counties",include = c("NY"),data = Z,values = "unionized") +
labs(title = "Percent Unionized",subtitle = "") +
scale_fill_continuous(low = "white",high = "red",na.value="light grey",name = "Unionization") + theme(legend.position = "right")
谢谢！
解决方法
这可以通过scale_fill_binned和guide_bins来实现。试试这个：
library(usmap)
library(ggplot2)
plot_usmap(regions = "counties",include = c("NY"),data = Z,values = "unionized") +
labs(title = "Percent Unionized",subtitle = "") +
scale_fill_binned(low = "white",high = "red",na.value="light grey",name = "Unionization",guide = guide_bins(axis = FALSE,show.limits = TRUE)) +
theme(legend.position = "right")
第二种选择是手动对变量进行装箱，然后使用scale_fill_manual来设置填充颜色，这使得设置标签变得容易，并具有自动添加NA的优点。对于色标，我使用colorRampPalette（默认情况下，colorRampPalette插值到rgb颜色空间中。要使用scale_fill_binned获得填充色，可以添加参数space = "Lab" ）。
library(usmap)
library(ggplot2)
Z$union_bin <- cut_interval(Z$unionized,n = 4,labels = c("0-25","26-50","51-75","76-100"))
plot_usmap(regions = "counties",values = "union_bin") +
labs(title = "Percent Unionized",subtitle = "") +
scale_fill_manual(values = colorRampPalette(c("white","red"))(5)[2:5],name = "Unionization") +
theme(legend.position = "right")
