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.
_______________________________________________
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]/
