Stochastic model of the spread of ebola

A few weeks ago I heard Oskar Hallatschek from the physics department give a talk about the spatial dynamics of evolutionary spread. He was modeling evolutionary drift, but the stochastic models can apply to epidemics too: the analogy is between mutations popping up in nearby populations and disease outbreaks. The talk was very timely considering the ebola crisis. In my mind I keep imagining a haunting series of animated simulations of spread events that he showed (you can see a static picture of the end result here).  I’m curious how the dynamics of ebola relate to these stochastic models.  It’s intuitively obvious that the way outbreaks spread depends on how many people are around, how contagious the disease is, etc, but it’s not so clear how to model and estimate these parameters.

The best study of the dynamics of this ebola outbreak that I’ve seen so far came out of PLOS.  People can speculate about what will happen all day, but I always appreciate when there’s some data backing up speculations.  Plus, their data visualizations were nice enough to be shown as evidence in Congress the other day.  Using standard stochastic models for epidemics, the authors estimated the basic reproductive number in West Africa using data from the beginning of July.  They found the reproductive number is around R0 = 1.8 which means that on average, one person with ebola will generate 1.8 more ebola cases.  In this model if R0 > 1, then the disease will tend to spread and become an epidemic.

Then, after the models for local outbreaks in West Africa were selected, the authors combined the models with airline data to estimate the spread of ebola internationally via simulations.  This is something that models like Hallatschek’s evolutionary spread dynamics and typical epidemic models aren’t able to capture because it incorporates “importation events”, big jumps across distances that wouldn’t be possible but for modern transportation.  The authors considered two cases, if the Nigeria outbreak would be contained during the timeline of their simulations and if it were not (fortunately it has since been contained).  They estimated both the probability of importation events (see the plots below) and the number of cases generated by an importation event in each country.  Judging by their results, it seems like the spread of ebola across the globe is quite likely, but vigilant public health efforts will be able to prevent a large number of cases.

This study came out more than 6 weeks ago so I’d be interested to see updated projections given the new data since then.