No description, website, or topics provided.
Branch: master
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Type Name Latest commit message Commit time
Failed to load latest commit information.
man BEAST and sampling-intensity helper functions. Jul 25, 2018
vignettes small correction to plot results May 28, 2017
NAMESPACE BEAST and sampling-intensity helper functions. Jul 25, 2018 Update Jul 13, 2017



The purpose of phylodyn is to facilitate phylodynamic inference and analysis in an approachable R package.


  1. Install (if necessary) package dependencies and helpers ape, spam and devtools using install.packages.

  2. Install INLA using install.packages("INLA", repos="")

  3. Load devtools using library(devtools).

  4. Install phylodyn using

    a. install_github("mdkarcher/phylodyn"), or

    b. install_github("mdkarcher/phylodyn", build_vignettes = TRUE) if you want some illustrative vignettes (note: using build_vignettes = TRUE will make the install take longer).


  1. SimpleBNPR: A short example showing how to use BNPR and BNPR-PS on simulated data, illustraring methodology in [2] and [5].

  2. NewYorkInfluenza: A case study analyzing influenza data from New York, reproducing analysis in [5] on data from [1].

  3. RegionalInfluenza: A case study analyzing influenza data from nine geographic regions, reproducing analsyis in [5] on data from [3].

  4. RegionalSeasonality: A case study analyzing influenza seasonality from nine geographic regions, reproducing analsyis in [5] on data from [3].

  5. SimplePhyloinfer: A short example comparing BNPR with a split HMC MCMC sampler approach, illustrating methodology in [4].

  6. LongPhyloinfer: A longer example comparing BNPR with multiple MCMC samplers, including split HMC as in SimplePhyloinfer, illustrating methodology in [4].

  7. LocalGenealogies: A short example of MCMC-based inference of effective population size trajectories from a sequence of local genealogies. Genealogies are assumed to be a realization of the Sequentially Markov Coalescent (SMC') model. The methodology is developed in [6]


Datasets below can be found at:

  1. New York influenza BEAST XML for inferring genealogy using sequence data from [1].

    • NewYork.xml
  2. Regional influenza BEAST XML for inferring genealogy using sequence data from [3].

    • Europe.xml
    • India.xml
    • JapanKorea.xml
    • NorthChina.xml
    • Oceania.xml
    • SouthAmerica.xml
    • SouthChina.xml
    • SoutheastAsia.xml
    • USACanada.xml


  1. A. Rambaut, O. G. Pybus, M. I. Nelson, C. Viboud, J. K. Taubenberger, E. C. Holmes The genomic and epidemiological dynamics of human influenza A virus. Nature, 453(7195): 615–619, 2008.

  2. J. A. Palacios and V. N. Minin. Integrated nested Laplace approximation for Bayesian nonparametric phylodynamics. In Proceedings of the Twenty-Eighth International Conference on Uncertainty in Artificial Intelligence, pages 726–735, 2012.

  3. D. Zinder, T. Bedford, E. B. Baskerville, R. J. Woods, M. Roy, M. Pascual. Seasonality in the migration and establishment of H3N2 Influenza lineages with epidemic growth and decline. BMC Evolutionary Biology, 14(1): 272, 2014.

  4. S. Lan, J. A. Palacios, M. Karcher, V. N. Minin, and B. Shahbaba. An Efficient Bayesian Inference Framework for Coalescent-Based Nonparametric Phylodynamics, Bioinformatics, 31(20): 3282-3289, 2015.

  5. M. D. Karcher, J. A. Palacios, T. Bedford, M. A. Suchard, and V. N. Minin. Quantifying and mitigating the effect of preferential sampling on phylodynamic inference. PLOS Computational Biology, 12:e1004789, 2016.

  6. J.A. Palacios, J. Wakeley, and S. Ramachandran. Bayesian nonparametric inference of population size changes from sequential genealogies. Genetics Vol. 201:281-304, 2015.

  7. M. Karcher M, J.A. Palacios, S. Lan, V.N. Minin. phylodyn: an R package for phylodynamic simulation and inference, Molecular Ecology Resources, 17, 96-100, 2017.