Learn R Programming

qtl (version 1.39-5)

replacemap.scanone: Replace the genetic map in QTL mapping results with an alternate map

Description

Replace the positions of LOD scores in output from scanone with values based on an alternative map (such as a physical map), with pseudomarker locations determined by linear interpolation.

Usage

"replacemap"(object, map)

Arguments

object
An object of class "scanone", as output by the function scanone.
map
A list containing the alternative genetic map. All chromosomes in object should have corresponding chromosomes in map, and markers must be in the same order in the two maps. There must be at least two markers on each chromosome in map that appear in object.

Value

object with the positions of LOD scores revised to match those in the input map.

Details

The positions of pseudomarkers are determined by linear interpolation between markers. In the case of pseudomarkers beyond the ends of the terminal markers on chromosomes, we use the overall lengths of the chromosome in object and map to determine the new spacing.

See Also

replacemap.cross, est.map, replacemap.scantwo

Examples

Run this code
data(fake.f2)

origmap <- pull.map(fake.f2)
newmap <- est.map(fake.f2)
fake.f2 <- replacemap(fake.f2, newmap)
fake.f2 <- calc.genoprob(fake.f2, step=2.5)
out <- scanone(fake.f2, method="hk")
out.rev <- replacemap(out, origmap)

Run the code above in your browser using DataLab