Pull genotype probabilities for an interval
Pull out the genotype probabilities for a given genomic interval
pull_genoprobint(genoprobs, map, chr, interval)
genoprobs |
Genotype probabilities as calculated by
|
map |
The marker map for the genotype probabilities |
chr |
Chromosome ID (single character sting) |
interval |
Interval (pair of numbers) |
A list containing a single 3d array of genotype probabilities, like the input genoprobs
but for the designated interval.
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))Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.