Skip to main content

Advertisement

Log in

Comparing Covariance Matrices by Relative Eigenanalysis, with Applications to Organismal Biology

  • Tools and Techniques
  • Published:
Evolutionary Biology Aims and scope Submit manuscript

Abstract

Most biologists are familiar with principal component analysis as an ordination tool for questions about within-group or between-group variation in systems of quantitative traits, and with multivariate analysis of variance as a tool for one useful description of the latter in the context of the former. Less familiar is the mathematical approach of relative eigenanalysis of which both of these are special cases: computing linear combinations for which two variance–covariance patterns have maximal ratios of variance. After reviewing this common algebraic–geometric core, we demonstrate the effectiveness of this exploratory approach in studies of developmental canalization and the identification of divergent and stabilizing selection. We further outline a strategy for statistical classification when group differences in variance dominate over differences in group averages.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10

Similar content being viewed by others

References

  • Anderson, T. W. (1958). An introduction to multivariate statistical analysis. New York: Wiley.

    Google Scholar 

  • Anderson, T. W. (1963). Asymptotic theory for principal component analysis. The Annals of Mathematical Statistics, 34, 122–148.

    Article  Google Scholar 

  • Atchley, W., & Rutledge, J. (1980). Genetic components of size and shape. I. Dynamics of components of phenotypic variability and covariablity during ontogeny in the laboratory rat. Evolution, 34, 1161–1173.

    Article  Google Scholar 

  • Badyaev, A. V., & Foresman, K. R. (2004). Evolution of morphological integration. I. Functional units channel stress-induced variation in shrew mandibles. The American Naturalist, 163, 868–879.

    Article  PubMed  Google Scholar 

  • Bookstein, F. L. (1991). Morphometric tools for landmark data: Geometry and biology. Cambridge: Cambridge University Press.

    Google Scholar 

  • Bookstein, F. L. (2014). Measuring and reasoning: Numerical inferences in the sciences. Cambridge: Cambridge University Press.

  • Bookstein, F.L., Connor, P.D., Huggins, J.E., Barr, H. M., Pimentel, K. D., & Streissguth, A. P. (2007). Many infants prenatally exposed to high levels of alcohol show one particular anomaly of the corpus callosum. Alcoholism: Clinical and Experimental Research, 31, 868–879.

    Article  CAS  Google Scholar 

  • Bookstein, F. L., Streissguth, A. P., Sampson, P. D., Connor, P. D., & Bar, H. M. (2002). Corpus callosum shape and neuropsychological deficits in adult males with heavy fetal alcohol exposure. Neuroimage, 15, 233–251.

    Article  PubMed  Google Scholar 

  • Cheverud, J. M. (1988). A comparison of genetic and phenotypic correlations. Evolution, 42, 958–968.

    Article  Google Scholar 

  • Coquerelle, M., Bookstein, F. L., Braga, J., Halazonetis, D. J., Weber, G. W., & Mitteroecker, P. (2011). Sexual dimorphism of the human mandible and its association with dental development. American Journal of Physical Anthropology, 145, 192–202.

    Article  PubMed  Google Scholar 

  • Debat, V., & David, P. (2001). Mapping phenotypes: Canalization, plasticity and developmental stability. Trends in Ecology & Evolution, 16, 555–561.

    Article  Google Scholar 

  • Falconer, D. S., & Mackay, T. F. C. (1996). Introduction to quantitative genetics. Essex: Longman.

    Google Scholar 

  • Felsenstein, J. (1985). Phylogenies and the comparative method. American Naturalist, 125, 1–15.

    Article  Google Scholar 

  • Felsenstein, J. (1988). Phylogenies and quantitative characters. Annual Review of Ecology, Evolution, and Systematics, 19, 445–471.

    Article  Google Scholar 

  • Flury, B. N. (1983). Some relations between the comparison of covariance matrices and principal component analysis. Computational Statistics & Data Analysis, 1, 97–109.

    Article  Google Scholar 

  • Flury, B. N. (1985). Analysis of linear combinations with extreme ratios of variance. Journal of the American Statistical Association, 80, 915–922.

    Article  Google Scholar 

  • Förstner W., & Moonen, B. (1999). A metric for covariance matrices. In: F. Krumm, V. S. Schwarze (Eds.), Quo vadis geodesia ...?, Festschrift for Erik W. Grafarend on the occasion of his 60th birthday. Stuttgart: Stuttgart University.

  • Gibson, G., & Wagner, G. (2000). Canalization in evolutionary genetics: A stabilizing theory? Bioessays, 22, 372–380.

  • Hallgrimsson, B., Brown, J. J., Ford-Hutchinson, A. F., Sheets, H. D., Zelditch, M. L., & Jirik, F. R. (2006). The brachymorph mouse and the developmental-genetic basis for canalization and morphological integration. Evolution & Development, 8, 61–73.

    Article  CAS  Google Scholar 

  • Hallgrimsson, B., & Hall, B. K. (2005). Variation: A central concept in biology. New York: Elsevier Academic Press.

    Google Scholar 

  • Hallgrimsson, B., Willmore, K., & Hall, B. K. (2002). Canalization, developmental stability, and morphological integration in primate limbs. American Journal of Physical Anthropology Supplement, 35, 131–158.

    Article  Google Scholar 

  • Hansen, T. F., & Houle, D. (2008). Measuring and comparing evolvability and constraint in multivariate characters. Journal of Evolutionary Biology, 21, 1201–1219.

    Article  CAS  PubMed  Google Scholar 

  • Houle, D., & Fierst, J. (2013). Properties of spontaneous mutational variance and covariance for wing size and shape in Drosophila melanogaster. Evolution, 67, 1116–1130.

    Article  PubMed  Google Scholar 

  • Howells, W. W. (1996). Howells craniometric data on the internet. American Journal of Physical Anthropology, 101, 441–442.

    Article  CAS  PubMed  Google Scholar 

  • Huttegger, S., & Mitteroecker, P. (2011). Invariance and meaningfulness in phenotype spaces. Evolutionary Biology, 38, 335–352.

    Article  Google Scholar 

  • Klingenberg, C. P., Debat, V., & Roff, D. A. (2010). Quantitative genetics of shape in cricket wings: Developmental integration in a functional structure. Evolution, 64, 2935–2951.

    PubMed  Google Scholar 

  • Koots, K. R., & Gibson, J. P. (1996). Realized sampling variances of estimates of genetic parameters and the difference between genetic and phenotypic correlations. Genetics, 143:1409–1416.

    CAS  PubMed Central  PubMed  Google Scholar 

  • Lande, R. (1979). Quantitative genetic analysis of multivariate evolution, applied to brain: Body size allometry. Evolution, 33, 402–416.

    Article  Google Scholar 

  • Manly, B. F. J., & Rayner, J. C. W. (1987). The comparison of sample covariance matrices using likelihood ratio tests. Biometrika, 74, 841–847.

    Article  Google Scholar 

  • Mardia, K. V., Kent, J. T., & Bibby, J. M. (1979). Multivariate analysis. London: Academic Press.

    Google Scholar 

  • Martin, G., Chapuis, E., & Goudet, J. (2008). Multivariate QST-FST comparisons: A neutrality test for the evolution of the g matrix in structured populations. Genetics, 180, 2135–2149.

    Article  PubMed Central  PubMed  Google Scholar 

  • Mitteroecker, P. (2009). The developmental basis of variational modularity: Insights from quantitative genetics, morphometrics, and developmental biology. Evolutionary Biology, 36, 377–385.

    Article  Google Scholar 

  • Mitteroecker, P., & Bookstein, F. L. (2009). The ontogenetic trajectory of the phenotypic covariance matrix, with examples from craniofacial shape in rats and humans. Evolution, 63, 727–737.

    Article  PubMed  Google Scholar 

  • Mitteroecker, P., & Bookstein, F. L. (2011). Classification, linear discrimination, and the visualization of selection gradients in modern morphometrics. Evolutionary Biology, 38, 100–114.

    Article  Google Scholar 

  • Mitteroecker, P., Gunz, P., Neubauer, S., & Müller, G. B. (2012). How to explore morphological integration in human evolution and development? Evolutionary Biology, 39, 536–553.

    Article  Google Scholar 

  • Mitteroecker, P., & Huttegger, S. (2009). The concept of morphospaces in evolutionary and developmental biology: Mathematics and metaphors. Biological Theory, 4, 54–67.

    Article  Google Scholar 

  • Morrison, D. F. (1976). Multivariate statistical methods. New York: McGraw-Hill.

    Google Scholar 

  • Nonaka, K., & Nakata, M. (1984). Genetic variation and craniofacial growth in inbred rats. Journal of Craniofacial Genetics and Developmental Biology, 4, 271–302.

    CAS  PubMed  Google Scholar 

  • Philipps, P. C., & Arnold, S. J. (1989). Visualizing multivariate selection. Evolution, 43, 1209–1222.

    Article  Google Scholar 

  • Rao, C. R. (1948). The utilization of multiple measurements in problems of biological classification. Journal of the Royal Statistical Society. Series B, 10, 159–203.

    Google Scholar 

  • Roff, D. (1995). The estimation of genetic correlations from phenotypic correlations: A test of Cheverud’s conjecture. Heredity, 74, 481–490.

    Article  Google Scholar 

  • Roff, D. (2000). The evolution of the G matrix: Selection or drift? Heredity (Edinb), 84, 135–142.

    Article  Google Scholar 

  • Smith, S. T. (2005). Covariance, subspace, and intrinsic Cramer-Rao bounds. IEEE Transactions on Signal Processing, 53, 1610–1630.

    Article  Google Scholar 

  • Tanner, J. M. (1963). Regulation of growth in size in mammals. Nature, 199, 845–850.

    Article  CAS  PubMed  Google Scholar 

  • Tyler, D. E., Critchley, F., Dümbgen, L., & Oja, H. (2009). Invariant co-ordinate selection. Journal of the Royal Statistical Society: Series B, 71, 549–592.

    Article  Google Scholar 

  • Zelditch, M. L., Bookstein, F. L., & Lundrigan, B. (1992). Ontogeny of integrated skull growth in the cotton rat Sigmodon fulviventer. Evolution, 46, 1164–1180.

    Article  Google Scholar 

  • Zelditch, M. L., Lundrigan, B. L., & Garland, T. (2004). Developmental regulation of skull morphology. I. Ontogenetic dynamics of variance. Evolution & Development, 6, 194–206.

    Article  Google Scholar 

  • Zelditch, M. L., Mezey, J. G., Sheets, H. D., Lundrigan, B. L., & Garland, T. (2006). Developmental regulation of skull morphology II: Ontogenetic dynamics of covariance. Evolution & Development, 8, 46–60.

    Article  Google Scholar 

