Why do we seek to quantify and understand disease dynamics? • Prevention and control of emerging infectious diseases (SARS, HIV/AIDS, H5N1 influenza “bird flu”) • Understanding the development and spread of drug resistant strains (malaria, tuberculosis, MRSA “the hospital super-bug”) Ed Ionides Infectious disease dynamics: a statistical perspective 4 Disease dynamics: epidemiology or ecology, or both? • Environmental host/pathogen dynamics are close to predator/prey relationships which are a central topic of ecology. • Analysis of diseases as ecosystems complements more traditional epidemiology (risk factors etc). • Ecologists typically seek to avoid extinctions, whereas epidemiologists typically seek the reverse. Things are not always this simple... – Helicobacter pylori bacteria used to live in the stomach of most humans. Some strains cause stomach ulcers and cancer. It is almost extinct in the developed world due to widespread use of

Inference for static parameters in state space models (i.e., unknown model parameters that do not vary in time) • Numerical issues have led to a considerable literature on this topic. • Iterated filtering maximizes the likelihood function via taking an average of filtered “local ” parameter estimates obtained by adding noise to the static parameters (which regularizes the numerical issues). This process is recursively repeated while reducing the added noise. • Iterated filtering, implemented by basic SMC, has a plug-and-play property: it requires simulation from the state process but not transition densities. Ed Ionides The theory and practice of iterated filtering 3 Plug-and-play methods for state space models • Statistical methods are plug-and-play if they require simulation from the dynamic model but not explicit likelihood ratios. • Bayesian plug-and-play: 1. Artificial parameter evolution (Liu and West, 2001)

Ed Ionides Time series analysis of infectious disease dynamics 1 Time series analysis of infectious disease dynamics: State of the art and future challenges Epidemic model hierarchies and model validation workshop

What is a “mechanistic ” approach to time series analysis? • Write down equations, based on scientific understanding of a dynamic system, which describe how it evolves with time. • Further equations describe the relationship of the state of the system to available observations on it. • Mechanistic time series analysis concerns drawing inferences from the available data about the hypothesized equations. • Questions of general interest: Are the data consistent with a particular model? If so, for what range of values of model parameters? Does one mechanistic model describe the data better than another? • A defining principle: the model structure should be chosen based on scientific considerations, rather than statistical convenience. Ed Ionides Time series analysis via mechanistic models 3 Example: why quantify biological population dynamics? • Conservation. Mankind is increasingly responsible for managing ecosystems. This requires a quantitative understanding of population behavior. • Public health. Pathogens are also biological populations. Despite successes of vaccination and medical treatment, new diseases are emerging (SARS, HIV/AIDS) and old ones re-emerging due to drug resistant strains (malaria, tuberculosis). Treating the pathogen as part of an ecosystem is one approach to understanding and controlling emergent and re-emergent diseases. • Basic scientific interest. Ed Ionides Time series analysis via mechanistic models 4 Time series data of sufficient quantity and quality to justify mechanistic modeling are increasingly available: Two recent examples • King, Ionides, Pascual and Bouma. Inapparent infections and cholera dynamics. To appear in Nature.

Abstract not found

Inference for partially observed Markov process models has been a longstanding methodological challenge with many scientific and engineering applications. Iterated filtering algorithms maximize the likelihood function for partially observed Markov process models by solving a recursive sequence of filtering problems. We present new theoretical results pertaining to the convergence of iterated filtering algorithms implemented via sequential Monte Carlo filters. This theory complements the growing body of empirical evidence that iterated filtering algorithms provide an effective inference strategy for scientific models of nonlinear dynamic systems. The first step in our theory involves studying a new recursive approach for maximizing the likelihood function of a latent variable model, when this likelihood is evaluated via importance sampling. This leads to the consideration of an iterated importance sampling algorithm which serves as a simple special case of iterated filtering, and may have applicability in its own right. 1. Introduction. Partially observed Markov process (POMP) models are

Inference for partially observed Markov process models has been a longstanding methodological challenge with many scientific and engineering applications. Iterated filtering algorithms maximize the likelihood function for partially observed Markov process models by solving a recursive sequence of filtering problems. We present new theoretical results pertaining to the convergence of iterated filtering algorithms implemented via sequential Monte Carlo filters. This theory complements the growing body of empirical evidence that iterated filtering algorithms provide an effective inference strategy for scientific models of nonlinear dynamic systems. The first step in our theory involves studying a new recursive approach for maximizing the likelihood function of a latent variable model, when this likelihood is evaluated via importance sampling. This leads to the consideration of an iterated importance sampling algorithm which serves as a simple special case of iterated filtering, and may have applicability in its own right. 1

modeling and inference for ecological systems

Time series analysis for nonlinear dynamical systems with applications to modeling of infectious diseases by

blowflies as a case study