numpy.corrcoef#
- numpy.corrcoef(x, y=None, rowvar=True, bias=<no value>, ddof=<no value>, *, dtype=None)[源]#
返回皮尔逊积矩相关系数。
有关更多详细信息,请参阅
cov
的文档。相关系数矩阵 R 和协方差矩阵 C 之间的关系为\[R_{ij} = \frac{ C_{ij} } { \sqrt{ C_{ii} C_{jj} } }\]R 的值在 -1 到 1 之间(包含 -1 和 1)。
- 参数:
- x类数组
包含多个变量和观测值的 1 维或 2 维数组。 x 的每一行代表一个变量,每一列代表所有这些变量的一个观测值。另请参阅下面的 rowvar。
- y类数组,可选
一组额外的变量和观测值。 y 的形状与 x 相同。
- rowvar布尔,可选
如果 rowvar 为 True(默认),则每行代表一个变量,列中包含观测值。否则,关系将转置:每列代表一个变量,而行中包含观测值。
- bias_无值,可选
无作用,请勿使用。
自 1.10.0 版本起已弃用。
- ddof_无值,可选
无作用,请勿使用。
自 1.10.0 版本起已弃用。
- dtype数据类型,可选
结果的数据类型。默认情况下,返回的数据类型将至少具有
numpy.float64
的精度。1.20 版本新增。
- 返回:
- Rndarray
变量的相关系数矩阵。
另请参阅
cov
协方差矩阵
备注
由于浮点四舍五入,生成的数组可能不是 Hermitian 矩阵,对角线元素可能不是 1,并且元素可能不满足不等式 abs(a) <= 1。实部和虚部被截断到区间 [-1, 1] 以试图改善这种情况,但这对于复数情况帮助不大。
此函数接受但会丢弃参数 bias 和 ddof。这是为了与此函数以前的版本向后兼容。这些参数对函数的返回值没有影响,在当前和以前的 numpy 版本中可以安全地忽略它们。
示例
>>> import numpy as np
在此示例中,我们生成两个随机数组
xarr
和yarr
,并计算行向和列向的皮尔逊相关系数R
。由于rowvar
默认为 true,我们首先计算xarr
变量之间的行向皮尔逊相关系数。>>> import numpy as np >>> rng = np.random.default_rng(seed=42) >>> xarr = rng.random((3, 3)) >>> xarr array([[0.77395605, 0.43887844, 0.85859792], [0.69736803, 0.09417735, 0.97562235], [0.7611397 , 0.78606431, 0.12811363]]) >>> R1 = np.corrcoef(xarr) >>> R1 array([[ 1. , 0.99256089, -0.68080986], [ 0.99256089, 1. , -0.76492172], [-0.68080986, -0.76492172, 1. ]])
如果我们添加另一组变量和观测值
yarr
,我们可以计算xarr
和yarr
中变量之间的行向皮尔逊相关系数。>>> yarr = rng.random((3, 3)) >>> yarr array([[0.45038594, 0.37079802, 0.92676499], [0.64386512, 0.82276161, 0.4434142 ], [0.22723872, 0.55458479, 0.06381726]]) >>> R2 = np.corrcoef(xarr, yarr) >>> R2 array([[ 1. , 0.99256089, -0.68080986, 0.75008178, -0.934284 , -0.99004057], [ 0.99256089, 1. , -0.76492172, 0.82502011, -0.97074098, -0.99981569], [-0.68080986, -0.76492172, 1. , -0.99507202, 0.89721355, 0.77714685], [ 0.75008178, 0.82502011, -0.99507202, 1. , -0.93657855, -0.83571711], [-0.934284 , -0.97074098, 0.89721355, -0.93657855, 1. , 0.97517215], [-0.99004057, -0.99981569, 0.77714685, -0.83571711, 0.97517215, 1. ]])
最后,如果我们使用选项
rowvar=False
,则列现在被视为变量,我们将找到xarr
和yarr
中变量之间的列向皮尔逊相关系数。>>> R3 = np.corrcoef(xarr, yarr, rowvar=False) >>> R3 array([[ 1. , 0.77598074, -0.47458546, -0.75078643, -0.9665554 , 0.22423734], [ 0.77598074, 1. , -0.92346708, -0.99923895, -0.58826587, -0.44069024], [-0.47458546, -0.92346708, 1. , 0.93773029, 0.23297648, 0.75137473], [-0.75078643, -0.99923895, 0.93773029, 1. , 0.55627469, 0.47536961], [-0.9665554 , -0.58826587, 0.23297648, 0.55627469, 1. , -0.46666491], [ 0.22423734, -0.44069024, 0.75137473, 0.47536961, -0.46666491, 1. ]])