1

我正在尝试使用我在网上找到的用于最大化期望的配方（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=int64m.dtype=float64

有没有人有任何想法？

4

1 回答 1

2

您似乎传递的是列表列表而不是数组。您可以执行以下操作：

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

在我的计算机上取平方根似乎在某些时候存在一些问题，但我认为这可能是因为由于某种原因该数据不适合该算法（但我不知道该算法试图做什么做，所以我不能肯定地说）。

于 2012-05-09T04:50:59.383 回答