原文:https://docs.scipy.org/doc/numpy/reference/generated/numpy.ma.var.html
校对:(虚位以待)
numpy.ma.
var
(self, axis=None, dtype=None, out=None, ddof=0, keepdims=<class numpy._globals._NoValue at 0x7faefc173ef0>) = <numpy.ma.core._frommethod instance>计算沿指定轴的方差。
返回数组元素的方差,衡量分布的扩展。默认情况下为展开的数组计算方差,否则在指定的轴上。
参数: | a:array_like
axis:无或int或tuple ints,可选
dtype:数据类型,可选
out:ndarray,可选
ddof:int,可选
keepdims:bool,可选 |
---|---|
返回: | 方差:ndarray,请参阅上面的dtype参数
|
笔记
The variance is the average of the squared deviations from the mean, i.e., var = mean(abs(x - x.mean())**2)
.
平均值通常计算为x.sum() / N
,其中N = len(x)
。但是,如果指定ddof,则使用除数N - ddof 代替。
在标准统计实践中,ddof=1
提供了假设无限总体方差的无偏估计量。ddof=0
提供正态分布变量的方差的最大似然估计。
请注意,对于复数,绝对值在平方之前进行,因此结果始终为实数和非负数。
对于浮点输入,使用输入具有的相同精度计算方差。根据输入数据,这可能导致结果不准确,特别是对于float32(请参见下面的示例)。使用dtype
关键字指定更高精度的累加器可以缓解此问题。
例子
>>> a = np.array([[1, 2], [3, 4]])
>>> np.var(a)
1.25
>>> np.var(a, axis=0)
array([ 1., 1.])
>>> np.var(a, axis=1)
array([ 0.25, 0.25])
在单精度中,var()可能不准确:
>>> a = np.zeros((2, 512*512), dtype=np.float32)
>>> a[0, :] = 1.0
>>> a[1, :] = 0.1
>>> np.var(a)
0.20250003
计算float64中的方差更准确:
>>> np.var(a, dtype=np.float64)
0.20249999932944759
>>> ((1-0.55)**2 + (0.1-0.55)**2)/2
0.2025