Skip to content
Back to formatted view

Raw Message

Message-ID: <AA89BFA18D20A8469B6E05BBC55E4471713F3034@EXCH10-MBX-06.home.ku.edu>
Date: 2011-11-05T00:38:20Z
From: Md Desa, Zairul Nor Deana Binti
Subject: set seed for random draws

Hello, all!
I need help on these two problems:

1) If I want to randomly draw numbers from standard normal (or other distributions) in loops e.g.:
 ty=0; ks=0
for (i in 1:5) {
        set.seed(14537+i)
        k<-rnorm(1)
        ks[i]<-.3*k+.9
        if (ty==0) {
            while ((ks<.2)||(ks>3)) {
            #set.seed(13237+i*100)
            k<-rnorm(1)
            ks[i]-.3*k+.9 }
        }
     }
....
....
....
     }

Question: Here I draw initial a, then if the drawn initial a satisfied 2 conditions I redraw a. I set.seed(13237) in the first draw of a, should I set.seed() in the redraw part?

2) I also have more loops after this i loop that also draw from normal(0,1). I want to randomly draws from normal(0,1) for loop j (inside loop j I draw another random numbers from N(0,1))
My question: Should I or shouldn't I set seed again and again for each loop? Why or why not.

I guess this problem concerned about setting seed as I want to have different number for each i.

Thanks!

Deana