Pull genotype probabilities for a particular position
Pull out the genotype probabilities for a particular position (by name)
pull_genoprobpos(genoprobs, map = NULL, chr = NULL, pos = NULL, marker = NULL)
genoprobs |
Genotype probabilities as calculated by
|
map |
A map object: a list (corresponding to chromosomes) of
vectors of marker positions. Can also be a snpinfo object (data
frame with columns |
chr |
A chromosome ID |
pos |
A numeric position |
marker |
A single character string with the name of the position to pull out. |
Provide either a marker/pseudomarker name (with the argument marker)
or all of map, chr, and pos.
A matrix of genotype probabilities for the specified position.
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)
pmar <- find_marker(gmap, 8, 40)
pr_8_40 <- pull_genoprobpos(pr, pmar)
pr_8_40_alt <- pull_genoprobpos(pr, gmap, 8, 40)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.