【计量】空间杜宾模型代码

1、空间杜宾模型 stata 代码clear allcd "F:stata14adomydoagglo" *调入权重矩阵use Wqiantaoinver.dta, clearkeep s*save "W2.dta" ,replace* 矩阵转换为 stata 应用格式spatwmat using W2, name(W)standardizematrix list W* 调入数据use data11.21.dta, clear*对数生成gen lnii=ln(ii)gen lnp=ln(p)gen lna=ln(a)gen lnt=ln(t)gen lnc=

2、ln(c)*设定面板数据格式xtset id year*变量描述性统计local x " ii c p a t" /5 个变量tabstat 'x', s(mea n sd min p25 p50 p75 max) /format(%6.4f) c(s)*(2)estimation*SDM estimation (个体固定 ind )/Spatial Durbin model (SDM)xsmle lnii lnc lnp lna lnt, wmat(W) model(sdm) fe type(ind) nsim(500) nolog effectsest

3、store M_sptial/不要加入 c三次方lna_2使用第三产业，加入 c_lng交乘项estat ic /AIC BIC test/Spatial Durbin model (时间固定 time )xsmle lnii lnc lnp lna lnt, wmat(W) model(sdm) fe type(time) nsim(500) nolog effects est store M_timeestat ic/Spatial Durbin model (双固定， both )xsmle lnii lnc lnp lna lnt , wmat(W) model(sdm) fe typ

4、e(both) nsim(500) nolog effects est store M_bothestat ic*结果汇总列表显示local m "M_sptial M_time M_both"esttab 'm', mtitle('m') nogap s(r2 ll aic bic N ) /star(* 0.1 * 0.05 * 0.01) b(%6.3f)esttab M_sptial M_time M_both using pane0000.rtf / 输出结果 rtf excel doc*hausman testxsmle lnii lnc lnp lna lnt , wmat(W) model(sdm) fe type(both) nsim(500) nolog estimates store sdm_fexsmle lnii lnc lnp lna lnt , wmat(W) model(sdm) re