Recology

R/etc.

Phylometa from R - UDPATE

    ggplot2 meta-analysis Phylogenetics R
 Source: .Rmd/.md

A while back I posted some messy code to run Phylometa from R, especially useful for processing the output data from Phylometa which is not easily done. The code is still quite messy, but it should work now. I have run the code with tens of different data sets and phylogenies so it should work.

I fixed errors when parentheses came up against numbers in the output, and other things. You can use the code for up to 4 levels of your grouping variable. In addition, there are some lines of code to plot the effect sizes with confidence intervals, comparing random and fixed effects models and phylogenetic and traditional models.

The code is in a Gist at https://gist.github.com/sckott/939970, and embedde below:

Use the first file to do the work, calling the second file using source().

This new code works with Marc’s new version of Phylometa, so please update: http://lajeunesse.myweb.usf.edu

Again, please let me know if it doesn’t work, if it’s worthless, what changes could make it better.

Some notes on tree formatting for Phylometa.

comments powered by Disqus