estimate_time {junctions} | R Documentation |
Estimate the time since the onset of hybridization, following equation 14 in Janzen et al. 2018
estimate_time(J = NA, N = Inf, R = Inf, H_0 = 0.5, C = 1)
J |
The number of junctions |
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) |
The number of generations passed since the onset of hybridization
Thijs Janzen
J <- number_of_junctions(N = 100, R = 1000, H_0 = 0.5, C = 1, t = 200) estimate_time(J = J, N = 100, R = 1000, H_0 = 0.5, C = 1) #should be 200 again