Discussion :

[AgroInst]

Bonjour et bonne année à tous,

J'ai un dataframe de la sorte :

A	Température1	Température 2	Temps	Lieu
25	10	-5	10	A
50	15	0	11	A
30	12	2	12	A
100	18	3	13	A
60	19	0	14	A
15	10	0	10	B
100	9	5	11	B
80	13	-1	12	B
85	14	1	13	B
15	16	1	14	B
etc.	etc.	etc.	etc.	etc.

Je voudrais représenter pour chaque Lieu, mes variables A, Température1, Température2 en fonction du temps sur le même graphique.
J'aimerais avoir deux axes Y : un pour A; un pour les températures.

quelqu’un pourrait-il m'aider ?

En tâtonnant pour l'instant j'ai écrit ça (en commençant avec 2 variables représentées sur les 3) :
Le code ne fonctionne pas + je n'arrive pas à le rendre itératif pour répéter le travail pour chaque lieu de façon automatique .

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
for (i in levels (DF$Lieu))
       par (mar=c(4,4,3,5))
>     plot (A~Temps, subset (DF, Lieu==i), col="green4", main = i, axes= F, ylim=c(0,150), type = "b",lwd=3)
>     axis(2, ylim=c(0,150))
>     mtext("A", side=2)
>     par(new = T)
>     plot(Temp1~Temps, subset (DF, Lieu==i), col="black", ylim(-5,25), type = "l")
>     axis(4, ylim=c(-5,25))
>     mtext("Temperature1, °C", side=4)
>     axis(1)
>     mtext("Temps", side=1)
>     legend (1,1, legend=c("A","Temp"), text.col =c("green4","black"))

Merci par avance pour votre aide,

[AgroInst]

Après avoir cherché plus, j'ai réussi à faire tourner mon script (je le joins pour qu'il puisse servir à d'autres) :

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
for (i in levels (DF$Lieu)){
+ par(mar=c(4,4,3,5))
+ plot(T1~Temps, subset(DF, Lieu==i), col="green4", main = i, axes=F,xlab="",ylab="", ylim=c(0,160), type="b", lwd=3) 
+ axis(2, ylim=c(0,160),col="green4",col.axis="green4",at=seq(0, 160, by=50))
+ mtext("Texte",side=2,line=2.5,col="green4")
+ 
+ par(new = T)
+ plot(T2 ~ Temps, subset(DF, Lieu==i), col="red",axes=F,xlab="",ylab="",ylim=c(0,30), type="l")
+ axis(4 ,col="red",col.axis="red",at=seq(0, 25, by=5))
+ mtext("Texte",side=4,line=2.5,col="red")
+ 
+ axis(1, ylim=c(10,45),col="black",col.axis="black",at=seq(10, 45, by=2))
+ mtext("Temps",side=1,line=2.5,col="black")
+ 
+ }