以下都用中国的例子
(一)得到数据库文件和坐标文件
1.之所以这一步是由于在实际分析中会删除个别省份,然后基于剩余省份制作shp文件后生成空间权重文件。找到中国地图省级、省会城市、城市、县的dbf和shp文件,注意:各级地图的dbf和shp文件的名称必须相同且是英文,如province.dbf和province.shp,provcap.dbf和provcap.shp等2.在stata中安装spmap、shp2dta和mif2dta命令是ssc install spmap\ shp2dta\mif2dta 3.在指定目录下,运行命令shp2dta using ProCap, data("ProCap_data") coor("ProCap_coordinates"),其中下划线部分分别是使用的数据库和坐标名称,转换后的数据库和坐标文件名称,然后用命令 merge ID using procap_coordinates,将文件复制到txt文档中编辑成符合geoda格式(第一行写上截面单位数和变量数,第二行是变量名称,用英文下双引号引起,用逗号隔开,剩下行是观测值也要用逗号隔开,当然可以先在excel中编辑再拷过来)然后打开geoda~tools~shape~points from ascii,主要设置位置,然后保存为shp文件,最后基于此文件生成空间权重矩阵,当然质点文件只能生成knearest和threshold distance距离权重。 (二)如果用geoda生成了权重数据,那么可以通过以下方法转换成stata文件 1.用geoda生成不同类型的权重文件,比如province_knear.gwt,打开另存为province_knear.txt 2.用记事本打开后复制出来或者在stata中输入命令insheet using province_knear.txt,删除关键词,保留位置变量和邻近数量变量,并重命名为v1和v2,并另存为spweight3.dta 3.在stata中安装spweightxt命令 4.输入命令
clear all
sysuse spweight3.dta, clear
list v1 v2
spweightxt v1 v2 , panel(34) time(8) matrix(W) eigw table