find_marker: Find markers by chromosome position

Description Usage Arguments Details Value Examples

Description

Find markers closest to specified set of positions, or within a specified interval.

Usage

1
find_marker(map, chr, pos, interval = NULL)

Arguments

map

A map object: a list (corresponding to chromosomes) of vectors of marker positions.

chr

A vector of chromosomes

pos

A vector of positions

interval

A pair of positions (provide either pos or interval but not both)

Details

If pos is provided, interval should not be, and vice versa.

If pos is provided, then chr and pos should either be the same length, or one of them should have length 1 (to be expanded to the length of the other).

If interval is provided, then chr should have length 1.

Value

A vector of marker names, either closest to the positions specified by pos, or within the interval defined by interval.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
iron <- read_cross2(system.file("extdata", "iron.zip", package="qtl2geno"))

# find markers by their genetic map positions
find_marker(iron$gmap, c(8, 11), c(37.7, 56.9))

# find markers by their physical map positions (two markers on chr 7)
find_marker(iron$pmap, 7, c(44.2, 108.9))

# find markers in an interval
find_marker(iron$pmap, 16, interval=c(35, 80))

rqtl/qtl2geno documentation built on May 28, 2019, 2:36 a.m.