Bonjour, je fais appel à vos services car je n’arrive pas à résoudre un problème de graphique.
Je dispose d’un tableau csv avec des colonnes : « modalite », « duree_retour », « session ».
Les données de « modalite » ont deux formes possibles : « temoin » et « contaminé » (évidement je n’ai pas mis d’accent dans le tableau). Les données de « session » sont de 4 types : « session1 » , « session2 » , « session3 », « session4 ». Et les données de « duree_retour » sont sous forme numérique :ex 13,5 ou 300,7 (exprimé en minutes – dans cette colonne certaines valeurs sont +Inf quand il n’y a pas de retour). J’aimerais réaliser un graphique qui me donne le taux de retour en fonction du temps et de la modalité et tout cela, par session. J’ai un script qui me donne déjà le taux de retour en fonction du temps par modalité. Mais je n’arrive pas à intégrer à ce graph les sessions. Voici le script :
(data_relacher_retour_marquage est le fichier csv)
Je pense que j’essaye de faire un graphique avec 2 variables qualitatives, je ne suis pas sûre de cela, mais je vous remercie d’avance si vous avez une idée pour modifier le script et me donner un sérieux coup de pouce. Bien cordialement.
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14 data_relacher_retour_marquage=merge(data_relacher_retour_unique,data_marquage,by.x="idtag_court",by.y="numtag_court",all.x=TRUE) data_relacher_retour_marquage$duree_retour[is.na(data_relacher_retour_marquage$duree_retour)]=+Inf duree_temoin=data_relacher_retour_marquage[data_relacher_retour_marquage$modalite=="temoin","duree_retour"] duree_contamine=data_relacher_retour_marquage[data_relacher_retour_marquage$modalite=="contamine","duree_retour"] dev.new() par(mfrow=c(1,2),mar=c(4.25,4.25,0.25,0.25)) plot(ecdf(duree_temoin),col="grey",xlab="Time after release (min)", ylab="Cumulative homing probability",xlim=c(0,500),ylim=c(0,0.25),main="",verticals = TRUE,do.points=FALSE,lwd=3) lines(ecdf(duree_contamine),verticals = TRUE,do.points=FALSE,lwd=3) legend("topleft",c("temoin","contamine"),col=c("grey","black"),lty=c(1,1),lwd=c(3,3),cex=0.75) par(mar=c(4.25,0.5,0.25,0.25),las=2) plot(ecdf(duree_temoin),col="grey",xlab="Time after release (min)", ylab="Cumulative homing probability",xlim=c(600,7000),ylim=c(0,0.25),main="",yaxt="n",verticals = TRUE,do.points=FALSE,lwd=3) lines(ecdf(duree_contamine),verticals = TRUE,do.points=FALSE,lwd=3)
Partager