- global M2 LY LX1 LX2 LX3 LX4 LX5 LX6 LX7 LX8
- ***创建spmat格式权重文件W
- use W, clear
- ***也可用该命令定义空间权重矩阵:spmat dta W a1-a25, normalize(row) replace
- spmat save W using W.spmat, replace
- spmat use W using W.spmat,replace
- ***生成空间加权解释变量W_LX1-W_LX8变量
- use data, clear
- spset N
- foreach var in $M2 {
- spmat lag double W_`var' W `var'
- }
- ***估计SAR/SEM/SDM模型
- spreg ml LY LX2 LX3 LX4 LX5 LX7 LX8, id(N) dlmat(W) robust //SEM
- spreg ml LY LX2 LX3 LX4 LX5 LX7 LX8, id(N) elmat(W) robust //SAR
- spreg ml LY LX2 LX3 LX4 LX5 LX7 LX8, id(N) dlmat(W) elmat(W) robust //SAC
- spreg ml LY LX2 LX3 LX4 LX5 LX7 LX8 W_LX2 W_LX3 W_LX4 W_LX5 W_LX7 W_LX8, id(N) elmat(W) robust //SDM,在SAR基础上自行加入空间加权解释变量W_LX1-W_LX8变量)
缺陷:没有现成的空间效应分解命令,只能手工计算。因此,做起来不如Matlab方便