我試圖使用C++ STD TechnicalReport1擴展,產生以下正態分佈的數字,但是這個代碼(改編自this article):只有 mt19937 eng;
eng.seed(SEED);
normal_distribution<double> dist;
// XXX if I use the one below it exits the for loop
// unifor
當我打電話給rnorm傳遞一個單值作爲均值時,很明顯會發生什麼:從Normal(10,1)生成一個值。 y <- rnorm(20, mean=10, sd=1)
但是,我看到被傳遞到rnorm整體載體的例子(或rcauchy,等等。);在這種情況下,我不確定R機器究竟做了什麼。例如: a = c(10,22,33,44,5,10,30,22,100,45,97)
y <- rnorm(a,