大家好，我是带我去滑雪！
地图相关图形绘制具有许多优点，这些优点使其在各种领域和应用中非常有用。例如：地图相关图形提供了一种直观的方式来可视化数据，使数据更容易理解和分析。通过地图，可以看到数据的空间分布，识别趋势和模式，并从中获取有关地理位置的见解。下面使用R语言绘制地图相关图。
（1）数据
使用美国50个州2005年的收入分配不均等指数，其中的6个收入不均等指数，除了泰尔熵指数之外，均介于0到1之间，数值越大代表收入分配越不均匀。详细数据在文末。
（2）收入不均等指数数据整理与点状图绘制
library(lattice)
ineq50=read.csv("ineq50.csv", header=F)
names(ineq50)=c("ATKIN05","GINI","RMEANDEV","THEIL","TOP1","TOP10")
save(ineq50, file="ineq50.RData")
load("ineq50.RData")
attach(ineq50)
vn=colnames(ineq50)
j=2
index=ineq50[,vn[j]]
data=data.frame(name = state.name,region = state.region,y=index)
region=state.region
dotplot(reorder(name, index)~index, ineq50, main=c("Inequality Index"),xlab=vn[j])
输出结果：
通过图像可以发现，康涅狄格州的收入分配最不平均，西弗吉尼亚州则最好。
（3）实现切段绘图
b=0.9 # Define cutting point
cuts =shingle(index,intervals = rbind(c(0, b),c(b, 2)))
dotplot(reorder(name, index) ~ index| cuts , ineq50,
strip = FALSE, layout = c(2, 1), levels.fos = 1:50,
scales = list(x = "free"), between = list(x = 0.5),
xlab = "Theil income inequality index",
par.settings = list(layout.widths = list(panel = c(2, 1))))
输出结果：
（4）按地理区域呈现收入不均等指数数据
ineq50$region=with(ineq50, reorder(region, index, median))
ineq50$name=with(ineq50,reorder(reorder(name, index), as.numeric(region)))
dotplot(name~ index | region, ineq50, strip = F, strip.left = T, layout = c(1, 4),
scales = list(y = list(relation = "free")),xlab = vn[j])
输出结果：
（5）在美国地图上呈现收入不均等指数
state.info = data.frame(name = state.name,long = state.centerx,lat = state.centery,y=index)
library(maps)
state.map = map("state", plot=FALSE, fill = FALSE)
panel.3dmap = function(..., rot.mat, distance, xlim, ylim, zlim,
xlim.scaled, ylim.scaled, zlim.scaled)
{
scaled.val = function(x, original, scaled) {
scaled[1] + (x - original[1]) * diff(scaled) / diff(original)
}
m = ltransform3dto3d(rbind(scaled.val(state.map$x, xlim, xlim.scaled),
scaled.val(state.map$y, ylim, ylim.scaled),zlim.scaled[1]),
rot.mat, distance)
panel.lines(m[1,], m[2,], col = "grey76")
}
cloud(index ~ long + lat, state.info,
subset = !(name %in% c("Alaska", "Hawaii")),
panel.3d.cloud = function(...) {
panel.3dmap(...)
panel.3dscatter(...)
},
type = "h", scales = list(draw = FALSE), zoom = 1.2,xlim = state.map$range[1:2],
ylim = state.map$range[3:4],xlab = NULL, ylab = NULL, zlab = NULL,
aspect = c(diff(state.maprange\[3:4\])/diff(state.maprange[1:2]), 0.3),
panel.aspect = 0.75, lwd = 2, screen = list(z = 30, x = -70),
par.settings = list(axis.line = list(col = "transparent"),
box.3d = list(col = "transparent", alpha = 0)))
输出结果：
可以发现，收入分配不均匀的州多集中在美国的东部和西部。
需要数据集的家人们可以去百度网盘（永久有效）获取：
链接：https://pan.baidu.com/s/1E59qYZuGhwlrx6gn4JJZTg?pwd=2138
提取码：2138
更多优质内容持续发布中，请移步主页查看。
点赞+关注,下次不迷路！