intensity {xoi}

Estimate intensity function

Description

Estimate intensity function for a chromosome.

Usage

intensity(cross, chr = NULL, window = 2.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 position and intensity. The input argument window is kept as an attribute.

Author(s)

Il youp Kwak

See Also

coincidence()

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 <- intensity(x)
plot(out, type="l", lwd=2, ylim=c(0, max(out[,2])))