A Math Model Is Predicting the Ebola Outbreak
motherboard.vice.com
motherboard.vice.com
The simplest ODE is something like:
du/dt = k(t) * u(t)
Whose solution locally looks like (for slow-varying k(t))
u(t) ∝ e^{ k(t) * t }
One possibility is k(t) > 0 sustained for 10 effective doubling times. One possibility is k(t) < 0. In this epidemic, either is plausible (?): the difference is only a small difference in some infection control parameters. A small uncertainty that translates into a factor of a thousand uncertainty in the outcome, because it gets blown up by a gigantic e^x.
And the estimates of k(t) seem to hover around the critical value of k(t) ~ 0. In the CDC model (their Excel spreadsheet is open source [0], and FYI won't import into LibreOffice), the shape of the epidemic is purely determined by the shape of their k(t) assumptions. Their defaults parameter exhibit, first, a fast growth phase, assuming k(t) > 0; then they assume a slight reduction to k(t) ~ 0, leading to slower growth; then slightly more reduction to k(t) < 0, causing the epidemic to halt. As far as I can tell k(t) is basically speculative, and completely determines the shape of u(t). So really you predict nothing.
A 'predictive model' should say 'if you do X then you will end up with Y' - and the X cannot be adjusting some number. The X has to be stuff like 'building ETU's in West Africa', or 'canceling all flights',... A predictive model should be able to say 'don't bother canceling flights, it's no use - instead do this...'.
Note: I'm going to be teaching a course in Erlang programming next month where the homework assignment is epidemic modeling - 100K's of concurrent processes moving around, getting exposed to each other - that sort of stuff. Rather ambitious, but I feel it's a good use case for Erlang - a pandemic is a 'viral' chat application.
This is just wrong. A predictive model does not necessarily have any "action" input. Example: weather forecast.
How does the model account for isolation capacity constraints?
In other words, hospitals have a finite capacity to deal with these kinds of patients.
How does the model account for infection and mortality within the medical community?
How does the model account for handling and or procedural errors?
How does the model account for the threshold beyond which N percent of medical professionals will refuse to treat patients?
Does that model estimate this threshold? This is pure conjecture on my part, I am guessing that if four or five nurses and doctors fall ill, or worst, die, it could trigger a really difficult scenario to deal with.
Does the model account for some of the tens of thousands of afflicted in West Africa travelling outside of the hot zone while not showing symptoms and then infecting general populations during the initial phase of their sickness?
Does the model assume a tolerance range for various parameters due to mutation or other factors, for example rate of reproduction, ease of communication, etc.?
the second model at least accounts for decay - gee - but is still essentially two terms with three parameters- the result is two competing decay terms.
This is, essentially, a trivial joke. it has no delay parameters for incubation, etc. I've only skimmed the article, but I can't find an actual attribution of these numbers to the model (the model itself is from a journal). It seriously appears like an editor said, "Let's get some scary-mathy-lookin stuff up there on Ebola).
if it bleeds...
That was funny. Yup, math is scary to a lot of people. Sad statement.
It doesn't account for finite medical capacity. It assumes a specified fraction of patients are isolated in hospitals at a given time, a fraction that increases, regardless of the size of the epidemic. I think the burden is on users to check if the results make any sense (it does give you the # of beds in use at a time, so you get that feedback).
The problem is, as the WHO has made very, very clear: "Problems with data gathering in Liberia continue. It should be emphasized that the reported fall in the number of new cases in Liberia over the past three weeks is unlikely to be genuine. Rather, it reflects a deterioration in the ability of overwhelmed responders to record accurate epidemiological data. It is clear from field reports and first responders that EVD cases are being under-reported from several key locations, and laboratory data that have not yet been integrated into official estimates indicate an increase in the number of new cases in Liberia. There is no evidence that the EVD epidemic in West Africa is being brought under control ...
Evidence obtained from responders and laboratory staff in [Liberia] indicates beyond doubt that there is widespread under-reporting of new cases, and that the situation in Liberia, and in Monrovia in particular, continues to deteriorate from week to week. Approximately 200 new probable and suspected cases, but very few confirmed cases, have been reported in the capital Monrovia in each of the past three weeks. A substantial proportion of these suspected cases are most probably genuine cases of EVD, and the reported fall in confirmed cases over the past three weeks reflects delays in matching laboratory results with clinical surveillance data. Efforts continue to urgently address problems with data acquisition in what is an extremely challenging environment, and it is likely that the figures will be revised upwards in due course."
In other words, awesome job building a model that fits garbage data.
http://www.plosone.org/article/fetchObject.action?uri=info%3...
Perhaps someone can help discuss the merits of this particular model? In my humble opinion, I suspect a lack of empirically gathered results and support for it's justifications. On the other hand, what qualifies as a representational event if we're supposed to prevent human death?
It all falls over when The Mule appears! (Sorry if this is a spoiler for you)