统计思维：将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)