如何使用 hidden_​​markov 修复广播状态和发射概率的错误？

Question

在其文档中的一个简单示例中应用轻微更改时，我遇到了 hidden_​​markov 包的问题。在下面的代码中，我尝试了 2 个状态和 3 个可能的观察结果（在文档的示例中，有 2 个可能的观察结果并且代码工作正常）：

states = ('s', 't')
possible_observation = ('A', 'B', 'C')

# Numpy arrays of the data

start_probability = np.matrix( '0.5 0.5')
transition_probability = np.matrix('0.6 0.4 ; 0.3 0.7')
emission_probability = np.matrix( '0.3 0.2 0.5 ; 0.3 0.1 0.6')

# Initialize class object

test = hmm(states,possible_observation,
           start_probability,
           transition_probability,
           emission_probability)

observations = ('A', 'B','B','A', 'C')
obs4 = ('B', 'C', 'A','B')
observation_tuple = []
observation_tuple.extend( [observations,obs4] )
quantities_observations = [18, 28]
num_iter=1000

e,t,s = test.train_hmm(observation_tuple,num_iter,quantities_observations)

运行代码后出现错误：

ValueError: operands could not be broadcast together with shapes (2,3) (1,2)

有趣的是，当我尝试 3 个状态和 3 个可能的观察结果（并根据此更改修改概率矩阵）时，代码运行良好。要么我遗漏了一些东西，要么状态的数量和可能的观察结果应该总是相等的，这是没有意义的。

Answer 1

根据您上面提供的矩阵，只有两个矩阵的形状为(2,3) 和(1,2)，分别是emission_probability 和start_probability，错误是矩阵维度不匹配。匹配矩阵维度以便它们可以用于计算点积是通过做的。

emission_probability = emission_probability.T
start_probability = start_probability.T

这一步必须在初始化类对象之前完成。