Discussion :

[Gloria_Im]

Bonjour à tous,
j'ai créé un hmm comme suit :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
 
hmm=initHMM(c("P1","P2","P3","P4"),
               c(".0008",".0007",".0004",".0007"),
               startProbs=c(1,.0,.0,.0),
               transProbs = matrix(c(NA,1,NA,NA,
                                     NA,NA,1,NA,
                                     NA,NA,NA,1,
                                     NA,NA,NA,1),nrow=4,ncol=4,byrow=TRUE),
 
               emissionProbs=matrix(c(.25,.25,.25,.25,
                                      .25,.25,.25,.25,
                                      .25,.25,.25,.25,
                                      .25,.25,.25,.25),nrow=4,ncol=4,byrow=TRUE))

Je veux apprendre à cet HMM en utilisant viterbitraining un simple ensemble d'observations

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
observations=c("6.405229e-002","0","2.745098e-002","0")

d'ailleurs voilà la commande d'apprentissage

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
viterbiT = viterbiTraining(hmm,observations,maxIterations=10, delta=1E-9, pseudoCount=1)

mais le problème que ça me génère une erreur indice hors limite mais je ne sais pas où est le problème.
Si quelqu'un peut m'aider je serais reconnaissante et merci d'avance

[A. D.]

Bonjour,

Tout d'abord, pensez à préciser les packages utilisés : il s'agit bien du package "HMM" ici?

Sinon, vous pouvez accéder au code des fonctions utilisées comme suit :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
> library(HMM)
 
> initHMM
> viterbiTraining

Et ensuite je vous conseille de tester ce code "morceau par morceau" afin de détecter l'endroit où se situe le problème d'indice.

Bonne continuation


Cordialement,

A.D.

[Gloria_Im]

Merci pour votre réponse.
Oui, c'est le package HMM et j'ai déjà fais tout ce que vous me suggérez avant même de poster ma question, mais toujours la même erreur !
C'est bizarre, je ne sais pas où est le problème

Je me demande que peut être le tableau des symboles émis par les états doit être des strings mais en fait ".000122" est une chaîne et .000122 et un float donc normalement ça devrait pas beuguer

Je n'ai pas encore trouvé de solution.

[Kevin.Rousseeuw]

Bonjour,

J'ai exactement le même problème. Pour essayer de le régler, j'ai affiché l'ensemble des matrices de transitions et d'émissions que la fonction viterbiTraining a générées.
Je me suis rendu compte qu'au bout d'un moment, le nombre d'états qu'il y avait dans le vecteur de sortie diminuait et donc, les matrices n'étaient plus les mêmes. Comme le nombre d'états qu'il génère et donc qu'il utilise est différent du nombre d'états initial, il se retrouve avec des valeurs manquantes qu'il est dans l'incapacité de traiter.

Par contre, je n'ai pas encore trouvé le moyen de palier le problème.

Kévin Rousseeuw