hmm baum welch算法hmm $ missionProbs中的错误[状态，观察[1]]：下标超出范围

Question

我正在尝试根据EEG信号运行Baum-Welch。我对这些值进行了四舍五入，因此我可以假装它是一个离散信号，只是为了简化。因此，我有一列带有四舍五入的心电图，另一列有医生对睡眠阶段的评分。像这样：

医生eeg

1唤醒10

2唤醒13

3唤醒16

4唤醒18

我的代码在这里：

y = data$eeg
states = c("NonREM1", "NonREM2", "NonREM3", "REM", "Wake")
start_p = c(NonREM1 = 0,NonREM2 = 0,NonREM3 = 0, REM = 0, Wake = 1)
library(markovchain)
trans_p <- markovchainFit(data$doctor)$estimate
sleep.values <- unique(data[2])
hmm = initHMM(states, Symbols = sleep.values, startProbs=start_p, 
transProbs=trans_p)

baumWelch(hmm, observation = y, maxIterations=100, delta=1E-9, 
pseudoCount=0)

我似乎也无法输入排放概率矩阵，因此将其留空。当我填充它时，我得到

E [，]中的错误

有人可以向我解释哪里出了问题吗？我整个下午都被困在这。 ：\

Answer 1

谁在卡车上向我开枪，伙计们需要帮助以证明是他或他尽快导致我严重的偏执狂，沮丧和吸毒