シミュレーションの乱数を描画する正しい方法

Question

シミュレーション用に乱数を生成しようとしています（以下の例では、簡単にするために一様分布を使用しています）。私のモデルで

set.seed(2) 

a <- runif(10000, 1, 999) 

draw <- function(x) { 
    runif(1, 1, 999) 
} 

b <- sapply(1:10000, draw) 
print(c(mean(a), mean(b))) 

、第一の方法のための乱数が[sim_numberを用いたシミュレーション内で参照されるであろう。なぜこれらの2つの方法は、異なる平均値（a: 503.2999, b: 497.5372）同じシード番号と10K回サンプリングするときを生成します]、第2の例ではrunif関数がシミュレーション関数自体の内部に配置されます。それを行う正しい方法がありますか？

Answer 1

答えが完全であるためには、無作為に描画する前にシードを設定する必要があります。