Slides and code from talk on 11 August 2014