From simardr@IRO.UMontreal.CA  Thu Jan  9 23:32:04 2003
From: simardr@IRO.UMontreal.CA (simardr@IRO.UMontreal.CA)
Date: Thu Jan  9 23:32:04 2003
Subject: [Rd] RNG.c: unif_rand, MARSAGLIA_MULTICARRY (PR#2437)
Message-ID: <200301092232.XAA12842@pubhealth.ku.dk>

Full_Name: Richard Simard
Version: 
OS: Linux
Submission from: (NULL) (132.204.25.139)


In the random number generator MARSAGLIA_MULTICARRY, your algorithm is
different than the one in the original reference that you give in the R manual:
Marsaglia in his post to the mailing list {\it sci.stat.math} on September 29,
1997.
The last line in the R program has a ^
while Marsaglia's algorithm has a +.
This gives very different random number generators.

The full line in RNG.c is:
    return ((I1 << 16)^(I2 & 0177777)) * i2_32m1; /* in [0,1) */



[Rd] RNG.c: unif_rand, MARSAGLIA_MULTICARRY (PR#2437)
simardr@IRO.UMontreal.CA