Are you calculating infections incorrectly?

[ErwinM]

Hi,

Your calculation for infections:
dI = a*E - gamma*I # where gamma = 1 / D_infectious
However, D_infectious is defined as the days someone is contagious. This seems unrelated to the number of people that recover from infection.

According to the wiki page for SEIR, gamma should represent the recovery rate:
gamma = 1 / days_to_recover

Am I understanding the code/wiki wrong?

Erwin

[ErwinM]

I misread:

You calculate R: REMOVED (not recovered), meaning no longer contributing to spread. And calculate recoveries separately. Gotcha.