View source: R/recomb_map_fix_ends_chr.R
recomb_map_fix_ends_chr | R Documentation |
Given an existing recombination map and a chromosome length in base pairs, extrapolates the map to ensure all positions are covered, and shifts to ensure position one in basepairs corresponds to position 0 in genetic position.
Recombination rates are extrapolated from the first and last 10Mb of data by default (separately per end).
Therefore fixes the fact that common maps start genetic position zero at base pair position >> 1
and do not extend to ends (some SNPs from modern projects fall out of range without fixes).
recomb_map_fix_ends_chr(map, pos_length, pos_delta = 10000000L)
map |
A tibble with two columns: |
pos_length |
The length of the chromosome in base pairs. |
pos_delta |
The size of the window used to extrapolate recombination rates. |
The extrapolated recombination map, shifted so the first non-trivial position maps to the genetic distance expected from the extrapolated rate at the beginning, then added a first trivial position (pos=1, posg=0
) and final basepair position at length of chromosome and expected genetic position from end extrapolation.
recomb_map_simplify_chr()
to simplify recombination maps to a desired numerical accuracy.
library(tibble) # create a toy recombination map with at least 10Mb at each end map <- tibble( pos = c( 3L, 15L, 100L, 120L ) * 1e6L, posg = c( 0, 10.4, 90.1, 110 ) ) # and length pos_length <- 150L * 1e6L # apply function! map_fixed <- recomb_map_fix_ends_chr( map, pos_length ) # inspect map_fixed
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.