Discussion :

[Haw27]

Bonjour,
Je cherche un moyen de créer une petite fonction qui décrit une dynamique de 2 populations "co-évoluant".
Voilà le script :

Code :

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

1
2
3
4
5
6
7
pop<-function(n1,n2,r1,r2){
d1<-n1*r1+m
d2<-n2*r2+m
m<-n1+n2=1
plot(d1,d2)
}
pop(n1=1,n2=2,r1=2,r2=3)

J'aimerais dire que la valeur de "m" vérifie que "n1+n2=1". Evidemment, le script ne marche pas !
Merci beaucoup de votre aide.

[TimoP]

Salut,

Si tu veux faire des dynamiques de pop, il vaut mieux regarder le package simecol et deSolve

Pour ton script…

Tu utilises la valeur de m sans la définir. Si tu veux juste tester m = (n2+n1==1), quelque chose comme

Code :

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

m <- ifelse((n2+n1)==1,1,0)

le fera.

Mais je ne vois pas ce que tu veux faire avec ta fonction. Si c'est avoir la dynamique de pop, il faut une boucle qui te donne le nombre d'itérations…

[Haw27]

Merci,
alors oui, je me suis aperçu qu'il faut créer une itération pour que le graphe ce fasse.
Alors voilà mon nouveau script :

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
pop<-function(n1=0.5,n2=0.5,r1=0.1,r2=2,times=100,m){
pop1<-rep(0,times)
pop2<-rep(0,times)
for (i in 1:times){
d1<-n1*r1+m
d2<-n2*r2+m
m<-ifelse((n1+n2)==1,1,0)
pop1[i]<-n1<-n1+d1
pop2[i]<-n2<-n2+d2 
}
plot(c(1,times),range(c(pop1,pop2)),types="n")
lines(1:times,pop1,col="blue")
lines(1:times,pop2,col="red")
}
pop()

je suis arrivé à le faire fonctionner sans ce "m". L'idée est que x1 et x2 soit des fréquences. Donc m devrait vérifier que x1+x2=1. m serait (d'après le bouquin) la fitness moyenne du système.
Alors bien sûr, ici je pourrais faire m = r1*x1 + r2*x2. Mais par la suite, si je fais une fonction avec beaucoup plus de population ou avec une équation de dynamique plus complexe, ce sera long (et difficile) de résoudre l'équation.
Merci beaucoup.