search for: seed2

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-DIF/R...

...verhead 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]]

...tive, which is very mind bending. The regression in out.ok is o. k. It checks by hand. For a more complete explanation (if more is wanted), including the printout from these summary commands on my machine and the check of out.ok "by hand", see http://www.stat.umn.edu/geyer/5931/mle/seed2.Rnw http://www.stat.umn.edu/geyer/5931/mle/seed2.pdf

...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 streams" approach, we are confide...