[R] RNG Cycle and Duplication

Duncan Murdoch murdoch at stats.uwo.ca
Thu Aug 14 23:04:14 CEST 2008


Shengqiao Li wrote:
> 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.
>   

This is no bug.  Cycle length and the set of possible values are 
different concepts. 

And please don't cross-post.

Duncan Murdoch
> 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
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



More information about the R-help mailing list