sim.drift: Simulate Genetic Drift

In rknx/drift: Monte-Carlo Random Drift

Description

This function allows you to simulate the change in allelic frequencies due to random drift in a Wright-Fisher model

Usage

sim.drift(n.pop,n,freq,n.gen,n.loci,n.time,s.plot,s.table)

Arguments

n.pop	Number of effective population [Default:25]
n	Ploidy: Haploid=1, Diploid=2... [Default:1]
freq	initial frequency of allele (0-1) [Default:0.5]
n.gen	Number of generations to simulate [Default:100]
n.loci	Number of genes/loci tested [Default:1]
n.time	Number of simulations to perform [Default:1]
s.plot	Save the image of plot? (TRUE/FALSE) [Default:FALSE]
s.table	Save value table as tab delimited file? (TRUE/FALSE) [Default:FALSE]

Value

Frequency Plots. If save is true, image and tabel text file are saved in working dierctory.

Examples

sim.drift(30,2,0.25,100,10,3,TRUE,FALSE)
sim.drift(75,1,0.5,250,8)
sim.drift(n.pop=25, n=2, freq=0.01, n.gen=50, n.loci=5, n.time=1, s.plot=TRUE, s.table=TRUE)
sim.drift()

rknx/drift documentation built on May 26, 2019, 5:43 p.m.