绘制美国热图与文字

Question

输入：

Percent.Turnout US.State 
70    CA 
80    NM 
76    RI 

我会为每个在美国50个州的数据。此外，US.State的州名缩写符合功能状态中的缩写.abb

我想要创建一个美国地图，其中Percent.Turnout印在每个州上。此外，使用ColorBrewer软件包，我想根据相对于其他状态的Percent.Turnout对每个状态进行着色。

我不是很熟悉ggplot语法，所以在基础R建议，将不胜感激（如果可行的话）

Answer 1

如果您想使用ggplot2，那么你所需要的重大的事情要做的是将状态缩写栏以小写形式映射到完整状态名称（为此，您可以使用state.name，但请确保在其上应用tolower()以获得正确的格式）。

从那里开始，只需将您的数据集加入州的地理空间信息并绘制数据即可。

# First, we need the ggplot2 library: 
> library(ggplot2) 
# We load the geospatial data for the states 
# (there are more options to the map_data function, 
# if you are intrested in taking a look). 
> states <- map_data("state") 
# Here I'm creating a sample dataset like yours. 
# The dataset will have 2 columns: The region (or state) 
# and a number that will represent the value that you 
# want to plot (here the value is just the numerical order of the states). 
> sim_data <- data.frame(region=unique(states$region), Percent.Turnout=match(unique(states$region), unique(states$region))) 
# Then we merge our dataset with the geospatial data: 
> sim_data_geo <- merge(states, sim_data, by="region") 
# The following should give us the plot without the numbers: 
> qplot(long, lat, data=sim_data_geo, geom="polygon", fill=Percent.Turnout, group=group) 

这是代码段的上面的输出：下面的代码段一步带你穿越这一步

现在，你说你想也将值Percent.Turnout添加到地图。在这里，我们需要找到各个州的中心点。您可以从我们上面检索到的地理空间数据（在states数据框中）计算出来，但结果看起来不会很令人印象深刻。值得庆幸的是，R具有在已经计算出我们国家的中心值，我们可以利用的是，如下所示：

# We'll use the state.center list to tell us where exactly 
# the center of the state is. 
> snames <- data.frame(region=tolower(state.name), long=state.center$x, lat=state.center$y) 
# Then again, we need to join our original dataset 
# to get the value that should be printed at the center. 
> snames <- merge(snames, sim_data, by="region") 
# And finally, to put everything together: 
> ggplot(sim_data_geo, aes(long, lat)) + geom_polygon(aes(group=group, fill=Percent.Turnout)) + geom_text(data=snames, aes(long, lat, label=Percent.Turnout)) 

这是最后一条语句的输出上面：