Learn R Programming

qtl (version 1.70)

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

# S3 method for scanone
replacemap(object, map)

Value

The input object with the positions of LOD scores revised to match those in the input 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.

Author

Karl W Broman, broman@wisc.edu

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)
fake.f2 <- subset(fake.f2, chr=18:19)
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