R做中国地图出错了

用R做中国分省的出生性别比地图，地图数据见下面附件
地图数据.zip (380.35 KB) 本附件包括：

• province2.shp
• province2.shx
• province2.dbf

2012-5-1 00:06:00 上传

代码如下：
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的代码，谢谢！

代码有错误，province$XINGBIEBI，这里要大写，你小写了

library(rgdal)
library(classInt)
province<-readOGR("E:\\","province2")
names(province)
province$XINGBIEBI
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")

呵呵，使者这两天你就在中国地图上纠结了

国家地图数据如何得到？谢谢！

kk22boy 发表于 2012-5-1 15:22
呵呵，使者这两天你就在中国地图上纠结了

还是kk22boy兄知我呀

吉林小王子 发表于 2012-5-1 00:24
代码有错误，province$XINGBIEBI，这里要大写，你小写了

library(rgdal)

一句话点醒梦中人。
多谢小王子兄！

dxystata 发表于 2012-5-1 16:48
国家地图数据如何得到？谢谢！

版主是问中国省级地图吗？我的是来自本论坛，感谢热心会员ywh19860616，详见下贴第16楼：
https://bbs.pinggu.org/thread-1410066-2-1.html

地图缺台湾省和钓鱼岛等，要加上。

dxystata 发表于 2012-5-1 20:41
地图缺台湾省和钓鱼岛等，要加上。

因为数据中没有台湾的数据，所以ywh19860616网友在绘制时删除了。
我也看着怪不舒服的，下来再向他要一下删除前的完整的shp文件。

耕耘使者 发表于 2012-5-1 20:57
因为数据中没有台湾的数据，所以ywh19860616网友在绘制时删除了。
我也看着怪不舒服的，下来再向他要一下 ...

记得给我一份。