Learn R Programming

phytools (version 0.2-40)

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. & A. S. Harrison. 2008. PCCA: A program for phylogenetic canonical correlation analysis. Bioinformatics, 24, 1018--1020.

See Also

phyl.pca