全部版块 我的主页
论坛 计量经济学与统计论坛 五区 数据交流中心 数据求助
1463 0
2017-04-19
R软件做空间计量分析时需要导入权重矩阵数据集,导进去后期用的时候说它不是空间近邻文件,请问空间近邻文件长什么样啊???应该怎样建立空间近邻的数据啊??代码如下:

###########读取权重矩阵并转化为权重对象
wgt=read.table("2015权重.txt",header=T)
mt=as.matrix(wgt)
row.names(mt)=c("长春市","吉林市","四平市","辽源市","通化市","白山市","松原市","白城市","延边朝鲜族自治州")
mt
结果为:
row.names(mt)=c("长春市","吉林市","四平市","辽源市","通化市","白山市","松原市","白城市","延边朝鲜族自治州")
> mt
                    长春市    吉林市    四平市    辽源市    通化市    白山市
长春市           0.0000000 0.3333333 0.3333333 0.0000000 0.0000000 0.0000000
吉林市           0.1666667 0.0000000 0.1666667 0.1666667 0.1666667 0.1666667
四平市           0.2500000 0.2500000 0.0000000 0.2500000 0.0000000 0.0000000
辽源市           0.0000000 0.3333333 0.3333333 0.0000000 0.3333333 0.0000000
通化市           0.0000000 0.3333333 0.0000000 0.3333333 0.0000000 0.3333333
白山市           0.0000000 0.3333333 0.0000000 0.0000000 0.3333333 0.0000000
松原市           0.3333333 0.0000000 0.3333333 0.0000000 0.0000000 0.0000000
白城市           0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

> ############将权重矩阵空间化
> crs=CRS(projargs=as.character(NA))
> crs
CRS arguments: NA
> s_wgt=SpatialPoints(mt,proj4string=crs)
> s_wgt
SpatialPoints:
                    长春市    吉林市    四平市    辽源市    通化市    白山市
长春市           0.0000000 0.3333333 0.3333333 0.0000000 0.0000000 0.0000000
吉林市           0.1666667 0.0000000 0.1666667 0.1666667 0.1666667 0.1666667
四平市           0.2500000 0.2500000 0.0000000 0.2500000 0.0000000 0.0000000
辽源市           0.0000000 0.3333333 0.3333333 0.0000000 0.3333333 0.0000000
通化市           0.0000000 0.3333333 0.0000000 0.3333333 0.0000000 0.3333333
白山市           0.0000000 0.3333333 0.0000000 0.0000

oordinate Reference System (CRS) arguments: NA

> ##########计算Moran'I系数
> mor_wgt=read.dat2listw(s_wgt)
Error in read.table(file) :
  'file' must be a character string or connection
> mor_wgt
错误: 找不到对象'mor_wgt'


这个应该怎么弄啊?



二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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