Bonjour!

j'ai un problème avec un exercice R que j'effectue. À noter que c'est un cours d'introduction, nous avons seulement appris des façons simple de faire du code. J'essaie de faire une fonction qui tire au hasard l'état au temps t+1 en fonction de l'état au temps t à partir d'un jeu de données.

Voici le début de ma matrice avec mes états de départ

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
> head(et.depart)
     borx bory  Etats
[1,] "0"  "0"   "M"  
[2,] "0"  "20"  "M"  
[3,] "0"  "40"  "T"  
[4,] "0"  "60"  "T"  
[5,] "0"  "80"  "T"  
[6,] "0"  "100" "T"
Et voici ma matrice d'états de transition:

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
             P(Bt+1|Xt) P(Mt+1|Xt) P(Tt+1|Xt) P(Rt+1|Xt)
P(Yt+1|Bt)     0.9059     0.0932     0.0000     0.0009
P(Yt+1|Mt)     0.0004     0.8454     0.1452     0.0090
P(Yt+1|Tt)     0.0000     0.0350     0.9594     0.0056
P(Yt+1|Rt)     0.0001     0.0020     0.0011     0.9968
Donc, je veux, par exemple, que si l'état de départ est "M", avoir un nouvel état dans une 4e colonne selon les probabilités de la matrice de transition. Je vous épargne tous les détails de ce que j'ai fait avant, mais voici où est mon problème:

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
#Creer la fonction de probabilite sachant l'etat a l'annee precedente
#Changement de colonnes dans la matrice d'etat
for(j in 4:ncol(ETATS)){
  #Changement de rangees dans la matrice d'etat
  for(i in 1:nrow(ETATS)){   
    #p.h pour savoir la probabilité
    p.h<-rmultinom(n=1,size=1, prob=?????
Comment faire pour s'assurer, sans écrire toutes les probabilités au long, que ma multinomiale ira chercher les bons chiffres dans ma matrice de transition? et de plus, si vous avez des pistes de comment construire ma fonction différemment, je ne suis même pas sûr que celle mentionnée plus haut fonctionne encore...