发帖
雷达卡
上篇帖子的关键词：
0、所有代码
use "G:\Program Files(x86)\Stata15\E1.dta"
xtset _ID year
spset
xtreg aqi job_house_01 log_conges_dayconges_job_house_01 log_pergdp log_buslog_gas fdi i.year, fe
xtreg aqi job_house_01 log_conges_dayconges_job_house_01 log_pergdp log_bus log_gasfdi i.year, re
spmatrix create idistance M if year == 2017
spxtregress aqi job_house_01 log_conges_dayconges_job_house_01 log_pergdp log_buslog_gas fdi i.year, fe dvarlag(M)
spxtregress aqi job_house_01 log_conges_dayconges_job_house_01 log_pergdp log_buslog_gas fdi i.year, re dvarlag(M)
1、引言，最终还是做出空间面板回归了，找到的回归有三类方法，我用的第二类。
第一类：陈强书上的 —— xlmle 优点是模型多，方法多，缺点是不知道怎么建立数据；
第三类：忘了，但确实存在。
第二类：stata网传的方法， spxtregressy x1 x2, fe dvarlag(M)
下面我对第二类方法进行讲解。
2、数据、资料的准备。
stata官网上的数据来源：
https://www.stata.com/gsearch.php?q=spatial&restrict=SP&site=stata&client=stata&proxystylesheet=stata&output=xml_no_dtd
https://www.stata.com/gsearch.php?q=spatial&restrict=SP&site=stata&client=stata&proxystylesheet=stata&output=xml_no_dtd
2-1 说明：右上图的download是整本书关于空间回归的实际案例教学；（不同于陈强），是第二个教学书籍。
2-2 说明2：这是对方案例中的数据
spxtregress 命令下的数据，就是我们要搞的空间面板回归。 2-1下载的整本pdf，也有对应的下载命令 pdf 199页
2-3 说明3：数据
总之，打开那个 犯罪率1960——1990.dta 是下面这个样子（如果你要学pdf的案例的话）（有时候直接点开那个dta文件，是另外一个错码的样子——你要把dta文件放到stata根目录下，然后用 use homicide_1960_1990 的命令打开。）
我们的目标是做空间面板，网传的方法都是 面板数据 + 自己做一个权重矩阵，我学到这一步，卡在了 不知道如何把权重矩阵与面板数据一一对应，然后做出空间面板回归。这个方法有待后来人去分享把。
我这里分享的是第二个做空间面板的方法。——在带shp，带地理坐标系的dta文件上，直接加上面板数据，然后用2-1下载的stata说明手册的spxtregress的方法做空间回归，我做成了面板回归，还tm是一一对应的，解决了我们的噩梦。
3、自己的空间面板数据加工
3-1 首先你得会用stata做空间权重矩阵，（我的另外一篇文章讲了，用arcgis把地理坐标系转成shp，用shp导入stata建立dta，此时的dta 带有空间数据的shp。然后用这个dta就可以建立距离倒序的 权重矩阵。 ——是的，我有了shp文件，也许用GEODA软件能做出contiguity的权重文件，后来人去试试。）
这个数据形式是另一篇文件的结果，也是我们这篇做空间面板回归的开始。
3-2 我们要做空间面板回归，数据加工就是——在后面加上我们的Y x1 x2 X...,
——我就手动复制的。！此时就完成了数据准备，跟2-3 的教学案例原始数据文件一样了。
4 回归。
按照pdf教学案例去做把。 P196-199的方法，就得到了结果。我的代码也贴到了本篇文章开头。
就这样，祝大家早日做出结果。
提升卡
置顶卡
沉默卡
变色卡
抢沙发
千斤顶
显身卡
京公网安备 11010802022788号
