Hello Everyone, I am fitting mixed models in which I have two phylogenies and I would like to understand them in the context of earlier work. In parafit (Legndre 2002) principal coordinates of the distance matrices (J-P) are used, where J is a matrix of ones and P the phylogenetic correlation matrix. For me, using unnormalised eigenvectors of solve(P) would result in nicer properties, and I wonder whether anyone could tell me whether there is some deep motivation for using the principal coordinates of J-A that I'm missing?
Kind Regards, Jarrod _______________________________________________ R-sig-phylo mailing list R-sig-phylo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-phylo