From f87c696a5927aba49ed0b9697d4c2597b6c7c5e8 Mon Sep 17 00:00:00 2001 From: Tanguy BARTHELEMY Date: Wed, 13 Sep 2023 17:03:06 +0200 Subject: [PATCH 1/8] =?UTF-8?q?Ajout=20d'une=20nouvelle=20fiche=20sur=20le?= =?UTF-8?q?s=20donn=C3=A9es=20temporelles?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .../Fiche_donnees_temporelles.qmd | 562 ++++++++++++++++++ DESCRIPTION | 6 +- _quarto.yml | 1 + 3 files changed, 568 insertions(+), 1 deletion(-) create mode 100644 03_Fiches_thematiques/Fiche_donnees_temporelles.qmd diff --git a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd new file mode 100644 index 00000000..73ce4d89 --- /dev/null +++ b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd @@ -0,0 +1,562 @@ +--- +editor: + markdown: + wrap: 72 +--- + +# Manipuler des données temporelles {#timeseries} + +## Tâches concernées et recommandations + +L'utilisateur souhaite manipuler des dates (import et export du format +`Date`, manipulation et formatage de la date, différentes fréquences) ou +des séries temporelles (création d'objets `ts` en R). + +::: {.callout-recommandation .icon} +## Information sur les packages utilisés + +- Les packages `base` et `stats` proposent déjà une bonne diversité de + fonctionnalités pour traiter les dates et séries temporelles ; +- Le package `lubridate` ajoute des fonctionnalités supplémentaires + sur le traitement des dates ; +- Le package `parsedate` permet de contrôler les normes de format de + date. +- Le package `zoo` permet de gérer les données journalières ou + infra-journalières ; +- Pour les graphiques, ce sont les packages `ggplot2` et `dygraphs` + qui seront utilisés. + +Pour installer les packages : + +```{r, eval=FALSE} +install.packages(c("lubridate", "parsedate", "zoo")) +install.packages(c("ggplot2", "dygraphs")) +``` +::: + +::: {.callout-remarque .icon} +- La fonction `as.Date()`{.r} existe dans les packages `base` et + `zoo`, mais n'a pas la même fonctionnalité. A chaque fois qu'elle + sera utilisée, on précisera son package d'origine au moyen de `::` + (`base::as.Date()`{.r} ou `zoo::as.Date()`{.r}). +::: + +## Définitions + +### Les dates + +Une date en R signifie généralement un jour d'une année. Elle peut +s'écrire sous différents formats : + +| Format | Type | Exemple (Noël 2023) | +|------------------|------------|---------------------| +| Format européen | JJ/MM/AAAA | 25/12/2023 | +| Format américain | AAAA-MM-JJ | 2023-12-25 | + +Une date peut aussi désigner un horaire ou un moment précis. On peut +alors spécifier les heures, les minutes, les secondes. + +### Les séries temporelles + +Une série temporelle est une données indexée sur le temps. Généralement +on parle de séries mensuelles lorsque le pas entre chaque date est un +mois mais on parle aussi de données annuelles, trimestrielles ou +journalières voire plus haute fréquence. + +## Création de date + +### Le format `Date` en R + +Pour créer une date en R, il suffit de faire appel à la fonction +`as.Date()`{.r} du package `base`. Par défaut, l'argument format +correspond au format américain : + +```{r} +M.Drucker_birth <- base::as.Date(x = "1942-09-12") +``` + +Il suffit de le changer pour créer une date à partir d'un autre format : + +```{r} +noel_2023 <- base::as.Date(x = "25/12/2023", format = "%d/%m/%Y") +``` + +La liste des formats de données est disponible sur la page d'aide de +`strptime()`{.r} (accessible via `help(strptime)`{.r}). + +### Les autres formats de date `POSIXct` et `POSIXlt` + +Pour les horaires ou pour désigner un instant, il faut créer un objet de +classe `POSIXt`. Pour cela, on peut utiliser les fonctions +`as.POSIXct()`{.r} et `as.POSIXlt()`{.r}. + +- La fonction `as.POSIXct()`{.r} (**ct** pour **Calendar Time**) va + stocker la date sous la forme d'un nombre de secondes depuis le 1er + janvier 1970. +- La fonction `as.POSIXlt()`{.r} (**lt** pour **List Time**) va + stocker la date sous la forme d'une liste contenant tous les + composants de la date (année, mois, jour, heure...). + +```{r} +pied_sur_la_lune <- as.POSIXct(x = "1969-07-21 02:56:20", + format = "%Y-%m-%d %H:%M:%S", tz = "UTC") +``` + +### Les différents fuseaux horaires + +Ici le fuseau horaire (argument `tz`) est fixé sur UTC (Coordinated +Universal Time = Temps Universel Coordonné). Mais il est possible +d'appliquer différents fuseaux horaires à un évènement. + +Par exemple : + +Neil Armstrong a posé le pied sur la lune le +`r as.POSIXct(pied_sur_la_lune, tz = "Europe/Paris")` à Paris mais le +`r as.POSIXct(pied_sur_la_lune, tz = "America/Los_Angeles")` à Los +Angeles. + +```{r} +heure_en_france <- as.POSIXct(x = pied_sur_la_lune, tz = "Europe/Paris") +heure_los_angeles <- as.POSIXct(x = pied_sur_la_lune, tz = "America/Los_Angeles") +``` + +Pour connaitre la liste des différents fuseaux horaires, il faut appeler +la fonction `OlsonNames()`{.r} (du nom de l'Olson Database). + +### Gérer les changements d'heure + +En France, depuis 1979, on avance d'une heure en mars pour l'heure +d'été. Ainsi le fuseau horaire européen est un peu plus compliqué que le +fuseau UTC (qui est linéaire tout au long de l'année). + +Il suit 2 fuseaux : + +- le CET (*Central European Time*) correspond à l'heure *d'hiver*. + Cela vaut UTC+1. +- En été, on utilise le CEST (*Central European Summer Time*), qui est + à UTC+2. + +On entend aussi parler de GMT (Greenwich Mean Time) qui correspond à +l'heure de Greenwich. Aujourd'hui, c'est l'heure UTC (même fuseau que +GMT) qui est la référence. + +R gère en interne **automatiquement** les changements d'horaire à partir +des fuseaux horaires. + +- Par exemple en hiver (CET = UTC+1) : + +```{r} +# définition en date et heure locale avec le bon fuseau horaire +chute_mur_berlin <- as.POSIXct(x = "1989-11-09 18:00", tz = "Europe/Berlin") + +# Heure locale +print(chute_mur_berlin) + +# Heure UTC (exemple en islande) +print(as.POSIXct(chute_mur_berlin, tz = "UTC")) +``` + +- En été (CEST = UTC+2) : + +```{r} +# définition en date et heure locale avec le bon fuseau horaire +victoire_fifa_1998 <- as.POSIXct(x = "1998-07-12 21:00", tz = "Europe/Paris") + +# Heure locale +print(victoire_fifa_1998) + +# Heure UTC (exemple au Burkina Faso) +print(as.POSIXct(victoire_fifa_1998, tz = "UTC")) +``` + +## Autres fonctions + +### Formater une date avec la fonction `format()`{.r} + +La fonction `format()`{.r} est utile avec les dates car elle permet de +formater une date selon *n'importe quelle* représentation. + +Par exemple pour des dates : + +```{r} +# On prend la date d'aujourd'hui +format(Sys.Date(), format = "Nous sommes le %A %d %B %Y.") +format(Sys.Date(), format = "Date du jour : %x") +``` + +Par exemple pour des temps : + +```{r} +# On prend la date d'aujourd'hui +format(Sys.time(), format = "Nous sommes le %d %B %Y et il est %Hh%M et %Ss.") +format(Sys.time(), format = "Il s'est écoulé %ss depuis le 1er janvier 1970.") +format(Sys.time(), format = "Heure : %X") +``` + +La liste des formats de données est disponible sur la page d'aide de +`strptime()`{.r} (accessible via `help(strptime)`{.r}). + +::: {.callout-remarque .icon} +### Les normes ISO8601 et RFC3339 + +Les normes ISO8601 et RFC3339 sont des conventions de représentation des +dates. Selon ces 2 normes, certains formats de dates sont acceptés ou +non. + +Par exemple, voici quelques format représentant la date du 24 mai 2023 à +8h43 (UTC) : + +- `"2023-05-24T08:43:00Z"`{.r} est un format accepté par ces 2 normes. +- `"2023-05-24t08:43:00z"`{.r} est un format accepté uniquement par la + norme RFC3339. +- `"2023-05-24T08:43Z"`{.r} est un format accepté uniquement par la + norme ISO8601. + +Pour savoir quels formats sont acceptés par ces normes, une infographie +est disponible [ici](https://ijmacd.github.io/rfc3339-iso8601/) : +https://ijmacd.github.io/rfc3339-iso8601/. + +On peut aussi utiliser le package `parsedate` qui permet de lire une +date au format ISO8601 + +```{r, warning = FALSE} +library("parsedate") + +parse_iso_8601("2023-05-24T08:43:00Z") # Accepté par ISO8601 +parse_iso_8601("2023-05-24t08:43:00z") # Refusé par ISO8601 +parse_iso_8601("2023-05-24T08:43Z") # Accepté par ISO8601 +``` +::: + +### Paramètres régionaux + +Pour obtenir des dates, d'autres fonctions existent comme : + +- `Sys.Date()`{.r} pour connaitre la date du jour, +- `Sys.time()`{.r} pour connaitre l'horaire actuel (date + heure), +- `Sys.timezone()`{.r} pour le fuseau horaire actuel. + +Il peut être utile de vouloir changer les paramètres régionaux sous R. +Pour cela, il faut faire appel à la fonction `Sys.setlocale()`{.r} + +```{r} +# Paramètres locaux en France +Sys.setlocale("LC_TIME", "fr_FR.UTF-8") +format(Sys.time(), format = "%c") + +# Paramètres locaux aux USA +Sys.setlocale("LC_TIME", "en_US") +format(Sys.time(), format = "%c") +``` + +::: {.callout-remarque .icon} +La fonction `Sys.setlocale()`{.r} a un impact sur l'affichage mais pas +sur la valeur de l'objet. Ainsi, cela ne change pas le fuseau horaire +(par exemple). +::: + +## Création de vecteur de date + +Pour créer un vecteur de date, il faut utiliser la fonction `seq()`{.r}. + +Les arguments de cette fonction sont : + +- `from` : une date de départ +- `to` : une date d'arrivée +- `by` : un pas +- `length.out` : une longueur + +`from` est obligatoire et il faut préciser au moins 2 autres arguments +parmi les 3 restants. + +Exemple : + +```{r} +date1 <- base::as.Date("2016-02-29") +date2 <- base::as.Date("2021-10-02") +date3 <- base::as.Date("2023-08-15") + +seq(from = date1, to = date3, by = "year") +seq(from = date2, to = date3, by = "quarter") + +time1 <- as.POSIXct(x = "2023-05-26 15:00") +time2 <- as.POSIXct(x = "2023-05-26 20:59:47") +time3 <- as.POSIXct(x = "2023-05-26 21:00") + +seq(from = time1, to = time3, by = "hour") +seq(from = time2, to = time3, by = "sec") +``` + +## Autres objets de la famille des dates + +### Les objets `difftime` + +Le objets de classe `difftime` représentent des durées. + +```{r} +age <- Sys.Date() - M.Drucker_birth +print(age) + +difftime(Sys.Date(), M.Drucker_birth) + +# Changement d'heure +heure_hiver <- as.POSIXct(x = "2023-03-26 01:00", tz = "Europe/Paris") +heure_ete <- as.POSIXct(x = "2023-03-26 03:00", tz = "Europe/Paris") +heure_ete - heure_hiver + +# Voyage sur différents fuseaux horaires +decollage_paris <- as.POSIXct(x = "2023-07-20 10:30", tz = "Europe/Paris") +arrivee_toronto <- as.POSIXct(x = "2023-07-20 01:00 PM", + format = "%Y-%m-%d %I:%M %p", + tz = "America/Toronto") +arrivee_toronto - decollage_paris +``` + +Avec la fonction `units()`{.r}, on peut changer l'unité de la durée. + +```{r} +# En minutes +units(age) <- "mins" +print(age) + +# En semaines +units(age) <- "weeks" +print(age) +``` + +Le package `lubridate` propose aussi d'autres formatages des durées : + +```{r, warning = FALSE} +library("lubridate") + +time_length(age, unit = "year") +``` + +### Fonctionnalités avancées avec `lubridate` + +Le package `lubridate` propose des fonctions de lecture de date +semblables aux fonctions (`base::as.Date()`{.r}, `format()`{.r}, +`as.POSIXct()`{.r}) + +On peut résumer les équivalences des principales fonctions dans le +tableau ci-dessous : + +| Fonction `lubridate` | Équivalent R `base` | Type | +|---------------------------|---------------------------|------------------| +| `today()`{.r}, `now()`{.r} | `Sys.Date()`{.r}, `Sys.time()`{.r} | Création de date | +| `years()`{.r}, `weeks()`{.r}, `hours()`{.r}, `seconds()`{.r} ... | \- | Création de durée | +| `with_tz()`{.r} | `base::as.Date()`{.r} ou `as.POSIXct()`{.r} avec l'argument `ts` | Modification d'une date | +| `ymd()`{.r}, `ydm_hm()`{.r}, ... | `base::as.Date()`{.r} ou `as.POSIXct()`{.r} | Lecture de date | + +Les fonction de création de durée permettent de modifier et décaler les +dates. + +Une date au format `Date` est stocké sous la forme d'un nombre de jours +ainsi : + +```{r} +# En R base +Sys.Date() + 5L # Date dans 5 jours + +# Avec lubridate +Sys.Date() + days(5L) # Date dans 5 jours +``` + +Par exemple, si je veux avoir la même date il y a 5 ou 4 ans : + +```{r} +fev_29 <- base::as.Date("2016-02-29") + +# Il y a 5 ans, 29 fevrier 2011 n'existe pas : +fev_29 - 365.2425 * 5 # Pas exactement la même date +fev_29 - years(5L) # Calcul impossible + +# Il y a 4 ans : +fev_29 - 365.2425 * 4 +fev_29 - years(4L) +``` + +## Les séries temporelles en R + +### Créer une série temporelle + +Pour créer une série temporelle, on utilise la fonction `ts()`{.r} du +package `stats`. L'argument `frequency` est le nombre de période en 1 +an. + +```{r} +# Pour une série annuelle, frequency = 1L +serie_annuelle <- ts(1:20, start = 2003L, frequency = 1L) + +# Pour une série trimestrielle, frequency = 4L +serie_trimestrielle <- ts(21:40, end = 2023L, frequency = 4L) + +# Pour une série mensuelle, frequency = 12L +serie_mensuelle <- ts(41:60, start = 2020L, frequency = 12L) +``` + +::: {.callout-remarque .icon} +⚠️ Attention ⚠️ + +La fonction `ts()`{.r} du package `stats` n'est pas faite pour +fonctionner avec des fréquences non entières. +::: + +### Autres formats de dates + +Les objets `ts` fonctionnent aussi avec des dates mais dans un format +différent. Pour délimiter notre série temporelle, il faut utiliser les +arguments `start` et/ou `end` ainsi qu'une fréquence (entière). + +Les arguments `start` et `end` doivent être au format `AAAA`, +`c(AAAA, PP)` ou `AAAA + (PP - 1)/T`, avec : + +- `AAAA` l'année en 4 caractère +- et `PP` le numéro de la période (dépendant de la fréquence) +- et `T` la fréquence de la série + +Plus d'information sur la page d'aide de la fonction `start()`{.r} +(accessible via `help(start)`{.r}). + +Ces formats ne sont lisibles qu'avec la fréquence. + +Exemple : + +- Pour `frequency = 4L`{.r}, `2021L`{.r} et `c(2021L, 1L)`{.r} + représentent la même date (1er trimestre 2021) +- Pour `frequency = 12L`{.r}, `2019 + 3/12 = 2019.25`{.r} et + `c(2019L, 4L)`{.r} représentent la même date + +Ainsi `c(2020L, 5L)`{.r} désigne mai 2020 pour des séries mensuelles +(`frequency = 12L`{.r}) mais le 1er trimestre 2021 pour des séries +trimestrielles (`frequency = 4L`{.r}). + +### Fonctions des objets `ts` + +Pour obtenir des informations sur un objet `ts`, on peut utiliser les +fonctions suivantes : + +- `start()`{.r} : retourne la date de début de la série, +- `end()`{.r} : retourne la date de fin de la série, +- `frequency()`{.r} : retourne la fréquence de la série, +- et `time()`{.r} : retourne le vecteur de date indexant la série + temporelle + +```{r} +start(serie_mensuelle) +end(serie_mensuelle) +frequency(serie_mensuelle) +time(serie_mensuelle) +``` + +Les dates en output de `time()`{.r} sont au format ts et non au format +`Date`. Pour les convertir, il faut utiliser la fonction +`zoo::as.Date()`{.r} : + +```{r, warning = FALSE} +library("zoo") + +zoo::as.Date(time(serie_mensuelle)) +``` + +### Séries hautes fréquences avec `zoo` + +Le package `zoo` a d'autres utilités, notamment gérer les séries +haute-fréquence. + +On peut construire des séries journalières : + +```{r} +#data : https://www.letour.fr/fr/parcours-general +date_tour_de_france <- seq(from = as.Date("2023-06-29"), + to = as.Date("2023-07-21"), by = "day") +kilometre_etape <- c(182, 209, 193.5, 182, 163, 145, 170, 201, + 182.5, 0, 167.5, 180, 169, 138, 152, 179, + 0, 22.4, 166, 185, 173, 133.5, 115.5) + +tour_de_france_ts <- zoo(x = kilometre_etape, + order.by = date_tour_de_france) +``` + +On peut construire des séries infra-journalières (heure par heure ou +encore plus haute-fréquence) : + +```{r} +#data : https://joint-research-centre.ec.europa.eu/photovoltaic-geographical-information-system-pvgis/pvgis-tools/hourly-radiation_en +heures_journee <- seq(from = as.POSIXct("2016-07-01 00:10:00"), + to = as.POSIXct("2016-07-01 23:10:00"), by = "hours") +temperature <- c(19.42, 19.21, 18.99, 18.78, 19.71, 20.64, 21.57, 22.82, + 24.06, 25.31, 26.25, 27.19, 28.12, 28.44, 28.75, 29.06, + 28.22, 27.38, 26.54, 24.81, 23.08, 21.35, 21.02, 20.69) + +temp_ts <- zoo(x = temperature, + order.by = heures_journee) +``` + +## Afficher des séries temporelles + +Plusieurs packages permettent l'affichage des séries temporelles. + +### Classique `plot()`{.r} + +La fonction `plot()`{.r} du package `graphics` (méthode `plot.ts()`{.r} +du package `stats`) permet d'e multiples visualisations'afficher +simplement des séries temporelles : + +```{r plot} +plot(x = tour_de_france_ts, col = "blue", + ylab = "Distance (en km)", xlab = "Time", + main = "Distance parcourue par jour (TDF 2023)") +plot(x = temp_ts, col = "red", + main = "Temperature") +``` + +### Package `ggplot2` + +Le package `ggplot2` propose une grande variété de graphes. Il est +nécessaire au préalable de convertir l'objet en `data.frame` pour +construire le graphique. + +```{r ggplot2, warning = FALSE} +library("ggplot2") + +sunspots_df <- data.frame(date = time(sunspots), value = sunspots) + +ggplot(sunspots_df, aes(x=date, y=value)) + + geom_line() +``` + +Pour plus d'information sur l'utilisation de `ggplot2`, voir la [fiche +sur les +graphiques](https://www.book.utilitr.org/03_fiches_thematiques/fiche_graphiques) + +### Package `dygraphs` + +Le package `dygraphs` propose aussi des graphiques pour séries +temporelles. L'avantage de ce package est l'interactivité et la +possibilité de zoomer dans les graphiques. + +```{r dygraphs, warning = FALSE} +library("dygraphs") + +dygraph(temp_ts) +``` + +On peut aussi afficher plusieurs courbes sur le même graphique : + +```{r dygraph_multiple_plot} +dygraph(cbind(fdeaths, mdeaths)) +``` + +## Pour en savoir plus + +- [Formats acceptés par les normes ISO8601 et + RFC3339](https://github.com/IJMacD/rfc3339-iso8601) ; +- [Gallerie de graphique avec + `dygraph`](https://r-graph-gallery.com/317-time-series-with-the-dygraphs-library.html) + ; +- [Manuel d'utilisation package + `dygraph`](https://rstudio.github.io/dygraphs/) ; +- [Les séries temporelles en + `ggplot2`](https://r-graph-gallery.com/279-plotting-time-series-with-ggplot2.html) + ; diff --git a/DESCRIPTION b/DESCRIPTION index 079eecc5..e91a38e9 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -67,7 +67,11 @@ Imports: renv, remotes, targets, - lintr + lintr, + parsedate, + lubridate, + zoo, + dygraphs Suggests: cranlogs, haven, diff --git a/_quarto.yml b/_quarto.yml index 82919af4..49b1deaa 100644 --- a/_quarto.yml +++ b/_quarto.yml @@ -68,6 +68,7 @@ book: - 03_Fiches_thematiques/Fiche_donnees_textuelles.qmd - 03_Fiches_thematiques/Fiche_survey.qmd - 03_Fiches_thematiques/Fiche_donnees_spatiales.qmd + - 03_Fiches_thematiques/Fiche_donnees_temporelles.qmd - 03_Fiches_thematiques/Fiche_analyse_de_donnees.qmd - part: "Produire des sorties avec R" chapters: From 354a9830e79261776489911ee062f3525bb2f133 Mon Sep 17 00:00:00 2001 From: Lino Galiana Date: Fri, 22 Sep 2023 17:31:10 +0200 Subject: [PATCH 2/8] Update 03_Fiches_thematiques/Fiche_donnees_temporelles.qmd Co-authored-by: Romain Lesur --- 03_Fiches_thematiques/Fiche_donnees_temporelles.qmd | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd index 73ce4d89..784abd3b 100644 --- a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd +++ b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd @@ -20,7 +20,7 @@ des séries temporelles (création d'objets `ts` en R). - Le package `lubridate` ajoute des fonctionnalités supplémentaires sur le traitement des dates ; - Le package `parsedate` permet de contrôler les normes de format de - date. + date ; - Le package `zoo` permet de gérer les données journalières ou infra-journalières ; - Pour les graphiques, ce sont les packages `ggplot2` et `dygraphs` From a956ad6a2f6fed0eb3303ef559d097858d9d8950 Mon Sep 17 00:00:00 2001 From: Tanguy Barthelemy <112711017+TanguyBarthelemy@users.noreply.github.com> Date: Tue, 26 Sep 2023 10:25:53 +0200 Subject: [PATCH 3/8] Update 03_Fiches_thematiques/Fiche_donnees_temporelles.qmd Co-authored-by: Romain Lesur --- 03_Fiches_thematiques/Fiche_donnees_temporelles.qmd | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd index 784abd3b..a885b83a 100644 --- a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd +++ b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd @@ -91,7 +91,7 @@ classe `POSIXt`. Pour cela, on peut utiliser les fonctions `as.POSIXct()`{.r} et `as.POSIXlt()`{.r}. - La fonction `as.POSIXct()`{.r} (**ct** pour **Calendar Time**) va - stocker la date sous la forme d'un nombre de secondes depuis le 1er + stocker la date sous la forme d'un nombre de secondes depuis le 1^er^ janvier 1970. - La fonction `as.POSIXlt()`{.r} (**lt** pour **List Time**) va stocker la date sous la forme d'une liste contenant tous les From abc98edbb2b4375810eb583286549c8649316618 Mon Sep 17 00:00:00 2001 From: Tanguy Barthelemy <112711017+TanguyBarthelemy@users.noreply.github.com> Date: Tue, 26 Sep 2023 10:28:14 +0200 Subject: [PATCH 4/8] Update 03_Fiches_thematiques/Fiche_donnees_temporelles.qmd Co-authored-by: Romain Lesur --- 03_Fiches_thematiques/Fiche_donnees_temporelles.qmd | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd index a885b83a..81c8c8ee 100644 --- a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd +++ b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd @@ -348,7 +348,7 @@ tableau ci-dessous : | `with_tz()`{.r} | `base::as.Date()`{.r} ou `as.POSIXct()`{.r} avec l'argument `ts` | Modification d'une date | | `ymd()`{.r}, `ydm_hm()`{.r}, ... | `base::as.Date()`{.r} ou `as.POSIXct()`{.r} | Lecture de date | -Les fonction de création de durée permettent de modifier et décaler les +Les fonctions de création de durée permettent de modifier et décaler les dates. Une date au format `Date` est stocké sous la forme d'un nombre de jours From 92275a7a3116370530d264e298cc9cef4dda82fe Mon Sep 17 00:00:00 2001 From: Tanguy Barthelemy <112711017+TanguyBarthelemy@users.noreply.github.com> Date: Tue, 26 Sep 2023 10:29:33 +0200 Subject: [PATCH 5/8] =?UTF-8?q?typo=20stock=C3=A9E?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Co-authored-by: Romain Lesur --- 03_Fiches_thematiques/Fiche_donnees_temporelles.qmd | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd index 81c8c8ee..dc968b9e 100644 --- a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd +++ b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd @@ -351,7 +351,7 @@ tableau ci-dessous : Les fonctions de création de durée permettent de modifier et décaler les dates. -Une date au format `Date` est stocké sous la forme d'un nombre de jours +Une date au format `Date` est stockée sous la forme d'un nombre de jours ainsi : ```{r} From 57e1eac8ee5797379603bbfd1815318e25a75dce Mon Sep 17 00:00:00 2001 From: Tanguy Barthelemy <112711017+TanguyBarthelemy@users.noreply.github.com> Date: Tue, 26 Sep 2023 10:30:06 +0200 Subject: [PATCH 6/8] exposant Co-authored-by: Romain Lesur --- 03_Fiches_thematiques/Fiche_donnees_temporelles.qmd | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd index dc968b9e..d8126c21 100644 --- a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd +++ b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd @@ -428,7 +428,7 @@ Exemple : `c(2019L, 4L)`{.r} représentent la même date Ainsi `c(2020L, 5L)`{.r} désigne mai 2020 pour des séries mensuelles -(`frequency = 12L`{.r}) mais le 1er trimestre 2021 pour des séries +(`frequency = 12L`{.r}) mais le 1^er^ trimestre 2021 pour des séries trimestrielles (`frequency = 4L`{.r}). ### Fonctions des objets `ts` From 2e238d93ea3b54ba374431cad021c65ad8be03fc Mon Sep 17 00:00:00 2001 From: Tanguy BARTHELEMY Date: Tue, 26 Sep 2023 11:35:06 +0200 Subject: [PATCH 7/8] correct some typos and accept review --- .../Fiche_donnees_temporelles.qmd | 33 +++++++++---------- 1 file changed, 16 insertions(+), 17 deletions(-) diff --git a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd index a885b83a..22868d6a 100644 --- a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd +++ b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd @@ -34,9 +34,9 @@ install.packages(c("ggplot2", "dygraphs")) ``` ::: -::: {.callout-remarque .icon} +::: {.callout-caution} - La fonction `as.Date()`{.r} existe dans les packages `base` et - `zoo`, mais n'a pas la même fonctionnalité. A chaque fois qu'elle + `zoo`, mais n'a pas la même fonctionnalité. À chaque fois qu'elle sera utilisée, on précisera son package d'origine au moyen de `::` (`base::as.Date()`{.r} ou `zoo::as.Date()`{.r}). ::: @@ -58,7 +58,7 @@ alors spécifier les heures, les minutes, les secondes. ### Les séries temporelles -Une série temporelle est une données indexée sur le temps. Généralement +Une série temporelle est une données indexées sur le temps. Généralement on parle de séries mensuelles lorsque le pas entre chaque date est un mois mais on parle aussi de données annuelles, trimestrielles ou journalières voire plus haute fréquence. @@ -72,7 +72,7 @@ Pour créer une date en R, il suffit de faire appel à la fonction correspond au format américain : ```{r} -M.Drucker_birth <- base::as.Date(x = "1942-09-12") +M_Drucker_birth <- base::as.Date(x = "1942-09-12") ``` Il suffit de le changer pour créer une date à partir d'un autre format : @@ -120,7 +120,7 @@ heure_en_france <- as.POSIXct(x = pied_sur_la_lune, tz = "Europe/Paris") heure_los_angeles <- as.POSIXct(x = pied_sur_la_lune, tz = "America/Los_Angeles") ``` -Pour connaitre la liste des différents fuseaux horaires, il faut appeler +Pour connaître la liste des différents fuseaux horaires, il faut appeler la fonction `OlsonNames()`{.r} (du nom de l'Olson Database). ### Gérer les changements d'heure @@ -203,10 +203,10 @@ Les normes ISO8601 et RFC3339 sont des conventions de représentation des dates. Selon ces 2 normes, certains formats de dates sont acceptés ou non. -Par exemple, voici quelques format représentant la date du 24 mai 2023 à +Par exemple, voici quelques formats représentant la date du 24 mai 2023 à 8h43 (UTC) : -- `"2023-05-24T08:43:00Z"`{.r} est un format accepté par ces 2 normes. +- `"2023-05-24T08:43:00Z"`{.r} ou `"2023-05-24T08:43:00+08:00"`{.r} sont des formats acceptés par ces 2 normes. - `"2023-05-24t08:43:00z"`{.r} est un format accepté uniquement par la norme RFC3339. - `"2023-05-24T08:43Z"`{.r} est un format accepté uniquement par la @@ -222,7 +222,7 @@ date au format ISO8601 ```{r, warning = FALSE} library("parsedate") -parse_iso_8601("2023-05-24T08:43:00Z") # Accepté par ISO8601 +parse_iso_8601("2023-05-24T08:43:00+08:00") # Accepté par ISO8601 parse_iso_8601("2023-05-24t08:43:00z") # Refusé par ISO8601 parse_iso_8601("2023-05-24T08:43Z") # Accepté par ISO8601 ``` @@ -294,10 +294,10 @@ seq(from = time2, to = time3, by = "sec") Le objets de classe `difftime` représentent des durées. ```{r} -age <- Sys.Date() - M.Drucker_birth +age <- Sys.Date() - M_Drucker_birth print(age) -difftime(Sys.Date(), M.Drucker_birth) +difftime(Sys.Date(), M_Drucker_birth) # Changement d'heure heure_hiver <- as.POSIXct(x = "2023-03-26 01:00", tz = "Europe/Paris") @@ -395,7 +395,7 @@ serie_trimestrielle <- ts(21:40, end = 2023L, frequency = 4L) serie_mensuelle <- ts(41:60, start = 2020L, frequency = 12L) ``` -::: {.callout-remarque .icon} +::: {.callout-caution} ⚠️ Attention ⚠️ La fonction `ts()`{.r} du package `stats` n'est pas faite pour @@ -415,7 +415,7 @@ Les arguments `start` et `end` doivent être au format `AAAA`, - et `PP` le numéro de la période (dépendant de la fréquence) - et `T` la fréquence de la série -Plus d'information sur la page d'aide de la fonction `start()`{.r} +Plus d'informations sur la page d'aide de la fonction `start()`{.r} (accessible via `help(start)`{.r}). Ces formats ne sont lisibles qu'avec la fréquence. @@ -500,8 +500,7 @@ Plusieurs packages permettent l'affichage des séries temporelles. ### Classique `plot()`{.r} La fonction `plot()`{.r} du package `graphics` (méthode `plot.ts()`{.r} -du package `stats`) permet d'e multiples visualisations'afficher -simplement des séries temporelles : +du package `stats`) permet d'afficher simplement des séries temporelles : ```{r plot} plot(x = tour_de_france_ts, col = "blue", @@ -513,7 +512,7 @@ plot(x = temp_ts, col = "red", ### Package `ggplot2` -Le package `ggplot2` propose une grande variété de graphes. Il est +Le package `ggplot2` propose une grande variété de graphiques. Il est nécessaire au préalable de convertir l'objet en `data.frame` pour construire le graphique. @@ -526,7 +525,7 @@ ggplot(sunspots_df, aes(x=date, y=value)) + geom_line() ``` -Pour plus d'information sur l'utilisation de `ggplot2`, voir la [fiche +Pour plus d'informations sur l'utilisation de `ggplot2`, voir la [fiche sur les graphiques](https://www.book.utilitr.org/03_fiches_thematiques/fiche_graphiques) @@ -552,7 +551,7 @@ dygraph(cbind(fdeaths, mdeaths)) - [Formats acceptés par les normes ISO8601 et RFC3339](https://github.com/IJMacD/rfc3339-iso8601) ; -- [Gallerie de graphique avec +- [Galerie de graphique avec `dygraph`](https://r-graph-gallery.com/317-time-series-with-the-dygraphs-library.html) ; - [Manuel d'utilisation package From e80bb5e465e93aa4139773a672c235d8e4815239 Mon Sep 17 00:00:00 2001 From: Lino Galiana Date: Mon, 4 Dec 2023 18:16:51 +0100 Subject: [PATCH 8/8] Update Fiche_donnees_temporelles.qmd --- .../Fiche_donnees_temporelles.qmd | 19 +++++++------------ 1 file changed, 7 insertions(+), 12 deletions(-) diff --git a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd index ac3a5b1e..23261247 100644 --- a/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd +++ b/03_Fiches_thematiques/Fiche_donnees_temporelles.qmd @@ -1,9 +1,3 @@ ---- -editor: - markdown: - wrap: 72 ---- - # Manipuler des données temporelles {#timeseries} ## Tâches concernées et recommandations @@ -12,7 +6,7 @@ L'utilisateur souhaite manipuler des dates (import et export du format `Date`, manipulation et formatage de la date, différentes fréquences) ou des séries temporelles (création d'objets `ts` en R). -::: {.callout-recommandation .icon} +::: {.callout-important} ## Information sur les packages utilisés - Les packages `base` et `stats` proposent déjà une bonne diversité de @@ -28,14 +22,15 @@ des séries temporelles (création d'objets `ts` en R). Pour installer les packages : -```{r, eval=FALSE} +```{r} +#| eval: False install.packages(c("lubridate", "parsedate", "zoo")) install.packages(c("ggplot2", "dygraphs")) ``` ::: ::: {.callout-caution} -- La fonction `as.Date()`{.r} existe dans les packages `base` et +La fonction `as.Date()`{.r} existe dans les packages `base` et `zoo`, mais n'a pas la même fonctionnalité. À chaque fois qu'elle sera utilisée, on précisera son package d'origine au moyen de `::` (`base::as.Date()`{.r} ou `zoo::as.Date()`{.r}). @@ -196,7 +191,7 @@ format(Sys.time(), format = "Heure : %X") La liste des formats de données est disponible sur la page d'aide de `strptime()`{.r} (accessible via `help(strptime)`{.r}). -::: {.callout-remarque .icon} +::: {.callout-note} ### Les normes ISO8601 et RFC3339 Les normes ISO8601 et RFC3339 sont des conventions de représentation des @@ -249,7 +244,7 @@ Sys.setlocale("LC_TIME", "en_US") format(Sys.time(), format = "%c") ``` -::: {.callout-remarque .icon} +::: {.callout-note} La fonction `Sys.setlocale()`{.r} a un impact sur l'affichage mais pas sur la valeur de l'objet. Ainsi, cela ne change pas le fuseau horaire (par exemple). @@ -396,7 +391,7 @@ serie_mensuelle <- ts(41:60, start = 2020L, frequency = 12L) ``` ::: {.callout-caution} -⚠️ Attention ⚠️ +## ⚠️ Attention ⚠️ La fonction `ts()`{.r} du package `stats` n'est pas faite pour fonctionner avec des fréquences non entières.