马尔可夫链概率计算- Python

Question

我有一个带有马尔可夫链模型状态转移概率的Python字典。

dict_m = {('E', 'F'): 0.29032258064516131, ('D', 'F'): 0.39726027397260272, ('D', 'D'): 0.30136986301369861, ('E', 'D'): 0.32258064516129031, ('E', 'E'): 0.38709677419354838, ('D', 'E'): 0.30136986301369861, ('F', 'F'): 0.68152866242038213, ('F', 'E'): 0.10191082802547771, ('F', 'D'): 0.19108280254777071, ('F', 'H'): 0.025477707006369428, ('H', 'F'): 1.0}

假设我有一个如下的数据序列。状态转换从D到E, E到F, F到E，等等

s = ['D','E','F','E','E','F','H','F']

现在我需要通过乘以dict_m中的概率来计算概率值。在这种情况下，所有状态转移的概率为

probability = 0.301370*0.290323*0.101911*0.387097*0.290323*0.025478*1.000000

谁能帮我解释一下这个?

print zip(s,s[1:])
[('D', 'E'), ('E', 'F'), ('F', 'E'), ('E', 'E'), ('E', 'F'), ('F', 'H'), ('H', 'F')]

Answer 1

您必须使用zip中的元组来索引字典，然后将所有数字相乘。

from operator import mul
print reduce(mul, (dict_m[t] for t in zip(s, s[1:])))