[R-sig-phylo] Standardising independent contrasts

2009-09-15 Thread Manabu Sakamoto
Dear list, I am trying to standardise independent contrasts following Garland et al. (1992), but I am unsure just how you would extract the sum of it's [contrast's] branch lenghts in R. I computed independent contrasts using the pic function in ape. I then took the contrasts and subtracted the

Re: [R-sig-phylo] Standardising independent contrasts

2009-09-15 Thread Liam J. Revell
The contrasts returned by pic() should already be standardized. Standardizing phylogenetic contrasts is performed by dividing each contrast (calculated between internal or terminal nodes) by a sum proportional to the square root of its expected variance. Assuming Brownian motion, for two