结果不一样 哭死了 咋办 要求结果一样 能复现 不然就是造假 第一张是刚做的  第二张是之前做的<br>
xtset city  year<br>
第一步：描述性分析<br>
ssc install asdoc<br>
asdoc sum y x1 x2 x3 x4 x5 <br>
图形分析<br>
graph twoway (scatter y x1)<br>
graph twoway (scatter y x1)(lfit y x1)<br>
graph twoway (scatter y x1)(lfit y x1)(qfit y x1)<br>
第二步相关性分析<br>
pwcorr y x1 x2 x3 x4 x5 <br>
第三步多重共线性分析<br>
regress y x1 x2 x3 x4 x5<br>
estat vif<br>
第四步<br>
F检验（判断是混合效应还是固定效应）<br>
xtreg y x1 x2 x3 x4 x5 ,fe<br>
LM检验（判断是混合效应还是随机效应）<br>
xtreg y x1 x2 x3 x4 x5,re<br>
xttest0<br>
第五步<br>
Hausman检验（判断固定效应还是随机效应）<br>
xtreg y x1 x2 x3 x4 x5,fe<br>
est store fe<br>
xtreg y x1 x2 x3 x4 x5 ,re<br>
est store re<br>
hausman fe re<br>
第六步基准回归<br>
选择固定模型<br>
xtreg y x1,fe r<br>
xtreg y x1 x2,fe r<br>
xtreg y x1 x2 x3,fe r<br>
xtreg y x1 x2 x3 x4,fe r<br>
xtreg y x1 x2 x3 x4 x5,fe r<br>
输出回归数据<br>
ssc install outreg2<br>
xtreg y x1,fe r<br>
outreg2 using ab,word replace<br>
xtreg y x1 x2,fe r<br>
outreg2 using ab,word<br>
xtreg y x1 x2 x3 ,fe r<br>
outreg2 using ab,word<br>
xtreg y x1 x2 x3 x4 ,fe r<br>
outreg2 using ab,word<br>
xtreg y x1 x2 x3 x4 x5 ,fe r<br>
outreg2 using ab,word<br>
第七步稳健性检验<br>
1.ols回归 <br>
reg y x1 x2 x3 x4 x5,robust （最小二乘法）<br>
est store reg<br>
esttab fe reg using reg2.rtf, replace b(%6.3f) se(%6.3f) se ar2(3) star(* 0.1 ** 0.05 *** 0.01) compress nogap mtitles("model1""model2" ) title("Table1")<br>
2.系统GMM<br>
ssc install xtabond2<br>
xtabond2 y x1 x2 x3 x4 x5,gmmstyle(L.y x5) ivstyle(x2) twostep robust<br>
输出结果<br>
xtreg y x1 x2 x3 x4 x5 ,fe r<br>
outreg2 using 123,word replace<br>
xtabond2 y x1 x2 x3 x4 x5,gmmstyle(L.y x5) ivstyle(x2) twostep robust<br>
outreg2 using 123,word<br>
       
第八步门槛效应/门限效应检验先从三门槛检验开始<br>
核心  科技创新  X1<br>
xthreg y,rx(x1) qx(x1) thnum(3) trim(0.01 0.01 0.01) grid(300) bs(300 300 300)<br>
双门槛<br>
xthreg y,rx(x1) qx(x1) thnum(2) trim(0.01 0.01 ) grid(300) bs(300 300 )<br>
单门槛<br>
xthreg y,rx(x1) qx(x1) thnum(1) trim(0.01 ) grid(300) bs(300 )<br>
固定资产投资<br>
三门槛<br>
xthreg y,rx(x1) qx(x2) thnum(3) trim(0.01 0.01 0.01) grid(300) bs(300 300 300)<br>
双门槛<br>
xthreg y,rx(x1) qx(x2) thnum(2) trim(0.01 0.01 ) grid(300) bs(300 300)<br>
单门槛<br>
xthreg y,rx(x1) qx(x2) thnum(1) trim(0.01 ) grid(300) bs(300)<br>
画出固定资产投资单门槛LR图<br>
_matplot e(LR), columns(1 2) yline(7.3523, lpattern(dash)) connect(direct) recast(line) ytitle("LR Statistics") xtitle("First Threshold") name(LR)<br>
输出门槛回归的结果<br>
xthreg y,rx(x1) qx(x2) thnum(1) trim(0.01 ) grid(300) bs(300)<br>
outreg2 using abd,word<br>
第九步中介效应<br>
中介变量x5<br>
xtreg y x1 x2 x3 x4,fe<br>
xtreg x5 x1 x2 x3 x4,fe<br>
xtreg y x5 x1 x2 x3 x4,fe<br>
输出<br>
xtreg y x1 x2 x3 x5,fe<br>
outreg2 using a,word replace<br>
xtreg x4 x1 x2 x3 x5,fe<br>
outreg2 using a,word<br>
xtreg y x4 x1 x2 x3  x5,fe<br>
outreg2 using a,word<br>
xtreg y x1 x2 x3 x4,fe<br>
outreg2 using a,word<br>
xtreg x5 x1 x2 x3 x4,fe