View source: R/MultivDriftTest.R
MultivDriftTest | R Documentation |
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.
MultivDriftTest(
population1,
population2,
G,
Ne,
generations,
iterations = 1000
)
population1 |
data.frame with original measurements for the ancestral population |
population2 |
data.frame with original measurements for the derived population |
G |
ancestral G matrix |
Ne |
effective population size estimated for the populations |
generations |
time in generations separating both populations |
iterations |
number of simulations to perform |
list with the 95 drift and the range of the observed magnitude of morphological change
Each trait is estimated independently.
Ana Paula Assis
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
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.