Question

这是我在StackOverflow上的第一篇文章，我可以使用一些帮助......如果我没有遵循正确的发布协议，请原谅我。

StackOverflow中还有另一个例子，我的工作量很大，但我无法弄清楚如何调整代码。最重要的是，我正在研究所提问题的解决方案。

这是链接：

Getting the next observation from a HMM gaussian mixture distribution

一些背景知识：

RHmm - 从R Forge下载的2.1.0版。

RStudio - 0.98.953

R - 3.0.2 32位

我正在尝试用我的代码弄清楚以下问题：

如何从上面的链接修改解决方案（预测下一个观察点）以使用我的Baum-Welch模型？防爆。 hm_model＆lt; - HMMFit（obs = TWII_Train，nStates = 5）
当我运行Bam-Welch版本的hm_model＆lt; - HMMFit时，R / RStudio会话中止（obs = TWII_Train，dis =＆＃34; MIXTURE＆＃34;，nStates = 5，nMixt = 4 ）。您可以重新创建错误并提出解决方法吗？

这是我的R代码：

library(quantmod) 

library(RHmm)

getSymbols("^TWII")

TWII_Subset <- window(TWII, start=as.Date("2012-01-01"), end = as.Date("2013-04-01"))

TWII_Train <- cbind(TWII_Subset$TWII.Close - TWII_Subset$TWII.Open,
                    TWII_Subset$TWII.Volume)


hm_model <- HMMFit(obs=TWII_Train, nStates=5)

VitPath <- viterbi(hm_model, TWII_Train)

Answer 1

我不是这个软件包的用户，这不是一个真正的答案，但评论会掩盖一些结构。似乎缺少模型的“比例”值（因此结构不同。“均值”值如下所示：

$ mean   :List of 5
  ..$ : num [1:2] 6.72 3.34e+06
  ..$ : num [1:2] -12.4 2420174.5
  ..$ : num [1:2] -2.4 1832546.5
  ..$ : num [1:2] -10.4 1432636.1
  ..$ : num [1:2] 5.02 1.96e+06

我还怀疑你应该使用2和5而不是4和5来表示m和n。看看模型的其余部分：

str(hm_model)

R中的隐马尔可夫模型 - 用RHmm预测下一次观测

1 个答案: