R中HMM的发射概率

Question

我们如何计算R中隐马尔可夫模型(HMM(的发射概率？

至于计算转移概率，我们使用函数

tr <- seqtrate(exampledata)

此函数返回一个转换矩阵。示例数据是顺序数据。

是否有一个函数可以返回我们的发射矩阵？

Answer 1

请查看 R 的 HMM 包 https://cran.r-project.org/web/packages/HMM/HMM.pdf

你可以在那里找到这样的例子

hmm = initHMM(c("A","B"), c("L","R"), transProbs=matrix(c(.8,.2,.2,.8),2),
emissionProbs=matrix(c(.6,.4,.4,.6),2))
print(hmm)
# Sequence of observations
observation = c("L","L","R","R")
baumWelch(hmm, observation, maxIterations=100, delta=1E-9, pseudoCount=0)

鲍姆韦尔奇算法返回更新的发射概率。