全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
4608 13
2012-05-01
悬赏 500 个论坛币 已解决
用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")


得图形如下:(略)

但这个图肯定是出了严重错误,由下面截取的部分数据表就可看出:
TT截图未命名2.bmp
错误一:最大值海南为135.64,而上图图例最大只到111.
错误二:西藏出生性别比仅为102.73,但地图依图例却在110.5之上。
作为对比,Geoda做出了正确的图


请高手修正上面R的代码,谢谢!
TT截图未命名3.bmp

原图尺寸 297.44 KB

TT截图未命名3.bmp

TT截图未命名.bmp

原图尺寸 752.18 KB

TT截图未命名.bmp

最佳答案

吉林小王子 查看完整内容

代码有错误,province$XINGBIEBI,这里要大写,你小写了 library(rgdal) library(classInt) province
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2012-5-1 00:24:51
代码有错误,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")
附件列表
tuxing.jpeg

原图尺寸 108.37 KB

tuxing.jpeg

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-5-1 15:22:06
呵呵,使者这两天你就在中国地图上纠结了
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-5-1 16:48:45
国家地图数据如何得到?谢谢!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-5-1 18:40:29
kk22boy 发表于 2012-5-1 15:22
呵呵,使者这两天你就在中国地图上纠结了
还是kk22boy兄知我呀
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

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

library(rgdal)
一句话点醒梦中人。
多谢小王子兄!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群