[R] RNG Cycle and Duplication

Shengqiao Li shli at stat.wvu.edu
Thu Aug 14 18:34:43 CEST 2008


Hello all,

I am generating large samples of random numbers. The RNG help page says: 
"All the supplied uniform generators return 32-bit integer values that are 
converted to doubles, so they take at most 2^32 distinct values and long 
runs will return duplicated values." But I find that the cycles are not 
the same as the 32-bit integer.

My test indicated that the cycles for Knuth's methods were 2^30 while 
Wichmann-Hill's cycle was larger than 2^32! No numbers were duplicated in 
10M numbers generated by runif using Wichmann-Hill. The other three 
methods had cycle length of 2^32.

So, anybody can explain this? And any improvement to the 
implementation can be made to increase the cycle length like the 
Wichmann-Hill method?


========================================
Shengqiao Li

Research Associate
The Department of Statistics
West Virginia University
Morgantown, WV 26506-6330



More information about the R-help mailing list