R语言空间可视化：绘制英国脱欧投票地图

library(sp)library(raster)download.file("http://biogeo.ucdavis.edu/data/gadm2.8/rds/GBR_adm2.rds","GBR_adm2.rds")UK=readRDS("GBR_adm2.rds")UK@data[159,"HASC_2"]="GB.NR"plot(UK, xlim = c(-4,-2), ylim = c(50, 59), main="UK areas")

plot(IRL,add=TRUE)

plot(FR,add=TRUE)

referendum=ddply(referendum,.(Region,HASC_code),summarise,Remain=sum(Remain),Leave=sum(Leave))

> sum(referendum$Leave)/(sum(referendum$Leave)+sum(referendum$Remain))[1] 0.5189184

plot(UK, col = map_data$Brexit, border = "gray1", xlim = c(-4,-2), ylim = c(50, 59), main="How the UK has voted?", bg="#A6CAE0")plot(IRL, col = "lightgrey", border = "gray1",add=TRUE)plot(FR, col = "lightgrey", border = "gray1",add=TRUE)

plot(UK, col = "grey", border = "gray1", xlim = c(-4,-2), ylim = c(50, 59),bg="#A6CAE0")plot(IRL, col = "lightgrey", border = "gray1",add=TRUE)plot(FR, col = "lightgrey", border = "gray1",add=TRUE)choroLayer(spdf = UK,df = map_data,var = "Percentage_Remain",breaks = seq(0,100,10),col = cols,legend.pos = "topright",legend.title.txt = "",legend.values.rnd = 2,add = TRUE)