使用depmix调整HMM时发生意外错误

时间:2019-06-01 23:32:34

标签: r hidden-markov-models

我正在尝试拟合泊松隐藏式马尔可夫模型,但出现以下错误。 注意:如果我使用family = gaussian(),它似乎可以工作。但是问题是我需要使用泊松分布,因为它们是计数。错误可能与table(y)有关,其中每个观察值只有一个计数。

library(depmixS4)
y <- c(13968L, 43557L, 18519L, 50253L, 10300L, 13250L, 8650L, 900L, 
       1013L, 11136L, 2877L, 7008L, 5542L, 4779L, 4900L, 5143L, 5464L, 
       5342L, 5728L, 5612L, 5607L, 6151L, 3908L, 4455L, 4885L, 5118L, 
       5175L, 5566L, 5788L, 6349L, 7014L, 7013L, 7748L, 8561L, 10798L, 
       10665L, 9800L, 9969L, 12922L, 15735L, 17447L, 21604L, 22772L, 
       24304L, 28280L, 28225L, 28048L, 26676L, 25318L, 26510L, 24306L, 
       23096L, 24358L, 26540L, 20067L, 28837L, 23523L, 20210L, 18111L, 
       17277L, 17198L, 16140L, 15817L, 15459L, 14763L, 14990L, 15187L, 
       12626L, 12460L, 12164L, 12077L, 12575L)

mod <- fit(depmix(y~1, nstates=2,family=poisson(), ntimes=length(y)))
summary(mod)
Error in glm.fit(x = object@x, y = object@y, weights = w, family = object@family,  : 
  NA/NaN/Inf in 'x'
In addition: Warning message:
step size truncated due to divergence

有关如何解决此错误的任何建议。

0 个答案:

没有答案