@@ -8,6 +8,7 @@ room_number <- #<COMPLETE> replace with 1/2/3/4
88
99# Combine interventions --------------------------------------------------
1010
11+ # step: complete the intervention or vaccination arguments
1112simulate_twointerventions <- epidemics :: model_default(
1213 population = population_object ,
1314 transmission_rate = transmission_rate ,
@@ -19,19 +20,28 @@ simulate_twointerventions <- epidemics::model_default(
1920 increment = 1.0
2021)
2122
23+ # step: paste table output in report
2224epidemics :: epidemic_peak(simulate_twointerventions )
2325
2426# Visualize effect --------------------------------------------------------
2527# Plot new infections
2628
29+ # step:
30+ # add intervention name
31+ # if your intervention is vaccination, then
32+ # activate the argument compartments_from_susceptible
33+ # run and paste plot output in report
34+
2735infections_baseline <- epidemics :: new_infections(
2836 data = simulate_baseline ,
29- by_group = FALSE
37+ # compartments_from_susceptible = "vaccinated", # if vaccination
38+ by_group = FALSE # if TRUE, then age-stratified output
3039)
3140
3241infections_twointerventions <- epidemics :: new_infections(
3342 data = simulate_twointerventions ,
34- by_group = FALSE
43+ # compartments_from_susceptible = "vaccinated", # if vaccination
44+ by_group = FALSE # if TRUE, then age-stratified output
3545)
3646
3747# Assign scenario names
0 commit comments