Я пытаюсь подобрать модель марков, скрытых пуассоном, но получаю следующую ошибку. Примечание: Если я использую 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
Любые предложения о том, как решить эту ошибку.