[R-sig-phylo] boxplots next to phylogeny

2015-07-14 Thread Martin Bulla
Dear list members, I wonder whether anyone of you has a code for phyldataplot (or any other plotting function) that allows for plotting the distributions of variable (e.g., boxplots) next to the phylogeny giving something like:

Re: [R-sig-phylo] boxplots next to phylogeny

2015-07-14 Thread Liam J. Revell
Hi Martin. The easiest way to do this is to create two side-by-side plots using par(mfcol=c(1,2)) and then put your tree your boxplot in the two plots; however one has to be very careful in that case to ensure that the tips align to the correct boxes. I posted a demo of this to my blog. Let

Re: [R-sig-phylo] Different estimate of Bloomberg's k in ape and caper in R

2015-07-14 Thread Slater, Graham
Hi Solomon, It’s because neither of the things you are estimating are Blomberg’s K and both are different from one another. In Caper, you are estimating Pagel’s Kappa (a branch length transformation that attempts to explain evolutionary change as punctuational [1] vs gradual [=1] ), while in

[R-sig-phylo] Different estimate of Bloomberg's k in ape and caper in R

2015-07-14 Thread Solomon Chak
Hi everyone, I'd like to get some opinions on why estimates of Bloomberg's k are different using the packages ape and caper in R. Please see the code below for example, in which I used the same data but got different estimates of k. In my own data, the estimated ks were completely opposite: