R: replicate()

apply系関数よりシンプルなreplicate()という関数があるのをいまごろ発見。サンプルに、一様乱数を足し合わせると正規乱数に近くなるというのをやってみた。

> r <- replicate(10000, sum(runif(12, 0, 1)))

[0,1]の一様乱数を12個足し合わせたものを10000個生成。

> r2 <- r - 6

6を引く。

> quantile(r2, c(0.005, 0.025, 0.5, 0.975, 0.995))
        0.5%         2.5%          50%        97.5%        99.5% 
-2.477040054 -1.974674143  0.002924659  1.956245504  2.547839763 

ヒストグラムを書いてみる。

>  hist(r2)

タグ：R