evolqg (version 0.3-4)

RelativeEigenanalysis: Relative Eigenanalysis

Description

Computes relative eigenvalues and eigenvectors between a pair of covariance matrices.

Usage

RelativeEigenanalysis(cov.x, cov.y, symmetric = FALSE)

Value

list with two objects: eigenvalues and eigenvectors

Arguments

cov.x

covariance matrix

cov.y

covariance matrix

symmetric

Logical. If TRUE, computes symmetric eigenanalysis.

Author

Guilherme Garcia, Diogo Melo

References

Bookstein, F. L., and P. Mitteroecker, P. "Comparing Covariance Matrices by Relative Eigenanalysis, with Applications to Organismal Biology." Evolutionary Biology 41, no. 2 (June 1, 2014): 336-350. doi:10.1007/s11692-013-9260-5.

Examples

data(dentus)
dentus.vcv <- dlply(dentus, .(species), function(df) var(df[, -5]))

dentus.eigrel <- RelativeEigenanalysis(dentus.vcv [[1]], dentus.vcv[[5]])