空间杜宾模型stata代码与案例数据更新

clear all

cd "F:\stata14\ado\mydo\agglo"

**调入权重矩阵

use Wqiantaoinver.dta, clear

keep s*

save "W2.dta" ,replace

**矩阵转换为stata应用格式

spatwmat using W2, name(W)standardize

matrix 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=ln(c)

*设定面板数据格式

xtset id year

**变量描述性统计

local x " ii c  p a  t"   //5个变量

          tabstat `x', s(mean 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 effects 

est 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_time

estat ic   

//Spatial Durbin model（双固定，both）

xsmle lnii lnc  lnp lna lnt , wmat(W) model(sdm) fe type(both)  nsim(500) nolog effects   

est store M_both  

estat ic

**结果汇总列表显示

local m "M_sptial M_time M_both"

          esttab `m', mtitle(`m') nogap s(r2 ll aicbic N ) ///

                star(* 0.1 ** 0.05 *** 0.01)b(%6.3f)    

        esttab M_sptial M_time M_both usingpane0000.rtf   //输出结果rtf excel doc

****hausman test

xsmle lnii lnc  lnp lna lnt , wmat(W) model(sdm) fe type(both) nsim(500) nolog 

estimates store sdm_fe

xsmle lnii lnc  lnp lna lnt , wmat(W) model(sdm) re type(both) nsim(500) nolog 

estimates store sdm_re

hausman sdm_fe sdm_re

stata空间计量模型案例数据集链接

链接：http://pan.baidu.com/s/1miFXiQ0 密码：yoo0

学术无界qq答疑群：660826713（遇到任何问题，请进群咨询）

学术无界公众号，欢迎您的关注和传播

欢迎有志于从事学术公益分享的小伙伴加入我们，联系qq：562344572.