Skip to content

Commit e55b4d9

Browse files
committed
record instructor notes for sections oen and two
1 parent 88177a8 commit e55b4d9

File tree

1 file changed

+95
-75
lines changed

1 file changed

+95
-75
lines changed

instructors/04-practical-tutors.qmd

Lines changed: 95 additions & 75 deletions
Original file line numberDiff line numberDiff line change
@@ -177,76 +177,76 @@ Write your answers to the questions above:
177177

178178
#### Outputs
179179

180-
| all compartments | new infections |
181-
|---|---|
182-
| ![image](https://hackmd.io/_uploads/r1o8OF_Rkx.png) | ![image](https://hackmd.io/_uploads/Syi1tY_R1x.png) |
180+
| country | all compartments | new infections |
181+
| -------- | -------- | -------- |
182+
| Italy | ![image](https://hackmd.io/_uploads/SyVnWdXvle.png) | ![image](https://hackmd.io/_uploads/Syua-OQPeg.png) |
183+
184+
```
185+
epidemics::epidemic_peak(data = simulate_baseline)
186+
#> demography_group compartment time value
187+
#> <char> <char> <num> <num>
188+
#> 1: [0,20) infectious 320 513985.5
189+
#> 2: [20,40) infectious 328 560947.3
190+
#> 3: 40+ infectious 329 932989.8
191+
```
192+
193+
| country | all compartments | new infections |
194+
| -------- | -------- | -------- |
195+
| Vietnam | ![image](https://hackmd.io/_uploads/BJFgQ_mvel.png) | ![image](https://hackmd.io/_uploads/HJ8-Qdmvlg.png) |
183196

184197
```
185-
# Get epidemic_peak
186198
epidemics::epidemic_peak(data = simulate_baseline)
187-
# Output:
188-
# demography_group compartment time value
189-
# <char> <char> <num> <num>
190-
# 1: [0,20) infectious 320 513985.6
191-
# 2: [20,40) infectious 328 560947.2
192-
# 3: 40+ infectious 329 932989.6
199+
#> demography_group compartment time value
200+
#> <char> <char> <num> <num>
201+
#> 1: [0,20) infectious 325 929142.4
202+
#> 2: [20,40) infectious 322 975411.0
203+
#> 3: 40+ infectious 317 1053391.8
204+
```
205+
206+
| country | all compartments | new infections |
207+
| -------- | -------- | -------- |
208+
| Zimbabwe | ![image](https://hackmd.io/_uploads/HJ_t7dXwgg.png) | ![image](https://hackmd.io/_uploads/Sk_9mdXDgx.png) |
209+
210+
```
211+
epidemics::epidemic_peak(data = simulate_baseline)
212+
#> demography_group compartment time value
213+
#> <char> <char> <num> <num>
214+
#> 1: [0,20) infectious 322 277709.4
215+
#> 2: [20,40) infectious 321 188967.9
216+
#> 3: 40+ infectious 317 111165.7
193217
```
194218

195219
#### Interpretation
196220

197221
Interpretation template:
198222

199-
+ In the population, the demographic group of `age from [0,20]` has a peak of infectious individuals at day `320` with a size of `513,986` individuals at that time.
223+
+ In an age-structured SEIR epidemic model for influenza transmission in Zimbabwe, the demographic group aged 0 to 20 years (`[0,20]`) reaches its peak number of *infectious individuals* on day `322`, with a peak size of `277,709` individuals.
200224

201225
Compare output types:
202226

203227
+ `epidemics::epidemic_peak(data = simulate_baseline)`
204-
+ The table output gives exact values for time and size of peak for *infectious individuals*.
228+
+ The table output gives exact values for time and size of peak for *infectious individuals* across age groups.
205229
+ `epidemics::new_infections(data = simulate_baseline)`
206-
+ We can not get exact value for time and size of peak of *new infections*.
207-
+ We can make qualitative comparisons.
230+
+ We can plot the trajectories of *new infections* across age groups, but not get exact value for time and size of peak directly.
231+
+ We can make qualitative comparisons between countries or scenarios.
208232
+ Comparing plots:
209-
+ The peak size of *new infections* is lower than the peak size of *infectious individuals*. New infections are defined as the daily outflow of individuals from the susceptible to the exposed compartment.
233+
+ The peak size of *new infections* is lower than the peak size of *infectious individuals*.
234+
+ New infections are defined as the daily outflow of individuals from the susceptible to the exposed compartment.
235+
+ Infectious are defined as the total number of cumulative amount of individuals in the infectious compartment at each time.
210236
+ The peak time may be similar in both outputs.
211-
+ Other packages that can estimate trend of new infections are `{EpiNow2}` and `{epichains}`.
237+
+ Other packages that can estimate the trend of new infections are `{EpiNow2}` and `{epichains}`.
212238

213239
Comparison between rooms:
214240

215-
| Vietnam | Zimbabwe |
216-
|---|---|
217-
| ![image](https://hackmd.io/_uploads/BkZGRY_Rkx.png) | ![image](https://hackmd.io/_uploads/SJVlkcORkl.png) |
218-
219-
- Population structure from Italy, Vietnam, and Zimbabwe influences the progression of the transmission in each population.
220-
221-
```
222-
# Italy
223-
contact_data$demography
224-
#> age.group population proportion year
225-
#> <char> <num> <num> <int>
226-
#> 1: [0,20) 11204261 0.1905212 2005
227-
#> 2: [20,40) 16305622 0.2772665 2005
228-
#> 3: 40+ 31298598 0.5322123 2005
229-
```
241+
- Population structure and age-specific social contact patterns in each country influence the progression of disease transmission.
242+
- Using `{socialmixr}`, the symmetric contact matrix contains the *mean number of contacts* that an individual in each age group (row) reports having with individuals of the same or another age group (column).
243+
+ Note: The contact matrix may look asymmetric, but it is *symmetric in total contacts*. That is, the total number of contacts from one group to another is the same in both directions — check this by multiplying the mean contacts by the population size for each group.
230244

231-
```
232-
# Vietnam
233-
contact_data$demography
234-
#> age.group population proportion year
235-
#> <char> <num> <num> <int>
236-
#> 1: [0,20) 31847968 0.3777536 2005
237-
#> 2: [20,40) 28759380 0.3411194 2005
238-
#> 3: 40+ 23701489 0.2811270 2005
239-
```
245+
| Italy | Vietnam | Zimbabwe |
246+
| --- | --- | --- |
247+
| ![image](https://hackmd.io/_uploads/HyYOZ57Dll.png) | ![image](https://hackmd.io/_uploads/HkTjb9QDxl.png) | ![image](https://hackmd.io/_uploads/Bkrhx5Xvlx.png) |
240248

241-
```
242-
# Zimbabwe
243-
contact_data$demography
244-
#> age.group population proportion year
245-
#> <char> <num> <num> <int>
246-
#> 1: [0,20) 8235388 0.5219721 2015
247-
#> 2: [20,40) 5179150 0.3282628 2015
248-
#> 3: 40+ 2362911 0.1497651 2015
249-
```
249+
Figures using `socialmixr::matrix_plot(contact_data$matrix * contact_data$demography$proportion)`
250250

251251
:::
252252

@@ -377,20 +377,59 @@ Write your answers to the questions above:
377377

378378
#### Outputs
379379

380-
| Intervention | Early start (day 100) | Late start (day 200) |
380+
| intervention | all compartments | new infections |
381381
|---|---|---|
382-
| **School Closure** | ![image](https://hackmd.io/_uploads/HyxrKNWDee.png) | ![image](https://hackmd.io/_uploads/B18MFNWPll.png) |
383-
| **Mask Mandate** | ![image](https://hackmd.io/_uploads/Hym4FEbDxl.png) | ![image](https://hackmd.io/_uploads/HknbYEbPlg.png) |
384-
| **Vaccination** | ![image](https://hackmd.io/_uploads/r12Xt4-Dll.png) | ![image](https://hackmd.io/_uploads/Hy6xKN-vex.png) |
382+
| School closure | ![image](https://hackmd.io/_uploads/SJcXBqmwxg.png) | ![image](https://hackmd.io/_uploads/B1G2HcXwlx.png) |
383+
384+
```
385+
epidemics::epidemic_peak(simulate_intervention)
386+
#> demography_group compartment time value
387+
#> 1: [0,20) infectious 568 190371.93
388+
#> 2: [20,40) infectious 566 121626.39
389+
#> 3: 40+ infectious 562 70255.15
390+
```
391+
392+
| intervention | all compartments | new infections |
393+
|---|---|---|
394+
| Mask mandate | ![image](https://hackmd.io/_uploads/S12cIcXPex.png) | ![image](https://hackmd.io/_uploads/BytgP5QPll.png) |
395+
396+
```
397+
epidemics::epidemic_peak(simulate_intervention)
398+
399+
#> demography_group compartment time value
400+
#> 1: [0,20) infectious 380 88058.78
401+
#> 2: [20,40) infectious 378 61155.76
402+
#> 3: 40+ infectious 374 38143.96
403+
404+
```
405+
406+
| intervention | all compartments | new infections |
407+
|---|---|---|
408+
| Vaccination | ![image](https://hackmd.io/_uploads/rJMVOcXPeg.png) | ![image](https://hackmd.io/_uploads/r1p_O9Qwgx.png) |
409+
410+
```
411+
epidemics::epidemic_peak(data = simulate_intervention)
412+
#> demography_group compartment time value
413+
#> 1: [0,20) infectious 318 155276.44
414+
#> 2: [20,40) infectious 317 107294.83
415+
#> 3: 40+ infectious 314 65867.72
416+
```
385417

386418
#### Interpretation
387419

388420
Interpretation Helpers:
389421

390-
+ School closure with short duration can delay the peak of new infections, but this will keep the same size.
391-
+ Mask mandate of 200 days during the time of the epidemic peak can delay and reduce the size of new infections.
392-
+ Vaccinations earlier in time will have a higher impact in reducing the size of the epidemic peak and extending its delay. Note that the effectiveness of vaccination can depend on various factors, including vaccine efficacy and timing relative to the outbreak.
422+
+ School closure starting on day 200 for a duration of 250 days can delay the peak of infectious across age groups by 240 days aprox., and reduce the total number of new infections in the whole population by 20,000 aprox.
423+
+ Mask mandate starting on day 200 for a duration of 250 days can delay the peak of infectious across age groups by 40 days aprox., and reduce the total number of new infections in the whole population by 50,000 aprox.
424+
+ Vaccinations starting on day 200 for a duration of 250 days can delay the peak of infectious across age groups by 4 days aprox., and reduce the total number of new infections in the whole population by 40,000 aprox.
425+
+ Note that the effectiveness of vaccination can depend on various factors, including **vaccine efficacy** and **timing relative to the outbreak**.
426+
427+
### Additional challenges
393428

429+
1. **How do the start time and duration of interventions influence the timing and size of the peak in new infections?**
430+
Try modifying the intervention start time from day 200 to day 100, or changing the duration from 250 days to 100 days, and observe the impact on the epidemic dynamics in Zimbabwe.
431+
432+
2. **How can interventions affect the timing and size of the peak in new infections across different countries?** Try changing the population from Zimbabwe to Vietnam or Italy to observe how country-specific factors like population structure and social contacts influence the epidemic curve.
394433

395434
:::
396435

@@ -539,25 +578,6 @@ Interpretation Helpers:
539578
540579
```
541580

542-
### Additional challenge
543-
544-
<!-- In Activity 1
545-
546-
- Modify the basic reproduction number (R₀) from 1.46 to 1.1 or 3. What changes do you observe in the outputs? (increase the total simulation time if needed) -->
547-
548-
<!-- Modify the basic reproduction number (R₀):
549-
550-
| R = 1.1 | R = 3 |
551-
|---|---|
552-
| ![image](https://hackmd.io/_uploads/H1UupFOAyl.png) | ![image](https://hackmd.io/_uploads/ryVoat_R1l.png) |
553-
554-
555-
- An epidemic with `R₀ = 1.1` has a days delayed and smaller outbreak based on number of infections (day 1200, 9000 new infections), compared with `R₀ = 3` with an earlier and higher peak than `R₀ = 1.46` (day 100, 1,000,000 new infections). -->
556-
557-
In Activity 2
558-
559-
- How can the start time or duration of interventions affect the timing and size of the peak in *new infections*? Try changing the start time from day 200 to day 100, or the duration from 250 days to 100 days.
560-
561581
# Continue your learning path
562582

563583
<!-- Suggest learners to Epiverse-TRACE documentation or external resources --->

0 commit comments

Comments
 (0)