tutorials

Simulation-based inference for epidemiological dynamics with pypomp

Aaron A. King, Edward L. Ionides and Kunyang He

Overview

This short course introduces statistical inference techniques and computational methods for dynamic models of epidemiological systems. The course explores deterministic and stochastic formulations of epidemiological dynamics and develop inference methods appropriate for a range of models. Special emphasis will be on exact and approximate likelihood as the key elements in parameter estimation, hypothesis testing, and model selection. Specifically, the course emphasizes sequential Monte Carlo techniques.

This is an adaptation of the previous SBIED short course from R-pomp to pypomp. As such, its objectives are:

  1. To introduce partially observed Markov process (POMP) models as tools for scientific investigation and public health policy.
  2. To give students the ability to formulate POMP models of their own.
  3. To teach efficient approaches for performing scientific inference using POMP models.
  4. To familiarize students with the pypomp package.
  5. To help students familiar with the R-pomp package to transfer that expertise to pypomp.
  6. To give students opportunities to work with inference methods, including methodologies that use automatic differentiation and GPUs which are not accessible via pomp.
  7. To provide documented examples for student re-use.

Lessons

0. Instructions for preparing your laptop for the course exercises. Not yet implemented.

1. Introduction: What is “Simulation-based Inference for Epidemiological Dynamics”? POMPs and pypomp.

2. Simulation of stochastic dynamic models.

3. Likelihood for POMPs: theory and practice.

4. Iterated filtering: theory and practice. Not yet implemented.

5. Case study: measles. Recurrent epidemics, long time series, covariates, extra-demographic stochasticity, interpretation of parameter estimates. Not yet implemented.

6. Case study: polio. Workflow for a real research problem. Not yet implemented.

7. Case study: Ebola. Model diagnostics and forecasting. Not yet implemented.

8. Case study: HIV and fluctuating sexual contact rates. Panel data. Not yet implemented.