我正在尝试使用r mvrnorm中的内置函数进行一些模拟。
我在同一代码中使用了3次此功能。
是否足以将set.seed放在开头，还是我每次必须在mvrnorm之前写它。
当然，我希望结果可简化，并将此模拟重复10,50和100次。
以下是我的代码
set.seed(30)
mu1<-cbind(matrix(0,nrow=1,ncol=5),t(c(1)))
mu2<-matrix(0,nrow=1,ncol=6)
Sigma<-matrix(1,nrow=p,ncol=6)+diag(6)*0.1
M1=mvrnorm(50,mu1,Sigma)
M2<-mvrnorm(50,mu2,Sigma)
M11<-mvrnorm(n,rep(0,10),diag(10)+matrix(0.3,nrow=10,ncol=10))
Data=data.frame(cbind(rbind(M1,M2),M11))