程序gwt2mat将GeoDa生成的gwt格式文件转成n*n矩阵，从而可以适用于matlab空间计量经济分析程序。

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

> setwd("…\\20180526将gwt权重文件转成n-by-n矩阵\\gwt2mat")

#用spdep包中的函数处理，偏离预期

> library(spdep)

> spmat1 <-   read.gwt2nb("province31.gwt")

Warning message:

In read.gwt2nb("province31.gwt") :   region.id not named POLY_ID

#用nb2mat转成矩阵，所有非零元素均为0.03，达不到预期效果

> spmat1 <- nb2mat(spmat1)

#用gwt2mat将gwt权重文件转为n*n矩阵

> source("gwt2mat.R")

> spmat2 <-   gwt2mat("province31.gwt")

Read 2790 items

#将spmat矩阵元素取平方倒数,再行标准化

> spmat2 <- apply(spmat2, 1,   function(x) {

+   x[x   != 0] <- 1 / (x[x != 0])^2

+     return(x / sum(x))

+ })

> spmat2 <- t(spmat2)

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

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

版权申明

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

网络链接

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

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

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

4、知乎账号：sdar

公众号功能

1、致力于打造中国R语言空间数据分析第壹号；

2、探讨R语言、ArcGIS、matlab、GeoDa、Stata、SPSS等软件在空间数据分析中的应用；

3、讨论产业集聚、空间溢出、区域创新网络、城市群发展等热门研究主题。