ywh19860616 发表于 2014-6-1 15:00 
程序给出来哈
sysuse np.dta,clear
tsset var1
reg np gf l.np l.gf
predict res,residual
keep res
save res1_data,replace
set matsize 1000
local reps=1000
mat B = J(`reps',1,.)
quietly{
forvalues i = 1(1)`reps'{
use res1_data,clear
bsample
gen id =_n
sort id
tempfile e_data
save e_data,replace
sysuse np,clear
regress np gf
predict np_hat
keep np_hat gf
gen id= _n
sort id
merge id using "`e_data'"
gen np_bs = np_hat + res
reg np_bs weight
mat b =e(b)
mat B[`i',1]=b[1,1]
}
}