R で、確率分布から乱数をサンプリングする runif や rnorm の引数はベクトル化されている。
> runif(10, min=c(0, 100), max=c(10, 110)) [1] 7.046310 103.526202 1.910480 103.118722 3.107876 102.799656 1.802613 109.406715 1.984414 [10] 102.594984 > rnorm(10, mean=c(0, 100)) [1] -0.7162998 98.9900047 0.2466831 99.2812851 0.8309554 99.0108599 -0.6513625 98.9389296 [9] 1.0515469 100.2959684
これは便利。