|
| 1 | +# Activity 5 |
| 2 | + |
| 3 | +# étape: donner le numéro du groupe |
| 4 | +room_number <- #<À COMPLETER> replace with 1/2/3/4 |
| 5 | + |
| 6 | +# charger les packages --------------------------------------------------------- |
| 7 | +library(cfr) |
| 8 | +library(epiparameter) |
| 9 | +library(tidyverse) |
| 10 | + |
| 11 | + |
| 12 | +# Importer les données d'entrées------------------------------------------------ |
| 13 | +# étape: copier et coller le lien fourni dans le fichier google doc. |
| 14 | +disease_dat <- readr::read_rds( |
| 15 | + #<À COMPLETER> |
| 16 | +) |
| 17 | + |
| 18 | +disease_dat |
| 19 | + |
| 20 | + |
| 21 | +# Créer un objet de type incidence --------------------------------------------- |
| 22 | +# étape: compléter les arguments nécessaires pour obtenir l'incidence des cas et |
| 23 | +# de décès. |
| 24 | +disease_incidence <- disease_dat %>% |
| 25 | + incidence2::incidence( |
| 26 | + #<À COMPLETER> |
| 27 | + ) |
| 28 | + |
| 29 | +plot(disease_incidence) |
| 30 | + |
| 31 | + |
| 32 | +# Confirmer le format des données d'entrées pour {cfr} ------------------------- |
| 33 | + |
| 34 | +# étape: vérifier si les noms de colonnes dans les données d'incidence sont |
| 35 | +# conformes aux noms requis par {cfr}: |
| 36 | +# date, cases, deaths |
| 37 | + |
| 38 | +disease_dat |
| 39 | + |
| 40 | +# Est-ce que les données d'entrées sont conformes au format attendu par {cfr} ? |
| 41 | +# Si oui, utiliser: |
| 42 | +disease_adapted <- disease_dat |
| 43 | +# OU |
| 44 | +# Si no, utiliser `cfr::prepare_data()` pour l'adapter: |
| 45 | +disease_adapted <- disease_incidence %>% |
| 46 | + cfr::prepare_data( |
| 47 | + #<À COMPLETER> |
| 48 | + ) |
| 49 | + |
| 50 | +disease_adapted |
| 51 | + |
| 52 | +# Accéder aux functions de distribution de délais ------------------------------ |
| 53 | +# étape: Accéder aux distributions de probability de délais entre apparition des |
| 54 | +# symptomes et le décès. |
| 55 | + |
| 56 | +# De quel délai a t-on besoin pour ajuster le RL? (dépend de la maladie) |
| 57 | +disease_delay <- epiparameter::#<À COMPLETER> |
| 58 | + |
| 59 | + |
| 60 | +# Estimer le RL naïf and ajusté aux délais ------------------------------------- |
| 61 | +# étape: Estimation du RL naïf and ajusté aux délais. |
| 62 | + |
| 63 | +# Estimer le RL naïf static |
| 64 | +disease_adapted %>% |
| 65 | + cfr::cfr_static() |
| 66 | + |
| 67 | +# Estimer le RL static ajusté au délai |
| 68 | +disease_adapted %>% |
| 69 | + cfr::cfr_static( |
| 70 | + delay_density = #<À COMPLETER> |
| 71 | + ) |
| 72 | + |
| 73 | +# étape: copier et coller les résultats obtenus, et répondre aux questions. |
| 74 | + |
| 75 | +# nolint end |
0 commit comments