What to do
A general note: Some of the tasks below are fairly open ended. You shouldn’t go through the tasks like a check-list. Instead, explore these simulations, play with them, query them. As you do, loop through iterations of formulating your expectations for a specific scenario/task, running the model and observing results, and if results don’t agree with your expectations, figure out why. Essentially, “do science/research”.
Task 1:
- Set the model parameters such that it corresponds to the following setting:
- A population size of 1000, 1 initially infected host, no pathogen in environment, simulation duration approximately 5 years.
- We’ll first look at direct transmission. Set transmission rate of bd = 0.001, environmental transmission to 0.
- Set environmental pathogen shedding and clearance to 0.
- Assume that the duration of the infectious period is about 2 months long.
- Turn off births and deaths for now.
- Run the simulation. You should get the usual direct transmission outbreak and SIR dynamics.
- Record the number and fraction of susceptible/infected/recovered remaining at the end of the outbreak.
Task 2:
- Now try various values for the rate of shedding into the environment and environmental clearance. Leave everything else unchanged.
- As you change those parameters away from 0, what do you expect to see?
- Think about your expectations, run the simulation and compare your expectations with the results.
- Anything surprising happening? Do you understand why you are seeing what you see?
Note: The ‘pathogen in the environment’ units are arbitrary here, and parameters are chosen such that everything fits (mostly) into a single plot. To compare with real data, where environmental pathogen load might be very high, you would need to show separate plots, or plot on a log scale.
Task 3:
- Now set p = 10 and c = 10 . Turn off direct transmission. Run to make sure you don’t get an outbreak.
- Turn on environmental transmission, set it to be = 0.001. What do you expect to see?
- Run simulations and compare with the direct transmission results.
Task 4:
- Now also turn on direct transmission by setting bd = be = 0.001. What do you expect to see?
- Run simulations and see how your expecations hold up.
Task 5:
- Now set initial number of infected to 0 and inital pathogen in the environment to 1. What do you expect to see?
- Run simulations and see how your expecations hold up.
Task 6:
- Try different combinations of direct and environmental transmission only and having only infected hosts or only pathogen in the environment at the start. What do you expect to see?
- Run simulations and see how your expecations hold up.
Task 7:
- Keep exploring by trying different parameters and transmission settings and see how they influence results.
- You can also turn on births/deaths and see what you get.
- As you continue your exploration, think about potential real infectious diseases that might be approximated by either one of the transmission types, and what settings for other model parameters would be for those ID.