用R做中国分省的出生性别比地图,地图数据见下面附件
地图数据.zip
大小:(380.35 KB)
马上下载
本附件包括:
- province2.shp
- province2.shx
- province2.dbf
代码如下:
library(rgdal)
library(classInt)
province<-readOGR("E:/statsoft/spss","province2")
pal<-terrain.colors(4)
q5<-classIntervals(province$xingbiebi,n=4,style="quantile")
q5Colours<-findColours(q5,pal)
plot(province,col=q5Colours)
legend("topleft",fill=attr(q5Colours,"palette"),legend=names(attr(q5Colours,"table")),bty="n")
得图形如下:(略)
但这个图肯定是出了严重错误,由下面截取的部分数据表就可看出:
错误一:最大值海南为135.64,而上图图例最大只到111.
错误二:西藏出生性别比仅为102.73,但地图依图例却在110.5之上。
作为对比,Geoda做出了正确的图
请高手修正上面R的代码,谢谢!