MultivDriftTest function

Multivariate genetic drift test for 2 populations

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
)

Arguments

• 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

Returns

list with the 95

drift and the range of the observed magnitude of morphological change

Note

Each trait is estimated independently.

Examples

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)

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

Author(s)

Ana Paula Assis