version originale

.gitignore

@@ -5,6 +5,7 @@
 
 # Session Data files
 .RData
+*.Rproj
 
 # Example code in package build process
 *-Ex.R
@@ -36,3 +37,5 @@ vignettes/*.pdf
 # Shiny token, see https://shiny.rstudio.com/articles/shinyapps.html
 rsconnect/
 
+data/
+

README.md

@@ -1,2 +0,0 @@
-# premiere_neige_montreal_2019
-

nomtreal_neige.Rmd

@@ -0,0 +1,70 @@
+---
+title: "Date de la première bordée"
+output:
+  pdf_document: 
+    fig_caption: yes
+    keep_tex: yes
+  html_notebook: default
+---
+
+Source des données:
+
+- [Rapport de données quotidiennes (2004-2012)](https://climat.meteo.gc.ca/climate_data/daily_data_f.html?hlyRange=1953-01-01%7C2013-02-14&dlyRange=1941-09-01%7C2013-02-13&mlyRange=1941-01-01%7C2013-02-01&StationID=5415&Prov=QC&urlExtension=_f.html&searchType=stnName&optLimit=specDate&StartYear=1840&EndYear=2019&selRowPerPage=25&Line=2&searchMethod=contains&txtStationName=Montreal&timeframe=2&Day=12&Year=2004&Month=11#)
+- [Rapport de données quotidiennes (2013-2019)](https://climat.meteo.gc.ca/climate_data/daily_data_f.html?hlyRange=2013-02-13%7C2019-11-11&dlyRange=2013-02-14%7C2019-11-11&mlyRange=%7C&StationID=51157&Prov=QC&urlExtension=_f.html&searchType=stnName&optLimit=yearRange&StartYear=1840&EndYear=2019&selRowPerPage=25&Line=4&searchMethod=contains&txtStationName=Montreal&timeframe=2&Day=12&Year=2013&Month=11#)
+
+```{r message=FALSE, warning=FALSE}
+library(tidyverse)
+```
+
+Liste des fichiers
+```{r}
+fichiers <- list.files(pattern = "*.csv")
+```
+
+Specification des colonnes
+```{r}
+col_spec <- cols(
+  .default = col_double(),
+  `Nom de la Station` = col_character(),
+  `Date/Heure` = col_date()
+)
+```
+
+Lecture des données
+```{r}
+data_neige <- fichiers %>%
+  lapply(
+    FUN = function(x)
+      x %>%
+      read_csv(
+        locale = locale(decimal_mark = ","),
+        col_types = col_spec,
+        na = c("T", "E", "M", "<31")
+      )
+  ) %>%
+  bind_rows()
+```
+
+Extraction des jours avec plus de 10 cm de neige entre le mois de septembre et décembre
+```{r}
+premiere_bordee <-
+  data_neige %>% select(`Date/Heure`, `Année`, Mois, Jour, `Neige tot. (cm)`) %>%
+  filter(`Neige tot. (cm)` >= 10) %>%
+  filter(Mois >= 9) %>%
+  arrange(`Année`, Mois, Jour) %>%
+  group_by(`Année`) %>%
+  filter(row_number() == 1) %>%
+  transmute(date = ISOdate(0, Mois, Jour), neige_cm = `Neige tot. (cm)`)
+premiere_bordee
+```
+
+Visualisation
+```{r}
+ggplot2::ggplot(data = premiere_bordee, aes(x = `Année`, y = date)) +
+  geom_point(aes(size = neige_cm)) +
+  ggtitle("Date de la première bordée (10 cm) à Montréal")
+```
+
+
+
+