我想在9对不同参数下，每对参数下取服从正态分布的随机数100个，构成一个9*100的矩阵，刚好每一行是一对参数下的100个随机数，请问怎么办？我写的循环有问题，一直返回的是最后一组参数下的随机数。u<-c(4,5,6)
sigma<-c(1,2,3)
m<-matrix(data=NA,nrow=9,ncol=100)
for (i in 1:3)
  for (j in 1:3)
    for (k in 1:9)
      m[k,]<-c(rnorm(100,u,sigma[j]))