sim_fin_chrom {junctions} | R Documentation |
Individual based simulation of the accumulation of junctions for a chromosome with an finite number of recombination sites.
sim_fin_chrom(pop_size, initial_heterozygosity, total_runtime, morgan, seed, R)
pop_size |
Population Size |
initial_heterozygosity |
Initial frequency of heterozygosity (e.g. in Hardy-Weinberg this equates to 2pq) |
total_runtime |
Maximum time after which the simulation is to be stopped |
morgan |
Size of the chromosome in Morgan |
seed |
Seed of the pseudo-random number generator |
R |
Number of genetic markers |
avgJunctions |
vector of the average number of junctions at time = [0, total_runtime] |
Thijs Janzen
N = 100 #population size R = 100 H_0 = 0.5 #initial heterozygosity maxT = 1000 #run time C = 1 #number of recombinations per meiosis seed = 42 v <- sim_fin_chrom(N, H_0, maxT, C, seed, R)