$DiscretePOMP.jl$ is a package for:
- Bayesian parameter inference, and
- Simulation of,
- Discrete-state-space Partially Observed Markov Processes, in Julia.
- It also includes automated tools for things like convergence diagnosis, model assessment and analysis.
Discrete-state-space (DSS) models are used throughout ecology and other scientific domains to represent systems comprised of interacting components (e.g. people or molecules.)
A well-known example is the Kermack-McKendrick susceptible-infectious-susceptible (SIR) model:
See the Simple example for a brief primer on DSS.
In applied situations (e.g. like a scientific study) such systems are often difficult to directly observe, and so they are referred to in context as being Partially Observed.
The dynamics (how the system changes over time) of the SIR, and other DSS models, can be represented in continuous time by [a set of coupled] Markov Processes. More specifically, we can define a probability density (a 'likelihood function' in Bayesian parlance) that governs the time-evolution of the system under study.
Combining these concepts, we have a general class of statistical model: Discrete-state-space Partially Observed Markov Processes, or Discrete POMP.
Furthermore, given some applicable [partially complete] scientific data, they yield a paradigm for (in this case, Bayesian) statistical inference based on that model class. That is, we can infer [the likely value of] unknown quantities, such as the unknown time of a known event (like the introduction of a pathogen,) or a model parameter that characterises the infectiousness of that pathogen.
To summarise, DPOMP models and associated methods allow us to learn about a given system of interest (e.g. an ecosystem, pandemic, chemical reaction, and so on,) even in when the available data is limited ['partial'].
Discrete POMP, and discrete-state-space models in general, have a wide range of applications including:
- Epidemiological modelling (e.g. SEIR models)
- Ecology (e.g. predator-prey dynamics)
- Many other potential use cases, e.g. physics; chemical reactions; social media.
The algorithms implemented by the package for simulation and inference include:
- The Gillespie direct method algorithm
- Data-augmented Markov chain Monte Carlo (MCMC)
- The model-based-proposal (MBP) algorithm
- Particle filters (i.e. Sequential Monte Carlo)
- SMC^2, or iterative-batch-importance sampling (IBIS)
A number of well-known models are provided as predefined examples:
- SIR, SEIR, and other epidemiological model
- The Lotka-Voltera predator-prey model
- Ross-MacDonald two-species malaria model
The package code was initially developed during the course of a postgraduate research project in infectious disease modelling at Biostatistics Scotland, and there is a heavy emphasis on epidemiology and epidemiological modelling throughout.
In practice though, this affects only the applied examples and naming conventions of the predefined models available with the package. Otherwise, the models and methods are applicable to many situations entirely outwith the field of ecology (such as chemical reactions.)
As a prerequisite, the package naturally requires a working installation of the Julia programming language. The package is not yet registered but can nonetheless must be added via the package manager Pkg in the usual way.
From the $Julia$ REPL, run:
using Pkg Pkg.add(url="https://github.com/ScottishCovidResponse/DataRegistryUtils.jl")
See the package code repository to inspect the source code.
- Simple example
- Package manual
- 1C. M. Pooley, S. C. Bishop, and G. Marion. Using model-based proposals for fast parameter inference on discrete state space, continuous-time Markov processes. Journal of The Royal Society Interface, 12(107):20150225–20150225, May 2015.
- 2N. Chopin, P. E. Jacob, and O. Papaspiliopoulos. SMC^2 : an efficient algorithm for sequential analysis of state space models: Sequential Analysis of State Space Models. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 75(3):397–426, June 2013.