向前、向后、和向前-向后过程的HMM r语言实现 & Viterbi algorithm vs Viterbi algorithm

forward, backward, forward-backward procedure by r

其实已经有r包可以实现了，自己写一遍过程便于理解。

• r包实现

exp(forward(hmm, observations))#forward
exp(backward(hmm, observations)) #backward
posterior(hmm, observations) #forward-backward

• 自己的代码实现

#set the probs
transProbs <- as.data.frame(matrix(c(0.7 , 0.4 , 0.3 , 0.6), nrow = 2))
States <- c("CpG","N")
rownames(transProbs) <- States
colnames(transProbs) <- States
startProbs <- data.frame('CpG' = 0.4, 'N' = 0.6)
emissionProbs <- as.