Description Usage Arguments Value Author(s) See Also Examples
This function simulates genetic drift at a biallelic genetic locus with no selection and no mutation in a sexually reproducing diploid population or set of populations. It is essentially redundant with drift.selection
, but in which there is no difference in relative fitness among genotypes; however, it also allows the user to visualize heterozygosity or genetic variation through time - options that are not yet implemented in drift.selection
.
1 | genetic.drift(p0=0.5, Ne=20, nrep=10, time=100, show="p", pause=0.1, ...)
|
p0 |
Starting frequency for the A allele. |
Ne |
Effective population size. |
nrep |
Number of replicate simulations. |
time |
Total time, in number of generations, for the simulation. |
show |
Various options for plotting. |
pause |
Pause between generations. |
... |
optional arguments. In |
The function creates one of several possible plots, depending on the value of show
.
The function also invisibly returns an object of class "genetic.drift"
that can be printed or re-plotted by the user using corresponding print
and plot
methods. (See examples.)
Liam Revell liam.revell@umb.edu
drift.selection
, founder.event
, selection
1 2 3 4 5 6 | ## Not run:
genetic.drift()
object<-genetic.drift(p0=0.7,show="heterozygosity")
plot(object,show="genotypes")
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.