Hidden Markov Model result(decoding )changes everytime

Question

I use HMM to predict the behavior of humans. Hidden states are Rest walk eat and observations are inside, outside and snack respectively.

hState = ['Rset', 'Walk', 'Eat']
Obs =  ['Inside 1', 'Outside 2', 'Snack 3']
order = [1, 2, 1, 2, 3 ,2, 2, 3, 1, 2 ,2, 2, 1, 2]
no = 3
sProb = np.array([0.1, 0.8, 0.1])
tProb = np.array([[0.2, 0.2, 0.6], 
                     [0.1, 0.8, 0.1], 
                     [0.6, 0.1, 0.3]])
eProb = np.array([[0.7, 0.2, 0.1],  
                    [0.1, 0.6, 0.3],
                    [0.2, 0.3, 0.5]])
h = hmm.MultinomialHMM(3, "full", sProb, tProb)
h.emissionprob_ = eProb
oreder1 = np.array([order]).T
result= hmm.GaussianHMM(n_components=no).fit(oreder1)
result.predict(oreder1)

But the result is changed whenever I run the code again. Sequences of the result are correct but the hidden state id is change The result of the first time run

array([2, 1, 2, 1, 0, 1, 1, 0, 2, 1, 1, 1, 2, 1])

The result of the second time run

array([0, 1, 0, 1, 2, 1, 1, 2, 0, 1, 1, 1, 0, 1])

What I want to know is how can I get constant value for each run or how to finalize the model. Thank you for your answer in advance.

Instead of sharing a screenshot of your code, include the code as a code block in your question. — Gordon Bean, Sep 29 '20 at 18:10

score 0 · Accepted Answer · answered Sep 29 '20 at 18:15

Hidden Markov Models (HMMs) are stochastic - they use a random-number generator to help make decisions. So if you run the code multiple times, each execution will be different because the random numbers used behind the scenes are different.

If you want to make each run use the same sequence of random numbers, you'll need to set the random number "seed". Look in the hmm documentation for how to set the seed. Setting it to a fixed value will ensure that each execution uses the same series of random numbers and should result in identical results for each run.

The question of finalizing the model is way beyond the scope of this site. Essentially, you'll have to learn about HMMs, how they work, what assumptions they make, and how your use-case fits (or doesn't!) those assumptions. Then you use your best judgement combined with rigorous testing/training to determine the final parameterization.

Thank you it works. Another small question do you know how to predict the next hidden state applying this same code — Yash, Sep 30 '20 at 04:34
I read, but many of them mentioned that it should calculate by multiplying the existing hidden state by the transition probability. But it seems there must have code level solution, like a method or something — Yash, Sep 30 '20 at 16:49

Hidden Markov Model result(decoding )changes everytime

1 Answers1