蒙特卡洛仿真实验,为什么只有第一次的参数的估计值?
4 个回复 - 3987 次查看
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 ...
2012-4-22 01:28 - nhtxy - Stata专版