вероятности эмиссии для HMM в R

Question

Как мы можем вычислить вероятности выбросов для скрытой марковской модели (HMM) в R?

Что касается вычисления вероятностей переходов, мы используем функцию

tr <- seqtrate(exampledata)

, и эта функция возвращает матрицу переходов,Пример данных - это последовательные данные.

Есть ли функция, которая возвращает нам матрицу выбросов?

Answer 1

Пожалуйста, посмотрите на пакет HMM от R 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)

Алгоритм baumWelch возвращает обновленные вероятности выбросов.