Bonsoir, je sollicite votre aide.
Au fait, je veux tracer la distribution intra-saisonnière des pluies et températures moyennes dans une station en comparant les observations aux modèles. Mon soucis est que je veux définir mes propres couleurs (les observations je veux les tracer en noir). Ci joint le script utilisé et le fichier de données.
Excusez moi du dérangement, mais je suis vraiment bloqué à ce niveau.
Cordialement.
la structure des données[
> dput(df)
structure(list(Mois = structure(c(5L, 4L, 9L, 2L, 8L, 7L, 6L,
1L, 12L, 11L, 10L, 3L), .Label = c("Août", "April", "Dec.", "Févr",
"Janv", "Jllt.", "Juin", "Mai", "Mars", "Nov.", "Oct.", "Sept."
), class = "factor"), Obs = c(0.04, 2.02, 0.03, 0.03, 1.26, 28.18,
79.8, 141.66, 106.62, 31.01, 0.25, 0.69), Ccma = c(0, 0, 0, 0,
6.5, 10.7, 44.6, 116.2, 63.8, 3.7, 0, 0), CNRM = c(0, 0.1, 0,
0.1, 3.5, 28.5, 90, 146.6, 123.1, 33.8, 1, 0), IPSL = c(0, 0,
0, 0, 0.1, 37.9, 54.2, 96.6, 58.9, 35.5, 0, 0), MOHC = c(0, 0.1,
0, 1.2, 5.3, 32.1, 63.3, 99.1, 79.3, 35.2, 0.4, 8)), .Names = c("Mois",
"Obs", "Ccma", "CNRM", "IPSL", "MOHC"), class = "data.frame", row.names = c(NA,
-12L))
Le code exécuté est ci dessous:
#le package GSODR permet de télécharger les données des synopes transmis par SMT à la NOAA
library("GSODR")
#plyr permet de faire un resumé mensuel
library(plyr)
library(snpar)
library(openxlsx)
library(ggplot2)
library(tidyr)
library(tidyverse)
par(mfrow=c(1,2))
#df <- read.table("/Users/ewengallic/Downloads/RCM_kaorr.txt", header = T, encoding = "latin1")
df<-read.table("D:/SOURCE/Datataylor/RCM_linrr.txt", header = TRUE, sep = "")
df <-
df %>%
mutate(mois_num = 1:12) %>%
mutate(mois_nom = factor(mois_num, levels = 1:12, labels = as.character(df$Mois)))
#df %>% gather(variable, value, Obs:MOHC) %>%
# ggplot(data = ., aes(x = mois_num, y = value)) +
# geom_point() +
# geom_line() +
# facet_wrap(~variable) +
# labs(x = "Mois", y = "Valeur") +
# scale_x_continuous(breaks = 1:12, labels = levels(df$mois_nom))
df %>% gather(Model, value, Obs:MOHC) %>%
ggplot(data = ., aes(x = mois_num, y = value, colour = Model)) +
geom_point(size=2) +
geom_line(size=2) +
# facet_wrap(~variable) +
labs(x = "Mois", y = "Pluie (mm)") +
scale_colour_discrete("Modèle") +
scale_x_continuous(breaks = 1:12, labels = levels(df$mois_nom))+
theme_gray(base_size = 15)+
theme(panel.background = element_rect(fill = "white",
colour = "white",
size = 0.5, linetype = "solid")) +
theme( axis.line = element_line(colour = "black",
size = 1, linetype = "solid"))