evolqg (version 0.3-4)

ExtendMatrix: Control Inverse matrix noise with Extension

Description

Calculates the extended covariance matrix estimation as described in Marroig et al. 2012

Usage

ExtendMatrix(cov.matrix, var.cut.off = 1e-04, ret.dim = NULL)

Value

Extended covariance matrix and second derivative variance

Arguments

cov.matrix

Covariance matrix

var.cut.off

Cut off for second derivative variance. Ignored if ret.dim is passed.

ret.dim

Number of retained eigenvalues

Author

Diogo Melo

References

Marroig, G., Melo, D. A. R., and Garcia, G. (2012). Modularity, noise, and natural selection. Evolution; international journal of organic evolution, 66(5), 1506-24. doi:10.1111/j.1558-5646.2011.01555.x

Examples

cov.matrix = RandomMatrix(11, 1, 1, 100)
ext.matrix = ExtendMatrix(cov.matrix, var.cut.off = 1e-6)
ext.matrix = ExtendMatrix(cov.matrix, ret.dim = 6)