r - 从均值为 5、标准差为 3 的正态分布中模拟 5000 个大小为 5 的样本

Question

我正在尝试从均值为 5 且标准差为 3 的正态分布中模拟 5000 个大小为 5 的样本。然后我想计算每个样本的均值并制作样本均值的直方图

我当前的代码没有给我错误，但我认为它不对:

nrSamples = 5000
e <- list(mode="vector",length=nrSamples)
for (i in 1:nrSamples) {
e[[i]] <- rnorm(n = 5, mean = 5, sd = 3)
}

sample_means <- matrix(NA, 5000,1)
for (i in 1:5000){
sample_means[i] <- mean(e[[i]])
}

关于如何解决这个问题的任何想法？我对 R 非常陌生!

Answer 1

在这种情况下您不需要列表。过度使用列表是新 R 用户的一个常见错误。

observations <- matrix(rnorm(25000, mean=5, sd=3), 5000, 5)
means <- rowMeans(observations)

现在 表示 是一个包含 5000 个元素的向量。