Re: [R] Random numbers

From: Duncan Murdoch <murdoch_at_stats.uwo.ca>
Date: Thu 22 Dec 2005 - 02:10:51 EST

On 12/21/2005 9:47 AM, Carl wrote:
> Hi All.
> I have R code whose functionality is being replicated within a C+
> program. The outputs are to be compared to validate the conversion
> somewhat - however (as is always the case) I have stuffed my code with
> random number calls.
>
> Random uniform numbers in C+ are being produced using the (Boost)
> mersenne-twister generators (mt11213b & mt19937) - which is the default
> type of generator in R (if I read things correctly). If it was all
> within R I would just set the seed for reproducibility.
>
> Basically - how do I specify in C+ for a set of random uniform numbers
> such that they are the same as from R? I have considered the possibility
> of storing/using the R generated random numbers in the C+ version for
> validation purposes - but there are a lot of them, and that strikes me
> as a generally ugly way of doing things.

I'd say the only reasonable way to do this is to call the R generators rather than trying to duplicate them. R tries hard to keep its generators consistent from version to version, but if you have an independent implementation of the same algorithm, it's going to be very hard to validate that you've really got things exactly identical.

The Writing R Extensions manual tells how to call the R generators from other programs. You can do it without going through interpreted R code, so there shouldn't be much in the way of a performance penalty.

Duncan Murdoch



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Dec 22 02:24:20 2005

This archive was generated by hypermail 2.1.8 : Thu 22 Dec 2005 - 04:28:54 EST