一)得到数据库文件和坐标文件
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
就可以生成中国34个省份观测单位8年的权重指标了。
注意:其中panel()项可以设置你要生成的权重观测单位数一般与原始文件一致,如果想减少几个省份,则先在原始文件减少,time()设置生成时间的维数,一般空间面板数据用的上。然后生成的权重数据分为横截面和面板类型分别保存在 Cross Section Weight Matrix File (Wcs.dta)和 Panel Weight Matrix File (Wxt.dta)中要查看只需输入 use Wcs.dta,clear 或者use Wxt.dta ,clear即可
(三)空间面板回归
主要说gmm估计
命令是 spgmmxt y x1 x2 , nc(7) wmfile(SPWxt) gmm(1) mfx(lin) test
暂时就说这么多,有空了再写详细点,主要想把在STATA中如何进行空间面板的操作写一下
评论前必须登录!
注册