*===============================
* 空间计量经济学与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第二版空间计量那章节。