vérification et inspection données

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)