统计思维：将for循环转换为重复函数并使用R存储输出

Question

我试图做一个简单的模拟，以便找到样本平均值的样本分布。 我希望能够严格使用repeat函数并完全删除for循环来编写它。 到目前为止，这是代码（对不起，对于空格，由于某些原因，stackoverflow对其进行了奇怪的格式化）：

X.samp <- sample(pop.1$height,100)
X.bar <- mean(X.samp)
X.bar <- rep(0,10^5)
for (i in 1:10^5) {
   X.samp <- sample(pop.1$height,100)
   X.bar[i] <- mean(X.samp)
}
hist(X.bar)

Answer 1

像这样吗

 dat <- rnorm(1000) # some data
 res <- replicate(1000, mean(sample(dat, 100)))
 hist(res)