search for: equidistribut

...7, 477-484. as a reference, but this algorithm does not discuss the generation of random values. help(RNG) indicates that the "Mersenne-Twister" is the default random number generator with reference: Matsumoto, M. and Nishimura, T. (1998) Mersenne Twister: A 623-dimensionally equidistributed uniform pseudo-random number generator, ACM Transactions on Modeling and Computer Simulation, 8, 3-30. I looked at the paper, but essentially I have no expertise in assessing whether this random number generator is "good" (which is probably a tricky concept in the first place...

...equal to zero in gsl_rng_set reproduces this. Later versions switched to 5489 as the default seed, you can choose this explicitly via gsl_rng_set instead if you require it. For more information see, Makoto Matsumoto and Takuji Nishimura, ?Mersenne Twister: A 623-dimensionally equidistributed uniform pseudorandom number generator?. ACM Transactions on Modeling and Computer Simulation, Vol. 8, No. 1 (Jan. 1998), Pages 3?30 The generator gsl_rng_mt19937 uses the second revision of the seeding procedure published by the two authors above in 2002. The original seed...

I seem to be able to use expected values that are decimal (e.g., 1.33) when using chisq.test but not when using fisher.test. This happens when using an array/matrix as input. Fisher.test returns: Error in sprintf(gettext(fmt, domain = domain), ...) : invalid format '%d'; use format %s for character objects. Thus, it appears fisher.test is looking for integers only. I tried putting the

Dear list useRs, I have to generate a random set of coordinates (x,y) in [-1 ; 1]^2 for say, N points. At each of these points is drawn a circle (later on, an ellipse) of random size, as in: > N <- 100 > > positions <- matrix(rnorm(2 * N, mean = 0 , sd= 0.5), nrow=N) > sizes<-rnorm(N, mean = 0 , sd= 1) > plot(positions,type="p",cex=sizes) My problem is to

A few simple questions from a novice R user. (I am running R version 0.62.3 Beta (Sept 8, 1998) under Windows NT.) 1) How do I time the execution of a function/program in R for Windows? Is there the equivalent of a dos.time function? 2) Can anyone send me details on the random number generator used by R (period, etc.). Also, why is it different from the SuperDuper RNG in S? 3) Will there

Whomever, I recently sent the "bug report" below toR-core at r-project.org and have just been asked to instead submit it to you. Although I am basically not an R user, I have installed version 3.3.1 and am also the author of a statistics program written in Visual Basic that contains a component which correctly implements the Mersenne Twister (MT) algorithm. I believe that it is