如果我使用mf和depmix()使用fit()包适合模型depmixS4，我想知道生成给定序列的对数可能性{{ 1}}，我该怎么办？
我知道在s包中我可以使用HiddenMarkov来获取对数可能性，但我发现forwardback(s, mf$Pi, mf$...)$LL中的forwardbackward()函数不起作用。< / p>
答案 0 :(得分：2)
最后，当我第一次使用S4时，我找到了自己的方式......
首先使用depmix()初始化模型（depmix obj），然后使用mf使用拟合setpars()中的参数更新它。之后forwardbackward()将起作用
init_mod <- depmix(respones, data, nstates) ## no solid value
mod <- setpars(init_mod, getpars(fm))
forwardbackward(mod)
任何人都知道更好的解决方案，请告诉我......