ryansmcgee / seirsplus

Models of SEIRS epidemic dynamics with extensions, including network-structured populations, testing, contact tracing, and social distancing.

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

The effects of checkpoint

cesargm2015 opened this issue · comments

commented

When I set a checkpoint between t=15 and t=75, I was expecting to see something like this:

Screen Shot 2020-04-11 at 14 11 53

The blue curve represents the number of infected people if no restrictions are applied. The green curve is during a simulated quarantine supposing a social contact reduction of 50% (that is, beta is cut in half) at time t=15 and goes back to the original value at the end of the checkpoint (I know this is different in the real world, just for the sake of simplicity).

But this is what I'm getting from seirsplus:

Screen Shot 2020-04-11 at 14 12 39

The checkpoint doesn't seem to take any effect. Below is a print of the settings I'm using:

Screen Shot 2020-04-11 at 14 19 44

What am I getting wrong?

The checkpoints are both being applied early relative to when the main wave takes off given your parameters, and what is happening during that phase is difficult to see due to the y axis being dominated by the larger peak. What does it look like if you only plot up to time 80, say? I suspect your reduction of beta is having an effect during that window, but it's just not visible in the plot. Also, your reduced beta corresponds to a R value of 1.92, which is still well above 1, so I wouldn't expect a drop in new cases, just a slowing of the rate of spread. If after looking into these aspects you still suspect something is wrong the checkpoints, please do let me know.