如果我使用 depmix()
和 fit()
使用 depmixS4
包拟合模型 mf
,我想知道生成给定序列 s
的对数似然,我应该怎么做?
我知道在 HiddenMarkov
包中我可以使用 forwardback(s, mf$Pi, mf$...)$LL
来获得对数似然,但是我发现 depmixS4
中的 forwardbackward()
函数的工作方式不同。
最佳答案
终于我自己找到了方法,因为第一次使用 S4...
首先用depmix()
初始化一个模型(depmix obj),然后用setpars()
在拟合的mf
中更新参数>。之后 forwardbackward()
将起作用
init_mod <- depmix(respones, data, nstates) ## no solid value
mod <- setpars(init_mod, getpars(fm))
forwardbackward(mod)
谁知道更好的解决方案请告诉我...
关于r - 使用 R 中的 depmixS4 使用拟合模型评估序列,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/31019622/