打印多个随机排列

Question

我正在尝试进行多个排列。从代码中：

# generate random Gaussian values
from numpy.random import seed
from numpy.random import randn
# seed random number generator
seed(1)
# generate some Gaussian values
values = randn(100)
print(values)

但是现在我想生成，例如，（值）20个排列。 使用代码：

import numpy as np
import random
from itertools import permutations
result = np.random.permutation(values)
print(result)

我只能观察到一个排列（或“手动”获得其他置换）。我希望我有很多排列（20或更多），因此自动计算每个置换量（来自值）的Durbin-Watson统计量。

from statsmodels.stats.stattools import durbin_watson
sm.stats.durbin_watson(np.random.permutation(values))

我该怎么办？

Answer 1

要从某些集合中获取20个排列，请Intial itertools.permutations iterator，然后使用next（）进行第一个20：

import numpy as np
import itertools as it

x = np.random.random(100)  # 100 random floats
p = it.permutations(x)  # an iterator which produces permutations (DON'T TRY TO CALCULATE ALL OF THEM)
first_twenty_permutations = [next(p) for _ in range(20)]

当然，这些不会是随机的排列（即它们以有组织的方式计算，尝试it.permutations（“ abcdef”），您会明白我的意思）。如果需要随机排列，则可以以相同的方式使用np.random.permunt.：

[np.random.permutation(x) for _ in range(20)]

然后计算Durbin Watson统计信息：

permutations = np.array([np.random.permutation(x) for _ in range(20)])
np.apply_along_axis(durbin_watson, axis=1, arr=permutations)