2013-04-15 54 views
8

我想用不同的種子生成scipy.stats.pareto.rvs(b,loc = 0,scale = 1,size = 1)。scipy.stats種子?

在numpy中,我們可以使用numpy.random.seed(seed = 233423)種子。

有沒有什麼辦法給scipy stats生成的隨機數播種。

注意:我不使用numpy pareto,因爲我想爲縮放賦予不同的值。

回答

12

scipy.stats只是使用numpy.random來生成它的隨機數,所以numpy.random.seed()也可以在這裏使用。例如,

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]兩次。

+0

謝謝,它的工作.. – ashok

+1

有沒有*沒有*的方式給scipy.stats一個numpy.random.RandomState對象?我真的不想碰默認的隨機流。 –

+2

顯然不是(好像現在),由Robert Kern在scipy用戶郵件列表中確認! –