空间计量全流程代码，do文档直接替换变量即可做出
1、描述性统计
2、计算莫兰指数，看莫兰指数显著性，画莫兰指数散点图(局部)
(1)全局莫兰指数
(2)局部莫兰指数散点图
3、选择模型
(1)LM检验(SEM or SAR)
(2)LR检验
检验 SDM 模型能否退化为 SEM、SAR 模型
(3)WALD检验
检验模型的适配性，一般与 LR 检验同时出现
4、固定效应与随机效应检验选择检验，Hausman 检验
5、检验地区固定效应、时间固定效应以及双固定效应，三种效应哪个最适合本文的研究
6、模型回归
(1)SDM 模型回归(空间杜宾模型)
(2)SAR 模型回归(空间滞后模型)
(3)SEM 模型回归(空间误差模型)