[面板数据求助] 用stata怎么做空间计量要具体命令 [推广有奖]

11楼

ning715 发表于 2016-4-5 16:54:42

学习了

12楼

supplywinter

发表于 2016-6-21 16:24:51

楼主可以考虑用MATLAB做啊，有专门的空间计量包，网上也有命令教程

已有 1 人评分	学术水平	收起理由
admin_kefu	+ 5	热心帮助其他会员

总评分: 学术水平 + 5 查看全部评分

13楼

huangli 发表于 2016-6-26 15:57:12

      *===============================
      * 空间计量经济学与STATA操作心得
      *===============================

         * 作者: XX
* 学校：xxxxxx
* 邮  箱：758938561@qq.com

*--------------------------------------------------

*    ==内容目录==
* 引言
* 1.*.gal格式空间权重的导入
* 2.Moran指数的计算及Moran I指数图
* 3.空间面板矩阵的管理与生成
* 4.空间面板回归的SEM和SAR

*引言
*————————————————————————————————
* pwd --设定stata默认工作位置将空间权重矩阵及全部数据定义在同一个路径下
* cd--默认路径的设定
* 使用ssc install及findit来安装外部命令
*ssc install spwmatrix
*ssc install sg162.pkg
*ssc install spatwmat

      *------------------------
*->1.将gal空间权重矩阵导入stata
*------------------------
   pwd
   cd D:\Stata12.0\fish
   spwmatrix import using Export_Output.gal, wname(W)
   matrix list W
   svmat W
   save "D:\Stata12.0\fish\W.dta",replace
   *-------------------------------
   *以上命令的意思为：定义工作路径，将*.gal文件调入stata，转换为空间权重矩阵（cross）
   *在显示过程中，会提示有位置不存在相邻，没关系，是海南岛。
   *矩阵是0-1 binary格式，无误就存起来
*——————————————————————————————————
*-----------------------
*2计算moran指数
*-----------------------
use 1.dta,clear
spatwmat using W.dta,name(W)
spatgsa _2003-_2014,weights(W) moran twotail
*未标准化的空间权重矩阵不能计算局部moran
spatwmat using W.dta, name(WW) standardize
spatlsa _2003, weights(WW) moran twotail
spatlsa _2003, weights(WW) moran graph(moran) symbol(n)

*------------------------
*以上命令的意义为：样本存在的情况下，调用空间权重矩阵W，双边检验Moran I的P值，如果P值小于5%，说明存在空间相关性。
*生成标准化的空间权重矩阵WW，计算剧本Moran I指数并画出局部Moran I的图
*——————————————————————————————————
*-----------------------
*3.面板空间矩阵的管理与生成spmat
*-----------------------
      use W.dta,clear
      spmat dta W W1-W31
      matrix list W
      spmat save W using WW.spmat,replace
*矩阵已经变成spmat格式，面板模型可以识别并调用，同时以WW名称存在pwd给出的工作路径下。此时，必需要存在W.dta和WW.spmat的corss权重和panel权重
*——————————————————————————————————
*-----------------------
*4.空间面板回归的SEM和SAR
*-----------------------
use fishpanel.dta,clear
xtset province year
spmat use WW using WW.spmat
xtreg y h l s inv agdp p dc,fe r
estimate store xtregfe
xsmle y h l s inv p agdp dc ,emat(WW) model(sem) fe r nolog
estimate store ssemfe
xsmle y h l s inv agdp p dc ,wmat(WW) model(sar) r nolog
estimate store sar
*=====================================================
*上述结果中未有空间杜宾模型，SDM。空间杜宾是SLM和SEM的爸爸，如果SDM检验不通过，就退化成SLM或者SEM。（利用test c1=c2=0检验）
*部分内容来自陈强的stata第二版空间计量那章节。