sim_fin_chrom {junctions}R Documentation

Individual Based Simulation of the accumulation of junctions

Description

Individual based simulation of the accumulation of junctions for a chromosome with an finite number of recombination sites.

Usage

sim_fin_chrom(pop_size, initial_heterozygosity, total_runtime, morgan, seed, R)

Arguments

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

Value

avgJunctions

vector of the average number of junctions at time = [0, total_runtime]

Author(s)

Thijs Janzen

Examples

  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)

[Package junctions version 1.1 Index]