Numpy __array_prepare__错误

Question

我正在尝试让我在网上找到的一个做期望最大化(http://code.activestate.com/recipes/577735-expectation-maximization/)的食谱起作用。我遇到了以下错误：

Traceback (most recent call last):
  File "./runem.py", line 7, in <module>
    print expectation_maximization([[1,2,3,4,5],[2,3,4,5,6],[9,8,7,4,1]], 2)
  File "/local/scratch-3/dk427/rp/em.py", line 83, in expectation_maximization
    Px[o,c] = pnorm(t[o,:], params[c]['mu'], params[c]['sigma'])
  File "/local/scratch-3/dk427/rp/em.py", line 18, in pnorm
    xmt = np.matrix(x-m).transpose()
TypeError: __array_prepare__ must return an ndarray or subclass thereof which is otherwise identical to its input

一定是算法中有一些缺陷，或者我给了它错误的输入，但我找不到哪里出了问题。我发现错误是由减法x-m引起的，而不是x.dtype=int64和m.dtype=float64，我认为这应该可以工作。

有谁有什么想法吗？

Answer 1

您似乎传入了一个列表列表，而不是一个数组。你可以这样做：

ts = np.array([[1,2,3,4,5],[2,3,4,5,6],[9,8,7,4,1]])
expectation_maximization(ts, 2)

在我的计算机上取平方根似乎有一些问题，但我认为这可能是因为出于某种原因，这些数据对这个算法并不好(但我不知道这个算法想要做什么，所以我不能确定)。