上篇帖子的关键词：

空间面板回归——用stata做出空间面板回归--距离倒数的权重矩阵

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的方法，就得到了结果。我的代码也贴到了本篇文章开头。

就这样，祝大家早日做出结果。