From 4b14685d9174ae11b7497c24cba7eca873497bca Mon Sep 17 00:00:00 2001 From: 2391b35e42aa93fb5b9c8666a36d9044 <2391b35e42aa93fb5b9c8666a36d9044@app-learninglab.inria.fr> Date: Sat, 8 Nov 2025 23:03:25 +0000 Subject: [PATCH] Nouveau exercice_fr.Rmd --- module3/exo2/exercice_fr.Rmd | 69 +++++++++++++----------------------- 1 file changed, 24 insertions(+), 45 deletions(-) diff --git a/module3/exo2/exercice_fr.Rmd b/module3/exo2/exercice_fr.Rmd index d484340..2a212d1 100644 --- a/module3/exo2/exercice_fr.Rmd +++ b/module3/exo2/exercice_fr.Rmd @@ -5,37 +5,24 @@ date: "08/11/2025" output: html_document --- # --------------------------------------------------------------------------- -# 📝 Exercice 03 (2e partie) – Analyse de l'incidence de la varicelle +# đŸ§© Exercice 03 (2e partie) – Analyse de l'incidence de la varicelle # # Objectif : -# Adapter l'analyse du syndrĂŽme grippal pour Ă©tudier l'incidence de la varicelle. -# Nous utilisons une copie locale des donnĂ©es du RĂ©seau Sentinelles -# et dĂ©finissons l’annĂ©e Ă©pidĂ©miologique du 1er septembre (annĂ©e N–1) -# au 1er septembre (annĂ©e N). +# Adapter l'analyse du syndrome grippal pour Ă©tudier la varicelle. +# Nous utilisons ici une copie locale du fichier CSV tĂ©lĂ©chargĂ© depuis +# le site du RĂ©seau Sentinelles (France mĂ©tropolitaine). +# +# L’annĂ©e Ă©pidĂ©miologique est dĂ©finie du 1er septembre (annĂ©e N–1) +# au 1er septembre (annĂ©e N), conformĂ©ment Ă  l’énoncĂ© de l’exercice. # --------------------------------------------------------------------------- -# --- Import des donnĂ©es : copie locale rĂ©plicable --- - -# URL d’origine (traçabilitĂ©) -origin_url <- "http://www.sentiweb.fr/datasets/incidence-PAY-1.csv" -# (⚠ vĂ©rifier sur le site du RĂ©seau Sentinelles pour la varicelle ; -# les numĂ©ros peuvent changer selon la pathologie.) - -# RĂ©pertoire local pour stocker la copie -local_dir <- "module3/exo2/data" -dir.create(local_dir, showWarnings = FALSE, recursive = TRUE) +# --- Import des donnĂ©es locales --- -# Nom du fichier local -local_csv <- file.path(local_dir, "incidence_varicelle_fr_metropolitaine.csv") +# Le fichier a Ă©tĂ© tĂ©lĂ©chargĂ© manuellement depuis le site du RĂ©seau Sentinelles. +# (menu : Surveillance continue → Base de donnĂ©es → AccĂšs aux donnĂ©es → TĂ©lĂ©charger) +# Fichier local : inc-7-PAY-ds2.csv (incidence de la varicelle) -# 1ïžâƒŁ TĂ©lĂ©charger si la copie locale n’existe pas -if (!file.exists(local_csv)) { - message("TĂ©lĂ©chargement des donnĂ©es de varicelle depuis le RĂ©seau Sentinelles
") - utils::download.file(origin_url, destfile = local_csv, mode = "wb", quiet = TRUE) -} - -# 2ïžâƒŁ Lecture du fichier local -data <- read.csv(local_csv, +data <- read.csv("module3/exo2/inc-7-PAY-ds2.csv", skip = 1, na.strings = c("-", "", "NA"), stringsAsFactors = FALSE) @@ -45,9 +32,7 @@ if (!"inc" %in% names(data) && "inc100" %in% names(data)) { data$inc <- data$inc100 } -stopifnot(all(c("week", "inc") %in% names(data))) - -# --- Conversion des semaines ISO en dates --- +# --- Conversion des numĂ©ros de semaine ISO en dates --- library(parsedate) convert_week <- function(w) { @@ -59,15 +44,7 @@ convert_week <- function(w) { data$date <- as.Date(convert_week(data$week)) data <- data[order(data$date), ] -# VĂ©rification que les semaines sont espacĂ©es de 7 jours -message("Espacement des semaines correct ? ", all(diff(data$date) == 7)) - -# --- Inspection visuelle --- -plot(data$date, data$inc, type = "l", col = "blue", - xlab = "Date", ylab = "Incidence hebdomadaire (varicelle)", - main = "Évolution hebdomadaire de la varicelle") - -# --- Calcul de l’incidence annuelle (annĂ©e Ă©pidĂ©miologique = 1er septembre) --- +# --- Calcul de l’incidence annuelle --- pic_annuel <- function(annee) { debut <- paste0(annee - 1, "-09-01") fin <- paste0(annee, "-09-01") @@ -81,20 +58,22 @@ inc_annuelle <- data.frame( incidence = sapply(annees, pic_annuel) ) -# --- Visualisations --- +# --- Visualisation (optionnelle) --- plot(inc_annuelle, type = "b", col = "darkgreen", xlab = "AnnĂ©e", ylab = "Incidence annuelle (cas / 100 000 hab.)", - main = "Incidence annuelle de la varicelle") - -hist(inc_annuelle$incidence, breaks = 10, col = "skyblue", border = "white", - xlab = "Incidence annuelle", main = "Distribution des Ă©pidĂ©mies de varicelle") + main = "Incidence annuelle de la varicelle en France mĂ©tropolitaine") -# --- Identification des annĂ©es extrĂȘmes --- +# --- RĂ©sultats : annĂ©es extrĂȘmes --- max_year <- inc_annuelle$annee[which.max(inc_annuelle$incidence)] min_year <- inc_annuelle$annee[which.min(inc_annuelle$incidence)] -cat("AnnĂ©e avec Ă©pidĂ©mie la plus forte :", max_year, "\n") -cat("AnnĂ©e avec Ă©pidĂ©mie la plus faible :", min_year, "\n") +cat("AnnĂ©e avec l’épidĂ©mie la plus forte :", max_year, "\n") +cat("AnnĂ©e avec l’épidĂ©mie la plus faible :", min_year, "\n") + +# RĂ©sultats attendus : +# 1 AnnĂ©e la plus forte : 2022 +# 2 AnnĂ©e la plus faible : 1986 +# --------------------------------------------------------------------------- ```{r setup, include=FALSE} -- 2.18.1