Skip to content

bsipos/phylosim

This branch is 17 commits ahead of sbotond/phylosim:master.

Folders and files

NameName
Last commit message
Last commit date

Latest commit

c110136 · Mar 15, 2020
Mar 1, 2020
Jan 6, 2014
Sep 6, 2016
Feb 6, 2012
Nov 4, 2017
Sep 6, 2016
Mar 15, 2020
Oct 27, 2010
Oct 8, 2010
Sep 4, 2016
Jan 6, 2014
Oct 8, 2010
Jun 12, 2013
May 27, 2009
May 27, 2009
Jan 6, 2014
Mar 20, 2013
Jan 6, 2014
Oct 8, 2010
Oct 8, 2010
Oct 8, 2010
Oct 8, 2010
Jan 6, 2014
Jun 12, 2013
Mar 15, 2020
Jun 12, 2013
Sep 4, 2016
Nov 19, 2019
Sep 4, 2016
Dec 3, 2012
Nov 18, 2019
Oct 8, 2010
Oct 11, 2010
Jan 31, 2012
Oct 8, 2010
Mar 15, 2020
Jun 12, 2013
Oct 13, 2010
Sep 4, 2016

Repository files navigation

Gillespie algorithm

PhyloSim

PhyloSim is an extensible object-oriented framework for the Monte Carlo simulation of sequence evolution written in 100 percent R. It is built on the top of the R.oo and ape packages and uses the Gillespie algorithm to simulate substitutions, insertions and deletions.

PhyloSim was brought to you by the Goldman group from EMBL-EBI.

Catalogued on GSR

Publication

Botond Sipos, Tim Massingham, Gregory E Jordan and Nick Goldman (2011) PhyloSim - Monte Carlo simulation of sequence evolution in the R statistical computing environment - BMC Bioinformatics 12:104 doi:10.1186/1471-2105-12-104

Download an install

The released packages are available from CRAN.

Key features

  • Simulation of the evolution of a set of discrete characters with arbitrary states evolving by a continuous-time Markov process with an arbitrary rate matrix.

  • Explicit implementations of the most popular substitution models (nucleotide, amino acid and codon substitution models).

  • Simulation under the popular models of among-sites rate variation, like the gamma (+G) and invariant sites plus gamma (+I+G) models.

  • The possibility to simulate under arbitrarily complex patterns of among-sites rate variation by setting the site specific rates according to any R expression.

  • Simulation of one or more separate insertion and/or deletion processes acting on the sequences and which sample the insertion/deletion length from an arbitrary discrete distribution or an R expression (so all the probability distributions implemented in R are readily available for this purpose).

  • Simulation of the effects of variable functional constraints over the sites by site-process specific insertion and deletion tolerance parameters which determine the rejection probability of a proposed insertion/deletion.

  • The possibility of having a different set of processes and site-process specific parameters for every site, which allows for an arbitrary number of partitions in the simulated data.

  • The possibility to evolve sites by a combination of substitution processes along a single branch.

  • Simulation of heterotachy and other cases of non-homogeneous evolution by allowing the user to set "node hook" functions altering the site properties at internal nodes.

  • The possibility to export the counts of various events ("branch statistics") as phylo objects (see the man page of exportStatTree.PhyloSim).

  • See the man page of the PhyloSim class and the package vignette for more features and examples.

Building from source

The package can be built from the source by issuing make pack on a *nix system. The building process need the standard unix tools, Perl and R with the ape, R.oo, and ggplot2 packages installed.

About

An R package for the Monte Carlo simulation of sequence evolution

Resources

License

Stars

Watchers

Forks

Packages

No packages published

Languages

  • R 99.9%
  • Other 0.1%