This function performs phylogenetic canonical correlation analysis (e.g., Revell & Harrison 2008; Bioinformatics).
Usage
phyl.cca(tree, X, Y, lambda=1.0, fixed=TRUE)
Arguments
tree
a phylogenetic tree in "phylo" format.
X
a data matrix with traits in columns.
Y
data matrix with traits in columns, to be correlated with X.
lambda
optionally, a (fixed) value for lambda.
fixed
optionally, a logical value indicating whether or not to estimate lambda using likelihood.
Value
A list with the following components:
corcanonical correlations.
xcoefcoefficients for the canonical variables for X.
ycoefcoefficients for the canonical variables for Y.
xscoresmatrix with the canonical scores for X.
yscoresmatrix with the canonical scores for Y.
chisqvector of $\chi^2$ values.
pP-values for the hypothesis test that the ith and all subsequent correlations are zero.
Details
(Optional) joint optimization of $\lambda$ is performed using optimize on the interval (0,1).
References
Revell, L. J., Harrison, A. S. (2008) PCCA: A program for phylogenetic canonical correlation analysis. Bioinformatics, 24, 1018-1020.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.