forvalue i=30/39{
clear
cd "E:\python_gis_spider\project_working\op_industry"
use CED_OP_ind.dta
set more off
set matsize 6000
xtset new_id year
drop if age>100

keep if cic2id==`i'
opreg ln_y, exit(exit) state(age lnkop) proxy(lnrinvop) cvars(soe export_firm) free(ln_l _Iyear* _Icity*) vce(bootstrap, reps(3))
predict TFP_OP_em
save CED_OP_ind`i'
}

运行结果


(running opreg on estimation sample)
an error occurred when bootstrap executed opreg  r（2000）
这时减去忽略的，我还有4万多观测值啊