統計思維：將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)