以全国高铁数据为例:
1.先在网上搜索全国的.shp 和.dbf地图格式文件
2.用geoda软件读入上述两个文件
3.在geoda中点击W生成空间权重矩阵,即.gal格式空间权重文件,filename.gal.
4.将生成的geoda中空间权重矩阵filename.gal文件读入stata中生成dat文件,该文件r及matalab可以直接用,但stata不行。代码如下:
spwmatrix import using D:\stata13\ado\personal\cjjjddjs.gal, wname(wcontig) rowstand xport(cjjjd,dat) cjjjd,dat为生成的dat文件。
5空间权重矩阵转换为.spmat格式,代码如下;
insheet using "E:\ZZZZ\sh\cjjjd,dat"
spmat dta weight v1-v149, replace normalize(row)
spmat save weight using weight.spmat
6.读入数据
use "E:\ZZZZ\she.dta"
7. 定义面板
xtset admincode year 个体 时间
8. 使用spmat文件
spmat use test using test.spmat
9. xsmle ca_3 tp pcgrp awesw aspcb psti denlh sd200 , wmat(test) emat(test) model(sac) fe type(ind) noeff nolog
ca_3为Y变量,tp pcgrp awesw aspcb psti denlh sd200为X变量,test为权重矩阵。