Download references

Acknowledgments

This research was supported by the Focus of Excellence grant “Biometrics of EvoDevo” from the Faculty of Life Sciences, University of Vienna, to Philipp Mitteroecker, and Grant DEB-1019583 to Fred Bookstein and Joseph Felsenstein from the National Sciences Foundation of the United States. We thank Katharina Puschnig for drawing the skull used in Figs. 6 and 8.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Philipp Mitteroecker.

Appendix

Appendix

Corresponding to the maneuvers here is a null model that is often helpful in deciding whether a relative eigenanalysis reveals anything worth thinking about. We’re not testing the ratio, only asking that it be at least as large as the value expected on the absence of any signal. The approach parallels the analogous decision regarding eigenvectors arising from successive eigenvalues of one single matrix (Coquerelle et al. 2011; Bookstein 2014). That approach, in turn, is a modification of the stepdown test for sphericity that is already in the advanced textbooks (cf. Morrison 1976, pp. 336–337). As this argument has not apparently been published before, we sketch it here for any interested reader. Our notation is borrowed from Anderson (1963), where several of the corresponding asymptotic likelihood ratio tests were published.

Suppose we are observing a covariance matrix on a sample of size N for a list of p variables that are all actually independent Gaussians of mean 0 and variance 1. Let the matrix \(\mathbf{U}\) be \(\sqrt{N}\) times the deviation of the empirically observed covariance matrix from the correct answer, which is the identity matrix of rank p. As samples grow large, \(\mathbf{U}\) becomes Gaussian with mean zero for every element and variance 2 along its diagonal, 1 elsewhere.

