私はPythonで乱数のいくつかの特殊な振る舞いを見つけました。特にnumpy.randomというモジュールを使用しました。私は二項分布を期待(興味のある読者のために、これはEhrenfest modelをシミュレート)および大n
正規分布のためであろう大N
の限界でnumpyのnumpy.histogram/random numbersの奇妙な振る舞いは?
n = 50
N = 1000
np.histogram(np.sum(np.random.randint(0, 2, size=(n, N)), axis=0), bins=n+1)[0]
:
は、以下の式を考えます。典型的な出力は、しかし、次のようになります
アレイ([
1、0、0、1、0、2、0、1、0、15、0、
12、0、18、 0、39、0、64、0、62、0、109、
0、97、0、107、0、114、0、102、0、92、0、
55、0、46、0、 35、0、10、0、9、0、4、
0、0、0、3、0、1、1
])
上記の文で、私は説明できませんヒストグラムのゼロの出現 - 私はここで何かを見逃していますか?
ああ、私はそれを逃した - おかげでたくさん! – rammelmueller
完全性のために: 'minlength = n + 1'を' np.bincount'に渡してスパン全体にわたってリストを得るとき、望ましい振る舞いが達成されます。 – rammelmueller