(03-31-2020, 09:59 AM)burger Wrote: (03-31-2020, 09:49 AM)BostonCard Wrote: (03-31-2020, 07:42 AM)burger Wrote: But there is little to no support for the bell-shaped curves they're using.
Again, Gaussian error function ≠ Gaussian function. It is not a bell-shaped curve; it's a sinusoidal curve. And the SEIR model would predict that a cumulative count of deaths would be a sinusoidal curve (though not, as far as I can tell, a Gaussian error function).
BC
The model uses the Gaussian error function for the cumulative number of deaths. The curve for the number of deaths by day is just the derivative of the error function, which is a gaussian curve. That's what they're showing in the results here http://covid19.healthdata.org/projections
I'm not sure why they're modeling the cumulative deaths rather than the individual counts by day, as that is surely a far easier equation to model with no integrals involved.
As I read it they used total deaths because that was a "hard" statistic, while the individual day counts are subject to testing limitations etc. Others on this thread have pointed out that deaths may be a better true indicator of cases. That would certainly be true if the cause of death was always properly assessed. In any case, that is why it was used (I think).
Integrals can actually be your friend. The integration smooths out the "noise" induced by non-disease related causes. Reports get delayed, like it did in Santa Clara County yesterday. We went low, low, high for reasons that had nothing to do with the disease, but rather reporting. That is noise. Integrating reduces noise, which can (usually does) make curve fitting easier, not harder.
They were just looking at "sigmoidal" (not sinusoidal, auto-correct in action?) functions. The rationale for this is something like the following. Deaths start out at zero and slowly slowly increase. Then, the number starts to increase faster, then still faster. Finally, the slope of total deaths starts to trend down, trends down more, and then there are no additional deaths (flat). That is an "S" shaped curve. All probability mass functions start at 0 and end at 1. The also have the property that they never decrease. So do other functions that are not probability mass functions. This curve fit has no relation to probabilities.
There are plenty of other ways to produce an estimator of the total number of deaths on day "n" out to infinity given the data from days "n-1" downto day 0. It is probable that an ARMA estimator will do so quite well for a small number of days, but its performance will go to hell as it's prediction goes further into the future. These guys wanted a closed form function that had the properties listed above. There is no way to know if it is predictive, and they didn't quantify how good the curve fit was on the existing data, as far as I could see. So, IMHO, it isn't worth much.
People hate uncertainty. The current trend is to put together some kind of mathematical model that will extrapolate into the future. The qualification of the model in terms of goodness of prediction is usually omitted because it wasn't done. The definition of a "good fit" is in the eyes of the beholder. Are you looking at RMS error, maximum error, average error, or some other norm? Oftentimes these issues are not critically examined, because we want answers, and a nonsense model is an answer. It is assumed to be "better than nothing". It may not be. It may be worse than nothing, because we think we know something that we don't.