空间计量最详尽代码及分析.pdfVIP

目录 1、描述性统计1 2、计算莫兰指数，看莫兰指数显著性，画莫兰指数散点图（局部）1 (1)全局莫兰指数1 (2)局部莫兰指数散点图2 3、选择模型2 （1）LM 检验（SEM or SAR ）2 （2 ）LR 检验————检验SDM 模型能否退化为SEM、SAR 模型3 （3 ）WALD 检验————检验模型的适配性，一般与LR 检验同时出现4 4 、固定效应与随机效应检验选择检验————Hausman 检验4 5、检验地区固定效应、时间固定效应以及双固定效应，三种效应哪个最适合本文的研究.5 6、模型回归5 （1）SDM 模型回归（空间杜宾模型）5 （2 ）SAR 模型回归（空间滞后模型）5 （3 ）SEM 模型回归（空间误差模型）6 （4 ）SAC 模型6 1、描述性统计 sum y x1 x2 x3 x4 logout, save(Descriptive) word replace: sum y x1 x2 x3 x4 //***输出描述性统计结果,点击Descriptive.rtf，即可出现word 版本的描述性统计表 2、计算莫兰指数，看莫兰指数显著性，画莫兰指数散点图（局部） 莫兰指数显著说明有必要使用空间计量模型。 cd F:\Stata14 //*定义路径 ssc install spmap //*安装spmap 命令 h spatreg //安装sptreg 命令, 以上三步非必要，视情况可以省略 use shuju //导入数据 spatwmat using F:\Stata14\W.dta, name(W) //调用空间权重矩阵，命名为W spatwmat using W.dta, name(W) standardize //将W 矩阵标准化 (1)全局莫兰指数 keep if year==2009 //*逐年计算全局莫兰指数，我计算的是2009-2015 年的。 spatgsa lny,weights(W) moran geary twotail keep if year==2010 spatgsa lny,weights(W) moran geary twotail keep if year==2011 spatgsa lny,weights(W) moran geary twotail keep if year==2012 spatgsa lny,weights(W) moran geary twotail keep if year==2013 spatgsa lny,weights(W) moran geary twotail keep if year==2014 spatgsa lny,weights(W) moran geary twotail keep if year==2015 spatgsa lny,weights(W) moran geary twotail (2)局部莫兰指数散点图 spatlsa a2015 ,weight(W) moran graph (moran) symbol(n) //*局部莫兰指数面板数据不能做，但是可以逐年做。第一列是省份，第二列变量命名为a2015. 是2015 年各省lny 的值。 //Moran 散点图第一象限至第四象限分别对应于某一区域与邻近区域之间的 4 种类型的 局部空间联系形势: 第一象限至第四象限分别表示 H－H 区域——高值地区被同样的高值 地区包围、L－H ——低值地区被低值地区包围、 L－L ——低值地区被低值地区包围， H －L ——高值地区被低值地区包围。 3、选择模型 （1）LM 检验（SEM or SAR ） 注意stata 的LM 检验必须先把空间权重矩阵转换成面板空间权重矩阵！ ssc install spcs2xt//先安装spcs2xt set matsize 5000 //这个命令我也不是很明白，写上去就对了。 use W //*使用矩阵W spcs2xt 北京-新疆,matrix(aaa)time(7) //把空间权重矩阵转换成面板空间权重矩阵，北京-新疆 是W 矩阵列变量；aaa 是面板空间权重矩阵的新名字；7 是你的研究年限，有几年就写几， 我的是2009-2015，因此是