本人在STATA中跑安慰剂检验,代码从别人那里借鉴来之后,自己却死活用不了,跑一部就出现错误 do "C:\Users\RRRRRR~1\AppData\Local\Temp\STD00000000.tmp"
. forvalue i=1/1000{
2.
unexpected end of file
r(612);
本人不知道到底是什么原因。希望会的大神能够帮帮我,我必将重金酬谢(价格可私聊)
先将代码复制如下:
forvalue i=1/1000{
sysuse aq,clear
g obs_id=_n
gen random_digit=runiform()
sort random_digit
g random_id=_n
preserve
keep random_id exp
rename exp random_exp
rename random_id id
label var id
save random_exp,replace
restore
drop random_digit random_id exp
rename obs_id id
label var id
save rawdata,replace
use rawdata, clear
merge 1:1 id using random_exp,nogen
reg risk random_exp gender age top1 zcfzl size roe board ddbl
g _b_random_exp= _b[random_exp] //提取x的回归系数
g _se_random_exp= _se[random_exp] //提取x的标准误
keep _b_random_exp _se_random_exp
duplicates drop _b_random_exp, force
save placebo`i', replace
use placebo1, clear
forvalue i=2/1000{
append using placebo`i' //纵向合并1000次回归的系数及标准误
}
gen tvalue= _b_random_exp/ _se_random_exp
kdensity tvalue, xtitle("t值") ytitle("分布") saving(placebo_test)
*-删除临时文件
forvalue i=1/1000{
erase placebo`i'.dta
}