search for: seed1

....Random.seed should be printed: test<-function(runs,seed){ .Random.seed<-seed for (i in 1:runs) { print(i) print(runif(1,0,1)) print(.Random.seed)} return(.Random.seed} Consider the following input/output: >RNGkind(kind="Marsaglia-Multicarry") > set.seed(20391) > seed1<-.Random.seed > seed1 [1] 401 -1607331462 -462081869 > test(2,seed1) [1] 1 [1] 0.4188851 [1] 401 -1607331462 -462081869 [1] 2 [1] 0.7713649 [1] 401 -1607331462 -462081869 [1] 401 -1607331462 -462081869 > seed1 [1] 401 -1607331462 -462081869 &...

Hello, Does anyone knows how to deal with zero subscript in R. I have this code: for (i in 1:nitems){ + for (j in 1:ncat-1) { + draw<-matrix(rnorm(nitems*(ncat-1),seed1,seed2),nitems,(ncat-1)) + d<-( sigma_d*draw ) + mu_d + draw<-matrix(rtnorm((nitems*(ncat-1)),mean = seed1, sd = seed2, lower = .1, upper = 1.5),nitems,(ncat-1)) + d<-(sigma_d*draw) + mu_d + write.matrix(cbind(b.i0,d), file = "F:/KU/MIRT group/MIMIC...

...out the SEXP overhead of course). But the seeds are not set accordingly and segmentation faults are not reproducible. More precisely, I need to replace the set_seed line in the following snippet with something else. #include <Rmath.h> int main(int argc,char* argv[]) { const unsigned int SEED1= 12345, SEED2= 67890; set_seed(SEED1, SEED2); ... } Thanks in advance, Oliver [[alternative HTML version deleted]]

...ch to create replicable simulations. It blends ideas about seed management from John M. Chambers Software for Data Analysis (2008) with ideas from the snowFT package by Hana Sevcikova and Tony R. Rossini. Here's my proposal. 1. Run a preliminary program to generate an array of seeds run1: seed1.1 seed1.2 seed1.3 run2: seed2.1 seed2.2 seed2.3 run3: seed3.1 seed3.2 seed3.3 ... ... ... run1000 seed1000.1 seed1000.2 seed1000.3 This example provides 3 separate streams of random numbers within each run. Because we will use the L'Ecuyer "many separate st...