由Python Numpy计算的协方差矩阵每次更改

用户名

我有一个1043 * 261矩阵,其数字在0和1之间非常小,我使用numpy.cov()函数计算了1043 * 1043协方差矩阵。我尝试运行了几次代码,并得到了相似(但不完全相同)的协方差矩阵,但是协方差矩阵中的元素因e-7的比例而略有不同。有时这会使协方差矩阵变为非PSD,这将对我造成严重的问题。

有谁知道为什么会存在差异以及如何解决差异?

附带的是我通过两次运行相同的代码获得的两个协方差矩阵。如果按元素对它们进行比较,将会看到一些细微的差异:

No. 1
[[  5.05639177e-06   2.44041401e-06   3.30187175e-06 ...,   1.66634014e-06
4.03972183e-06   1.18433575e-06]
[  2.44041401e-06   9.67277658e-06   9.04356309e-06 ...,   2.50668884e-06
5.43371939e-06   4.74297546e-06]
[  3.30187175e-06   9.04356309e-06   2.09334309e-05 ...,   3.13977728e-06
8.69946165e-06   6.15981652e-06]
..., 
[  1.66634014e-06   2.50668884e-06   3.13977728e-06 ...,   4.20175297e-06
4.16076781e-06   1.59827406e-06]
[  4.03972183e-06   5.43371939e-06   8.69946165e-06 ...,   4.16076781e-06
2.58010941e-05   3.02797946e-06]
[  1.18433575e-06   4.74297546e-06   6.15981652e-06 ...,   1.59827406e-06
3.02797946e-06   6.60805238e-06]]

No.2
[[  5.05997030e-06   2.42187179e-06   3.30788097e-06 ...,   1.66495376e-06
4.03676937e-06   1.17413702e-06]
[  2.42187179e-06   9.60677140e-06   9.05219266e-06 ...,   2.50338648e-06
5.42679569e-06   4.75547515e-06]
[  3.30788097e-06   9.05219266e-06   2.04172017e-05 ...,   3.13058624e-06
8.67976701e-06   6.28137859e-06]
..., 
[  1.66495376e-06   2.50338648e-06   3.13058624e-06 ...,   4.20175297e-06
4.16076781e-06   1.59827884e-06]
[  4.03676937e-06   5.42679569e-06   8.67976701e-06 ...,   4.16076781e-06
2.58010941e-05   3.02810307e-06]
[  1.17413702e-06   4.75547515e-06   6.28137859e-06 ...,   1.59827884e-06
3.02810307e-06   6.63834973e-06]]

非常感谢你!

韦德拉克

numpy.cov 似乎是确定性的:

import numpy

randoms = numpy.random.random((1043, 261))

covs = [numpy.cov(randoms) for _ in range(10)]
all((c==covs[0]).all() for c in covs)
#>>> True

我以为问题出在其他地方。

另外请注意,这个结果与数字1000持有大小

本文收集自互联网,转载请注明来源。

如有侵权,请联系[email protected] 删除。

编辑于
0

我来说两句

0条评论
登录后参与评论

相关文章

来自分类Dev

Python Numpy计算的协方差矩阵每次更改

来自分类Dev

Python:手动协方差矩阵

来自分类Dev

椭圆的协方差矩阵

来自分类Dev

numpy polyfit-协方差矩阵

来自分类Dev

Python:从列表创建协方差矩阵

来自分类Dev

如何在Python中获取n个数组的nxn个协方差矩阵?

来自分类Dev

在Python中建立协方差矩阵

来自分类Dev

numpy协方差和公式的协方差矩阵产生不同的结果

来自分类Dev

numpy中的加权协方差矩阵

来自分类Dev

Python溢出的原因是许多概率的乘积,大协方差矩阵的求逆以及其行列式的计算

来自分类Dev

如何计算3d numpy数组的协方差矩阵?

来自分类Dev

不使用numpy计算协方差矩阵

来自分类Dev

Python中的不确定性包:使用给定的协方差矩阵获取数据不确定性

来自分类Dev

在python中跨矩阵计算单元内的方差

来自分类Dev

从OpenCV Expectation Max中的协方差矩阵计算方差

来自分类Dev

Python:手动协方差矩阵

来自分类Dev

相关协方差矩阵到方差协方差矩阵

来自分类Dev

numpy polyfit-协方差矩阵

来自分类Dev

NUMPY:嵌套列表的协方差矩阵

来自分类Dev

Python:从列表创建协方差矩阵

来自分类Dev

Numpy协方差矩阵非归一化

来自分类Dev

如何计算块中的协方差矩阵?

来自分类Dev

协方差矩阵不是正定的

来自分类Dev

计算改进的kppm模型的方差协方差矩阵

来自分类Dev

从出现列表计算协方差矩阵

来自分类Dev

如何在 Python 中生成随机协方差矩阵?

来自分类Dev

Python中尺度矩阵(协方差矩阵)的高效计算

来自分类Dev

如何计算数据框的协方差矩阵

来自分类Dev

从协方差的 Pandas DataFrame 创建一个 numpy 协方差矩阵

Related 相关文章

热门标签

归档