Python:Numpy 标准差错误

Question

这是一个简单的测试

import numpy as np
data = np.array([-1,0,1])
print data.std()

>> 0.816496580928

我不明白这个结果是如何产生的？很明显:

( (1^0.5 + 1^0.5 + 0^0.5)/(3-1) )^0.5 = 1

在 matlab 中它给了我 std([-1,0,1]) = 1。你能帮我了解一下 numpy.std() 的工作原理吗？

Answer 1

这个问题的症结在于你需要除以N(3)，而不是N-1(2)。正如 Iarsmans 指出的，numpy 将使用总体方差，而不是样本方差。

所以真正的答案是 sqrt(2/3) 这正是:0.8164965...

如果您碰巧尝试故意为自由度使用不同的值(默认值 0)，请使用关键字参数 ddof 使用非 0 的正值:

np.std(data, ddof=1)

...但是这样做 here 会重新引入您原来的问题，因为 numpy 会除以 N - ddof。