pull_genoprobint: Pull genotype probabilities for an interval

In qtl2: Quantitative Trait Locus Mapping in Experimental Crosses

Pull genotype probabilities for an interval

Description

Pull out the genotype probabilities for a given genomic interval

Usage

pull_genoprobint(genoprobs, map, chr, 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)

Value

A list containing a single 3d array of genotype probabilities, like the input genoprobs but for the designated interval.

See Also

find_marker(), pull_genoprobpos()

Examples

iron <- read_cross2(system.file("extdata", "iron.zip", package="qtl2"))

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

qtl2 documentation built on May 29, 2024, 11:46 a.m.