2013-04-15 11 views
8

Ich versuche, scipy.stats.pareto.rvs (b, loc = 0, Maßstab = 1, Größe = 1) mit unterschiedlichen Samen zu generieren.scipy.stats Samen?

In numpy können wir mit numpy.random.seed (seed = 233423) säen.

Gibt es eine Möglichkeit, die von scipy stats generierte Zufallszahl zu seed.

Hinweis: Ich benutze nicht numpy Pareto, weil ich verschiedene Werte für die Skalierung geben möchte.

Antwort

12

scipy.stats verwendet nur numpy.random, um seine Zufallszahlen zu generieren, so dass numpy.random.seed() auch hier funktioniert. Z. B. wird

import numpy as np 
from scipy.stats import pareto 
b = 0.9 
np.random.seed(seed=233423) 
print pareto.rvs(b, loc=0, scale=1, size=5) 
np.random.seed(seed=233423) 
print pareto.rvs(b, loc=0, scale=1, size=5) 

[ 9.7758784 10.78405752 4.19704602 1.19256849 1.02750628] zweimal drucken.

+0

Vielen Dank, es funktioniert .. – ashok

+1

Gibt es * keine * Möglichkeit, scipy.stats ein numpy.random.RandomState-Objekt zu geben? Ich würde den standardmäßigen zufälligen Stream lieber nicht anfassen. –

+2

Anscheinend nicht (wie jetzt), bestätigt von Robert Kern auf scipy-users mailing list! –

Verwandte Themen