Bonjour tout le monde
j’essai d'écrire un script en R sur l'algorithme progressif (forward) de Baum Welch mais je ne suis vraiment pas doué en R.
on a ces données
la loi initial µ1=(0.3, 0.6, 0.1)
la matrice de transition T1
0.3 0.2 0.5
0.3 0.4 0.3
0.3 0.3 0.4
la matrice d'émission E1
0.3 0.2 0.2 0.3
0.3 0.1 0.2 0.4
0.2 0.3 0.4 0.1
le nombre d'observation n = 10
les valeurs observées X(n) =(1, 3, 2, 4, 1, 4, 1, 3, 2, 4)
je veux calculer F1 , . . . , Fn
en gros : F1(1) F1(2) F1(3) F2(1) F2(2) F2(3) ... Fn(1) Fn(2) Fn(3)
sachant que
Fk(u) = P(X1=i1 X2=i2 , ...Xk=ik, Sk=u)
j'ai réussi à introduire quelques données que voici dans R
µ1=t(c(0.3,0.6,0.1)) #la loi initial
(T1=matrix(c(0.3,0.2,0.5,0.3,0 .4,0.3,0.3,0.3,0.4), nrow=3,byrow=TRUE)) # la matrice transitoire à n=1
(E1=matrix(c(0.3,0.2,0.2,0.3,0 .3,0.1,0.2,0.4,0.2,0.3,0.4,0.1 ), ncol=4,byrow=TRUE)) # la matrice d'émission à n=1
mais à présent je ne sais pas trop comment faire pour la suite
si vous connaissez comment faire ou avez besoin de plus d'informations n'hésitez pas!
-----