熊猫中的矩阵乘法

时间:2013-05-09 23:18:13

标签: python pandas

我将数字数据存储在两个DataFrame x和y中。来自numpy的内部产品起作用,但来自熊猫的点产品不起作用。

In [63]: x.shape
Out[63]: (1062, 36)

In [64]: y.shape
Out[64]: (36, 36)

In [65]: np.inner(x, y).shape
Out[65]: (1062L, 36L)

In [66]: x.dot(y)
---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-66-76c015be254b> in <module>()
----> 1 x.dot(y)

C:\Programs\WinPython-64bit-2.7.3.3\python-2.7.3.amd64\lib\site-packages\pandas\core\frame.pyc in dot(self, other)
    888             if (len(common) > len(self.columns) or
    889                     len(common) > len(other.index)):
--> 890                 raise ValueError('matrices are not aligned')
    891 
    892             left = self.reindex(columns=common, copy=False)

ValueError: matrices are not aligned

这是一个错误还是我使用熊猫错了?

1 个答案:

答案 0 :(得分:27)

xy的形状不仅必须正确,而且还必须正确 x的列名必须与y的索引名称匹配。除此以外 pandas/core/frame.py中的此代码将引发ValueError:

if isinstance(other, (Series, DataFrame)):
    common = self.columns.union(other.index)
    if (len(common) > len(self.columns) or
        len(common) > len(other.index)):
        raise ValueError('matrices are not aligned')

如果您只想计算矩阵产品而不使x的列名与y的索引名匹配,则使用NumPy点函数:

np.dot(x, y)

x的列名必须与y的索引名称匹配的原因是因为pandas dot方法会重新索引xy所以如果x的列顺序与y的索引顺序不自然匹配,则在执行矩阵产品之前,它们将匹配:

left = self.reindex(columns=common, copy=False)
right = other.reindex(index=common, copy=False)

NumPy dot函数不会这样做。它只会根据底层数组中的值计算矩阵产品。


以下是重现错误的示例:

import pandas as pd
import numpy as np

columns = ['col{}'.format(i) for i in range(36)]
x = pd.DataFrame(np.random.random((1062, 36)), columns=columns)
y = pd.DataFrame(np.random.random((36, 36)))

print(np.dot(x, y).shape)
# (1062, 36)

print(x.dot(y).shape)
# ValueError: matrices are not aligned