*-SDM估计
xsmle y x1 x2 x3, wmat(W) model(sdm) fe type(both) robust nolog
eatat ic
***Test for SAR
test [Wx]x1=[Wx]x2=[Wx]x3=0
***Test for SEM
testnl ([Wx]x1=-[Spatial]rho*[Main]x1) ([Wx]x2=-[Spatial]rho*[Main]x2) ([Wx]x3=-[Spatial]rho*[Main]x3)
***Test for SAC
qui xsmle y x1 x2 x3, wmat(W) emat(W) model(sac) fe type(both) robust nolog
estat ic
LM检验
use quanzhong.dta(权重文件),clear
save "101w.dta",replace
spwmatrix import using 101w.dta,wname(w2b) dta xtw(n)n为观测对象数量
reg y v1 v2 v2
spatdiag,weights(w2b)