clear
set mem 10m
forvalues i=1(1)1000 {
set obs 500
gen y1=.
gen y2=.
replace y1=0
replace y2=0
matrix m = (0,0)
matrix sd =(sqrt(2),sqrt(2))
drawnorm e1 e2, n(500) means(m) sds(sd)
replace y1=.6*y1[_n-1]+.4*y2[_n-1]+e1 if _n>2
replace y2=.3*y1[_n-1]+.7*y2[_n-1]+e2 if _n>2
regress y1 y2
matrix c=e(b)
svmat c, names(beta)
}
这是我现在写的一个基本程序，里面有些命令实在不知道怎么弄了。计划是要做1000次，得出1000个beta回归值再对这些回归值作样本分布的分析，但是我写这个出来的结果只有最新得到的一个数据。应该要怎么改啊？
拜谢先！