[This article was first published on Recology
, and kindly contributed to R-bloggers
]. (You can report issue about the content on this page here
Want to share your content on R-bloggers? click here
if you have a blog, or here
if you don't.
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.
Get the code at my website:
(FYI: there are emticons in the code on the website,
but it pastes fine into a text editor for me)
-This new code works with Marc’s new version of Phylometa: http://lajeunesse.myweb.usf.edu/publications.html
Again, please let me know if it doesn’t work, if it’s worthless, what changes could make it better.