On 3 October 2013 at 15:39, Renaud Gaujoux wrote:
| Last time I looked at RNG.c (or filename like that) in R source code,
| you could not get anything else than Marsaglia.

The Marsaglia generator to seed all distributions.  

So if you want the 20+ distributions rlnorm(), rweibull(), ... you are good.
If you to comparese different U(0,1) generators you are not. I though you
wanted the former--sorry if I misunderstood.

| That's the thing, set.seed is the only exported function that can
| touch the RNG and only enables to set the 2-length seed for Marsaglia,
| and the kind cannot be changed, see (which also interestingly tells
| you something about the order of evaluation of printf arguments...):
| 
| Rscript -e "RNGkind('Marsaglia'); .Random.seed[2:3] <- c(123L, 456L); 
rnorm(2)"
| [1] -0.2934974 -0.3343770

Nice. There is the (partial) equivalence (subject to ordering).

Dirk 

-- 
Dirk Eddelbuettel | e...@debian.org | http://dirk.eddelbuettel.com
_______________________________________________
Rcpp-devel mailing list
Rcpp-devel@lists.r-forge.r-project.org
https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel

Reply via email to