爱笑的傻子 发表于 2019-12-1 16:12 
请问您解决了吗,stata可以做LR、LM和ward检验吗,有没有具体命令呀,非常感谢!
LM可以把空间权重矩阵拓展到面板,然后用 spatdiag命令,拓展矩阵 spcs2xt s1-s30,matrix(W) time(13) ,time里面是年份;wald和LR都是在论坛里看到的,具体帖子忘记了。。。
wald检验-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
LR检验:
xsmle ratioef gdp cpi fdi_service_number population income ,fe model(sdm) wmat(ivdw) type(both) nolog effects nsim(500) durbin(gdp population fdi_service_number) ///DSM
est store sdm_fdiservicenumber
xsmle ratioef gdp cpi fdi_service_number population income ,fe model(sar) wmat(ivdw) type(both) nolog effects nsim(500) ///SAR
est store sar_fdiservicenumber
xsmle ratioef gdp cpi fdi_service_number population income ,fe model(sem) emat(ivdw) type(both) nolog effects nsim(500) //SEM
est store sem_fdiservicenumber
lrtest sdm_fdiservicenumber sar_fdiservicenumber ///LR test spatial lag
lrtest sdm_fdiservicenumber sem_fdiservicenumber ///LR test spatial error