qtl2 (version 0.36)

pull_genoprobint: Pull genotype probabilities for an interval

Description

Pull out the genotype probabilities for a given genomic interval

Usage

pull_genoprobint(genoprobs, map, chr, interval)

Value

A list containing a single 3d array of genotype probabilities, like the input genoprobs

but for the designated interval.

Arguments

genoprobs

Genotype probabilities as calculated by calc_genoprob().

map

The marker map for the genotype probabilities

chr

Chromosome ID (single character sting)

interval

Interval (pair of numbers)

See Also

find_marker(), pull_genoprobpos()

Examples

iron <- read_cross2(system.file("extdata", "iron.zip", package="qtl2"))
iron <- iron[,"8"]
gmap <- insert_pseudomarkers(iron$gmap, step=1)
pr <- calc_genoprob(iron, gmap, error_prob=0.002)

pr_sub <- pull_genoprobint(pr, gmap, "8", c(25, 35))