J'aimerais tracer sur le même graphe un barplot représentant la pluviométrie et plusieurs courbes (représentant l'humidité du sol à différentes profondeurs). Toutes les données sont journalières. Voilà mes deux data frames :
- Humidité du sol :
Code : Tout sélectionner
head(umidade_data)
date S_10cm.1.1_pol S_20cm.1.1_pol S_30cm.1.1_pol S_40cm.1.1_pol S_60cm.1.1_pol S_100cm.1.1_pol
1 2014-08-14 0.1383 0.1731 0.2059 0.1880 0.1932 0.3341
2 2014-08-15 0.1439 0.1813 0.2095 0.1915 0.2009 0.3429
3 2014-08-16 0.1481 0.1832 0.2115 0.1943 0.2012 0.3450
4 2014-08-17 0.2731 0.3508 0.3240 0.2584 0.2085 0.3456
5 2014-08-18 0.2700 0.3463 0.3211 0.2623 0.2156 0.3463
6 2014-08-19 0.2661 0.3397 0.3168 0.2607 0.2143 0.3462
- Pluviométrie :
Code : Tout sélectionner
> head(rain_data)
date rain_data...2.
1 2014-10-03 0
2 2014-10-04 0
3 2014-10-07 0
4 2014-10-08 0
5 2014-10-09 0
6 2014-10-10 0
J'ai essayé avec les fonctions graphiques de base avec un premier barplot suivi des autres courbes superposées, les courbes d'affichent toutes mais les échelles des courbes sont différentes et il est donc impossible de comparer les valeurs :
Code : Tout sélectionner
plot_pluvio <- barplot(height = rain_data$rain_data...2., names.arg = rain_data$date, xlab="Date", ylab="Daily precipitation (mm)", col = "black")
#Plotar Umidade no mesmo grafico
par(new = TRUE)
plot(x = umidade_data$date, y = umidade_data$S_10cm.1.1_pol, xlab = "", ylab = "", type = "l", axes = FALSE, col = "green", lty=1)
par(new = TRUE)
plot(x = umidade_data$date, y = umidade_data$S_20cm.1.1_pol, xlab = "", ylab = "", type = "l", axes = FALSE, col = "blue", lty=1)
par(new = TRUE)
plot(x = umidade_data$date, y = umidade_data$S_30cm.1.1_pol, xlab = "", ylab = "", type = "l", axes = FALSE, col = "red", lty=1)
par(new = TRUE)
plot(x = umidade_data$date, y = umidade_data$S_40cm.1.1_pol, xlab = "", ylab = "", type = "l", axes = FALSE, col = "darkorchid1", lty=1)
par(new = TRUE)
plot(x = umidade_data$date, y = umidade_data$S_60cm.1.1_pol, xlab = "", ylab = "", type = "l", axes = FALSE, col = "darkorange", lty=1)
par(new = TRUE)
plot(x = umidade_data$date, y = umidade_data$S_100cm.1.1_pol, xlab = "", ylab = "", type = "l", axes = FALSE, col = "gold4", lty=1)
Je pense qu'il y a une possibilité avec lattice ou ggplot2 mais je n'arrive pas à les utiliser... Quelqu'un aurait-il une idée ?
Merci d'avance !
Pierre