Leverhulme Research Project Grant (RPG-2015-049): 2018-2019.
Paul Blackwell
Dylan Childs (Primary Investigator)
Mark Rees
An important class of physiologically structured population models is the dynamic energy budget (DEB), which outlines the relationship between an individual’s resource acquisition, growth, maintenance and reproduction. By aggregating individuals through their competition for a shared resource, the DEB can be scaled up to a population, creating individual variation in life histories and allowing the resulting population dynamics to be modelled.
DEBs are often evaluated against data in an informal way by selecting parameters from the literature, but here we aim to infer the parameters in a robust way. We use a DEB tailored to the Soay sheep of St Kilda and use recapture data spanning over 30 years to calibrate this model.
This inverse problem in challenging, and we approach it likelihood-free, implementing approximate Bayesian computation (ABC). This problem is high dimensional, not only in parameter-space, but in the DEB output, it being a whole population’s individual life histories. We implement efficient ways to deal with this high dimension, including sequential inference methods to handle high parameter dimension and auxiliary modelling to produce low dimension summary statistics that retain high informative value for the original DEB.
PhD thesis: 2014-2018.
Paul Blackwell (supervisor)
This PhD thesis concerns the statistical modelling of animal movement paths given observed GPS locations. With observations being in discrete time, mechanistic models of movement are often formulated as such. This popularity remains despite an inability to compare analyses through scale invariance and common problems handling irregularly timed observations. A natural solution is to formulate in continuous time, yet uptake of this has been slow, often excused by a difficulty in interpreting the `instantaneous’ parameters associated with a continuous-time model.
The aim here was to bolster usage by developing a continuous-time model with interpretable parameters, similar to those of popular discrete-time models that use turning angles and step lengths to describe the movement process. Movement is defined by a continuous-time, joint bearing and speed process, the parameters of which are dependent on a continuous-time behavioural switching process, thus creating a flexible class of movement models. Further, we allow for the observed locations derived from this process to have unknown error. Markov chain Monte Carlo inference is presented for parameters given irregular, noisy observations. The approach involves augmenting the observed locations with a reconstruction of the underlying continuous-time process.
Example implementations showcasing this method are given featuring simulated and real datasets. Data from elk (Cervus elaphus), which have previously been modelled in discrete time, demonstrate the interpretable nature of the model, finding clear differences in behaviour over time and insights into short-term behaviour that could not have been obtained in discrete time. Observations from reindeer (Rangifer tarandus) reveal the effect observation error has on the identification of large turning angles—a feature often inferred in discrete-time modelling. Scalability to realistically large datasets is shown for lesser black-backed gull (Larus fuscus) data.