做门限回归找门槛值时,总是出现错误,出现的错误是file C:\Users\VR\AppData\Local\Temp\ST_00000008.tmp could not be opened
同样的命令,在别人电脑上没有这个情况。我觉得是电脑问题,但不知道怎么解决。求大家帮忙
我的原始命令是:
use D:\stata11\ado\personal\work\lner.dta,clear
tsset id year
xtdes
drop if dq == "西藏"
cap drop group
gen group = 2
replace group = 1 if ///
dq=="陕西" | dq=="甘肃" | dq=="宁夏" | dq=="青海" | ///
dq=="新疆" | dq=="重庆" | dq=="广西" | dq=="内蒙" | ///
dq=="四川" | dq=="贵州" | dq=="云南"
replace group = 3 if ///
dq=="北京" | dq=="天津" | dq=="河北" | dq=="福建" | ///
dq=="上海" | dq=="江苏" | dq=="浙江" | dq=="山东" | ///
dq=="广东" | dq=="海南"
replace group = 4 if ///
dq=="黑龙江" | dq=="吉林" | dq=="辽宁"
label define group_lab 1 "西部" 2 "中部" 3 "东部" 4 "东北部"
label value group group_lab
*东北部
drop if (group==1 | group==2 | group==3 )
*门限回归 //将er作为自身的门限变量
*-0-Summary Statistics
set seed 1357
local x "er il il2 is ics tp did lo ia hc"
tabstat `x', s(mean sd min p25 p50 p75 max) format(%6.3f) c(s)
xtdes
*-1-Estimating
xtthres te er il il2 is ics tp did lo ia hc,thres(er) dthres(er) ///
min(3) bs1(500) bs2(500) bs3(300)
*-2-Graphing
xttr_graph
graph work01.wmf
xttr_graph,m(22)
graph work02.wmf
xttr_graph,m(21)
graph work03.wmf
*-3- 估计结果
dis e(rhat21)
dis e(rhat22)
global q1=min(e(rhat21),e(rhat22)) //取出门槛值
global q2=max(e(rhat21),e(rhat22))
dis "$q1" //较小的门限值
dis "$q2" //较大的门限值
dropvars d1 d2 er_x_* //生成虚拟变量
gen d1=(er<=$q1)
gen d3=(er>$q2)
gen er_x_er1=er*d1 //交乘项
gen er_x_er3=er*d3
local x "il il2 is ics tp did lo ia hc" //解释变量
xtreg te `x' er_x_er* er,fe //常规标准误
est store fe
xtreg te `x' er_x_er* er,fe robust //稳健型标准误
est store fe_robust
local m "fe fe_robust"
esttab `m',mtitle(`m') nogap s(r2 r2_w N F) ///
star(* 0.1 ** 0.05 *** 0.01) ///
order(`x' er_x_er1 er er_x_er3)