calculate_mat {junctions} | R Documentation |
Function that calculates the maximum time after hybridization after which the number of junctions can still be reliably used to estimate the onset of hybridization. This is following equation 15 in Janzen et al. 2018.
calculate_mat(N = Inf, R = Inf, H_0 = 0.5, C = 1)
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 maximum accurate time
Thijs Janzen
calculate_mat(N = Inf, R = 1000, H_0 = 0.5, C = 1)