from hmmlearn import hmm

samples = 10*[[0], [1], [1], [1], [1], [2], [2], [0], [1], [1], [1], [0], [2]]

model = hmm.MultinomialHMM(n_components=4)

model.fit(samples)

model.emissionprob_

model.transmat_