File tree Expand file tree Collapse file tree 1 file changed +3
-10
lines changed Expand file tree Collapse file tree 1 file changed +3
-10
lines changed Original file line number Diff line number Diff line change @@ -66,8 +66,7 @@ serial_interval <- epiparameter_db(
66
66
disease = "mers",
67
67
epi_name = "serial",
68
68
single_epiparameter = TRUE
69
- ) #%>%
70
- #epiparameter::discretise()
69
+ )
71
70
72
71
# Set seed for random number generator
73
72
set.seed(33)
@@ -711,12 +710,10 @@ We can count how many chains reached the 100-case threshold using `{dplyr}` func
711
710
summary_chains %>%
712
711
dplyr::arrange(desc(cases_total)) %>%
713
712
dplyr::filter(cases_total > 100)
714
-
715
- # similar to printing the summary of <epichains>
716
- # the printed vector have the size of each chain
717
- #summary(multiple_epichains)
718
713
```
719
714
715
+ This output should give us equivalent results to ` summary(multiple_epichains) ` .
716
+
720
717
::::::::::::::::::::::::::::::::::: spoiler
721
718
722
719
### Observed cases vs simulated chains
@@ -872,10 +869,6 @@ simulated_chains_mpox <- epichains::simulate_chains(
872
869
generation_time = function(x) generate(x = mpox_serialint, times = x)
873
870
)
874
871
875
- # Print the summary of <epichains>
876
- # the printed vector have the size of each chain
877
- #summary(simulated_chains_mpox)
878
-
879
872
# visualize ---------------------------------------------------------------
880
873
881
874
# daily aggregate of cases
You can’t perform that action at this time.
0 commit comments