A note on simulating null distributions for G matrix comparisons
Autor: | Keyne Monro, Sandra Hangartner, Michael B. Morrissey |
---|---|
Přispěvatelé: | The Royal Society, University of St Andrews. School of Biology, University of St Andrews. Centre for Biological Diversity |
Rok vydání: | 2019 |
Předmět: |
0106 biological sciences
0301 basic medicine QH301 Biology Value (computer science) Biology 010603 evolutionary biology 01 natural sciences QH301 03 medical and health sciences Matrix (mathematics) Null distribution Simple (abstract algebra) Genetics Humans Applied mathematics Computer Simulation QA Mathematics QA Ecology Evolution Behavior and Systematics Tensor analysis Models Genetic Quantitiatve genetics Null (mathematics) Genetic Variation 3rd-DAS Covariance Biological Evolution 030104 developmental biology Differentiation G matrix Variety (universal algebra) General Agricultural and Biological Sciences Adaptive evolution |
Zdroj: | Evolution. 73:2512-2517 |
ISSN: | 1558-5646 0014-3820 |
DOI: | 10.1111/evo.13842 |
Popis: | MBM is supported by a University Research Fellowship from the Royal Society (London). KM is supported by a Future Fellowship from the Australian Research Council. Genetic variances and covariances, summarised in G matrices, are key determinants of the course of adaptive evolution. Consequently, understanding how G matrices vary among populations is critical to answering a variety of questions in evolutionary biology. A method has recently been proposed for generating null distributions of statistics pertaining to differences in G matrices among populations. The general approach facilitated by this method is likely to prove to be very important in studies of the evolution of G . We have identified an issue in the method that will cause it to create null distributions of differences in G matrices that are likely to be far too narrow. The issue arises from the fact that the method as currently used generates null distributions of statistics pertaining to differences in G matrices across populations by simulating breeding value vectors based on G matrices estimated from data, randomising these vectors across populations, and then calculating null values of statistics from G matrices that are calculated directly from the variances and covariances among randomised vectors. This calculation treats breeding values as quantities that are directly measurable, instead of predicted from G matrices that are themselves estimated from patterns of covariance among kin. The existing method thus neglects a major source of uncertainty in G matrices, which renders it anticonservative. We first suggest a correction to the method. We then apply the original and modified methods to a very simple instructive scenario. Finally, we demonstrate the use of both methods in the analysis of a real data set. Postprint |
Databáze: | OpenAIRE |
Externí odkaz: |