Discussion :

[Kev_ser]

Bonjour,

Je chercher à simuler un modèle de marche aléatoire renforcée par arrête.

L'idée est que nous avons une bifurcation à deux branche dans laquelle des fourmis doivent choisir entre la branche droite et gauche.

La probabilité conditionnelle des Xi est donné par :

Xi = 1 si la fourmi choisie la branche droite et -1 si elle choisie la branche gauche

Zn est le nombre de fois que la branche droite a été choisie après n passages de fourmis.

alpha et c sont des paramètres qui, selon leurs valeurs, indiquent si la sélection (ou non ) d'une branche au bout d'un certain temps est lente ou rapide.

Voici ce que j'ai essayé :

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
16
17
18
19
20
21
 
 
K = 20 #nombre de trajectoires
N = 500 # durée d'une simulation
 
renforcement = function(alpha,c,Zn){
  return((1+(Zn/c))^alpha)
}
 
alpha = 10
c = 0.1
for (i in 1:K){
  Passages_droite= c()
  nbr_droite = 0
  for (j in 1:N){
    proba = renforcement(alpha,c,nbr_droite)/(renforcement(alpha,c,nbr_droite)+renforcement(alpha,c,j-nbr_droite))
    onpasse_droite = runif(1)<proba
    Passages_droite[j]=onpasse_droite
    nbr_droite=nbr_droite+onpasse_droite
  }
}

Mon code ne fonctionne pas ... avez-vous des idées sur comment l'améliorer?
Aussi, avez-vous une idée de comment représenter un graphe avec les Xi en fonction de n (en considérant la probabilité conditionnelle)?

Je vous remercie