1+ # nolint start
2+
3+ # Pratique 2
4+ # Activité 1
5+
6+ # étape : indiquez votre numéro de chambre
7+ room_number <- # <COMPLETE> remplacer par 1/2/3/4
8+
9+ # Charger les paquets -----------------------------------------------------------
10+ library(epiparameter )
11+ library(EpiNow2 )
12+ library(tidyverse )
13+
14+
15+ # Lire les cas signalés -----------------------------------------------------
16+ # étape : Coller les liens URL sous forme de chaîne pour lire les données d'entrée.
17+
18+ # pour covid
19+ dat <- read_rds(
20+ " paste/link/url/here/covid" # <COMPLETE>
21+ ) %> %
22+ dplyr :: select(date , confirm )
23+ # ou
24+ # for ebola
25+ dat <- read_rds(
26+ " paste/link/url/here/ebola" # <COMPLETE>
27+ ) %> %
28+ dplyr :: select(date , confirm = cases )
29+
30+
31+ # Vérifier le format des données d'incidence ----------------------------------------
32+ # étape : Vérifier les noms des colonnes dans les données d'incidence
33+ # correspondre aux exigences d'EpiNow2 concernant les noms des colonnes : date et confirmer
34+
35+ dat
36+
37+ # Définir un temps de génération de {epiparameter} à {EpiNow2} ---------------
38+
39+ # étape : accéder à un intervalle sériel
40+ dat_serialint <- epiparameter :: epiparameter_db(
41+ # <COMPLETE>
42+ )
43+
44+ # étape : extraire les paramètres de l'objet {epiparameter}
45+ dat_serialint_params <- epiparameter :: # <COMPLETE>
46+
47+ # étape : adapter {epiparamètre} à l'interface de distribution {EpiNow2}
48+ dat_generationtime <- EpiNow2 :: # <COMPLETE>
49+
50+
51+ # Définir les délais entre l'infection et la déclaration du cas pour {EpiNow2} -----------
52+
53+ # étape : définir le délai entre l'apparition des symptômes et la déclaration du cas
54+ # Vous devez interprétez la description fournie dans le tableau Entrées
55+
56+ dat_reportdelay <- EpiNow2 :: # <COMPLETE>
57+
58+ # étape : définir un délai entre l'infection et l'apparition des symptômes
59+ dat_incubationtime <- epiparameter :: epiparameter_db(
60+ # <COMPLETE>
61+ )
62+
63+ # étape : période d'incubation : extraire les paramètres de distribution
64+ dat_incubationtime_params <- epiparameter :: # <COMPLETE>
65+
66+ # étape : période d'incubation : discrétiser et extraire la valeur maximale (p = 99 %)
67+ dat_incubationtime_max <- dat_incubationtime %> % # <COMPLETE>
68+
69+ # étape : période d'incubation : adaptation à l'interface de distribution {EpiNow2}
70+ dat_incubationtime_epinow <- EpiNow2 :: # <COMPLETE>
71+
72+ # étape : impression des données requises
73+ dat_generationtime
74+ dat_reportdelay
75+ dat_incubationtime_epinow
76+
77+
78+ # Définir le nombre de cœurs parallèles pour {EpiNow2} --------------------------
79+ # étape : exécuter cette étape de configuration
80+
81+ withr :: local_options(list (mc.cores = parallel :: detectCores() - 1 ))
82+
83+
84+ # Estimer la transmission à l'aide de EpiNow2::epinow() ---------------------------
85+ # étape : exécuter epinow() en utilisant les cas d'incidence et les délais requis
86+ # avec les fonctions EpiNow2::*_opts() pour le temps de génération, les délais et stan.
87+
88+ estimates <- EpiNow2 :: epinow(
89+ data = dat ,
90+ # <COMPLETE>
91+ stan = EpiNow2 :: stan_opts(samples = 1000 , chains = 3 )
92+ )
93+
94+
95+ # Imprimer les résultats du graphique et du tableau récapitulatif ------------------------------------
96+ # étape : coller les résultats du graphique et du tableau dans le document partagé
97+
98+ summary(estimates )
99+ plot(estimates )
100+
101+ # étape : répondre aux questions dans le document
102+
103+ # nolint end
0 commit comments