使用MASS包中的mvrnorm

我需要使用

seed（12346）

生成一个多变量正态分布的随机样本，包含100列和5000行

到目前为止，我已经做到了：

set.seed（12346）
Preg1您可以使用
mu=rep（01100）
。
rep
功能用于重复值

set.seed(12346)
ncol = 100
Preg1<-data.frame(mvrnorm(n = 5000, mu = rep(0, ncol), Sigma = diag(ncol)))

dim(Preg1)
# [1] 5000  100


在这种特殊情况下，因为
Sigma
是一个单位矩阵，所以每一列实际上是独立的。因此，这相当于独立生成每个列（甚至每个绘图），我们可以采用以下任一方法：

x = replicate(n = ncol, rnorm(5000))
dim(x)
# [1] 5000  100

z = matrix(rnorm(5000 * ncol), ncol = ncol)
dim(z)
# [1] 5000  100

您可以使用
mu=rep（01100）
。
rep
功能用于重复值

set.seed(12346)
ncol = 100
Preg1<-data.frame(mvrnorm(n = 5000, mu = rep(0, ncol), Sigma = diag(ncol)))

dim(Preg1)
# [1] 5000  100


在这种特殊情况下，因为
Sigma
是一个单位矩阵，所以每一列实际上是独立的。因此，这相当于独立生成每个列（甚至每个绘图），我们可以采用以下任一方法：

x = replicate(n = ncol, rnorm(5000))
dim(x)
# [1] 5000  100

z = matrix(rnorm(5000 * ncol), ncol = ncol)
dim(z)
# [1] 5000  100

也许
mu=rep（01100）
？也许
mu=rep（01100）
？