Rで確率・統計の基礎２ 〜 中心極限定理

前回「Rで確率・統計の基礎１　〜　母集団と標本」の続き。

前回とは異なり、常に10個の抽出を行って平均を取得し、それらをリストmean_listに入れる。
そのヒストグラムを描き、Xの平均を縦線で追加で描く。

> mean_list<-list()
> for(i in 1:10000){
+  mean_list[[i]] <- mean(sample(X, 10, replace = TRUE))
+ }
> hist(unlist(mean_list), breaks=500, xlab = "Mean of 10 samples from X", main = "Conv of sample dist", cex.main = 0.8)
> abline(v = mean(X), lwd=3, col="white", lty=2)
> 

次は、50%の確率で0か1を選択したデータで同じことを行う。

> population<-sample(c(0,1), 100000, replace=TRUE)
> hist(population, main="Non-normal", cex.main=0.8)
> abline(v=mean(population), lwd=3, lty=3)

> mean_list <- list()
> for(i in 1:10000){
+  mean_list[[i]] <- mean(sample(population, 10, replace = TRUE))
+ }
> hist(unlist(mean_list), xlab = "Avg of 10 samples", main = "Dist of Avg", cex.main = 0.8)
> abline(v=mean(population), lwd=3)