Bayesian experimental design for control and surveillance in epidemiology

Effective public health interventions must balance an array of interconnected challenges, and decisions must be made based on scientific evidence from existing information. Building evidence requires extrapolating from limited data using models. But when data are insufficient, it is important to recognize the limitations of model predictions and diagnose how they can be improved. This dissertation shows how principles from Bayesian experimental design can be applied to surveillance and control efforts to allow researchers to get more out of their data and direct limited resources to best effect. We argue a Bayesian perspective on data gathering, where design decisions are made to maximize utility on average over a joint distribution of beliefs and outcomes, is better suited to the epidemiological setting where observational studies are the norm. We illustrate these ideas using a range of models and topics across epidemiology. We focus first on Chagas disease, where in Guatemala an endemic vector continues to cause a high rate of domiciliary infestation in rural communities, and shortages of insecticides and resources for critical house improvements hamper control efforts. Using an adaptive sampling and geospatial modeling framework, we show that interpolating from a traditional design goal of minimizing prediction uncertainty to targeting houses of high risk can satisfy competing objectives, namely, to efficiently identify houses in need of treatment while mitigating sampling bias. We next focus on tick surveillance in the southeastern United States. By framing tick collection surveys as a design problem over time and space, we show optimal survey design can yield greater information compared to random or convenience sampling. Finally, we shift attention from experimental design to the closely related concept of practical identifiability. We propose a novel method to quantify practical identifiability which reflects the average amount of posterior shrinkage that would occur in a Bayesian analysis, without requiring computationally expensive techniques like Markov Chain Monte Carlo. With this method, we demonstrate the limits of using epidemiological models to derive standard statistics such as the basic reproductive number early in an outbreak

Global parameter identification of stochastic reaction networks from single trajectories

We consider the problem of inferring the unknown parameters of a stochastic biochemical network model from a single measured time-course of the concentration of some of the involved species. Such measurements are available, e.g., from live-cell fluorescence microscopy in image-based systems biology. In addition, fluctuation time-courses from, e.g., fluorescence correlation spectroscopy provide additional information about the system dynamics that can be used to more robustly infer parameters than when considering only mean concentrations. Estimating model parameters from a single experimental trajectory enables single-cell measurements and quantification of cell--cell variability. We propose a novel combination of an adaptive Monte Carlo sampler, called Gaussian Adaptation, and efficient exact stochastic simulation algorithms that allows parameter identification from single stochastic trajectories. We benchmark the proposed method on a linear and a non-linear reaction network at steady state and during transient phases. In addition, we demonstrate that the present method also provides an ellipsoidal volume estimate of the viable part of parameter space and is able to estimate the physical volume of the compartment in which the observed reactions take place.Comment: Article in print as a book chapter in Springer's "Advances in Systems Biology

Bayesian analysis of community dynamics

Elucidating the mechanisms responsible for the patterns of species abundance, diversity, and distribution within and across ecological systems is a fundamental research focus in ecology. Species abundance patterns are shaped in a convoluted way by interplays between inter-/intra-specific interactions, environmental forcing, demographic stochasticity, and dispersal. Comprehensive models and suitable inferential and computational tools for teasing out these different factors are quite limited, even though such tools are critically needed to guide the implementation of management and conservation strategies, the efficacy of which rests on a realistic evaluation of the underlying mechanisms. This is even more so in the prevailing context of concerns over climate change progress and its potential impacts on ecosystems. This thesis utilized the flexible hierarchical Bayesian modelling framework in combination with the computer intensive methods known as Markov chain Monte Carlo, to develop methodologies for identifying and evaluating the factors that control the structure and dynamics of ecological communities. These methodologies were used to analyze data from a range of taxa: macro-moths (Lepidoptera), fish, crustaceans, birds, and rodents. Environmental stochasticity emerged as the most important driver of community dynamics, followed by density dependent regulation; the influence of inter-specific interactions on community-level variances was broadly minor. This thesis contributes to the understanding of the mechanisms underlying the structure and dynamics of ecological communities, by showing directly that environmental fluctuations rather than inter-specific competition dominate the dynamics of several systems. This finding emphasizes the need to better understand how species are affected by the environment and acknowledge species differences in their responses to environmental heterogeneity, if we are to effectively model and predict their dynamics (e.g. for management and conservation purposes). The thesis also proposes a model-based approach to integrating the niche and neutral perspectives on community structure and dynamics, making it possible for the relative importance of each category of factors to be evaluated in light of field data

Modeling and Optimization of Dynamical Systems in Epidemiology using Sparse Grid Interpolation

Infectious diseases pose a perpetual threat across the globe, devastating communities, and straining public health resources to their limit. The ease and speed of modern communications and transportation networks means policy makers are often playing catch-up to nascent epidemics, formulating critical, yet hasty, responses with insufficient, possibly inaccurate, information. In light of these difficulties, it is crucial to first understand the causes of a disease, then to predict its course, and finally to develop ways of controlling it. Mathematical modeling provides a methodical, in silico solution to all of these challenges, as we explore in this work. We accomplish these tasks with the aid of a surrogate modeling technique known as sparse grid interpolation, which approximates dynamical systems using a compact polynomial representation. Our contributions to the disease modeling community are encapsulated in the following endeavors. We first explore transmission and recovery mechanisms for disease eradication, identifying a relationship between the reproductive potential of a disease and the maximum allowable disease burden. We then conduct a comparative computational study to improve simulation fits to existing case data by exploiting the approximation properties of sparse grid interpolants both on the global and local levels. Finally, we solve a joint optimization problem of periodically selecting field sensors and deploying public health interventions to progressively enhance the understanding of a metapopulation-based infectious disease system using a robust model predictive control scheme