evolqg: MultivDriftTest – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

MultivDriftTest

Multivariate genetic drift test for 2 populations

Description

This function estimates populations evolving through drift from an ancestral population, given an effective population size, number of generations separating them and the ancestral G-matrix. It calculates the magnitude of morphological divergence expected and compare it to the observed magnitude of morphological change.

Usage

MultivDriftTest(
  population1,
  population2,
  G,
  Ne,
  generations,
  iterations = 1000
)

Arguments

`population1`	dataframe with original measurements for the ancestral population
`population2`	dataframe with original measurements for the derived population
`G`	ancestral G matrix
`Ne`	effective population size esitmated for the populations
`generations`	time in generations separating both populations
`iterations`	number of simulations to perform

Value

list with the 95 drift and the range of the observed magnitude of morphological change

Note

Each trait is estimated independently.

Author(s)

Ana Paula Assis

References

Hohenlohe, P.A ; Arnold, S.J. (2008). MIPod: a hypothesis testing framework for microevolutionary inference from patterns of divergence. American Naturalist, 171(3), 366-385. doi: 10.1086/527498

Examples

## Not run: 
data(dentus)
A <- dentus[dentus$species== "A",-5]
B <- dentus[dentus$species== "B",-5]
G <- cov(A)
MultivDriftTest(A, B, G, Ne = 1000, generations = 250)

## End(Not run)

evolqg

Tools for Evolutionary Quantitative Genetics

v0.2-8

MIT + file LICENSE

Authors

Ana Paula Assis, Diogo Melo, Edgar Zanella, Fabio Andrade Machado, Guilherme Garcia

Initial release

2020-11-14