coincidence {xoi}

Estimate coincidence function

Description

Estimate coincidence function for a chromosome.

Usage

coincidence(cross, chr = NULL, window = 5, ncalc = 500)

Arguments

cross

Cross object; must be a backcross. See qtl::read.cross() for format details.

chr

Chromosome to consider (only one is allowed). If NULL, the first chromosome is considered.

window

Window size

ncalc

Total number of points for calculations.

Value

Data frame with columns distance and coincidence. The input argument window is kept as an attribute.

Author(s)

Il youp Kwak

See Also

intensity(), est.coi()

Examples

map1 <- sim.map(103, n.mar=104, anchor=TRUE, include.x=FALSE, eq=TRUE)
x <- sim.cross(map1, n.ind=2000, m=6, type="bc")

out <- coincidence(x, ncalc=101)
plot(out, type="l", lwd=2, ylim=c(0, max(out[,2])))