#设置为自己的工作文件夹

> setwd("…\\gal2mat")

# 用nb2mat将gal文件转成矩阵

> library(spdep)

> mynb <-   read.gal("province31.gal")

> spmat1 <- nb2mat(mynb,   zero.policy=TRUE)

#用自编函数将gal文件转成矩阵

#载入函数gal2mat

> source("gal2mat.R")

> spmat2 <-   gal2mat("province31.gal")

Read 198 items

Warning messages:

1: In gal2mat("province31.gal") :   Unit 20 has no neighbor.

2: In FUN(newX[, i], ...) : The sum of row is   0.

#比较spmat1和spmat2是否一致，结果显示为一致

> all(spmat1 == spmat2)

[1] TRUE

#写出n*n空间权重矩阵

> write.csv(spmat2, "spmat.csv")

版权申明

本号所有图文资料，除特别说明外，其版权归浙江工业大学王庆喜领衔的“思达工作室”所有。转发、引用请注明原始出处。

网络链接

1、科学网博客：http://blog.sciencenet.cn/u/workshopofsdar

2、网易博客：http://wqx1976.blog.163.com/

3、人大经济论坛账号：R语言区域经济

4、知乎账号：sdar