求助！在R的循环里怎么对矩阵进行编号呢？？ [推广有奖]

R刚学了几天，在这有个弱弱的问题请教一下各位大侠
   写了一个循环，在这个循环里要再一个条件下随机抽取不同组数据放到矩阵里

   每完成一个循环都要重新抽取数据，对数据进行分析处理得到一个返回向量，简单点说就像下面那段，随机抽一组正态分布数据，LS得到beta的估计
然后重复这个过程，得到10组beta的估计值……不过最后那个beta[i]直接对矩阵里的数值进行赋值，不是矩阵的循环编号，所以总出问题……
哪个高手帮帮忙解决一下，如果是算法问题，那怎么才能重复得到一组beta估计的过程呢？
for(i in 1:10){
Sigma <- matrix(c(1,0,0,1),2,2)
a<-mvrnorm(n=10, rep(0,2),Sigma)
x1<-a[,1]
x2<-a[,2]
e<-rnorm(1000,mean=0,sd=1)   
b<-matrix(c(1,1),2,1)         
y<-a%*%b+e
r1<-lm(y~x1+x2)               

m1<-r1$coefficients
m2<-matrix(c(0,1,0,0,0,1),3,2)
beta[i]<-m1%*%m2               

}

论坛新人……所以没有报酬啦……恳请热心的高手~~