diff --git a/.gitignore b/.gitignore index dc63652..388f14a 100644 --- a/.gitignore +++ b/.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/ + diff --git a/README.md b/README.md deleted file mode 100644 index 55de052..0000000 --- a/README.md +++ /dev/null @@ -1,2 +0,0 @@ -# premiere_neige_montreal_2019 - diff --git a/nomtreal_neige.Rmd b/nomtreal_neige.Rmd new file mode 100644 index 0000000..2b69ab2 --- /dev/null +++ b/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") +``` + + + +