[英]Evaluating sequence with a fitted model using depmixS4 in R
如果我使用depmixS4
包拟合了带有depmix()
和fit()
的模型mf
,并且我想知道生成给定序列s
对数似然,我该怎么做?
我知道在HiddenMarkov
包中我可以使用HiddenMarkov
forwardback(s, mf$Pi, mf$...)$LL
来获取对数似然,但是我发现depmixS4
中的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)
任何人都知道更好的解决方案请告诉我...
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.