RelativeEigenanalysis function

Relative Eigenanalysis

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

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

Arguments

• cov.x: covariance matrix
• cov.y: covariance matrix
• symmetric: Logical. If TRUE, computes symmetric eigenanalysis.

Returns

list with two objects: eigenvalues and eigenvectors

Examples

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

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

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.

Author(s)

Guilherme Garcia, Diogo Melo