Dear all, I was trying to find the pylogenetic signal for a trait and I used phylosig():
phylosig(tree, trait, method="lambda", test=TRUE, nsim=999) If I use, for instance, a lambda correlation in a gls, am I expected to obtain the same value? gls(trait~var1+var2, corPagel(value=1,phy=tree,fixed=F)...) Is the lambda estimated for a model the same as the one calculated for a trait alone? If I have a different result, what's it due to? Thanks in advance. -- Com os melhores cumprimentos, Sérgio Ferreira Cardoso. -------------------- Best regards, Sérgio Ferreira Cardoso MSc. Paleontology Museu da Lourinhã, GEAL. LATR/IST/CTN - Campus Tecnológico e Nuclear. Lisboa, Portugal ᐧ [[alternative HTML version deleted]] _______________________________________________ R-sig-phylo mailing list - R-sig-phylo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-phylo Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/