[R-sig-phylo] brownian motion, residuals and PCA

2013-08-21 Thread Miriam Santos-Herrera
Dear all, Im an newbie in this world of comparative methods and excuse me in advance if I am asking something very obvious. I have performed phylogenetic size correction and a phylogenetic pca to my data using phyl.resid and phyl.pca. My concern is that these functions assume that my traits

Re: [R-sig-phylo] brownian motion, residuals and PCA

2013-08-21 Thread Liam J. Revell
Hi Miriam. It is the correlation structure (phylogenetic signal) of the residual error of y|x that matters, not the structure of the individual variables. You can use an OU model for the correlation structure of the residual error in PGLS. This would look something like this (if 'yourdata'