At 17:49 14/02/00 +0100, you wrote:
>Dear BUGS users.
>
>Does anybody know how to get information about internal specifications in
>WinBUGS (or BUGS)?
>
>I have had this problem before, but right now I am looking for information
>on the random number generator WinBUGS uses. What kind is it? What is the
>period?
>
>Thanks in advance
>
>Anders.
>
>
>
Hi Anders,
(Win)BUGS uses a linear congurence generator, its period is 2^32
Andrew
PROCEDURE Rand* (): REAL;
CONST
a = 16807; m = 2147483647; q = m DIV a; r = m MOD a; realm = 2147483647.0;
VAR
gamma: INTEGER;
BEGIN
gamma := a * (iRand MOD q) - r * (iRand DIV q);
IF gamma > 0 THEN iRand := gamma ELSE iRand := gamma + m END;
RETURN iRand / realm
END Rand;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|