number_of_junctions {junctions}R Documentation

Calculate the average number of junctions

Description

Calculate the average number of junctions after t generations, provided information on the initial heterozygosity, population size and the number of generations.

Usage

  number_of_junctions(N = Inf, R = Inf, H_0 = 0.5, C = 1, t = 100)

Arguments

N

Population size

R

Number of genetic markers

H_0

Initial frequency of heterozygosity (e.g. in Hardy-Weinberg this equates to 2pq)

C

Mean number of crossovers per meiosis (e.g. size in Morgan of the chromosome)

t

A single time point or a vector of time points for which the number of junctions is estimated.

Value

The number of junctions for at time = t

Author(s)

Thijs Janzen

Examples

  jt <-  number_of_junctions(N = 100, R = 1000, H_0 = 0.5, C = 1, t = 1000)
  jt2 <- number_of_junctions(N = 100, R = 1000, H_0 = 0.5, C = 1, t = 0:1000)

[Package junctions version 1.1 Index]