使用R语言进行美国失业率时空分析（包括绘图）

今天写一篇利用R语言，针对面板数据的简单分析与绘图。让我们直接开始把。

一、数据准备

这次的示例数据非常简单，只有一个shp格式的美国区县矢量数据，我们在QGIS中打开数据查看一下它的属性表。事实上我们需要的数据都在属性表的字段中。

二、开始分析

1、引入必要的包，将本地的数据导入，查看数据是否正确。进行简单可视化。

library(sf)
library(tmap)

shp <- st_read("G:/Maping/demo/data/data/USA_data/ncovr/NAT.shp")

head(shp)

数据当中"UE60", "UE70", "UE80", "UE90"字段分别表示1960，1970，1980，1990年的美国各区县的失业率。下面我们别将各年份的数据可视化出来。

ue60 = tm_shape(shp) +
  tm_polygons("UE60", palette = "Reds", title = "1960 Unemployment Rate") +
  tm_compass(position = c(0.8, 0.86)) +
  tm_scale_bar(position = c(0.77,0.06), breaks = c(0,500,1000)) +
  tm_layout(legend.title.size = 0.8)
ue60

#如果想保存到本地，运行下面这句话就行了
tmap_save(ue60, "G:/Maping/ue60.png", dpi = 300)