我想用矩估计法来估计正态分布参数。若已知总体均值为2，方差为4，我如何用R生成1000组rnorm(100,2,4)的样本，再利用这1000组样本运用矩估计的方法估计出样本的参数？最后运用MSE比较总体参数与样本参数的偏差度。请问如何生成1000组？我的一段小程序是：
for(i in 1:1000){a=paste('x',i,sep='')
+ b=rnorm(100,2,4);assign(a,b)}
但是我只能敲一个x1生成100个数、敲一个x2生成100个数、...我如何将这1000组数放到一个数组里呢？