44# Activity 3
55
66# step: fill in your room number
7- room_number <- 1
7+ room_number <- 1 # valid for all, account by specific changes
88
99# Combine interventions --------------------------------------------------
1010
@@ -16,10 +16,10 @@ simulate_twointerventions <- epidemics::model_default(
1616 recovery_rate = recovery_rate ,
1717 # Intervention
1818 intervention = list (
19- # transmission_rate = intervention_mask_mandate#,
20- contacts = intervention_schoolclosure
19+ # transmission_rate = intervention_mask_mandate#, #<CHANGE-BY-ROOM>
20+ contacts = intervention_schoolclosure # <CHANGE-BY-ROOM>
2121 ),
22- vaccination = intervention_vaccinate ,
22+ vaccination = intervention_vaccinate , # <CHANGE-BY-ROOM>
2323 time_end = 1000 ,
2424 increment = 1.0
2525)
@@ -38,30 +38,30 @@ epidemics::epidemic_peak(simulate_twointerventions)
3838
3939infections_baseline <- epidemics :: new_infections(
4040 data = simulate_baseline ,
41- compartments_from_susceptible = " vaccinated" , # if vaccination
41+ compartments_from_susceptible = " vaccinated" , # if vaccination #<CHANGE-BY-ROOM>
4242 by_group = FALSE # if TRUE, then age-stratified output
4343)
4444
4545infections_twointerventions <- epidemics :: new_infections(
4646 data = simulate_twointerventions ,
47- compartments_from_susceptible = " vaccinated" , # if vaccination
47+ compartments_from_susceptible = " vaccinated" , # if vaccination #<CHANGE-BY-ROOM>
4848 by_group = FALSE # if TRUE, then age-stratified output
4949)
5050
5151# Assign scenario names
5252infections_baseline $ scenario <- " Baseline"
53- infections_twointerventions $ scenario <- " School closure + Vaccination"
54- # infections_twointerventions$scenario <- "Mask mandate + School closure"
55- # infections_twointerventions$scenario <- "Mask mandate + Vaccination"
53+ infections_twointerventions $ scenario <- " School closure + Vaccination" # <CHANGE-BY-ROOM>
54+ # infections_twointerventions$scenario <- "Mask mandate + School closure" #<CHANGE-BY-ROOM>
55+ # infections_twointerventions$scenario <- "Mask mandate + Vaccination" #<CHANGE-BY-ROOM>
5656
5757# Compare interventions --------------------------------------------------
5858
5959# Combine the data from all scenarios
6060compare_interventions <- dplyr :: bind_rows(
6161 infections_baseline ,
62- infections_schoolclosure ,
63- # infections_mask_mandate,
64- # infections_vaccinate,
62+ infections_schoolclosure , # <CHANGE-BY-ROOM>
63+ # infections_mask_mandate, #<CHANGE-BY-ROOM>
64+ # infections_vaccinate, #<CHANGE-BY-ROOM>
6565 infections_twointerventions
6666)
6767
0 commit comments