---
title: "Exemple de journal"
author: "Nathan Drogue"
date: '21 août 2024'
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

## Données d'entrée
Je n'ai actuellement pas de données particulièrement intéressantes à analyser. J'ai donc fait appel aux données météo France pour le département 63 comme base pour cet exercice. Elle peuvent être retrouvée [ici](https://www.data.gouv.fr/fr/datasets/r/d834495a-86d7-40e5-b663-db0b6cc34d84).

Voici un aperçu du jeu de données source :
```{r data, prompt=F, warning=F}
library(data.table)
raw_data <- data.table::fread("./Q_63_1871-1949_RR-T-Vent.csv")
head(raw_data)
```

## Température moyenne
Je m'intéresse en particulier à confronter la pluviométrie et l'altitude de ces différents sites. En plus, il faut que j'effectue une conversion sur la colonne temps qui n'est pas très lisible, et je ne veux que les données du XXème siècle et après...
```{r temp, prompt=F, warning=F}
library(dplyr)
rr_20eme <- raw_data %>%
  dplyr::select(NOM_USUEL, AAAAMMJJ, RR, ALTI) %>%
  dplyr::mutate(JOUR = as.Date.character(AAAAMMJJ, format = "%Y%m%d")) %>%
  dplyr::select(-AAAAMMJJ) %>%
  dplyr::filter(year(JOUR)>="1900")
head(rr_20eme)
```

Maintenant, j'aimerai savoir toutes les villes qui sont à ma disposition, ainsi que leurs cumul moyennes respectives par année.
```{r temp, prompt=F, warning=F}
print(unique(rr_20eme$NOM_USUEL))
rr_grouped <- rr_20eme %>%
  dplyr::mutate(ANNEE = year(JOUR)) %>% 
  dplyr::group_by(NOM_USUEL, ANNEE) %>%
  dplyr::summarise(RR_YCUM = sum(RR, na.rm = T),
                   ALTI = unique(ALTI))
head(rr_grouped)
```

Très bien, essayons maintenant de voir s'il pourrait exister un lien entre l'altitude d'un site et la pluviométrie qu'il a reçu.
```{r corr, warning=F, prompt=F}
library(tidyverse)
ggplot(rr_grouped) +
  aes(x = ALTI, y = RR_YCUM, color = ANNEE) +
  geom_point(alpha = .45) +
  scale_color_gradient(low="blue", high="red")
  labs(x = "Altitude", y = "Pluviométrie annuelle cumulée")
```
Bon, la répartition temporelle est très hétérogène, mais il semblerait qu'il y ait quelque chose à creuser. Je verrai ça demain... ;)