20170 发表于 2018-10-31 10:34
你好,可否把您的工具箱具体讲一下呢?被这个问题困扰好久了
双空间权重矩阵的回归可以通过stata中的spm命令来实现,当然也可以通过西佛吉尼亚大学区域研究所的某位教授的matlab命令来实现。
下面是我以前写的双空间权重矩阵的命令,供你参考,但是其中那个mmat文件是如何创建的,我也忘记了,你可以在stata中看一下spm这个命令的帮助。祝好运!
cd "J:\Users\yu-xu\Desktop\map\BS_rail\BS_rail"
use BSbefore2010.dta, clear
xtset i year
mata mata matuse W1.mmat, replace
mata st_matrix("W1",W1)
mata mata matuse W2.mmat, replace
mata st_matrix("W2",W2)
spm mean lnk lnex lnden, model(sar) sarwmat(W1) sarw2mat(W2)
est store t1
use BSafter2010.dta, clear
xtset i year
spm mean lnk lnex lnden, model(sar) sarwmat(W1) sarw2mat(W2)
est store t2
local gg "t1 t2"
esttab `gg' using ".\panel3.rtf", mtitle(`gg') p(%6.3f) b(%6.3f) nogap s(r2 N F chi2) obslast star(* 0.10 ** 0.05 *** 0.01) compress replace