Learn R Programming

phytools (version 0.4-45)

phyl.cca: Phylogenetic canonical correlation analysis

Description

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.

See Also

phyl.pca