valeriupredoi / COVID-19_LINEAR

Simple script to plot daily number of cases
7 stars 0 forks source link

Plotting the doppeling time #2

Closed TakeItAndRun closed 4 years ago

TakeItAndRun commented 4 years ago

Would it be possible to plot the doubling time as a function of time. Using always the last X days for the exponential regression. This would really show the efficiency of containment effects. How to choose X? In this paper:

https://annals.org/aim/fullarticle/2762808/incubation-period-coronavirus-disease-2019-covid-19-from-publicly-reported

they plot the Cumulative distribution function of the COVID-19 incubation:

Incubation period

"The estimated median incubation period of COVID-19 was 5.1 days (CI, 4.5 to 5.8 days). We estimated that fewer than 2.5% of infected persons will display symptoms within 2.2 days (CI, 1.8 to 2.9 days) of exposure, whereas symptom onset will occur within 11.5 days (CI, 8.2 to 15.6 days) for 97.5% of infected persons."

So I would initially choose X=5.

I think this data also means that any change of the containment of the corona virus would only been seen after about the median incubation period later in the data?

valeriupredoi commented 4 years ago

yeah good call! I suspect the containment will show up like you say, on average, after x=incubation time - it's quite visible in China 2019-ncov_lin_02-02-2020 Quarantine started January 23rd and you can see a visible slowdown January 27. I'll have a stab at plugging in the doubling time tomorrow, cheers for the suggestion :beer:

valeriupredoi commented 4 years ago

I have now started the automated recording of the daily mean (over 12 countries) doubling time in this file - I am interested in seeing a decrease over time suggesting the effect of local quarantines (this does not include the UK since there is no actual quarantine instated here just yet sigh). Hopefully we'll see a noticeable change in a couple weeks :+1:

TakeItAndRun commented 4 years ago

Are you starting all your plots on the March 1th? As time goes on your range gets longer and longer. If there are changes in the rise over time, than the quality of your fits should be getting poorer and poorer.


There might be a way to find the best range to do the regression over: Do the regression for the last 2 days. Then do the regression for one day longer each. Pick the regression with the best fit (R closest to one)

valeriupredoi commented 4 years ago

yeah good call! I have now changed to two regressions for select countries that show slowdowns in the past few days; I can always run an MCMC to find the curve parameters but I must do my day job (from home) so not much time to implement it :cry:

TakeItAndRun commented 4 years ago

Nice work. It would be so valuable to get your work out there to make people see that/which measures are effective and necessary.

valeriupredoi commented 4 years ago

cheers! am gonna close this since all is now documented, have a look at https://valeriupredoi.github.io/ and stay safe and have a good weekend! :beer: