Message Boards Message Boards

Coronavirus logistic growth model: China

Posted 5 years ago

MODERATOR NOTE: coronavirus resources & updates: https://wolfr.am/coronavirus


POSTED BY: Robert Rimmer
3 Replies

enter image description here -- you have earned Featured Contributor Badge enter image description here

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!

POSTED BY: EDITORIAL BOARD

Thanks Robert! I really enjoyed reading your post. I have a question and I appreciate if you let me know your thoughts. If I use only some of the data of, e.g., Hubei China, the logistics model seems to give very different results. For example, let's say that we are on the 20th day of outbreak. If we use the available data of the first 20 days, the model fails to reproduce the data of the following 16 days.

POSTED BY: Mohammad Bahrami
Posted 5 years ago

The data at the beginning of the epidemic counts are not likely to be very accurate. See some of the discussion on later posts. There are some basic numerical reasons also for this: there are only 2 or 3 significant digits of information in the early numbers, where as the later numbers have higher precision. But the main reason is likely that the early data is also probably not accurate, especially in Hubei where the laboratories probably were quickly overwhelmed. They stopped using laboratory to count cases somewhere in the middle, and the information at WHO suggests that the death counts were only for lab confirmed cases. Then later they added the other death counts in a lump, causing the jump in the accumulated data. What you should try is using the recent data which contains more information and is likely to be more accurate. Also try looking at sequential fits as data accumulates.

nlmd = NonlinearModelFit[#, 
     logistic, {{k, 0.3}, {L, 3000}, {x0, 20}}, x] & /@ 
   Table[Take[deathCases, i], {i, 5, Length@deathCases}];

Then look at the sequential parameters with ListPlot to see how the parameters stabilize with time. It is a good idea to give NonlinearModelFit some reasonable starting parameters.

POSTED BY: Robert Rimmer
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract