probs_to_grid: Subset genotype probability array to pseudomarkers on a grid

Description Usage Arguments Details Value See Also Examples

Description

Subset genotype probability array (from calc_genoprob() to a grid of pseudomarkers along each chromosome.

Usage

1

Arguments

probs

Genotype probabilities as output from calc_genoprob() with stepwidth="fixed".

grid

List of logical vectors that indicate which positions are on the grid and should be retained.

Details

This only works if calc_genoprob() was run with stepwidth="fixed", so that the genotype probabilities were calculated at a grid of markers/pseudomarkers. When this is the case, we omit all but the probabilities on this grid. Use calc_grid() to find the grid positions.

Value

Same list as input, but subset to just include pseudomarkers along a grid. The map attribute is similarly subset.

See Also

calc_grid(), map_to_grid()

Examples

1
2
3
4
5
6
7
grav2 <- read_cross2(system.file("extdata", "grav2.zip", package="qtl2geno"))
map_w_pmar <- insert_pseudomarkers(grav2$gmap, step=1)
probs <- calc_genoprob(grav2, map_w_pmar, error_prob=0.002)
sapply(probs, dim)
grid <- calc_grid(grav2$gmap, step=1)
probs_sub <- probs_to_grid(probs, grid)
sapply(probs_sub, dim)

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