coal: Simulate and visualize a coalescent process

In popgenr: Accompaniment to Population Genetics with R: An Introduction for Life Scientists

Description

A function that provides a simple starting off point to simulate a coalescent process.

Usage

coal(length, number, muscale, reps, prnt)

Arguments

length	Length of sequence to simulate.
number	Number of starting lineages.
muscale	Mutation rate, scaled by 4Neu.
reps	Number of replicates to run.
prnt	Whether to print calculation output to the terminal (prnt=1), or not (prnt=0).

Examples

#assign genotype counts
length=250

number=8

muscale=6.25 # theta = 4Neu = 25, theta/4 = 6.25

reps=100

prnt=0

#run function
coal(length, number, muscale, reps, prnt)

popgenr documentation built on Feb. 25, 2021, 5:10 p.m.