reproducing Box-Muller numbers
Brian
I'll make a call to set.seed clear the saved value. That's all we intend to be reproducible, setting the state by set.seed.
I'm not sure what you mean by this. Are you suggesting that setting .Random.seed directly is not intended to give reproducible results? I presume not, but if so, then set.seed needs better control over the seed (i.e. vector seeds).
Also, could "user-supplied" be added as an option for normal.kind in RNGkind.
I'm sure the Box-Muller in R is better than my own attempt, but I would like
to
be able to reproduce results I've obtained in the past.
Not very easily, and we are in feature freeze for 1.0.0 now, so not before 1.1.0. It would be easy to patch in your own version, but the difficulty is to add all the protection needed for user-supplied code.
When are you expecting 1.1.0? I can easily patch the version used by my library, but for obvious reasons I would prefer not to post a library that does this. A very large number of my code tests do not work if I cannot reproduce previous rnorm random numbers. Paul Gilbert -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._