Discussion :

[aymenbech]

Bonjour,

je travaille sur le HMM
en fait je donne comme donnée ( data ) un vecteur [3 0 0 1 0 2 0 1 6 8] comme entrée pour la fonction principale de HMM qui est dhmm_em (

Code :

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

[LL, prior2, transmat2, obsmat2] = dhmm_em(data, prior1, transmat1, obsmat1, 'max_iter', 10);

et ça me genere ces fautes :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
??? Attempted to access obsmat(:,0); index must be a positive integer or logical.
 
Error in ==> multinomial_prob at 19
  B(:,t) = obsmat(:, data(t));
 
Error in ==> dhmm_em>compute_ess_dhmm at 103
 obslik = multinomial_prob(obs, obsmat);
 
Error in ==> dhmm_em at 47
 [loglik, exp_num_trans, exp_num_visits1, exp_num_emit] = ...
 
Error in ==> dhmm_em_demo at 105
[LL, prior2, transmat2, obsmat2] = dhmm_em(data, prior1, transmat1, obsmat1, 'max_iter', 10);

est ce c'est concernant data ou quoi , si quelqu'un peut m'aider

[Jerome Briot]

Je ne connais pas cette Toolbox mais il semble bien que le vecteur data ne doivent contenir que des valeurs entières strictement positives.

On retrouve cette condition au début du code contenu dans le fichier multinomial_sample.m (par exemple) :

Code :

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

1
2
3
if any(X(:)==0)
  error('data must only contain positive integers')
end

[shaiHulud]

les états cachés doivent être numérotés 1,...,n. La raison est surement une instruction du genre transition= ProbMatrix(state,:) cachée dans le code

[aymenbech]

on a data =[3 0 0 1 0 2 0 1 6 8]

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
% initial guess of parameters
prior1 = normalise(rand(Q,1));
transmat1 = mk_stochastic(rand(Q,Q));
obsmat1 = mk_stochastic(rand(Q,O));
 
% improve guess of parameters using EM
[LL, prior2, transmat2, obsmat2] = dhmm_em(data, prior1, transmat1, obsmat1, 'max_iter', 10);
LL;
 
% use model to compute log likelihood
loglik = dhmm_logprob(data, prior2, transmat2, obsmat2);
% log lik is slightly different than LL(end), since it is computed after the final M step
%plot (prior1,'displayName');figure(gcf);
transmat1;
obsmat1;

j'ai pas compris ou est le probleme ???

[Invité]

Bonjour,

As-tu pris le temps de lire les réponses?

Je ne connais pas cette Toolbox mais il semble bien que le vecteur data ne doivent contenir que des valeurs entières strictement positives.

les états cachés doivent être numérotés 1,...,n.

Voir aussi la Subscript indices must either be real positive integers or logicals.
Il est là ton problème:

on a data =[3 0 0 1 0 2 0 1 6 8]