generateUniformGeneticMap: Generates a uniform genetic map.

Description Usage Examples

View source: R/recombination.R

Description

Generates a uniform genetic map by approximating 1 cm / Mbp. Only used for examples.

Usage

1

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
library("sim1000G")

examples_dir = system.file("examples", package = "sim1000G")
vcf_file = sprintf("%s/region.vcf.gz", examples_dir)
vcf = readVCF( vcf_file, maxNumberOfVariants = 100 ,
               min_maf = 0.12 ,max_maf = NA)

# For realistic data use the function readGeneticMap
generateUniformGeneticMap()

pdf(file=tempfile())
plotRegionalGeneticMap(seq(1e6,100e6,by=1e6/2))
dev.off()

adimitromanolakis/sim1000G documentation built on Sept. 29, 2021, 3:44 a.m.