plot_age_posterior_distributions.R

2012-10-20T00:52:02Z (GMT) by Carlos Peña
<p>This script plots density distributions of age nodes of phylogenetic trees in R.</p> <p>It appears that different genes might estimate <strong>very</strong> different age estimates for your phylogenetic trees. This seems to be the case with the COI gene that tends to pull your timings towards the past. Saturation in the 3rd codon position might be accused for responsibility  (you can inspect the saturation level with <a href="http://nsg-databases.blogspot.fi/2012/10/plots-de-saturacion-para-secuencias.html">some cool plots</a>).</p> <p><br><a href="http://dx.doi.org/10.1093/sysbio/syq045">Brandley et al (2011</a>) use a nice way to figure out whether some of your genes might be giving very different time estimates for your trees. You can select 1000 random trees from your BEAST run and plot the distributions of the ages for the crown group of different genes, different codon positions and the combined analyses.<br><br>Something like this plot consisting on a simulation of a gen1 estimating a crown age of 30Mya, gen2 estimating an age of 50Mya and the combined analysis giving an age of 40Mya.</p>