Learn R Programming

evolqg (version 0.3-4)

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
)

Value

list with the 95

drift and the range of the observed magnitude of morphological change

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

Author

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

Run this code
# \donttest{
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)
# }

Run the code above in your browser using DataLab