Code and software
R
- I am the lead developer of treess, a package for diagnosing the convergence of MCMC samples of phylogenies and other objects beyond the reach of standard approaches. Provides Fréchet generalizations of the effective sample size and visual diagnostics.
- I am the lead developer of nbbp, a package for inference of the effective reproduction number from final size outbreak data under a negative binomial branching process model.
- R code for two computational labs I helped to develop is available here. One lab explores the Luria-Delbruck experiment and modifications thereof. This lab is designed to work in conjunction with a wet lab recapitulating the Luria-Delbruck experiment, but can be done with any data from such an experiment (including the original data from Luria and Delbruck). The other lab explores models of population cycling, and includes all necessary datasets.
python
- I am the lead developer of cladecombiner, a package leveraging phylogenetic taxonomy to streamline the process of aggreagating viral lineages for nowcasting and short-term forecasting of variant dynamics.
- I have contributed to the development of tooling to test the efficacy of variant nowcasting models via retrospective analyses.
- I helped create a simulation based python streamlit app for exploring the behavior of branching processes in the presence of multiple observation processes.
Other
RevBayes
RevBayes is a probabilistic programming language written in C++. It offers a friendly scripting-based interface for model specification, a large number of phylogenetic models, and highly customizable MCMC. In order to streamline the user experience, the original repository and its history was archived and the rebooted version is now in use. Among other things, I have contributed
- A general time-varying model of a partially-observed branching process (now here, originally here). I also wrote a tutorial on using this model.
- An implementation of the elliptical slice sampler (now here, originally here)
BEAST X
BEAST X is an extraordinarily long-lived scientific software project for phylogenetic inference, in service since the early aughts. Users can specify a wide variety of phylogenetic models through an XML interface, while the source code itself is written in java and many of the more computationally intensive tasks are offloaded to C++ via the BEAGLE library. Among other things, I have
- Assisted in the implementation of efficient and differentiable likelihoods for partially-observed branching processes.
- Written code for tracking the implied total number of infected individuals from a partially-observed branching process.
- Implemented the ExpGamma distribution.