intensity: Estimate intensity function

In kbroman/xoi: Tools for Analyzing Crossover Interference

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])))

kbroman/xoi documentation built on May 1, 2023, 9:35 p.m.