That theorem corresponds to a \(\mathbf{U}\) that was generated to describe some multivariate Gaussian distribution’s ordinary eigenvalues. The relative eigenvalues, which we have been depicting here as the eigenvalues of \(\mathbf{T}^{-1/2}\mathbf{S}\mathbf{T}^{-1/2}\), are also the eigenvalues of \(\mathbf{S}\mathbf{T}^{-1}\). For large samples, the distribution of the deviation of \(\mathbf{T}^{-1}\) from the identity is the same as the distribution of the deviation of \(\mathbf{T}\). The effect of this additional factor of \(\mathbf{T}^{-1}\) on what was already the deviation of the sample described by \(\mathbf{S}\) is to alter \(\mathbf{U}\) by another term, additive in this metric, of exactly the same distribution. Their sum, scaled by \(\sqrt{N}\), thus is in the limit a set of Gaussians with terms of mean zero and variances 4 down the diagonal, 2 elsewhere.

At equation 3.9, pages 132–133 of this same paper, Anderson shows, by expressing the eigenvalues themselves in terms of the elements of \(\mathbf{U}\), that for the ordinary eigenproblem the log likelihood ratio for the null hypothesis of sphericity for any q consecutive eigenvalues (not necessarily the full set of all p corresponding to the p original variables) is the quantity log a/g, where a is the arithmetic mean of the eigenvalues and g is their geometric mean. Using the theorem about \(\mathbf{U}\) for ordinary eigenvalues, he shows that in the limit of large samples this log likelihood ratio is distributed approximately as 1/Nq times a χ2 on (q − 1)(q + 2)/2 degrees of freedom, where q is the number of eigenvalues being compared (in our applications, usually, q = 2). It follows, then, that in the relative eigenanalysis application the same quantity log a/g is distributed as 2/Nq times the same χ2 (Fig. 11).

Fig. 11
figure 11

From the expected value of the corresponding χ2 distributions comes a permissive criterion for narrative validity of relative eigenvectors: restrict the text only to those for which the relative eigenvalue λ i has a ratio to its successor λ i+1 as least as large as the value given in this diagram. Solid line curve for ordinary eigenvalues. Dotted line curve for relative eigenvalues

Rights and permissions

Reprints and permissions

About this article

Cite this article

Bookstein, F.L., Mitteroecker, P. Comparing Covariance Matrices by Relative Eigenanalysis, with Applications to Organismal Biology. Evol Biol 41, 336–350 (2014). https://doi.org/10.1007/s11692-013-9260-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s11692-013-9260-5

Keywords

Navigation