qtl2 (version 0.36)

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

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.

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.

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)