number_of_junctions {junctions} | R Documentation |
Calculate the average number of junctions after t generations, provided information on the initial heterozygosity, population size and the number of generations.
number_of_junctions(N = Inf, R = Inf, H_0 = 0.5, C = 1, t = 100)
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. |
The number of junctions for at time = t
Thijs Janzen
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)