scipy.stats.

circvar#

scipy.stats.circvar(samples, high=6.283185307179586, low=0, axis=None, nan_policy='propagate', *, keepdims=False)[source]#

计算角度观测样本的圆形方差。

给定 \(n\) 个以弧度为单位测量的角度观测值 \(x_1, \cdots, x_n\)，它们的圆形方差定义为 ([2]，方程 2.3.3)

\[1 - \left| \frac{1}{n} \sum_{k=1}^n e^{i x_k} \right|\]

其中 \(i\) 是虚数单位，\(|z|\) 表示复数 \(z\) 的长度。上式中的 \(|z|\) 被称为 平均结果长度。

参数:

samplesarray_like

角度观测值的输入数组。完整角度的值等于 (high - low)。

highfloat, 可选

角度主值的边界。默认为 2*pi。

lowfloat, 可选

角度主值的边界。默认为 0。

axisint 或 None, 默认为 None

如果为整数，则为输入沿其计算统计量的轴。输入的每个轴切片（例如，行）的统计量将出现在输出的相应元素中。如果为 None，则在计算统计量之前将对输入进行扁平化。

nan_policy{‘propagate’, ‘omit’, ‘raise’}

定义如何处理输入 NaN。

propagate: 如果在沿其计算统计量的轴切片（例如，行）中存在 NaN，则输出的相应条目将为 NaN。
omit: 执行计算时将忽略 NaN。如果在沿其计算统计量的轴切片中剩余的数据不足，则输出的相应条目将为 NaN。
raise: 如果存在 NaN，则会引发 ValueError。

keepdimsbool, 默认为 False

如果将其设置为 True，则缩减的轴将保留在结果中，作为大小为一的维度。使用此选项，结果将与输入数组正确广播。

返回值:

circvarfloat

圆形方差。返回值在 [0, 1] 范围内，其中 0 表示无方差，1 表示大方差。

如果输入数组为空，则返回 np.nan。

参见

circmean: 圆形平均值。
circstd: 圆形标准差。

备注

在小角度的极限情况下，如果以弧度为单位测量，圆形方差接近“线性”方差的一半。

从 SciPy 1.9 开始，np.matrix 输入（不推荐用于新代码）在执行计算之前将转换为 np.ndarray。在这种情况下，输出将是适当形状的标量或 np.ndarray，而不是 2D np.matrix。同样，虽然会忽略掩码数组的掩码元素，但输出将是标量或 np.ndarray，而不是具有 mask=False 的掩码数组。

参考文献

[1]

Fisher, N.I. Circular data 的统计分析. 剑桥大学出版社, 1993.

[2]

Mardia, K. V. and Jupp, P. E. 方向统计. 约翰·威利父子公司, 1999.

示例

>>> import numpy as np
>>> from scipy.stats import circvar
>>> import matplotlib.pyplot as plt
>>> samples_1 = np.array([0.072, -0.158, 0.077, 0.108, 0.286,
...                       0.133, -0.473, -0.001, -0.348, 0.131])
>>> samples_2 = np.array([0.111, -0.879, 0.078, 0.733, 0.421,
...                       0.104, -0.136, -0.867,  0.012,  0.105])
>>> circvar_1 = circvar(samples_1)
>>> circvar_2 = circvar(samples_2)

绘制样本。

>>> fig, (left, right) = plt.subplots(ncols=2)
>>> for image in (left, right):
...     image.plot(np.cos(np.linspace(0, 2*np.pi, 500)),
...                np.sin(np.linspace(0, 2*np.pi, 500)),
...                c='k')
...     image.axis('equal')
...     image.axis('off')
>>> left.scatter(np.cos(samples_1), np.sin(samples_1), c='k', s=15)
>>> left.set_title(f"circular variance: {np.round(circvar_1, 2)!r}")
>>> right.scatter(np.cos(samples_2), np.sin(samples_2), c='k', s=15)
>>> right.set_title(f"circular variance: {np.round(circvar_2, 2)!r}")
>>> plt.show()