0

运行时出现以下错误train.hmm()

ValueError: operands could not be broadcast together with shapes (7,20) (1,7)

我知道什么是广播错误并浏览了文档。但我无法更改矩阵的尺寸,因为它会在另一个地方弹出错误。代码如下所示。

start_probability = np.matrix( '0.5 0.02 0.18 0.05 0.01 0.2 0.04 ')
transition_probability = np.matrix('0.9 0.012 0.012 0.012 0.002 0.05 0.012 \
;  0.03 0.35 0.05 0.05 0.02 0.45 0.05 \
;  0.03 0.05 0.45 0.05 0.02 0.35 0.05 \
;  0.1 0.1 0.1 0.4 0.1 0.1 0.1 \
;  0.1 0.1 0.1 0.1 0.4 0.1 0.1 \
;  0.2 0.05 0.05 0.05 0.05 0.4 0.2 \
;  0.12 0.12 0.12 0.12 0.01 0.12 0.39')
emission_probability = np.matrix(np.ones((7, 20)) * 0.05)

test = hmm(states,possible_observation,start_probability,transition_probability,emission_probability)
observations = ['A', 'S','T','A']
obs4 = ['C', 'A','G']
observation_tuple = []
observation_tuple.extend( [observations,obs4] )
quantities_observations = [10, 20]
num_iter=1000
e,t,s = test.train_hmm(observation_tuple,num_iter,quantities_observations)

错误对应于最后一行。如果我转置start_prbabilityoremission_probability矩阵,我会test = hmm()在行本身中得到一个错误。

我在这里做错了什么?

4

1 回答 1

2

将第 349 行更改hmm_class.py为,

emProbNew = emProbNew/ np.reshape(em_norm.transpose(),[-1,1])

不是我所知道的最好的解决方案,但我认为作者忽略了一个事实,即Python并不总是能够知道他在这里所说的除法是什么意思。

于 2017-06-20T19:30:06.127 回答