This repository contains the scripts used to run the ChromoSSE phylogenetic models of chromosome number evolution. These models allow one to infer the mode of chromosome number evolution; is chromosome evolution occurring primarily within lineages, primarily at lineage splitting, or in clade-specific combinations of both? Furthermore, they estimate the location and timing of possible chromosome speciation events over the phylogeny as well as ancestral chromosome numbers.
ChromoSSE was implemented in a Bayesian statistical framework, specifically in the software RevBayes, to accommodate uncertainty in parameter estimates while leveraging the full power of likelihood based methods. Provided here are the Rev scripts needed to run ChromoSSE in RevBayes. Please modify these scripts to use your own data, and feel free to contact me with any questions.
Please make sure you are using RevBayes v1.0.7 or greater. Some Rev language functions have changed and earlier versions will not work.
I've written a tutorial on how to use RevBayes for a wide range of chromosome evolution analyses using ChromEvol, BiChroM, ChromoSSE and other related models:
Freyman, W.A. and S. Höhna. 2017. Cladogenetic and anagenetic models of chromosome number evolution: a Bayesian model averaging approach. Systematic Biology syx065 [preprint] [journal] [pdf]
cd examples/aristolochia/
rb aristolochia_analysis.Rev
Use the RevGadgets R package and this example R script to generate figures like this:
Will Freyman, UC Berkeley
willfreyman@gmail.com
http://willfreyman.org
Copyright 2016 Will Freyman
License: GNU GPLv3 http://www.gnu.org/licenses/gpl.html