用stata生成空间权重的步骤和操作详解

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

以下都用中国的例子
（一）得到数据库文件和坐标文件
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中如何进行空间面板的操作写一下！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享4 收藏548 回帖

关键词：Stata 空间权重 tata Coordinates coordinate 中国地图 数据库文件 下划线 英文 空间

非常感谢楼主！期待您写出运用stata分析空间面板数据模型的详细步骤，如原始数据的导入、几类空间面板权重的生成、moran\G指数的计算，真是大功德一件了！

首先非常感谢楼主提供的宝贵经验！其次我仔细揣摩你提供的方法，有两点疑问，希望你能够帮忙解答。其一，你在第一步转换得到的数据库和坐标文件的作用是什么，或者说它与第二步的关系是什么。其二，我们知道SPWEIGHT命令的V1和V2只能是0或者1，那么请问若是构建地理距离权重矩阵或者经济距离权重矩阵，V1和V2又该如何设定呢，若是不能设定，那么如何得出面板数据的空间权重矩阵呢，这些问题一直困扰我，所以希望你能帮忙解答一下。总结来说，就是我该如何通过STATA构建地理/经济距离权重矩阵，详细步骤是什么，麻烦你了。

zhichi  ge  xiexie

it is goooooooooooooood

不错，期待LZ接下来的详解

楼主 我在stata中输入spweightxt v1 v2 , panel(34)  time(5) matrix(W) ptable，得不出结果。
显示matrix __00000K11 would have missing values
楼主帮我看看。geoda生成的权重文件在附件中。谢谢！

geoda生成的权重文件（按v1和v2整理的文件）。

谢谢好经验~！！！！！！！！！！

非常感谢您的分享，但是我想请问， spweightxt v1 v2 , panel(34)  time(8) matrix(W) eigw table中的“ eigw table”是什么意思呢？