moment#
- scipy.stats.moment(a, order=1, axis=0, nan_policy='propagate', *, center=None, keepdims=False)[源代码]#
计算样本关于均值的第 n 阶矩。
矩是衡量一组点形状的特定定量指标。它通常用于计算偏度和峰度的系数,因为它与它们密切相关。
- 参数:
- aarray_like
输入数组。
- orderint 或 int 的一维 array_like,可选
返回中心矩的阶数。默认值为 1。
- axisint 或 None,默认值:0
如果为 int,则为计算统计量的输入轴。输入的每个轴切片(例如,行)的统计量将出现在输出的相应元素中。如果
None
,则在计算统计量之前将输入展平。- nan_policy{‘propagate’, ‘omit’, ‘raise’}
定义如何处理输入 NaN。
propagate
:如果计算统计量的轴切片(例如,行)中存在 NaN,则输出的相应条目将为 NaN。omit
:执行计算时将省略 NaN。如果计算统计量的轴切片中剩余的数据不足,则输出的相应条目将为 NaN。raise
:如果存在 NaN,则会引发ValueError
。
- centerfloat 或 None,可选
计算矩的中心点。这可以是样本均值、原点或任何其他点。如果 None(默认),则将中心计算为样本均值。
- keepdimsbool,默认值:False
如果将其设置为 True,则缩减的轴将作为大小为 1 的维度保留在结果中。使用此选项,结果将与输入数组正确广播。
- 返回:
- 关于 `center` 的第 n 阶矩ndarray 或 float
如果 axis 为 None,则为给定轴或所有值上的适当矩。矩计算的分母是观测值的数量,不进行自由度校正。
注释
数据样本的 k 阶矩为
\[m_k = \frac{1}{n} \sum_{i = 1}^n (x_i - c)^k\]其中 n 是样本数,而 c 是计算矩的中心。此函数使用平方求幂 [1] 来提高效率。
请注意,如果 a 是一个空数组 (
a.size == 0
),则具有一个元素的数组moment
(moment.size == 1) 的处理方式与标量moment
(np.isscalar(moment)
) 相同。这可能会产生意外形状的数组。从 SciPy 1.9 开始,在执行计算之前,
np.matrix
输入(不建议用于新代码)将转换为np.ndarray
。在这种情况下,输出将是标量或适当形状的np.ndarray
,而不是二维np.matrix
。类似地,虽然会忽略掩码数组的掩码元素,但输出将是标量或np.ndarray
,而不是mask=False
的掩码数组。参考文献
示例
>>> from scipy.stats import moment >>> moment([1, 2, 3, 4, 5], order=1) 0.0 >>> moment([1, 2, 3, 4, 5], order=2) 2.0