diff --git a/module3/exo1/mon-analyse.Rmd b/module3/exo1/mon-analyse.Rmd index 9c8d0d66ee85da8c8f4bf459a21292dde530dbb8..c61a98003b37b9dbcc1e69cf47b5e5e9e9d47269 100644 --- a/module3/exo1/mon-analyse.Rmd +++ b/module3/exo1/mon-analyse.Rmd @@ -9,6 +9,7 @@ output: html_document knitr::opts_chunk$set(echo = TRUE) ``` +# importer les données Je télécharge les données du site sentinelles en ignorant la première ligne qui contient les commentaires @@ -34,3 +35,56 @@ class(data$week) class(data$inc) ``` Ce sont bien des entiers, tout va bien + +# vérification et inspection des données + +## conversion des dates +conversion des dates grâce à la bibliothèque parsedate : 199501 en 1995-01-02 UT (date du lundi) + +```{r} +library(parsedate) +convert_week<-function(w){ + ws<-paste(w)#convertir en chaine de caracteres + iso<-paste0(substring(ws,1,4),"-W",substring(ws,5,6))#1995-W01 + as.character(parse_iso_8601(iso))#1995-01-02 UT +} +``` + +mettre les dates à la norme iso dans une nouvelle colonne date +```{r} +data$date<-as.Date(convert_week(data$week)) +``` + +on vérifie que c'est bien de classe date + +```{r} +class(data$date) +``` + +mettre les données dans l'ordre chronologique +```{r} +data<-data[order(data$date),] +``` + +vérification +```{r} +head(data) +``` + +vérification que les données ont bien été faites toutes les semaines (toutes les 7 jours) + +```{r} +all(diff(data$date)==7) +``` +## inspection + +regardons à quoi ressemblent les données +```{r} +with(data,plot(date,inc,type="l")) +``` +zoom sur les 200 derniers points (4 ans) +```{r} +with(tail(data,200),plot(date,inc,type="l")) +``` + +on voit un pic en début d'année (hiver)