R呈现权重结构初步，更加直观的看到邻接关系

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. options(digits=3)
   
2. ##设置工作路径，指定文件夹
   
3. setwd("F:\\china31")
   
4. getwd()
   
5. 
   
6. library(spdep)
   
7. library(maptools)
   
8. library(RColorBrewer)
   
9. library(classInt)
   
10. 
    
11. ###读取shp文件，并显示数据
    
12. poly<-readShapePoly("region.shp")
    
13. data<-slot(poly,"data")
    
14. class(poly)
    
15. 
    
16. ###生成rook标准邻接对象并生成权重对象
    
17. nb_rook<-read.gal("region.gal")
    
18. listw_rook<-nb2listw(nb_rook,style="B")
    
19. listw_rook
    
20. mat_rook<-listw2mat(listw = listw_rook)
    
21. mat_rook[1:11,1:11]
    
22. listw_rook$weights
    
23. 
    
24. ###查看邻接文件的属性，包括类别、id、是否对称等
    
25. attributes(nb_rook)
    
26. class(nb_rook)
    
27. ###汇总权重文件的邻接情况
    
28. summary(nb_rook)
    
29. cards<-card(nb_rook)
    
30. table(cards)
    
31. plot(cards,unlist(lapply(listw_rook$weights,sum)),xlab="number of links", ylab="row sums of weights")
    
32. attach(data)
    
33. colnames(data)
    
34. 
    
35. ###queen近邻标准在地图上的呈现
    
36. plot(poly,col="white",border="grey60")
    
37. coords<-coordinates(poly)
    
38. plot(nb_rook,coords,add=T)

复制代码

china31.zip (3.68 MB, 需要: 10 个论坛币)

2015-7-18 23:11:32 上传

shapefile
需要: 10 个论坛币  [购买]

扫码加我 拉你入群

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

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

分享0 收藏4 回帖

关键词：Library options Brewer Option BRARY

1. ##如果要展示高阶权重的话，运行下面的命令
   
2. nb6_rook<-nblag(nb_rook,6)
   
3. plot(nb6_rook[[5]],coords,add=T)
   
4. 
   
5. ##接着添加标记要素
   
6. attr<-c("浙江","云南","新疆","西藏","四川","陕西","山西","山东",
   
7.         "青海","宁夏","内蒙古","辽宁","江西","吉林","湖南","湖北",
   
8.         "黑龙江","河南","北京","天津","海南","贵州","广西","甘肃",
   
9.         "福建","安徽","上海","重庆","江苏","广东","河北")
   
10. text(coords,labels=attr,pos=4,col="black",cex=0.75,offset=0.05)
    
11. locator<-locator()
    
12. text(locator,labels = "任建辉制作",pos=4,col="dark red")

复制代码

虽不能下载 先感谢楼主