unfair-casino.py

General Hidden Markov Model Library 一个通用的隐马尔科夫模型的C代码库

#/usr/bin/env python2.3
from ghmm import *

throws = IntegerRange(1,6)
dist = DiscreteDistribution(throws)

m = HMMFromMatrices(throws,dist,
                    [[0.9,0.1],[0.2,0.8]],
                    [[0.16, 0.16, 0.16, 0.16, 0.16, 0.16],[0.16, 0.9, 0.16, 0.16, 0.16, 0.16]],
                    [0.8,0.2])
# How do we want to do that ...


observations = [...] # Or require EmissionSequence(throws,[...])
states = m.viterbi(EmissionSequence(throws, observations)) 

# Alternatively: lambda knows EmissionDomain: Can call EmissionSequence
states = m.viterbi(observations)

print #sequence + viterbi

# posterior decoding
#forward = lamdba.forward(observations)

### How can we know, that the state with index 0 is 'fair' ???

#posterior_fair =