time_error {junctions} | R Documentation |
Calculate the error in the estimate of the onset of hybridisation, following Equations 3 & 4 in the Supplementary information of Janzen et al. 2018.
time_error(J = NA, N = Inf, R = Inf, H_0 = 0.5, C = 1, t = 1, relative = TRUE)
J |
The number of junctions at time t |
N |
Population Size |
R |
Number of genetic markers |
H_0 |
Frequency of heterozygosity at t = 0 |
C |
Mean number of crossovers per meiosis (e.g. size in Morgan of the chromosome) |
t |
Inferred time |
relative |
Boolean flag, if TRUE: return the relative error, if FALSE: return error in generations |
the error in the time estimate
Thijs Janzen
time_error(J = 100, N = Inf, R = 1000, H_0 = 0.5, C = 1)