@@ -11,51 +11,55 @@ library(tidyverse)
1111
1212# Read reported cases -----------------------------------------------------
1313# for covid
14- dat <- read_rds(" paste/link/url/here/covid" ) %> %
14+ dat <- read_rds(
15+ " paste/link/url/here/covid" # <COMPLETE>
16+ ) %> %
1517 dplyr :: select(date , confirm )
1618# or
1719# for ebola
18- dat <- read_rds(" paste/link/url/here/ebola" ) %> %
20+ dat <- read_rds(
21+ " paste/link/url/here/ebola" # <COMPLETE>
22+ ) %> %
1923 dplyr :: select(date , confirm = cases )
2024
2125
2226# Define a generation time from {epiparameter} to {EpiNow2} ---------------
2327
2428# access a serial interval
25- ebola_serialint <- epiparameter :: epiparameter_db(
29+ dat_serialint <- epiparameter :: epiparameter_db(
2630 # <COMPLETE>
2731)
2832
2933# extract parameters from {epiparameter} object
30- ebola_serialint_params <- epiparameter :: # <COMPLETE>
34+ dat_serialint_params <- epiparameter :: # <COMPLETE>
3135
3236# adapt {epiparameter} to {EpiNow2} distribution inferfase
33- ebola_generationtime <- EpiNow2 :: # <COMPLETE>
37+ dat_generationtime <- EpiNow2 :: # <COMPLETE>
3438
3539
3640# Define the delays from infection to case report for {EpiNow2} -----------
3741
3842# define delay from symptom onset to case report
39- ebola_reportdelay <- EpiNow2 :: # <COMPLETE>
43+ dat_reportdelay <- EpiNow2 :: # <COMPLETE>
4044
4145# define a delay from infection to symptom onset
42- ebola_incubationtime <- epiparameter :: epiparameter_db(
46+ dat_incubationtime <- epiparameter :: epiparameter_db(
4347 # <COMPLETE>
4448)
4549
4650# incubation period: extract distribution parameters
47- ebola_incubationtime_params <- epiparameter :: # <COMPLETE>
51+ dat_incubationtime_params <- epiparameter :: # <COMPLETE>
4852
4953# incubation period: discretize and extract maximum value (p = 99%)
50- ebola_incubationtime_max <- ebola_incubationtime %> % # <COMPLETE>
54+ dat_incubationtime_max <- dat_incubationtime %> % # <COMPLETE>
5155
5256# incubation period: adapt to {EpiNow2} distribution interface
53- ebola_incubationtime_epinow <- EpiNow2 :: # <COMPLETE>
57+ dat_incubationtime_epinow <- EpiNow2 :: # <COMPLETE>
5458
5559# print required input
56- ebola_generationtime
57- ebola_reportdelay
58- ebola_incubationtime_epinow
60+ dat_generationtime
61+ dat_reportdelay
62+ dat_incubationtime_epinow
5963
6064
6165# Set the number of parallel cores for {EpiNow2} --------------------------
0 commit comments