R构建空间距离权重的说明，抛砖引玉

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

在做空间计量的时候，经常碰到如何构建距离权重的问题，那么spdep包是如何实现的，跟你想的是不是一样呢，那么我试着根据我的理解展示一下，欢迎拍砖指正。

1. > ##加载软件包
   
2. > library(maptools)
   
3. > library(spdep)
   
4. > ##读取shp文件和权重文件（.gal)
   
5. > columbus <-readShapePoly(system.file("etc/shapes/columbus.shp",package="spdep")[1])
   
6. > col.gal.nb <- read.gal(system.file("etc/weights/columbus.gal",package="spdep")[1])
   
7. > summary(col.gal.nb)
   
8. Neighbour list object:
   
9. Number of regions: 49
   
10. Number of nonzero links: 230
    
11. Percentage nonzero weights: 9.579342
    
12. Average number of links: 4.693878
    
13. Link number distribution:
    
14. 
    
15. 2  3  4  5  6  7  8  9 10
    
16. 7  7 13  4  9  6  1  1  1
    
17. 7 least connected regions:
    
18. 1 6 31 39 42 46 47 with 2 links
    
19. 1 most connected region:
    
20. 20 with 10 links
    
21. > ##生成点对应的横纵坐标
    
22. > coords <- coordinates(columbus)
    
23. > ##nbdists计算的是两点间的欧式距离。
    
24. > dlist <- nbdists(col.gal.nb, coords)
    
25. > plot(col.gal.nb,coords)
    
26. > text(coords,labels=attr(col.gal.nb,"region.id"),pos=4)
    
27. > ##最后给出邻接关系中，两点距离的倒数
    
28. > dlist1 <- lapply(dlist, function(x) 1/x)

复制代码

   值得说明的是：1、nbdists命令实际算的仅是邻居的距离，最后也不是矩阵的形式。
                2、需要学会编写简单的函数和循环条件语句
     

扫码加我 拉你入群

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

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

分享0 收藏8 回帖

关键词：抛砖引玉 distribution percentage connected neighbour 空间

看看再说

怎么你前面两个包安装 不了呢

Nelsh--Deng 发表于 2015-10-16 19:39
怎么你前面两个包安装 不了呢

install.packages(c("maptools","spdep")

jianhui80 发表于 2015-10-16 23:13
install.packages(c("maptools","spdep")

你好，你上面的程序都是对的，我刚刚已经演示了，请问你会做我国各个省市之间的空间距离倒数 的空间权重矩阵吗？希望如果你知道，可以教我一下，我现在有点急用，谢谢了

jianhui80 发表于 2015-10-16 23:13
install.packages(c("maptools","spdep")

你好，你上面的程序都是对的，我刚刚已经演示了，请问你会做我国各个省市之间的空间距离倒数的空间权重矩阵吗？希望如果你知道，可以教我一下，我现在有点急用，谢谢了

能不能加一下你QQ ，谢谢了

加上一条，通过nb2listw(col.gal.nb, glist=dlist1)转换，这样才能运用到空间计量模型中

求问大神，R可以允许我自己用excel写空间权重矩阵吗？

是不是必须要像geoda一样，必须有个shp文件、然后软件自己做的空间权重矩阵才能使用？

其实用excel自己算空间矩阵，更灵活，可是不知如何导入R，在使用moran.test时就遇到了问题，提示
Error: unexpected input in "moran.test(morandata,nb2listw?

原程序如下：
data<-read.csv("F:/neizidiqu.csv",header=F)
juzhen<-read.csv("F:/linjie.csv",header=F)
morandata<-data[,1]
juzhenone<-as.matrix(juzhen)
moran.test(morandata,nb2listw（juzhenone）)

edsioncjlu 发表于 2017-2-19 08:42
求问大神，R可以允许我自己用excel写空间权重矩阵吗？

是不是必须要像geoda一样，必须有个shp文件、然后 ...

nb listw matrix三者可以相互转换