RNG R/Splus compatibility

Paul Gilbert pgilbert@bank-banque-canada.ca
Thu, 06 May 1999 10:15:44 -0400

Here is a correction to one of the functions I posted yesterday, which did not
work correctly when kind was a list. With a null argument this function returns
current RNG and seed information. Otherwise, it sets the values and returns the
new settings. Thus, in a function one might

  old <- set.RNG()
  effective <- set.RNG(new)   #(= new)

Another approach would be that used by par() (i.e. return the old settings) so
one would

  old <- set.RNG(new)
  effective <- set.RNG()  #(= new)

Does anyone have thoughts on this?

Paul Gilbert

    set.RNG <- function(kind=NULL, seed=NULL, rnorm=NULL)
      {# with a null argument this also serves as get.RNG so it is possible to
       # old <- set.RNG()
       # on.exit(set.RNG(old))
       # set.RNG("default")  can affect both kind and seed, but not rnorm
       if (is.list(kind))
         list(kind=RNGkind(kind$kind), seed=set.seed(kind$seed),
              rnorm= set.RNORMtransform(kind$rnorm) )
         list(kind=RNGkind(kind), seed=set.seed(seed),
              rnorm=set.RNORMtransform(rnorm) )

r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch