Re: [R] Initialising Mersenne-Twister with one integer

From: Prof Brian Ripley <ripley_at_stats.ox.ac.uk>
Date: Mon 25 Sep 2006 - 08:48:16 GMT

On Mon, 25 Sep 2006, Gad Abraham wrote:

> Hi,
>
> It seems to me that the Mersenne-Twister PRNG can be initialised using
> one integer instead of 624 integers, since inside RNG.c code there's a
> function defined as MT_sgenrand(Int32).
>
> How do I actually set this seed within R?

set.seed(), on the help page for ?.Random.seed.

> I've tried:
>
> > .Random.seed <- c(3, 1)
> > runif(1)
> Error in runif(1) : .Random.seed has wrong length

>From the help page

      '.Random.seed' is an integer vector, containing the random number
      generator (RNG) *state* for random number generation in R.  It can
      be saved and restored, but should not be altered by the user.

> In addition, is '3' actually the correct rng.kind for the Mersenne-Twister?
>
> I'm using R version 2.2.1, 2005-12-20 on Ubuntu Dapper Linux 686.

Not current, but I suspect the help page is the same in that version.

-- 
Brian D. Ripley,                  ripley@stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________
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
and provide commented, minimal, self-contained, reproducible code.
Received on Mon Sep 25 19:05:57 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 27 Sep 2006 - 00:30:07 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.