Re: [R-sig-phylo] Doubts and problems (bug?) when using comparative.data and pgls functions in caper‏

2011-12-14 Thread Orme, David
Hi Edwin, That is a bug - I've fixed it in the developmental version on R-forge but haven't yet rolled out a new version to CRAN. Many thanks for reporting it! Phylogenetic independent contrasts, where comparisons need to be made between pairs of branches, require methods adjustments for

[R-sig-phylo] split data set and topology on PGLS

2011-12-14 Thread Renata Brandt
Hello everybody. I have some questions for you. When performing a phylogenetic regression between two continuous traits, the dependent trait clearly divides my data set in two distinct morphologies. This division overlaps with what I previously thought were two different ecomorphs. The

Re: [R-sig-phylo] split data set and topology on PGLS

2011-12-14 Thread Liam J. Revell
Hi Renata. I agree with Ted - you should be able to fit an ANCOVA model with gls(y~x+ecomorph,...,correlation=corBrownian(phy=tree)), or something like this, where ecomorph is a factor. Alternatively, you could use the method we devised in Revell Collar (2009; Evolution) in which we first