[SciPy-user] Generating random variables in a joint normal distribution?
Parvel Gu
parvel.gu@gmail....
Thu Nov 1 01:08:38 CDT 2007
Hi,
Thanks so much for the reply. I could go furthur in this problem.
And here is the following issue:
In my understanding, the continues calls to random.multivariate_normal
would get a serials of random vairables which are all following the
joint normal distribution. Then what about how to reset the random
generating in each iteration?
Assuming that I have to simulate something for 50 iterations then to
get the average value. In each iteration I need a serial of random
pairs which follows the joint normal distrubition. The interations are
expected to be independent to each other.
Thus:
In iteration 1:
serial = []
in sub loop
serial.append(random.multivariate_normal(m, cov))
...
And in iteration 2, I want a fresh serial which should not be affected
by the previous one s. Is there any problem to empty the serial
then still call random.multivariate_normal(m, cov) ?
In my understanding if there is no refresh or something performed to
the internal random generator, the random.multivariate_normal would
continue to generate variables in the context of the previous ones.
Br,
Parvel
More information about the SciPy-user
mailing list