RelativeEigenanalysis {evolqg}

Relative Eigenanalysis

Description

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

Usage

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

Arguments

cov.x

covariance matrix

cov.y

covariance matrix

symmetric

Logical. If TRUE, computes symmetric eigenanalysis.

Value

list with two objects: eigenvalues and eigenvectors

Author(s)

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]])