Hello, I'm sorry for not specifying that. But yes, that's exactly it. Thanks for the help.
Cheers, Sérgio. 2015-06-15 12:47 GMT+01:00 Liam J. Revell <[email protected]>: > Hi Sergio. > > You didn't specify what BM & vf are, but if: > > BM<-corBrownian(1,phy) > ## and > vf<-diag(vcv(phy)) > > then, yes, this is correct. There is a function in phytools, phy.resid, > that also does this calculation. > > All the best, Liam > > Liam J. Revell, Assistant Professor of Biology > University of Massachusetts Boston > web: http://faculty.umb.edu/liam.revell/ > email: [email protected] > blog: http://blog.phytools.org > > > On 6/15/2015 7:41 AM, Sergio Ferreira Cardoso wrote: > >> Dear all, >> >> I just read this paper: Revell, L. J. (2009). Size‐correction and >> principal >> components for interspecific comparative studies. Evolution, 63(12), >> 3258-3268. There is no coding for R, but it says it can be made using APE. >> >>> From what I could understand this is as simples as: >>> >> >> ->bm.gls<-gls(trait~Bodymass, >> correlation=BM,data=DF,weights=varFixed(~vf)) >> ->residuals(bm.gls) >> >> Is this correct? Or do I need to ask for some special kind of residuals? >> >> Best regards, >> Sérgio. >> >> >> -- Com os melhores cumprimentos, Sérgio Ferreira Cardoso. -------------------- Best regards, Sérgio Ferreira Cardoso MSc. Paleontology candidate Departamento de Ciências da Terra - FCT /Universidade Nova de Lisboa Geociências - Universidade de Évora Lisboa, Portugal [[alternative HTML version deleted]] _______________________________________________ R-sig-phylo mailing list - [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-phylo Searchable archive at http://www.mail-archive.com/[email protected]/
