unsmooth_gmap: Unsmooth genetic map

View source: R/unsmooth_gmap.R

unsmooth_gmapR Documentation

Unsmooth genetic map

Description

Performs the reverse operation of smooth_gmap(), in case one wants to go back to the original genetic map.

Usage

unsmooth_gmap(gmap, pmap, alpha = 0.02)

Arguments

gmap

Genetic map, as a list of numeric vectors; each vector gives marker positions for a single chromosome.

pmap

Physical map, as a list of numeric vectors; each vector gives marker positions for a single chromosome, with the same chromosomes and markers as gmap.

alpha

Proportion of mixture to take from constant recombination.

Details

An interval of genetic length d_g and physical length d_p is changed to have length (d_g - \alpha d_p r)/(1-\alpha) where r = L_g / L_p is the chromosome-specific recombination rate.

Value

A genetic map like the input gmap, but with the reverse operation of smooth_gmap() applied, provided that exactly the same physical map and alpha are used.

See Also

smooth_gmap()

Examples

iron <- read_cross2(system.file("extdata", "iron.zip", package="qtl2"))
gmap_adj <- smooth_gmap(iron$gmap, iron$pmap)
gmap_back <- unsmooth_gmap(gmap_adj, iron$pmap)


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