用bootstrap做?
sysuse auto, clear
capture program drop bsdelta
program bsdelta, rclass
reg price foreign weight
local reg1_b = _b[foreign]
reg price foreign weight length
local reg2_b = _b[foreign]
return scalar delta = `reg1_b' - `reg2_b'
end
bs delta = r(delta), reps(1000): bsdelta