freqdep: Numerical analysis of a frequency dependent seletion model

In KlausVigo/popgenSim: Population genetic simulations & numerical analysis

Description

This function performs numerical analysis of a frequency dependent selection model based on Rice (Evolutionary Theory). The fitnesses of the three genotypes in the model are as follows, where f(XX) denotes the frequency of the XX genotype: w(AA)=1-3*f(Aa)+3*f(aa); w(Aa)=1-s*f(Aa); w(aa)=1-3*f(Aa)+3*f(AA).

Usage

	freqdep(p0=0.01,s=0,time=100,show="p",pause=0)

Arguments

p0	Starting frequency for the A allele.
s	Parameter that determines the strength of selection against heterozygotes when they are common.
time	Number of generations to run the analysis.
show	Various options for plotting. "p" shows the frequency of A through time; "surface" plots the mean fitness as a function of p; "deltap" shows the change in p as a function of p; "cobweb" creates a cobweb plot showing p(t) by p(t+1). The default is show="p".
pause	Pause between generations. pause=0.01 (for instance) might smooth animation.

Value

The function creates one of several possible plots, depending on the value of show.

Author(s)

Liam Revell [e-mail]

See Also

sexratio

Examples

	freqdep(time=100)
	freqdep(s=1.5,time=100)
	freqdep(s=2,show="cobweb",time=100)

KlausVigo/popgenSim documentation built on May 8, 2019, 5:38 p.m.