Predecir emisiones futuras del modelo HMM ajustado
He instalado unHMM
modelar a mis datos usandohmm.discnp
paquete en R de la siguiente manera:
library(hmm.discnp)
zs <- hmm(y=lis,K=5)
Ahora quiero predecir las futuras observaciones de K (emisiones) de este modelo. Pero solo puedo obtener la secuencia de estado más probable para las observaciones que ya tengo a través deViterbi
algoritmo.
yo tengot
emisiones ya, es decir(y(1),...,y(t))
.
Quiero el futuro mas probableK
emisiones de la equipadaHMM
objeto, es decir(y(t+1),...y(t+k))
.
¿Hay una función para calcular esto? si no, ¿cómo lo calculo manualmente?