global ylist y
global xlist x1 x2 x3 x4 x5 x6 x7 x8
xtset id year
**SDM模型
xsmle $ylist $xlist,wmat(w) model(sdm) re vce(cluster id) nolog
est store m1
xsmle $ylist $xlist,wmat(w) model(sdm) fe type(ind) vce(cluster id) nolog
est store m2
xsmle $ylist $xlist,wmat(w) model(sdm) fe type(time) vce(cluster id) nolog
est store m3
xsmle $ylist $xlist,wmat(w) model(sdm) fe type(both) vce(cluster id) nolog
est store m4
xsmle $ylist $xlist,wmat(w) model(sdm) noeffects nolog
est store m5
**导出回归结果
esttab m1 m2 m3 m4 m5 using regn1000.rtf,r2 se star(* 0.1 ** 0.05 *** 0.01)replace nogap
**豪斯曼检验
xsmle $ylist $xlist,wmat(W) model(sdm) hausman nolog