MODERATOR NOTE: coronavirus resources & updates: https://wolfr.am/coronavirus
The notebook here is also (without output) placed in GitHub: https://github.com/antononcube/SystemModeling/blob/master/Projects/Coronavirus-propagation-dynamics/WL-notebooks/Basic-experiments-workflow-for-simple-epidemiological-models.nb .
Many thanks for the effort you put into this. I also much enjoy your website and learn a lot from it.
Thank your kind words!
I put a fix into the notebook -- I forgot to run it with the most recent package updates.
Nice live-stream today!
-- you have earned Featured Contributor Badge
Your exceptional post has been selected for our editorial column Staff Picks http://wolfr.am/StaffPicks and Your Profile is now distinguished by a Featured Contributor Badge and is displayed on the Featured Contributor Board. Thank you, keep it coming, and consider contributing to the The Notebook Archive!
One possible issue is from an oversimplification. In actuality the beta parameters are going to be functions of time and space. To some extent quarantining can account for this, at least the temporal aspect, if you can model it by some plausible kinetics.
Another issue is in the data itself. There are believed to be many mild cases. For these, recovery=infected (up to a time lag) and if both go unreported, this will have the effect of making recovery seem slower than it actually is.
That stated, both posts are really nice.
Probably refined models can be made that have more population compartments and have time-dependent rates that reflect the medical observations and enacted policies.
My work is posted now: [Notebook] Epidemiological Models for Influenza and COVID-19
I'll get a a look at your repo later today.
Our models are very similar, right down to some of the assumptions! I see that your fits suffer from the same problems as mine: a long exponential tail on the confirmed cases (aka, infected) and the recovered rises too slowly.
Thank you for your comment, Bob!
The purpose of that notebook is to show the general workflow with a simple model.
There are multiple ways to make the SEI2R model more comprehensive. The directions I have in mind can be found in the repository Coronavirus propagation dynamics.