请选择 进入手机版 | 继续访问电脑版
楼主: fei355
110416 222

[学习心得] 用stata生成空间权重的步骤和操作详解   [分享]

副教授

22%

还不是VIP/贵宾

-

威望
0
论坛币
19691 个
通用积分
659.2204
学术水平
78 点
热心指数
80 点
信用等级
66 点
经验
15012 点
帖子
327
精华
1
在线时间
855 小时
注册时间
2008-12-4
最后登录
2021-3-14

fei355 发表于 2013-8-19 14:38:54 |显示全部楼层
以下都用中国的例子
(一)得到数据库文件和坐标文件
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中如何进行空间面板的操作写一下!
关键词:Stata 空间权重 tata Coordinates coordinate 中国地图 数据库文件 下划线 英文 空间

已有 17 人评分经验 论坛币 学术水平 热心指数 信用等级 收起 理由
李氏玉 + 1 + 1 精彩帖子
xujingtang + 160 精彩帖子
np84 + 100 精彩帖子
米霍克 + 2 + 4 + 3 + 3 精彩帖子
红色沸鱼 + 1 + 1 精彩帖子
liangshuang522 + 3 + 3 热心帮助其他会员
xddlovejiao1314 + 80 + 60 + 5 + 5 + 5 精彩帖子
百变星空1227 + 1 + 1 + 1 精彩帖子
DaphneDawn + 1 + 1 + 1 对论坛有贡献
日新少年 + 1 精彩帖子

总评分: 经验 + 440  论坛币 + 92  学术水平 + 24  热心指数 + 22  信用等级 + 17   查看全部评分

本帖被以下文库推荐

stata SPSS
笑笑笑 发表于 2013-9-7 17:23:56 |显示全部楼层
非常感谢楼主!期待您写出运用stata分析空间面板数据模型的详细步骤,如原始数据的导入、几类空间面板权重的生成、moran\G指数的计算,真是大功德一件了!
纵一苇之所如,凌万顷之茫然。
回复

使用道具 举报

PX0706 发表于 2013-9-8 12:23:19 |显示全部楼层
首先非常感谢楼主提供的宝贵经验!其次我仔细揣摩你提供的方法,有两点疑问,希望你能够帮忙解答。其一,你在第一步转换得到的数据库和坐标文件的作用是什么,或者说它与第二步的关系是什么。其二,我们知道SPWEIGHT命令的V1和V2只能是0或者1,那么请问若是构建地理距离权重矩阵或者经济距离权重矩阵,V1和V2又该如何设定呢,若是不能设定,那么如何得出面板数据的空间权重矩阵呢,这些问题一直困扰我,所以希望你能帮忙解答一下。总结来说,就是我该如何通过STATA构建地理/经济距离权重矩阵,详细步骤是什么,麻烦你了。
已有 2 人评分论坛币 学术水平 热心指数 信用等级 收起 理由
小小小拖鞋 + 1 + 1 + 1 精彩帖子
crystal8832 + 10 + 1 鼓励积极发帖讨论

总评分: 论坛币 + 10  学术水平 + 1  热心指数 + 2  信用等级 + 1   查看全部评分

回复

使用道具 举报

gdczlhd 发表于 2013-9-8 22:00:58 |显示全部楼层
zhichi  ge  xiexie
寻梦?撑一支长篙,向青草更青处漫溯;满载一船星辉,在星辉斑斓里放歌。
回复

使用道具 举报

xge2000 发表于 2013-9-10 07:56:11 |显示全部楼层
it is goooooooooooooood
回复

使用道具 举报

evarei2 发表于 2013-9-20 13:23:50 |显示全部楼层
不错,期待LZ接下来的详解
回复

使用道具 举报

xiaoping2006 发表于 2013-11-7 16:59:17 |显示全部楼层
楼主 我在stata中输入spweightxt v1 v2 , panel(34)  time(5) matrix(W) ptable,得不出结果。
显示matrix __00000K11 would have missing values
楼主帮我看看。geoda生成的权重文件在附件中。谢谢!

sj_regiondist.txt

10.34 KB

回复

使用道具 举报

xiaoping2006 发表于 2013-11-7 17:04:07 |显示全部楼层
geoda生成的权重文件(按v1和v2整理的文件)。

spweight3.xls

30.5 KB

已有 1 人评分论坛币 学术水平 收起 理由
我的素质低 + 20 + 3 热心帮助其他会员

总评分: 论坛币 + 20  学术水平 + 3   查看全部评分

回复

使用道具 举报

大同爱跻 在职认证  发表于 2013-11-27 10:31:54 |显示全部楼层
谢谢好经验~!!!!!!!!!!
回复

使用道具 举报

brave-man 发表于 2014-1-3 18:07:45 |显示全部楼层
非常感谢您的分享,但是我想请问, spweightxt v1 v2 , panel(34)  time(8) matrix(W) eigw table中的“ eigw table”是什么意思呢?
已有 1 人评分经验 热心指数 收起 理由
crystal8832 + 10 + 1 鼓励积极发帖讨论

总评分: 经验 + 10  热心指数 + 1   查看全部评分

回复

使用道具 举报

您需要登录后才可以回帖 登录 | 我要注册

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2021-4-17 19:13