stata命令见下(生成变量的过程略去了,因为xtivreg能执行,应该变量生成不会有问题吧):
- xtset prvcode sgnyea
- gen lngdp = ln(gdp)
- gen lnpergdp =ln(pergdp)
- gen lnk= ln(k)
- gen lnepprv =ln(epprv)
- gen lnpltcha =ln(pltcha)
- gen pltcha1 =L.pltcha
- gen pltcha2 =L2.pltcha
- gen ipltcha1 =ln(pltcha1)
- gen ipltcha2 =ln(pltcha2)
- gen K_L =k/epprv
- gen lnK_L =ln(K_L)
- gen lnK_L_squ= lnK_L * lnK_L
- gen K_P =k/pltcha
- . gen lnK_P =ln(K_P)
- . gen lnK_P_squ= lnK_P * lnK_P
- xtset prvcode sgnyea
- xtreg lnpergdp lnk lnepprv lnpltcha lnK_L_squ lnK_P_squ i.sgnyea, fe
- est store fe
- xtreg lngdp lnk lnepprv lnpltcha lnK_L_squ lnK_P_squ i.sgnyea, re
- est store re
- hausman fe re, constant sigmamore
- sum sgnyea prvcode pergdp ifa k pltcha epprv
- sort prvcode
- by prvcode: reg lnpergdp lnk lnepprv lnpltcha lnK_L_squ lnK_P_squ
- by prvcode: reg lnpergdp lnk lnepprv lnpltcha lnK_L_squ lnK_P_squ
- xtivreg lnpergdp lnk lnepprv lnK_L_squ lnK_P_squ i.sgnyea (lnpltcha=ipltcha1 ipltcha2), fe
- est store fe_iv
- dmexogxt
- local models "fe re fe_iv"
- esttab `models', b(%6.3f) se(%6.3f) mtitle(`models') sca(r2 r2_w)
- xtivreg2 lnpergdp lnk lnepprv lnK_L_squ lnK_P_squ i.sgnyea (lnpltcha = ipltcha1 ipltcha2), fe small
- xtoverid
三个命令执行时stata窗口中显示的内容:
. dmexogxt
2003b: operator invalid
r(198);
. xtivreg2 lnpergdp lnk lnepprv lnK_L_squ lnK_P_squ i.sgnyea (lnpltcha = ipltcha1 ipltcha2), fe
factor variables not allowed
r(101);
. xtoverid
2003b: operator invalid
r(198);
stata初学者,还望高人指点


雷达卡







京公网安备 11010802022788号







