time_error {junctions}R Documentation

Estimate the error in the time estimate

Description

Calculate the error in the estimate of the onset of hybridisation, following Equations 3 & 4 in the Supplementary information of Janzen et al. 2018.

Usage

time_error(J = NA, N = Inf, R = Inf, H_0 = 0.5, C = 1, t = 1, relative = TRUE)

Arguments

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

Value

the error in the time estimate

Author(s)

Thijs Janzen

Examples

  time_error(J = 100, N = Inf, R = 1000, H_0 = 0.5, C = 1)

[Package junctions version 1.1 Index]