Discussion :

[justinquem]

Bonjour,

Je viens de représenter des boxplots avec en abscisse des variables qualitatives de type Culture (canne, maïs et blé) et en ordonnée des valeurs quantitatives d'émissions de CO2.
Pour chaque culture j'ai affiché les points représentant les différents résultats donnés dans différentes études.
Je voudrais maintenant afficher le nom de l'étude (en utilisant la fonction text) pour les maximum et minimum de chaque type de culture.
(je ne peux pas afficher tous les noms des études car cela rendrait le graphique très brouillon)

Comment faire? J'arrive à donner les références en comptant à la main les numéros des lignes à afficher mais je voudrais automatiser tout ça!!

Merci de votre aide,

Justin

[A. D.]

Bonjour,

Je ne suis pas sûre d'avoir bien compris ce que vous souhaitez faire...
Dans l'idéal, il est préférable de poster un court mais représentatif extrait de son jeu de données ainsi que le code utilisé (entre balises codes, bien sûr), cela facilité souvent la compréhension du problème

Sinon, si ce que vous souhaitez est automatiser la récupération des noms des études correspondant au minimum et au maximum de taux de CO2 pour chaque type de Culture, vous pouvez procéder comme ci-dessous (en adaptant à votre cas bien sûr) :

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
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
 
# creation d'un jeu de donnees pour tester le code
 
x<-100*rnorm(20)
y<-sample(c("canne","mais","blé"),20,replace=TRUE)
z<-paste("Etude",1:20,sep="")
 
M<-data.frame(CO2=x,Culture=y,Etude=z)  # si les donnees ne sont pas dans un data.frame, soit les passer en data.frame soit adapter un peu le code
 
 
# code utilise
 
tab<-table(M$Culture)
Cult<-names(tab)
 
lvs<-levels(M$Etude)
 
 
 
MinMax<-rep(0,6*2)
dim(MinMax)<-c(2,6)
 
colnames(MinMax)<-paste(c("Min","Max"),sort(rep(Cult,2)),sep="-")
 
k<-0
for(i in 1:length(tab))
{
	S<-M[which(M$Culture==Cult[i]),]
 
	MinMax[2,i+k]<-lvs[S$Etude[which(S$CO2==min(S[,1]))]]
	MinMax[2,i+k+1]<-lvs[S$Etude[which(S$CO2==max(S[,1]))]]
 
	MinMax[1,i+k]<-min(S[,1])
	MinMax[1,i+k+1]<-max(S[,1])
 
	k<-k+1
}

En espérant que cela vous aide à avancer...
Bonne continuation


Cordialement,

A.D.

[justinquem]

Voilà les données de mon graphique:

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
Ref                                          
Dumortier et al. 2009: 5    
EPA 2010             : 4   
JRC 2010 IMPACT      : 4   
CARB 2009            : 3                    
Tyner et al. 2010    : 3                           
EPA 2009             : 2                        
(Other)              :26
 
Culture
Ble  :11            
Canne: 7   
Mais :29
 
EmTot      
Min.   : 32.0  
1st Qu.: 74.5  
Median : 92.0  
Mean   :101.8  
3rd Qu.:127.5  
Max.   :214.0

Voilà mon code:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
 
grille <- read.csv(file="Decomposition_Ethanol.csv", sep=";", dec=",",na.strings = "NA")
 
plot(grille$Culture,grille$EmTot,las=2, main="Emissions totales en g CO2e / MJ (CAS+ACV)")
 
points(grille$Culture,grille$EmTot,col=c("red"),pch=19)
 
#lignes <- c(3,1,12,27,16,47,8,38,23)
 
text(grille[lignes,"Culture"],grille[lignes,"EmTot"],grille[lignes,"Ref"],pos=4,cex=0.9)

Il faut que je remplace la ligne "lignes", par le code que vous m'avez donné. Je vais essayer de le décrypter!!

[justinquem]

Ca y est j'ai réussi, merci A.D!

Bonne journée,

Justin.