Question

我已经使用rnorm创建了一个均值设置为1000的代码，现在我想创建一个由10000个元素组成的矢量，每个元素都是我创建的这种正态分布的均值。我该怎么办？

我尝试使用

a = replicate(10000, (mean(rnorm(4, mean=1000, sd=4))))

但是随后我得到了大约1004的值，这使我相信这不是我得到的分布的均值。有什么我想念的吗？

Answer 1

如评论中所指出：您在样本中仅使用四个随机数的事实可能导致经验均值与真实均值（即1000）相差甚远。

随着每次迭代中随机数的增加，经验均值将接近真实均值（当谈论 unbiased 估计量时，这正是这个意思）。

例如参见以下

set.seed(123)
replicate(5, (mean(rnorm(4, mean=1000, sd=4))))
[1] 1000.8386 1001.0402 1000.4514 1001.7425  998.7598
replicate(5, (mean(rnorm(1000, mean=1000, sd=4))))
[1] 1000.0734 1000.1480  999.9374  999.9468  999.8768
replicate(5, (mean(rnorm(100000, mean=1000, sd=4))))
[1]  999.9981  999.9939 1000.0004 1000.0001  999.9966

如何使用向量中的每个值作为平均值重复R中的rmrm？

1 个答案: