如何使用 hmmlearn 在 Python 中运行隐藏的马尔可夫模型?

问题描述

我尝试使用 hmmlearn from GitHub 来运行二进制隐藏马尔可夫模型.这不起作用:

I tried to use hmmlearn from GitHub to run a binary hidden markov model. This does not work:

import hmmlearn.hmm as hmm
transmat = np.array([[0.7, 0.3],
                      [0.3, 0.7]])
emitmat = np.array([[0.9, 0.1],
                    [0.2, 0.8]])
obs = np.array([0, 0, 1, 0, 0])
startprob = np.array([0.5, 0.5])
h = hmm.MultinomialHMM(n_components=2, startprob=startprob,
                       transmat=transmat)
h.emissionprob_ = emitmat
# fails
h.fit([0, 0, 1, 0, 0])
# fails
h.decode([0, 0, 1, 0, 0])
print h

我得到这个错误:

ValueError:零维数组不能串联

使用这个模块的正确方法是什么?请注意，我使用的是与 sklearn 分离的 hmmlearn 版本，因为显然 sklearn 不再维护 hmmlearn.

What is the right way to use this module? Note I am using the version of hmmlearn that was separated from sklearn, because apparently sklearn doesn't maintain hmmlearn anymore.

推荐答案

Fit 接受 序列列表 而不是单个序列(因为通常你可以从不同的运行中观察到多个独立的序列你的实验/观察).因此，只需将您的列表放入另一个列表中

Fit accepts list of sequences and not a single sequence (as in general you can have multiple, independent sequences observed from different runs of your experiments/observations). Thus simply put your list inside another list

import hmmlearn.hmm as hmm
import numpy as np

transmat = np.array([[0.7, 0.3],
                      [0.3, 0.7]])
emitmat = np.array([[0.9, 0.1],
                    [0.2, 0.8]])

startprob = np.array([0.5, 0.5])
h = hmm.MultinomialHMM(n_components=2, startprob=startprob,
                       transmat=transmat)
h.emissionprob_ = emitmat
# works fine
h.fit([[0, 0, 1, 0, 0]])
# h.fit([[0, 0, 1, 0, 0], [0, 0], [1,1,1]]) # this is the reason for such
                                            # syntax, you can fit to multiple
                                            # sequences
print h.decode([0, 0, 1, 0, 0])
print h

给予

(-4.125363362578882, array([1, 1, 1, 1, 1]))
MultinomialHMM(algorithm='viterbi',
        init_params='abcdefghijklmnopqrstuvwxyzABCDEFGHIJKLMNOPQRSTUVWXYZ',
        n_components=2, n_iter=10,
        params='abcdefghijklmnopqrstuvwxyzABCDEFGHIJKLMNOPQRSTUVWXYZ',
        random_state=<mtrand.RandomState object at 0x7fe245ac7510>,
        startprob=None, startprob_prior=1.0, thresh=0.01, transmat=None,
        transmat_prior=1.0)

这篇关于如何使用 hmmlearn 在 Python 中运行隐藏的马尔可夫模型?的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持！