RelativeEigenanalysis: Relative Eigenanalysis

In evolqg: Evolutionary Quantitative Genetics

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

evolqg documentation built on Aug. 8, 2023, 5:12 p.m.