Publications
You can run, you can hide: The epidemiology and statistical mechanics of zombies
We use a popular fictional disease, zombies, in order to introduce techniques used in modern epidemiology modeling, and ideas and techniques used in the numerical study of critical phenomena. We consider variants of zombie models, from fully connected continuous time dynamics to a full scale exact stochastic dynamic simulation of a zombie outbreak on the continental United States.
Overshoot during phenotypic switching of cancer cell populations
The dynamics of tumor cell populations is hotly debated: do populations derive hierarchically from a subpopulation of cancer stem cells (CSCs), or are stochastic transitions that mutate differentiated cancer cells to CSCs important? Here we argue that regulation must also be important. We sort human melanoma cells using three distinct cancer stem cell (CSC) markers-CXCR6, CD271 and ABCG2-and observe that the fraction of non-CSC-marked cells first overshoots to a higher level and then returns to the level of unsorted cells.
Perspective: Sloppiness and emergent theories in physics, biology, and beyond
Large scale models of physical phenomena demand the development of new statistical and computational tools in order to be effective. Many such models are "sloppy," i.e., exhibit behavior controlled by a relatively small number of parameter combinations. We review an information theoretic framework for analyzing sloppy models. This formalism is based on the Fisher information matrix, which is interpreted as a Riemannian metric on a parameterized space of models. Distance in this space is a measure of how distinguishable two models are based on their predictions.
Variational method for estimating the rate of convergence of Markov-chain Monte Carlo algorithms
We demonstrate the use of a variational method to determine a quantitative lower bound on the rate of convergence of Markov chain Monte Carlo (MCMC) algorithms as a function of the target density and proposal density. The bound relies on approximating the second largest eigenvalue in the spectrum of the MCMC operator using a variational principle and the approach is applicable to problems with continuous state spaces.
Sloppiness, robustness, and evolvability in systems biology
The functioning of many biochemical networks is often robust - remarkably stable under changes in external conditions and internal reaction parameters. Much recent work on robustness and evolvability has focused on the structure of neutral spaces, in which system behavior remains invariant to mutations. Recently we have shown that the collective behavior of multiparameter models is most often sloppy: insensitive to changes except along a few 'stiff' combinations of parameters, with an enormous sloppy neutral subspace.