set obs 100
set seed 500//这两个有什么用的?
prog monte
gen x1=5*runiform()
gen x2=rnormal(5,2)
gen err=invnormal(uniform())
gen y=5+x1+0.5*x2+err
reg y x1 x2
end
simulate _b _cons, reps(1000): monte//
cap pr drop monte
pr monte clear
set ob 100
g x1=5*runiform()
g x2=rnormal(5, 2)
g err=rnormal()
g y=5+x1+0.5*x2+err
reg y x1 x2
end
simulate, r(1000) nodots seed(500): monte